我想在顶部用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)
答案 0 :(得分:2)
我找到了答案,我会把它放在可能遇到同样问题的其他人身上。
解决方案很简单:(只要raster和shapefile在同一个CRS中)
plot(r)
plot(shape, add = TRUE)