Укладка нескольких участков в R (Растр + линия высот) - PullRequest
0 голосов
/ 15 мая 2018

В основном я пытаюсь построить два графика на одном графике, и я хочу, чтобы оба масштаба (в х) были идентичны.Первый график - это растр, второй - высота каждой ячейки, когда растр обрезан горизонтальной линией.

Пока это код, который я написал.

library(raster)
library(rgdal)

T1_Domain <- raster("C:/Users/rm053078/Desktop/Report/2018/Results/EA benchmarking/Test1/GIS/test1DEM.asc")
T1_Section <- readOGR("C:/Users/rm053078/Desktop/Report/2018/Results/EA benchmarking/Test1/GIS/Section.shp")
T1_lat_profile <- extract(T1_Domain, T1_Section)


par(mfrow=c(2,1), xaxs = "i")

plot(T1_Domain,
     ylab = "y",
     xlab = "x",
     xlim = c(T1_Domain@extent@xmin,T1_Domain@extent@xmax))

points(400,50, pch = 4)
text(400,50,labels = "1", pos = 4)
points(600,50, pch = 4)
text(600,50,labels = "1", pos = 4)

plot(x=seq(from = T1_Domain@extent@xmin, to = T1_Domain@extent@xmax, length.out = length(T1_lat_profile[[1]])),
     y=T1_lat_profile[[1]],
     type = "l",
     xlab = "x",
     ylab = "z",
     xlim = c(T1_Domain@extent@xmin,T1_Domain@extent@xmax))

Хотякогда код выполняется, кажется, что я добился того, чего хочу. Изображение при выполнении .

Моя проблема заключается в том, что когда я запускаю это как блок rmarkdown, получаемый мной график отличается Изображение при выполнении в rmarkdown .

Есть идеи, как этого добиться?

...