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

Question

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

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

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)

Answer 1

我找到了答案，我将其放在这里供可能遇到同样问题的其他人使用。

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

plot(r)

plot(shape, add = TRUE)