将坐标参考系统设置为SpatialPolygon列表

Question

我想将坐标参考系统设置为名为SpP

的SpatialPolygon列表

library(sp)
Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Sr2 = Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))
Sr3 = Polygon(cbind(c(4,4,5,10,4),c(5,3,2,5,5)))
Sr4 = Polygon(cbind(c(5,6,6,5,5),c(4,4,3,3,4)), hole = TRUE)
Srs1 = Polygons(list(Sr1), "s1")
Srs2 = Polygons(list(Sr2), "s2")
Srs3 = Polygons(list(Sr3, Sr4), "s3/4")
SpP = SpatialPolygons(list(Srs1,Srs2,Srs3), 1:3)

我想为列表中的每个功能提供一个坐标参考系统。 CRS("+proj=longlat +datum=WGS84"

Answer 1

对于任何Spatial*对象，您可以使用proj4string函数设置CRS。

对于你的例子：

proj4string(SpP) <- CRS("+proj=longlat +datum=WGS84")

有关详细信息，请参阅?proj4string和?CRS。