如何将NAD 83与RGDAL包转换为纬度和经度？

Question

我有坐标，所有这些都应该位于DC中，但我无法弄清楚如何将它们从NAD 83转换为R中的纬度和经度。我正在使用RGDAL包中的spTransform()函数并获得错误关于非符合数据。

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的其他组合产生相同的错误。我相信错误是因为坐标太大了，但我不确定为什么会。坐标的文档如下：

值在马里兰州州平面仪表NAD 83地图投影中。

我是非常新的映射和预测，并赞赏任何帮助。

Answer 1

查找espg：4269：

http://spatialreference.org/ref/>

及其一个LOL长系统。所以你的大数字（米）太大了。

如果您在这些坐标中的数据中有一个shapefile，那么您可能会有一个.prj文件，它将具有投影规范，否则您必须在spatialReference.org上追逐：

NAD83的各种变化，这里还有“州平面”。我不是太确定哪个是哪个。 EPSG：代码更标准，然后有一堆ESRI：代码。 SR-ORG：其中的用户提供用户在网站上提供。

ESRI代码看起来足够最接近您给出的文本。让我们尝试：

> 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

.

在DC附近的任何地方？实际上，epsg：2804和epsg：3559给出相同的答案，可能更多的'标准'...