我正在尝试使用以下代码在多个栅格图中绘制相同的shapefile:
ras <- list.files("filepath", pattern = "\\.tif$", full=TRUE)
s <- stack(ras)
plot(s)
plot(world, add=TRUE)
[不幸的是,使用下面的代码添加shapefile似乎无效。
plot(world, add=TRUE)
我也在下面尝试过代码:
levelplot(cropped, legend.args = list(text = 'Elevation (m)'))
任何帮助将不胜感激。
谢谢
我曾经绘制许多栅格,并绘制数据点,并覆盖了世界地图以可视化我的结果。
我会先绘制我的点,然后添加地图。这听起来像您想要做的
尽管我没有两次使用plot命令。我会尝试
plot(s)
map(region="world", add=T)
因为我无权访问您的数据,所以我无法自己对此进行测试