Konvertieren zwischen wgs84 und nad83

Dies ist eine Anfängerfrage, die sich mit der Konvertierung zwischen wgs84 und nad83 befasst. Ich habe Luftbilder, die ich so genau wie möglich berechnen möchte. Es gibt mehrere Faktoren, die die üblichen Arten von Korrekturen schwierig machen, z. B. fragwürdige Parameter.

Daher brauche ich eine Möglichkeit, meinen Fortschritt zu überprüfen, wenn ich Korrekturen vornehme. Anfangs war es in Ordnung, nur gegen Google Earth zu prüfen, aber zu diesem Zeitpunkt ist der etwa 3 m große Unterschied zwischen meinen Daten (nad83) und Google Earth (wgs84) ein begrenzender Faktor geworden. Also muss ich nad83 in wgs84 umwandeln.

Ich habe rgdal in R und dotspatial in C # verwendet, um Projektionen durchzuführen, aber ich denke, der Wechsel zwischen nad83 und wgs84 wird nicht als "Projektion" angesehen (habe ich recht?). Der Versuch, ein Projekt mit "+ proj = NAD83 ellps = WGS84" aufzurufen, löst nur einen Fehler aus.

Es ist in Ordnung, eine 7-Parameter-Koordinatenänderung selbst zu schreiben, außer ich bin nicht sicher, wie ich eine der beiden Koordinaten in (X, Y, Z) -Koordinaten umwandeln soll (ich gehe davon aus, dass dies bedeutet, dass von Lat Long auf die Oberfläche der Ellipsoide konvertiert wird, auf denen sie basieren) ). Bei den bestehenden Fragen zu nad83 und wgs84 scheint es sich eher um die zu verwendende Methode zu handeln, als um die Beantwortung einer sehr grundlegenden Frage wie meiner, dh wie macht man das überhaupt?

Typischerweise liegen NAD83 und WGS84 innerhalb eines Meters voneinander. Ihre Bedenken in Bezug auf Unterschiede von 2,5 Fuß, die weniger als einen Meter betragen, deuten darauf hin, dass Sie diese Datumstransformation durchführen müssen .

Kurz gesagt, für diese Berechnung muss bekannt sein, wann die Koordinaten erfasst wurden, damit ihre Bewegung im Laufe der Zeit berücksichtigt werden kann (hauptsächlich aufgrund der Bewegungen der darunter liegenden tektonischen Platten und teilweise aufgrund der "mit Erdbeben verbundenen episodischen Bewegung"). Sie konvertieren zuerst die Rohdaten in äquivalente Koordinaten, die auf bestimmte Standarddaten wie 1997.0 oder 2002.0 bezogen sind. Anschließend schlagen Sie 14 Parameter nach, die (a) den Versatz, die Drehung und die Neuskalierung (dh eine Transformation mit sieben Parametern) beschreiben, die zum Ändern von Daten zu einem bestimmten Datum erforderlich sind, sowie (b) die ersten Ableitungen dieser Parameter in Bezug auf die Zeit. Das Multiplizieren von (b) mit der verstrichenen Zeit und das Addieren von (a) ergibt eine Annäherung erster Ordnung an die erforderliche Sieben-Parameter-Transformation, die Sie dann anwenden.

Obwohl das Verfahren nicht mathematisch schwierig ist, ist es eher aufwändig und von Einzelheiten abhängig, wie Ihre Daten ursprünglich gesammelt und verarbeitet wurden. Eine NOAA-Veröffentlichung, Transformationen zwischen NAD83 und WGS84 , bietet eine klare, lesbare Erklärung. Es wird empfohlen, die interaktive Software HTDP (Horizontal Time-Dependent Position) zu verwenden, die auf der NOAA-Website unter http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml verfügbar ist . Der Fortran-Quellcode besteht hauptsächlich aus den Koeffizienten des plattentektonischen Modells plus etwas Schnittstellencode, aber ich schätze, dass ungefähr 3.000 Zeilen tatsächlich Berechnungen durchführen. Es könnte ein erheblicher Schmerz sein, zu portieren.

Ich wäre mir nicht sicher, ob die Genauigkeit von Google Earth-Luftbildern bei einem Versatz von 2,5 Fuß liegt. Es könnte sich nur um eine nicht ganz orthogonale Fotografie handeln.

Was Sie vielleicht brauchen, ist nur eine affine Transformation, um Ihre Fotos auf ihre abzustimmen.