我在R中使用以下代码将我的UTM坐标转换为十进制度
utmcoor<-SpatialPoints(cbind(starling2$Easting,starling2$Northing), proj4string=CRS("+proj=utm +zone=51"))
longlatcoor<-spTransform(utmcoor,CRS("+proj=longlat"))
如何将结果('latlongcoor'输出)添加为原始数据框'starling2'中的两个新列?
答案 0 :(得分:1)
在SpatialPointsDataFrame中,坐标的行为类似于数据的“rownames”:
starling2$X <- coordinates(longlatcoor)[,1]
starling2$y <- coordinates(longlatcoor)[,2]