在 R 中的栅格图层上绘制 shapefile答案

【问题标题】：Plotting a shapefile on a raster layer in R在 R 中的栅格图层上绘制 shapefile
【发布时间】：2017-05-05 15:03:54
【问题描述】：

我想在顶部绘制一个带有 shapefile 点的栅格图层。我已经检查了以前的答案，但我仍然有问题。我可以毫无问题地分别绘制点 shapefile 和栅格图层，但不能一起绘制。

据我所知，它们应该在相同的投影和位置。

require(maptools) 

myproj <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
shape <- readShapeSpatial("directory/mypoints.shp", proj4string = CRS(myproj))

plot(r <- raster(listVI[200]))
plot(shape)

【问题讨论】：

标签： r plot raster shapefile

【解决方案1】：

我找到了答案，我会把它放在这里给其他可能遇到同样问题的人。

解决方法很简单：（只要 raster 和 shapefile 在同一个 CRS 中）

plot(r)

plot(shape, add = TRUE)

【讨论】：

也可以使用layer 和rasterVis。请参阅?rasterVis::levelplot 中的示例。