2017-07-21 46 views
0

我使用R中下面的代碼到我的UTM座標轉換爲十進制度轉換UTM爲十進制度R中

utmcoor<-SpatialPoints(cbind(starling2$Easting,starling2$Northing), proj4string=CRS("+proj=utm +zone=51")) 
longlatcoor<-spTransform(utmcoor,CRS("+proj=longlat")) 

如何我添加的結果(「latlongcoor」輸出)作爲兩個新列我的原始數據框'starling2'?

回答

0

在SpatialPointsDataFrame座標表現像 「rownames」 爲您的數據:

starling2$X <- coordinates(longlatcoor)[,1] 
starling2$y <- coordinates(longlatcoor)[,2]