Wie konvertiere ich NAD 83-Koordinaten mit dem rgdal-Paket in Breiten- und Längengrade?

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

  •  11-12-2019
  •  | 
  •  

Frage

Ich habe Koordinaten, die alle in DC liegen sollten, aber ich kann nicht herausfinden, wie ich sie von NAD 83 in Breiten- und Längengrade in R umrechnen kann.Ich verwende das spTransform() Funktion im rgdal-Paket und erhalten Sie eine Fehlermeldung über nicht konforme Daten.

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

Andere Kombinationen von proj4strings führen zum gleichen Fehler.Ich glaube, der Fehler liegt daran, dass die Koordinaten zu groß sind, aber ich bin mir nicht sicher, warum das so ist.Die Dokumentation für die Koordinaten finden Sie unten:

Die Werte sind in der Kartenprojektion Maryland State Plane Meter NAD 83 angegeben.

Ich bin sehr neu im Bereich Kartierung und Projektionen und bin für jede Hilfe dankbar.

War es hilfreich?

Lösung

Suchen Sie nach espg:4269:

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

und es ist ein Lat-Long-System.Ihre großen Zahlen (die Meter sind) sind also viel zu groß.

Wenn Sie irgendwo ein Shapefile mit Daten in diesen Koordinaten haben, verfügen Sie möglicherweise über eine .prj-Datei mit der Projektionsspezifikation. Andernfalls müssen Sie diese auf spatialreference.org suchen:

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

Es gibt verschiedene Variationen von NAD83 und hier und da gibt es auch „State Plane“.Ich bin mir nicht ganz sicher, welches genau welches ist.Das EPSG:Codes sind eher Standard, dann gibt es eine Reihe von Esri:Codes.Die sr-org:Diese werden vom Benutzer auf der Website bereitgestellt.

Der Esri-Code kommt dem von Ihnen angegebenen Text am nächsten.Lass es uns versuchen:

> 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

Irgendwo in der Nähe von DC?Tatsächlich geben epsg:2804 und epsg:3559 die gleichen Antworten und sind wahrscheinlich eher „Standard“ ...

Lizenziert unter: CC-BY-SA mit Zuschreibung
Nicht verbunden mit StackOverflow
scroll top