基本上我试图在同一个图中绘制两个图,我希望两个尺度(在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不同。
知道如何实现这个目标吗?