Konwersja między wgs84 i nad83

To pytanie dla początkujących dotyczy konwersji między wgs84 i nad83. Mam zdjęcia lotnicze, które próbuję obliczyć najdłuższe współrzędne jak najdokładniej. Istnieje kilka czynników, które utrudniają zwykłe typy poprawek, np. Wątpliwe parametry ... więc wiąże się to z pewnymi próbami i błędami.

Dlatego potrzebuję sposobu, aby sprawdzić swoje postępy podczas wprowadzania poprawek. Na początku samo sprawdzanie z Google Earth było w porządku, ale w tym momencie różnica około 2,5 stóp między moimi danymi (nad83) i Google Earth (wgs84) stała się czynnikiem ograniczającym. Więc muszę przekształcić nad83 na wgs84.

Używam rgdal w R i dotspatial w C # do robienia projekcji, ale myślę, że zmiana pomiędzy nad83 a wgs84 nie jest uważana za „projekcję” (mam rację?). Próba wywołania projektu za pomocą „+ proj = NAD83 ellps = WGS84” po prostu generuje błąd.

Nie przeszkadza mi pisanie 7-parametrowej zmiany współrzędnych, z tym wyjątkiem, że nie jestem pewien, jak przekonwertować jedną ze współrzędnych na (X, Y, Z) współrzędnych (co, jak zakładam, oznacza konwersję z długości długiej na powierzchnię elipsoid, na których są oparte ). Istniejące pytania dotyczące nad83 i wgs84 wydają się dotyczyć, której metody użyć, zamiast odpowiedzi na bardzo podstawowe pytanie, takie jak moje, tj. Jak w ogóle to zrobić?

Zazwyczaj NAD83 i WGS84 znajdują się w odległości jednego metra od siebie. Twoje obawy związane z różnicami 2,5 stóp, które są mniejsze niż metr, wskazują, że musisz wykonać tę transformację układu odniesienia .

W skrócie, obliczenia te wymagają wiedzy o tym, kiedy współrzędne zostały zebrane, aby można było uwzględnić ich ruch w czasie (przede wszystkim z powodu ruchów leżących pod nimi płyt tektonicznych, a częściowo z powodu „ruchu epizodycznego związanego z trzęsieniami ziemi”). Najpierw konwertujesz nieprzetworzone dane na równoważne współrzędne odniesione do określonych standardowych dat, takich jak 1997.0 lub 2002.0. Następnie przeglądasz 14 parametrów opisujących (a) przesunięcie, obrót i przeskalowanie (tj. Transformację siedmioparametrową) potrzebne do zmiany punktów odniesienia w danym dniu wraz z (b) pierwszymi pochodnymi tych parametrów w odniesieniu do czasu. Pomnożenie (b) przez upływający czas i dodanie do (a) daje przybliżenie pierwszego rzędu wymaganej transformacji siedmioparametrowej, którą następnie zastosujesz.

Chociaż procedura nie jest trudna matematycznie, jest raczej zaangażowana i uzależniona od szczegółów, w jaki sposób Twoje dane były pierwotnie gromadzone i przetwarzane. Artykuł NOAA, Transformacje między NAD83 i WGS84 , zawiera jasne, czytelne wyjaśnienie. Zaleca się korzystanie z interaktywnego oprogramowania HTDP, dostępnego na stronie internetowej NOAA pod adresem http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml . Kod źródłowy Fortran składa się przede wszystkim ze współczynników modelu płytowego tektonicznego oraz jakiegoś kodu interfejsu, ale szacuję, że około 3000 wierszy faktycznie wykonuje obliczenia. Przeniesienie może być znaczącym bólem.

Nie byłbym pewien co do dokładności zdjęć lotniczych Google Earth, gdyby był przesunięty do 2,5 stopy. Może to być po prostu kwestia nie dokładnie ortogonalnej fotografii.

To, czego możesz potrzebować, to po prostu transformacja afiniczna, aby dopasować swoje zdjęcia do ich.