NAD 83座標をrgdalパッケージで緯度と経度に変換するにはどうすればよいですか？

Question

私は座標系を持っていますが、そのすべてが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 

他のproj4ストリングの組み合わせでも同じエラーが発生します。私は、座標が大きすぎるのでエラーがあると信じていますが、なぜそうなるのか分かりません。座標のドキュメントは次のとおりです。

メリーランド州平面のメーター値はNAD 83マップ投影です。

私は非常にマッピングと予測に新しいです、そして、どんな助けもありがとうございます。

Answer 1

見上げてespg：4269：

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

とその緯度・長期のシステム。だからあなたの大きな数字（メートル）は大きすぎます。

：あなたはこれらの座標のデータとどこでもシェープファイルを持っていれば

あなたはそうあなたがspatialreference.orgでそれを追いかける必要があるでしょう、投影スペックを持つことになり、それとともに.prjファイルがある可能性があります

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

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は同じ回答ですが、おそらくもっと標準的なものでしょう...