Question

基本上我试图在同一个图中绘制两个图，我希望两个尺度（在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不同。

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

在R中堆叠多个图（光栅+高程线）

0 个答案: