将SpatialPolygonsDataFrame从地理坐标投影到UTM时出错

Question

我想将2012年秘鲁人口普查的坐标投影到另一个具有以下坐标参考系统的栅格图层：

  

＆＃34; + proj = utm + zone = 18 + south + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + units = m + no_defs＆＃34;。

帖子https://gis.stackexchange.com/questions/31743/projecting-sp-objects-in-r的回答在我的案例中不起作用。

crs(census)

  

CRS论点：    + proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0,0,0

newCRScensus=CRS("+proj=utm +zone=18 +south +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs ") 
projcensus=spTransform(census,crs=newCRScensus)

  

spTransform出错（人口普查，crs = CRS（＆＃34; + proj = utm + zone = 18 + south + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + units = m + no_defs＆＃34;））：     第二个参数需要是CRS类

我不理解错误信息，因为第二个参数确实是CRS类。

欢迎任何有关投射地理坐标的其他方式的建议。

Answer 1

spTransform的第二个参数称为CRSobj，而不是crs，因此正确的调用应为

projcensus=spTransform(census, newCRScensus)

或

projcensus=spTransform(census, CRSobj = newCRScensus)