Conversione tra wgs84 e nad83

Questa è una domanda da principiante che riguarda la conversione tra wgs84 e nad83. Ho delle foto aeree che sto cercando di calcolare le coordinate lat lunghe il più vicino possibile. Ci sono diversi fattori che rendono difficili i soliti tipi di correzioni, ad esempio parametri discutibili, ... quindi sono coinvolti alcuni tentativi ed errori.

Quindi ho bisogno di un modo per controllare i miei progressi mentre apporto correzioni. All'inizio solo il controllo su google earth andava bene, ma a questo punto la differenza di circa 2,5 piedi tra i miei dati (nad83) e google earth (wgs84) è diventata un fattore limitante. Quindi ho bisogno di trasformare nad83 in wgs84.

Ho usato rgdal in R e dotspatial in C # per fare proiezioni, ma penso che cambiare tra nad83 e wgs84 non sia considerato una "proiezione" (ho ragione?). Il tentativo di chiamare il progetto con "+ proj = NAD83 ellps = WGS84" genera un errore.

Sto bene scrivendo io stesso un cambio di coordinate a 7 parametri, tranne per il fatto che non sono sicuro di come convertire una in (X, Y, Z) coordinate (che presumo significhi convertire da lat long alla superficie degli ellissoidi su cui si basano ). Le domande esistenti su nad83 e wgs84 sembrano riguardare quale metodo usare, piuttosto che rispondere a una domanda molto semplice come la mia, ovvero come si fa a farlo?

In genere, NAD83 e WGS84 si trovano a un metro l'uno dall'altro. Le vostre preoccupazioni circa le differenze di 2,5 piedi, che sono meno di un metro, indicano che è necessario eseguire questa trasformazione di riferimento .

In breve, questo calcolo richiede la conoscenza di quando sono state raccolte le coordinate in modo che il loro movimento nel tempo possa essere giustificato (principalmente a causa dei movimenti delle placche tettoniche sottostanti e in parte a causa del "movimento episodico associato ai terremoti"). Per prima cosa converti i dati grezzi in coordinate equivalenti a cui fanno riferimento specifiche date standard, come 1997.0 o 2002.0. Quindi si cercano 14 parametri che descrivono (a) l'offset, la rotazione e il riscaling (ovvero una trasformazione a sette parametri) necessari per modificare i riferimenti in una determinata data insieme a (b) le prime derivate di tali parametri rispetto al tempo. Moltiplicando (b) per il tempo trascorso e aggiungendo a (a) si ottiene un'approssimazione del primo ordine alla trasformazione a sette parametri necessaria, che si applica quindi.

Sebbene la procedura non sia matematicamente difficile, è piuttosto coinvolta e dipende dai dettagli su come i dati sono stati originariamente raccolti ed elaborati. Un documento NOAA, Transformations Between NAD83 e WGS84 , fornisce una spiegazione chiara e leggibile. Ti consiglia di utilizzare il software interattivo "Posizione dipendente dal tempo orizzontale", HTDP, disponibile sul sito Web NOAA all'indirizzo http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml . Il codice sorgente di Fortran è costituito principalmente dai coefficienti del modello tettonico a placche più un certo codice di interfaccia, ma stima che circa 3.000 righe eseguano effettivamente calcoli. Potrebbe essere un dolore significativo port.

Non sarei sicuro della precisione delle foto aeree di Google Earth se si tratta di un offset di 2,5 piedi. Potrebbe essere solo una questione di fotografia non esattamente ortogonale.

Ciò di cui potresti aver bisogno è solo una trasformazione affine per allineare le tue foto alle loro.