R-将空间 - 螺旋形框架从地理坐标投射到UTM的错误

Question

我想投射2012年秘鲁人口普查的坐标，以与具有以下坐标参考系统的另一个栅格层一起工作：

" proj = utm Zone = 18 South Ellps = WGS84 TOWGS84 = 0,0,0,0,0,0,0,0,0 单位= m no_defs"。

帖子的答案https://gis.stackexchange.com/questions/31743/projecting-projecting-sp-bobjects-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(" proj = utm Zone = 18 South ellps = wgs84 towgs84 = 0,0,0,0,0,0,0,0,0,0 单位= m no_defs"((： 第二个参数需要是CRS类

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

欢迎对另一种投射地理坐标方式的任何建议。

Answer 1

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

projcensus=spTransform(census, newCRScensus)

或

projcensus=spTransform(census, CRSobj = newCRScensus)