我有一个坐标列表(纬度和经度)和两个栅格(一个有250x250平方米,另一个有1000x1000平方米网格);两种栅格具有相同的CRS(etrs89)。当我在250x250平方米的光栅上绘制点(纬度/经度)时,它可以工作,但对于1000x1000米的光栅,它不能。请帮忙!!
这是我正在使用的代码 -
library(rgdal)
library(raster)
points1 <- cbind(ds$lon, ds$lat)
points2 <- SpatialPoints(points1, proj4string=CRS("+proj=longlat +datum=WGS84"))
points3 <- spTransform(points2, projection(raster250))
#am just replacing this raster while trying to plot stuff on the 1000 m sq raster
plot(raster250)
points(points3, col="red")
两种栅格的CRS是 -
+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs