Como converter coordenadas NAD 83 em latitude e longitude com o pacote rgdal?

Question

Tenho coordenadas, todas elas deveriam estar localizadas em DC, mas não consigo descobrir como convertê-las de NAD 83 para latitude e longitude em R.estou usando o spTransform() função no pacote rgdal e obter um erro sobre dados não 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

Outras combinações de proj4strings geram o mesmo erro.Acredito que o erro seja porque as coordenadas são muito grandes, mas não sei por que isso aconteceria.A documentação das coordenadas está abaixo:

Os valores estão na projeção do mapa NAD 83 dos medidores do plano do estado de Maryland.

Sou muito novo em mapeamento e projeções e qualquer ajuda será apreciada.

Answer 1

Procure espg:4269:

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

e é um sistema lat-long.Portanto, seus números grandes (que são metros) são grandes demais.

Se você tiver um shapefile em qualquer lugar com dados nessas coordenadas, poderá ter um arquivo .prj com as especificações de projeção, caso contrário, você terá que persegui-lo em espacialreference.org:

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

Existem diversas variações no NAD83, e também há 'plano de estado' aqui e ali.Não tenho muita certeza de qual é qual.O EPSG:os códigos são mais padronizados, então há um monte de esri:códigos.O sr-org:aqueles são fornecidos pelo usuário no site.

O código esri parece mais próximo do texto que você forneceu.Vamos tentar:

> 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

Em algum lugar perto de DC?Na verdade, epsg:2804 e epsg:3559 dão as mesmas respostas e são provavelmente mais 'padrão'...