Convert UTM à degré décimal en R

Question

J'utilise le code suivant dans R pour convertir mes coordonnées UTM en degrés décimaux

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

comment ajouter les résultats (sortie « latlongcoor ») que deux nouvelles colonnes mon dataframe original 'starling2'?

Answer 1

Dans un SpatialPointsDataFrame les coordonnées se comportent comme les "rownames" pour vos données:

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