Convertendo entre wgs84 e nad83

Esta é uma pergunta para iniciantes que trata da conversão entre wgs84 e nad83. Tenho fotos aéreas que estou tentando calcular as coordenadas mais longas o mais próximo possível. Existem vários fatores que tornam complicados os tipos usuais de correções, por exemplo, parâmetros questionáveis, ... portanto, há alguma tentativa e erro envolvido.

Portanto, preciso de uma maneira de verificar meu progresso enquanto faço as correções. No começo, apenas a verificação no google earth era boa, mas nesse momento a diferença de 2,5 pés entre meus dados (nad83) e o google earth (wgs84) se tornou um fator limitante. Então, eu preciso transformar nad83 em wgs84.

Eu tenho usado rgdal em R e dotspatial em C # para fazer projeções, mas acho que mudar entre nad83 e wgs84 não é considerado uma "projeção" (estou certo?). Tentar chamar o projeto com "+ proj = NAD83 ellps = WGS84" apenas gera um erro.

Eu estou bem em escrever uma alteração de coordenada de 7 parâmetros, exceto por não ter certeza de como converter uma para coordenadas (X, Y, Z) (que eu suponho que significa converter de lat muito tempo para a superfície dos elipsóides em que eles se baseiam). ) As perguntas existentes sobre nad83 e wgs84 parecem ser sobre qual método usar, em vez de responder a uma pergunta muito básica como a minha, ou seja, como alguém faz isso?

Normalmente, o NAD83 e o WGS84 estão a um metro um do outro. Suas preocupações com diferenças de 2,5 pés, inferiores a um metro, indicam que você precisa executar essa transformação de dados .

Resumidamente, esse cálculo requer conhecimento de quando as coordenadas foram coletadas para que seus movimentos ao longo do tempo possam ser contabilizados (principalmente devido aos movimentos das placas tectônicas subjacentes e parcialmente devido ao "movimento episódico associado aos terremotos"). Você primeiro converte os dados brutos em coordenadas equivalentes referenciadas em datas padrão específicas, como 1997.0 ou 2002.0. Em seguida, você consulta 14 parâmetros que descrevem (a) o deslocamento, rotação e redimensionamento (ou seja, uma transformação de sete parâmetros) necessários para alterar dados em uma determinada data, juntamente com (b) as primeiras derivadas desses parâmetros em relação ao tempo. Multiplicar (b) pelo tempo decorrido e adicionar a (a) fornece uma aproximação de primeira ordem à transformação de sete parâmetros necessária, que você aplica.

Embora o procedimento não seja matematicamente difícil, ele é bastante envolvido e depende de detalhes de como seus dados foram originalmente coletados e processados. Um artigo da NOAA, Transformations Between NAD83 and WGS84 , fornece uma explicação clara e legível. Ele recomenda o uso do software interativo "Posição Horizontal Dependente do Tempo", HTDP, disponível no site da NOAA em http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml . O código fonte do Fortran consiste principalmente nos coeficientes do modelo tectônico de placas mais algum código de interface, mas eu estimo que cerca de 3.000 linhas efetivamente fazem cálculos. Pode ser uma dor significativa para o porto.

Eu não teria certeza sobre a precisão das fotos aéreas do Google Earth se o deslocamento fosse de 2,5 pés. Pode ser apenas uma questão de não ser exatamente uma fotografia ortogonal.

O que você pode precisar é apenas de uma transformação afim para alinhar suas fotos às delas.