NAD 83座標をRGDALパッケージで緯度と経度に変換する方法

StackOverflow https://stackoverflow.com//questions/11649580

  •  11-12-2019
  •  | 
  •  

質問

座標を持っていますが、すべてDCに配置する必要がありますが、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のマッピングプロジェクションです。

マッピングと予測には非常に新しい、そして助けが高く評価されています。

役に立ちましたか？

解決

ルックアップESPG：4269：

href="http://spatialReference.org/ref/epsg/4269/,refeg/4269/" REL="NOREFERRER"> http://spatialreeference.org/ref/epsg/4269/

とそのLAT-LONDシステム。それで、あなたの大きな数字（メートル）は大きすぎるようです。

これらの座標にデータを持つShapeFileを持っている場合は、投影スペックを持つ.prjファイルがある可能性があります。

http://spatialReference.org/ref/?search=nad83 +maryland&通り; srtext=検索

NAD83にバリエーションの盛り合わせがあり、ここに「州平面」もあります。私は正確にそれが間違いないものではありません。 EPSG：コードはより標準的なものであり、次にESRIの束があります。 SR-ORG：1つはサイトでユーザーが提供されています。

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は同じ答えを与え、おそらくより「標準」...

ライセンス： CC-BY-SA と 帰属
所属していません StackOverflow
scroll top