在R中堆叠多个图(光栅+高程线)

时间:2018-05-15 16:25:04

标签: r plot r-markdown raster

基本上我试图在同一个图中绘制两个图,我希望两个尺度(在x中)是相同的。第一个图是栅格,第二个是栅格被水平线切割时每个单元的高程。

到目前为止,这是我编写的代码。

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))

虽然代码执行时,似乎我实现了我想要的。 Image when executed

我的问题是当我将其作为rmarkdown chunk运行时,我得到的图表是Image when executed in rmarkdown不同。

知道如何实现这个目标吗?

0 个答案:

没有答案