¿Cómo convertir coordenadas NAD 83 a latitud y longitud con el paquete rgdal?

Question

Tengo coordenadas, todas las cuales deberían estar ubicadas en DC, pero no puedo entender cómo convertirlas de NAD 83 a latitud y longitud en R.estoy usando el spTransform() funciona en el paquete rgdal y obtiene un error sobre datos no conformes.

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

Otras combinaciones de proj4strings producen el mismo error.Creo que el error se debe a que las coordenadas son demasiado grandes, pero no estoy seguro de por qué.La documentación de las coordenadas se encuentra a continuación:

Los valores están en la proyección cartográfica NAD 83 de los metros del plano del estado de Maryland.

Soy muy nuevo en mapeo y proyecciones, y se agradece cualquier ayuda.

Answer 1

Busque espg:4269:

http://spatialreference.org/ref/epsg/4269/

y es un sistema de latitud larga.Entonces tus números grandes (que son metros) son demasiado grandes.

Si tiene un archivo de forma en algún lugar con datos en estas coordenadas, es posible que tenga un archivo .prj que tendrá la especificación de proyección; de lo contrario, tendrá que buscarlo en espacialreference.org:

http://spatialreference.org/ref/?search=nad83+maryland&srtext=Buscar

Hay una variedad de variaciones en NAD83, y también hay 'avión estatal' aquí y allá.No estoy muy seguro exactamente de cuál es cuál.El epsg:Los códigos son más estándar, luego hay un montón de Esri:códigos.La sr-org:Estos son proporcionados por el usuario en el sitio.

El código de Esri parece lo suficientemente cercano al texto que proporcionó.Intentemos:

> proj4string(nad83_coords)=CRS("+init=esri:102285")
> spTransform(nad83_coords,CRS("+init=epsg:4326"))
SpatialPoints:
             x        y
[1,] -77.03642 38.93586
[2,] -77.02437 38.92650
[3,] -77.01942 38.91656
[4,] -77.02128 38.92368
[5,] -77.02294 38.91614

¿En algún lugar cerca de DC?En realidad, epsg:2804 y epsg:3559 dan las mismas respuestas y probablemente sean más "estándar"...