Conversion entre wgs84 et nad83

Il s'agit d'une question pour les débutants qui concerne la conversion entre wgs84 et nad83. J'ai des photos aériennes pour lesquelles j'essaie de calculer les coordonnées latales aussi près que possible. Il existe plusieurs facteurs qui compliquent les types de corrections habituels, par exemple les paramètres discutables, ... il y a donc quelques essais et erreurs impliqués.

J'ai donc besoin d'un moyen de vérifier mes progrès lorsque je fais des corrections. Au début, la simple vérification par rapport à Google Earth était correcte, mais à ce stade, la différence d'environ 2,5 pieds entre mes données (nad83) et Google Earth (WGS84) est devenue un facteur limitant. J'ai donc besoin de transformer nad83 en wgs84.

J'ai utilisé rgdal en R et dotspatial en C # pour faire des projections, mais je pense que changer entre nad83 et wgs84 n'est pas considéré comme une "projection" (ai-je raison?). Tenter d'appeler le projet avec "+ proj = NAD83 ellps = WGS84" renvoie simplement une erreur.

Je suis d'accord pour écrire moi-même un changement de coordonnées à 7 paramètres, sauf que je ne sais pas comment convertir l'une ou l'autre en coordonnées (X, Y, Z) (ce qui, je suppose, signifie convertir de lat long à la surface des ellipsoïdes sur lesquels ils sont basés. ). Les questions existantes sur nad83 et wgs84 semblent concerner la méthode à utiliser, plutôt que de répondre à une question très basique comme la mienne, c'est-à-dire comment procéder?

En règle générale, NAD83 et WGS84 sont à moins d'un mètre l'un de l'autre. Vos préoccupations concernant les différences de 2,5 pieds, qui sont inférieures à un mètre, indiquent que vous devez effectuer cette opération transformation de référence .

En bref, ce calcul nécessite de savoir quand les coordonnées ont été collectées afin que leur mouvement dans le temps puisse être pris en compte (principalement en raison des mouvements des plaques tectoniques sous-jacentes et en partie en raison du "mouvement épisodique associé aux tremblements de terre"). Vous convertissez d'abord les données brutes en coordonnées équivalentes référencées à des dates standard spécifiques, telles que 1997.0 ou 2002.0. Ensuite, vous recherchez 14 paramètres décrivant (a) le décalage, la rotation et la mise à l'échelle (c'est-à-dire une transformation à sept paramètres) nécessaires pour changer les références à une date donnée avec (b) les premières dérivées de ces paramètres par rapport au temps. La multiplication (b) par le temps écoulé et l'ajout à (a) donne une approximation de premier ordre à la transformation à sept paramètres nécessaire, que vous appliquez ensuite.

Bien que la procédure ne soit pas mathématiquement difficile, elle est plutôt impliquée et dépendante des détails de la façon dont vos données ont été initialement collectées et traitées. Un article de la NOAA, Transformations entre NAD83 et WGS84 , fournit une explication claire et lisible. Il vous recommande d'utiliser le logiciel interactif "Horizontal Time-Dependent Position", HTDP, disponible sur le site Web de la NOAA à l' adresse http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml . Le code source de Fortran se compose principalement des coefficients du modèle tectonique des plaques et d'un code d'interface, mais j'estime qu'environ 3 000 lignes effectuent des calculs. Il pourrait être très difficile de mettre en communication.

Je ne serais pas sûr de la précision des photos aériennes de Google Earth si elle est inférieure à 2,5 pieds. Il pourrait simplement s'agir d'une photographie pas exactement orthogonale.

Ce dont vous pourriez avoir besoin est juste une transformation affine pour aligner vos photos avec les leurs.