Мне нужно преобразовать координаты WGS84 (широта / долгота) в NAD83 (широта / долгота) в R. Как мне сказали, должна быть очень небольшая разница в значении, учитывая достаточно десятичных знаков. Мой код выглядит следующим образом:
points <- data.frame(Lat=44.8579141141,Lon=-98.5289141141)
coordinates(points) <- ~Lon+Lat
proj4string(points) <- CRS("+proj=longlat")
wgs84 <- spTransform(points, CRS("+init=epsg:4326"))
nad83 <- spTransform(wgs84, CRS("+init=epsg:4269"))
Я считаю, что это правильный код, но мои значения не отличаются между wgs84@coords и nad83@coords. Я прошу прощения за общий характер этого вопроса, но я ' Мне интересно, является ли это правильным кодом для выполнения работы, даже если нет очевидной разницы, поскольку окончательные координаты только go из 5 цифр.