为什么不在光栅上绘制点?

时间:2018-03-30 21:40:27

标签: r spatial raster r-raster

我有一个坐标列表(纬度和经度)和两个栅格(一个有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

0 个答案:

没有答案