Преобразование UTM в десятичную степень в R

Я использую следующий код в R для преобразования моих координат UTM в десятичные градусы

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

как добавить результаты (вывод «latlongcoor») в виде двух новых столбцов в исходный фрейм данных «starling2»?


person TAW    schedule 21.07.2017    source источник


Ответы (1)


В SpatialPointsDataFrame координаты ведут себя как «имена строк» ​​для ваших данных:

starling2$X <- coordinates(longlatcoor)[,1]
starling2$y <- coordinates(longlatcoor)[,2]
person Elio Diaz    schedule 19.09.2017