Как преобразовать координаты NAD 83 в широту и долготу с помощью пакета rgdal?

У меня есть координаты, все они должны быть расположены в DC, но я не могу понять, как преобразовать их из NAD 83 в широту и долготу в R. Я использую функцию spTransform()в пакете rgdal и получаю сообщение об ошибке данные, не соответствующие -.

library(rgdal)
nad83_coords <- data.frame(x=c(396842.6, 397886.9, 398315.5, 398154.3, 398010.3), y=c(140887.1, 139847.0, 138743.9, 139534.5, 138697.3))
coordinates(nad83_coords) <- c('x', 'y')
proj4string(nad83_coords) <- CRS("+init=epsg:4269")
Error in `proj4string<-`(`*tmp*`, value = <S4 object of class "CRS">) : 
  Geographical CRS given to non-conformant data: 398315.5 140887.1

Другие комбинации proj4strings приводят к той же ошибке. Я считаю, что ошибка связана с тем, что координаты слишком велики, но я не уверен, почему это может быть. Документация по координатам ниже:

Values are in the Maryland State Plane meters NAD 83 map projection.

Я очень новичок в картографировании и проекциях, и любая помощь приветствуется.

7
задан Erik Shilts 25 July 2012 в 12:20
поделиться