Suite

Convertir les coordonnées géographiques en zone UTM 29N avec Javascript

Convertir les coordonnées géographiques en zone UTM 29N avec Javascript


J'ai un problème pour convertir les lat longs en UTM avec l'API javscript ArcGIS.

Disons que j'ai une carte, basée sur la zone UTM 29N, également connue sous le nom de WKID 32629.

J'ai un ensemble de coordonnées connu d'un GPS, dans ce format :

Latitude : 62,220596 Longitude : -6,564331

Les coordonnées géographiques sont dans WKID 4326.

Pour convertir à l'aide de javascript, j'ai essayé d'utiliser la méthode de projet du service de géométrie.

Si nous regardons ce post dans le forum ArcGIS, il s'agit du même problème. Ici, j'ai essayé d'utiliser le même code et d'ajuster le wkid de la sortie spatialReference à32629, et vous voyez qu'il ne renvoie queNaN, essayez dans le JSfiddle de changer le wkid en 4326, et vous verrez les coordonnées géographiques.

Modifier 1

En regardant leHTTPdemandes faites au serveur, je vois que cette url est demandée (en utilisant maintenant la carte des îles Féroé)

http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Geometry/GeometryServer/project?f=json&outSR=32629&inSR=4326&geometries={"geometryType":"esriGeometryPoint","geometries":[{"x":62.01387 ,"y":-6.7866,"spatialReference":{"wkid":4326}}]}&callback=dojo.io.script.jsonp_dojoIoScript2._jsonpCallback

Cela ressemble à une URL correcte, contenant outSR, inSR, les géométries et le rappel. Cependant, simuler des requêtes vers cette URL me donne soit un500 Erreur de serveur interne, ou cette chaîne :

dojo.io.script.jsonp_dojoIoScript2._jsonpCallback({"error":{"code":400,"message":"Impossible de terminer l'opération du projet.","details":["Error projecting geometries"]}});

Des idées?


Je n'ai pas regardé votre code, mais dans l'URL que vous avez fournie, vos X et Y sont inversés. Pour les îles Féroé, X doit être de -6,7866 (longitude) tandis que Y doit être de 62,01387 (latitude). Voici une URL fonctionnelle de ce que vous essayez de faire :

http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Geometry/GeometryServer/project?f=json&outSR=32629&inSR=4326&geometries={"geometryType":"esriGeometryPoint","geometries":[{"x":- 6.7866,"y":62.01387,"spatialReference":{"wkid":4326}}]}&callback=dojo.io.script.jsonp_dojoIoScript2._jsonpCallback

Cela renvoie correctement :

dojo.io.script.jsonp_dojoIoScript2._jsonpCallback({"geometryType":"esriGeometryPoint","geometries":[{"x":615863.053581134,"y":6877701.89119242}]});

En ce qui concerne les 500, il me semble que j'en reçois assez fréquemment avec sampleserver3 ce soir également. Essayez plutôt d'utiliser sampleserver1, celui-ci semble être plus stable.


En modifiant le WKID dans le code sur JSFiddle, vous demandez au service de géométrie de se projeter d'une zone UTM à une autre. Cela n'a pas vraiment de sens...

Essayez ceci : http://jsfiddle.net/swingley/Yab7E/

Le code là-bas prend un point de clic de carte (qui est Web Mercator), le convertit en lat, long, envoie le point lat, long au service de géométrie et dit le projeter sur WGS84 UTM Zone 29N.

Edit : Voici une URL générée à partir de ce violon qui renvoie un point dans la zone UTM 29N :

http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Geometry/GeometryServer/project?f=json&outSR=32629&inSR=4326&geometries={%22geometryType%22%3A%22esriGeometryPoint%22%2C%22geometries%22%3A% 5B{%22x%22%3A-7.080688421874865%2C%22y%22%3A62.19107532141271%2C%22spatialReference%22%3A{%22wkid%22%3A4326}}%5D}&callback=dojo.io.script.jsonp_dojoIoScript3. _jsonpRappel

J'ai mis en sac le formulaire de code ici et cela fonctionne très bien. Côté client intégré, vous n'avez pas réellement besoin d'un appel au serveur. Le côté serveur intégré dans un script, ce que j'ai également fait, garde tout encapsulé.


Vous demandez une coordonnée UTM Zone 29N pour un point au Michigan.
La zone UTM 29N n'existe pas dans le Michigan, vous êtes en dehors des limites du système de coordonnées.
Pour démontrer, sur votre exemple d'origine, effectuez un zoom arrière 14 fois, puis cliquez sur environ un demi-pouce sous le M dans "Mercator". Vous obtiendrez un résultat car vous êtes maintenant dans les limites du système de coordonnées UTM Zone 29N.
Ou essayez de passer à la zone 16N, wkid 32616, qui existe dans le Michigan et vous verrez que votre exemple de code renvoie un point sans aucun problème.
En plus de cela, ESRI doit vraiment corriger la documentation de l'erreur. Maintenant que je sais quel était le problème, l'erreur a du sens, mais pouvez-vous vraiment comprendre à partir de cette erreur que vous essayez de projeter en dehors des limites du système de coordonnées ?


La solution de 10 secondes pour celui-ci consiste à copier/coller votre latitude/longueur sur http://www.hamstermap.com - vous les convertirez instantanément en UTM.


1. Vue d'ensemble

Le chapitre 1 a décrit plusieurs des propriétés distinctives des données géographiques. La première est que les données géographiques sont nécessairement généralisées et que cette généralisation a tendance à varier avec l'échelle. Une deuxième propriété distinctive est que la forme complexe et presque sphérique de la Terre complique les efforts pour spécifier des positions exactes sur la surface de la Terre. Ce chapitre explore les implications de ces propriétés en éclairant les concepts d'échelle, de géométrie de la Terre, de systèmes de coordonnées, les "références horizontales" qui définissent la relation entre les systèmes de coordonnées et la forme de la Terre, et les différentes méthodes de transformation des données de coordonnées entre les grilles 3D et 2D, et d'une donnée à une autre.

Par rapport au chapitre 1, le chapitre 2 peut sembler long, technique et abstrait, en particulier pour ceux pour qui ces concepts sont nouveaux.

Objectifs

Les étudiants qui terminent avec succès le chapitre 2 devraient être capables de :

  1. démontrer votre capacité à spécifier des emplacements géospatiaux à l'aide de coordonnées géographiques
  2. convertir les coordonnées géographiques entre deux formats différents
  3. expliquer le concept de référence horizontale
  4. calculer le changement d'un emplacement de coordonnées dû à un changement d'un système de référence horizontal à un autre
  5. estimer l'ampleur du "décalage de référence" associé à l'ajustement de NAD 27 à NAD 83
  6. reconnaître le type de transformation approprié pour géo-enregistrer deux ensembles de données ou plus
  7. décrire les caractéristiques du système de coordonnées UTM, y compris sa base dans la projection cartographique Transverse Mercator
  8. tracer les coordonnées UTM sur une carte
  9. décrire les caractéristiques du système SPC, y compris la projection cartographique sur laquelle il est basé
  10. convertir les coordonnées géographiques en coordonnées SPC
  11. interpréter des diagrammes de distorsion pour identifier les propriétés géométriques de la sphère qui sont préservées par une projection particulière et
  12. classer les graticules projetés par famille de projection.

"Essaye ça!" Activités

Prenez une minute pour terminer l'une des activités Essayez ceci que vous rencontrez tout au long du chapitre. Ce sont des exercices amusants et stimulants pour vous aider à mieux comprendre les idées présentées dans le chapitre.


Pourquoi utiliser les coordonnées UTM

UTM Fournit une relation de distance constante n'importe où sur la carte. Dans les systèmes de coordonnées angulaires comme la latitude et la longitude, la distance couverte par un degré de longitude diffère à mesure que vous vous rapprochez des pôles et n'est égale qu'à la distance couverte par un degré de latitude à l'équateur. Étant donné que la navigation terrestre se fait dans une très petite partie du monde à un moment donné à l'aide de cartes à grande échelle. Le système UTM permet de lier directement le système de numérotation des coordonnées à un système de mesure de distance.


Militaire

Le BUCS est le principal moyen de calculs de conversion. Le formulaire DA 5600-R (Calcul - Conversion des coordonnées géographiques en coordonnées UTM (BUCS)) est utilisé pour les calculs de conversion. (Voir Figure 11-1). (Une copie reproductible de ce formulaire est incluse dans la section Formulaires vierges de ce livre.) Le sphéroïde correct doit être déterminé.

une. Le formulaire DA 5600-R suit le même format de base que les autres formulaires BUCS. Le haut du formulaire est destiné à l'enregistrement des données administratives. Les données dans ce domaine comprennent les suivantes :

b. La partie suivante du formulaire fournit des notes pour les opérations spécifiques du programme et d'autres notes nécessaires pour remplir ce formulaire.

c. La partie suivante du formulaire est divisée en deux colonnes principales. La colonne de gauche fournit des instructions pour que l'ordinateur utilise le programme. La colonne de droite est pour l'enregistrement des données. Deux conversions peuvent être calculées sur chaque formulaire.

ré. La colonne de gauche est divisée en trois colonnes plus petites : ÉTAPE, PROMPT et PROCÉDURE. La colonne STEP est la séquence numérique que l'ordinateur utilise au fur et à mesure qu'il descend le formulaire. La colonne PROMPT indique à l'ordinateur ce qui apparaîtra sur l'écran BUCS à chaque étape particulière. La colonne PROCÉDURE indique à l'ordinateur l'action qu'il entreprendra à chaque étape ou invite particulière.

e. Les procédures de calcul du formulaire DA 5600-R sont présentées dans le tableau 11-1.

11-2. CONVERSION DES COORDONNÉES UTM EN CALCUL DE COORDONNÉES GÉOGRAPHIQUES

une. Les coordonnées géographiques sont nécessaires pour le calcul de certaines méthodes d'observation astro. Souvent, les coordonnées UTM seront les seules données disponibles. Dans ces situations, l'arpenteur doit convertir les coordonnées de la grille UTM en coordonnées géographiques.

b. Le calcul de la conversion des coordonnées UTM en coordonnées géographiques est effectué sur le formulaire DA 5601-R (Calcul - Conversion des coordonnées UTM en coordonnées géographiques (BUCS)) à l'aide du BUCS. (Voir Figure 11-2.) Les procédures de calcul du formulaire DA 5601-R sont présentées dans le Tableau 11-2.

une. Le point d'origine des coordonnées ordonnées est l'équateur. Le point d'origine des coordonnées est est le méridien central (CM) de chaque zone de grille. Le CM reçoit une fausse valeur de 500 000. Les coordonnées est augmentent lorsque vous vous déplacez vers l'est et diminuent lorsque vous vous déplacez vers l'ouest.

b. La figure 11-3 montre deux zones UTM adjacentes - 14 et 15. Visualisez le point P situé dans la zone 14. Les coordonnées indiquées pour le point P sont (800 000 et 3 700 000). Si les coordonnées du point P devaient être transformées dans la grille adjacente (zone 15), l'action entreprise équivaudrait à superposer la zone de grille 15 sur la zone de grille 14 comme indiqué sur la figure 11-3. En fait, la transformation des coordonnées n'implique que la continuation mathématique de la grille adjacente à la grille utilisée et les corrections ultérieures de localisation causées par le changement de référence nord de la grille. Dans la figure 11-3, bien que l'emplacement du point P au sol soit inchangé, les valeurs des coordonnées changeront avec la transformation. Les valeurs au niveau de la zone 15 seraient inférieures à 100 000 mètres en direction est et supérieures à 3 700 000 mètres en direction nord. A partir de cette figure, il est évident que lors de la transformation des coordonnées, les coordonnées est changeront sensiblement en raison du système de numérotation des coordonnées de chaque zone et du changement de référence nord de la grille. Il n'y aura qu'un petit changement dans les coordonnées nord basé uniquement sur la différence dans la référence nord de la grille.

11-4. CALCULS DE TRANSFORMATION

Le BUCS est le principal moyen de transformation informatique. Le formulaire DA 5602-R (calcul - transformation de zone à zone - coordonnées de grille UTM et azimut de grille UTM (BUCS)) est utilisé pour les calculs. (Voir Figure 11-4.) (Une copie reproductible de ce formulaire est incluse dans la section Formulaires vierges de ce manuel.) Ce formulaire suit le même format de base que les formulaires BUCS décrits précédemment. Le tableau 11-3 montre les instructions pour le calcul du formulaire DA 5602-R.

une. La terre n'est pas une sphère, mais un ellipsoïde (ou sphéroïde), légèrement aplati aux pôles et légèrement bombé à l'équateur. Un ellipsoïde est une figure mathématique générée par la révolution d'une ellipse autour d'un de ses axes. (Voir Figure 11-5.) L'ellipsoïde est utilisé comme surface de référence pour la réduction mathématique des données géodésiques et cartographiques.

b. Une donnée est un ensemble de grandeurs numériques ou géométriques spécifiant le système de coordonnées de référence utilisé pour le contrôle géodésique dans le calcul des coordonnées de points sur la terre. Les datums peuvent avoir une étendue globale ou locale. Une référence locale définit un système de coordonnées qui n'est utilisé que sur une région d'étendue limitée. Une donnée globale spécifie le centre de l'ellipsoïde de référence qui doit être situé au centre de masse de la Terre.

c. La Defense Mapping Agency (DMA) TM 8358.1 identifie cinq ellipsoïdes préférés aux fins de la production de nouvelles cartes. Le DMA a identifié 21 référentiels géodésiques préférés dans les différents ellipsoïdes. À l'heure actuelle, plus de 1 000 références sont utilisées dans le monde. Actuellement, le TM 5-241-1 est utilisé pour déterminer les sphéroïdes et les datums, mais le DMA TM 8358.1 le remplace.

ré. Les coordonnées d'un point basé sur un datum montreront une différence par rapport au même point basé sur un autre datum. Dans les cas extrêmes, la différence dépassera 1 000 mètres d'erreur radiale. La Figure 11-6 montre un exemple de la différence entre les datums pour la même zone cartographique.

e. Le problème de la transformation des informations d'arpentage entre les datums existe parce que chacun des ellipsoïdes prend une dimension légèrement différente pour la Terre (le demi-grand axe, le demi-petit axe et l'aplatissement de l'ellipsoïde). Ce problème est encore aggravé par les références dans chacun des ellipsoïdes. Chaque système de référence a son propre point d'origine sur lequel le réseau d'arpentage principal dans une zone est ajusté.

Le module BUCS DDCT comprend trois programmes : le programme 14 (références répertoriées), le programme 15 (Gauss-Kruger) et le programme 16 (défini par l'utilisateur).

une. Programme 14--Références répertoriées. Transforme les coordonnées UTM et/ou géographiques (gee) entre les références répertoriées dans DMA TR 8350.2. Le programme effectue également des transformations de zone à zone lorsqu'il est sélectionné.

b. Programme 15--Gauss-Kruger. Transforme les coordonnées Gauss-Kruger (GK) (UTM et/ou géo) en coordonnées UTM. Les références GK sont principalement utilisées dans les pays du Pacte de Varsovie.

c. Programme 16--Défini par l'utilisateur. Transforme les coordonnées UTM et/ou géo entre un datum défini par l'utilisateur et des datums répertoriés ou un autre datum défini par l'utilisateur. Le DMA ou un QG supérieur diffusera les paramètres de référence définis par l'utilisateur si nécessaire.

Noter. Les décalages de référence effectués entre WGS-72 et toute autre référence à l'aide du programme 14 ou 15 ne produisent pas les mêmes solutions que le programme 16 lors de la saisie des paramètres de transformation en tant que références définies par l'utilisateur. Cela se produit parce qu'une formule plus précise pour WGS-72 est utilisée dans les programmes 14 et 15. Le programme 16 utilise les formules abrégées de Molodensky et les décalages de référence pour WGS-72. Le programme 16 est utilisé pour les transformations de référence uniquement lorsqu'il est dirigé par un QG supérieur et des paramètres spécifiques sont fournis.

Noter. Le programme DDCT imprime un formulaire vierge pour chaque programme lorsqu'il est connecté à l'imprimante BUCS.

La fonction principale du module BUCS DDCT est le calcul des transformations de coordonnées entre les références désignées par l'utilisateur. Les programmes utilisent la méthode abrégée de Molodensky pour transformer les systèmes géodésiques locaux en WGS-72 et en WGS-84.

11-8. FONCTIONNALITÉS DE SAUVEGARDE, DE HAUT DE FICHIER, DE CORRECTION D'ERREUR ET D'ANNULATION

Les capacités SAUVEGARDE, HAUT DE FICHIER, CORRECTION D'ERREUR et ABANDON utilisées dans le module d'enquête se trouvent également dans le module DDCT. (Voir chapitre 12.) Ces fonctions sont identiques à une exception près. L'opérateur ne peut pas revenir à l'invite MODULE DÉSIRÉ : ## et sélectionner un module différent. Si un module différent est souhaité, utilisez la fonction ABORT.

Les programmes BUCS DDCT utilisent les messages d'erreur ou d'avertissement indiqués dans le Tableau 11-4.

11-10. PROGRAMMES DDCT

Les programmes BUCS DDCT sont conçus pour être conviviaux. L'utilisateur est guidé à travers chaque programme par une série d'invites simples et compréhensibles. Les invites rencontrées à l'écran doivent être suivies étape par étape pour déterminer les données ou la solution correctes dans ce programme. Les données et les calculs DDCT peuvent être enregistrés sur les formulaires de calcul appropriés développés pour être utilisés avec le BUCS ou enregistrés directement par l'imprimante lors de l'utilisation du BUCS Special. Les abréviations des invites sont expliquées dans les instructions du programme où elles apparaissent.

une. Suivez les étapes du Tableau 11-5 pour exécuter les programmes DDCT lorsque le numéro de module est connu.

b. Pour exécuter les programmes DDCT lorsque le numéro de module est inconnu, suivez la procédure du Tableau 11-6.

Noter. Lors de l'exécution des programmes DDCT, le BUCS demande MODULE DÉSIRÉ : 00. Ce que le BUCS demande est le numéro du programme.

une. Objectif. Ce programme transforme les coordonnées entre deux références quelconques dont les paramètres de transformation ellipsoïde et abrégé de Molodensky sont stockés en permanence dans la ROM BUCS DDCT. Une liste de ces références figure à l'annexe E. Le programme dispose de quatre fonctions distinctes disponibles pour les transformations de coordonnées.

(1) UTM à UTM. Cette fonction convertit les coordonnées UTM d'un datum en coordonnées UTM dans le même datum ou dans un autre. Les transformations de zone à zone peuvent être incluses dans le cadre d'une transformation UTM à UTM. Dans de tels cas, la zone de grille désignée pour la donnée de sortie doit être adjacente à la zone de grille désignée pour la donnée d'entrée.

(2) UTM vers GEO. Cette fonction convertit les coordonnées UTM d'un datum en coordonnées géographiques du même datum ou d'un autre.

(3) GEO vers UTM. Cette fonction convertit les coordonnées géographiques d'un datum en coordonnées UTM du même datum ou d'un autre.

(4) GEO à GEO. Cette fonction convertit les coordonnées géographiques d'un datum en coordonnées géographiques du même datum ou d'un autre.

b. Capacités. Ce programme convertira les coordonnées d'un datum dans le même ou un autre datum. Les coordonnées d'entrée et de sortie peuvent être n'importe quelle combinaison de coordonnées UTM et géographiques. De plus, ce programme peut être utilisé pour calculer des transformations de zone à zone lorsque les coordonnées d'entrée et de sortie sont données dans la grille UTM. La zone de grille désignée dans la donnée de sortie doit être adjacente à la zone de grille désignée dans la donnée d'entrée. Voir le tableau des zones de grille à l'annexe F pour déterminer les zones de grille adjacentes. Lors de l'exécution de transformations de zone à zone, plusieurs conditions d'erreur ou d'avertissement existent.

(1) Les transformations de coordonnées de zone à zone peuvent produire des abscisses négatives. Si la transformation se fait vers une zone de quadrillage à l'ouest et que le point est trop à l'ouest de la limite de la zone, une abscisse négative est produite. Les distances qui causent des abscisses négatives sont fonction de la latitude (nord ou sud) et sont répertoriées dans le tableau 11-7.

(2) Les transformations de coordonnées de zone à zone peuvent produire le message d'avertissement d'erreur IMAGE OVFL pour les coordonnées est. Si la transformation se fait vers une zone de quadrillage à l'ouest et que le point est trop à l'est de la limite de la zone, un message d'erreur ou d'avertissement IMAGE OVFL s'affiche. Le message d'erreur s'affiche si la coordonnée Est dépasse ou égale la limite légale d'un million de mètres. Les distances à l'origine du message d'erreur ou d'avertissement sont fonction de la latitude (nord ou sud) et sont répertoriées dans le tableau 11-8.

(3) BUCS DDCT ne reconnaît pas les zones de grille de taille anormale 31, 33, 35 et 37 entre 72 & 176 et 84 & 176 latitude nord. De plus, il ne reconnaît pas l'absence de zones de grille 32, 34 et 36 entre les mêmes latitudes. Il n'y a pas de grandes masses continentales dans cette région à l'exception du groupe d'îles norvégiennes de Svalbard.

Noter. Lorsque la coordonnée à transformer est proche de l'équateur, une ordonnée négative ou un message d'erreur IMAGE OVFL peut s'afficher. Cela se produit lorsque le décalage de référence place la coordonnée dans un hémisphère différent, par exemple, l'équateur se déplace du nord au sud du point et vice versa. Pour les zones de quadrillage au nord de l'équateur, l'équateur a une valeur de zéro. Si le décalage du datum de transformation fait que le point de coordonnées au nord de l'équateur se retrouve au sud de l'équateur, soustrayez un million de la coordonnée d'ordonnée d'entrée. Une fois la transformation calculée, ajoutez 1 million à la coordonnée d'ordonnée et soustrayez 10 millions si la somme est supérieure à 10 millions. Le point est maintenant dans l'hémisphère nord.

c. Données requises. Voir le tableau 11-9 pour les données requises pour utiliser le programme 14.

ré. Instructions pour le programme 14--Références répertoriées et formulaire DA 7287-R. Les instructions d'utilisation du programme 14 pour calculer le formulaire DA 7287-R (Calcul de la transformation de coordonnées Datum-to-Datum (Program 14 - Listed Datums) (BUCS)) (Figures 11-7 à 11-11) sont décrites dans le Tableau 11- dix.

11-12. PROGRAMME 15--GAUSS-KRUGER

une. Objectif. Ce programme transforme les coordonnées entre les références GK et UTM dont les paramètres de transformation ellipsoïde et abrégé de Molodensky se trouvent dans le module BUCS-DDCT. Une liste de ces références figure à l'annexe E. Le programme dispose de quatre fonctions distinctes disponibles pour les transformations de coordonnées.

(1) Krassovski à UTM. Cette fonction convertit les coordonnées GK dans un système de référence basé sur Krassovsky en coordonnées UTM.

(2) UTM vers Krassovski. Cette fonction convertit les coordonnées UTM en coordonnées GK dans un système de référence basé sur Krassovsky.

(3) Bessel à UTM. Cette fonction convertit les coordonnées GK dans un système de référence basé sur Bessel en coordonnées UTM.

(4) UTM vers Bessel. Cette fonction convertit les coordonnées UTM en coordonnées GK dans un système de référence basé sur Bessel.

b. Capacités. Ce programme convertit les coordonnées des datums GK en datums UTM et vice versa. Les coordonnées d'entrée et de sortie peuvent être n'importe laquelle des combinaisons ci-dessus de références de coordonnées GK et UTM.

Le programme 15 calcule incorrectement les transformations de coordonnées Gauss-Kruger dans l'hémisphère sud. Ajoutez 4 000 mètres à la coordonnée nord lors du calcul des transformations UTM vers GK et GK vers UTM. Le seul pays utilisant la projection GK qui est touché est la Somalie.

c. Données requises. Voir le tableau 11-11 pour les données requises pour utiliser le programme 15.

ré. Instructions pour le programme 15 et le formulaire DA 7288-R. Les instructions pour le calcul du formulaire DA 7288-R (Calcul de la transformation de coordonnées Datum-to-Datum (Program 15--Gauss-Kruger [GK]) (BUCS)) (Figures 11-12 à 11-15) sont décrites dans le Tableau 11 -12. (Une copie reproductible de ce formulaire est incluse dans la section Formulaires vierges de ce livre.)

11-13. PROGRAMME 16--DÉFINI PAR L'UTILISATEUR

une. Objectif. Ce programme transforme les coordonnées entre deux références définies par l'utilisateur, entre une référence définie par l'utilisateur et une référence répertoriée, et entre une référence répertoriée et une référence définie par l'utilisateur. Voir l'annexe E pour les références répertoriées. Le DMA ou un QG supérieur diffusera les paramètres de référence définis par l'utilisateur si nécessaire.

(1) Référence définie par l'utilisateur à référence définie par l'utilisateur. Cette fonction convertit les coordonnées UTM d'un datum défini par l'utilisateur en coordonnées UTM dans le même ou un autre datum défini par l'utilisateur.

(2) Référence définie par l'utilisateur à référence répertoriée. Cette fonction convertit les coordonnées UTM d'un datum défini par l'utilisateur en coordonnées UTM dans un datum répertorié.

(3) Référence répertoriée en référence définie par l'utilisateur. Cette fonction convertit les coordonnées UTM d'un datum répertorié en coordonnées UTM dans un datum défini par l'utilisateur.

Noter. Les paramètres de transformation de la donnée Adindan (code 1) ont été mis à jour par la Defense Mapping Agency. Pour effectuer des transformations impliquant la référence Adindan, utilisez les paramètres mis à jour dans le Tableau 11-13.

Noter. Les décalages de référence effectués entre WGS-72 et toute autre référence à l'aide du programme 14 ou 15 ne produisent pas les mêmes solutions que le programme 16 lors de la saisie des paramètres de transformation en tant que références définies par l'utilisateur. Cela se produit parce qu'une formule plus précise pour WGS-72 est utilisée dans les programmes 14 et 15. Le programme 16 s'appuie sur les formules abrégées de Molodensky et les décalages de référence pour WGS-72. Utilisez le programme 16 pour les transformations de données uniquement lorsqu'elles sont dirigées par un QG supérieur et lorsque des paramètres spécifiques sont fournis. Les transformations de zone à zone peuvent être incluses dans l'une quelconque des transformations ci-dessus. Dans de tels cas, les zones de grille doivent être adjacentes. Voir le tableau des zones de quadrillage à l'annexe F.

b. Capacités. Ce programme convertit les coordonnées entre des références définies par l'utilisateur ou entre une référence définie par l'utilisateur et une référence répertoriée et calcule également les transformations de zone à zone.

Noter. Le demi-grand axe et le demi-petit axe ne peuvent pas être tous les deux 6360000.000, ou le programme peut s'arrêter.

c. Données requises. Voir le tableau 11-14 pour les données requises pour utiliser le programme 16.

ré. Instructions pour le programme 16 et le formulaire DA 7289-R. Les instructions d'utilisation du programme 16 pour calculer le formulaire DA 7289-R (Calcul de la transformation de coordonnées Datum-to-Datum (Program 16 - User-Defined) (BUCS) (Figures 11-16 à 11-18) sont décrites dans le Tableau 11 -15. (Une copie reproductible de ce formulaire est incluse dans la section Formulaires vierges de ce livre.)


Contenu

L'invention d'un système de coordonnées géographiques est généralement attribuée à Eratosthène de Cyrène, qui a composé son désormais perdu La géographie à la Bibliothèque d'Alexandrie au IIIe siècle av. [3] Un siècle plus tard, Hipparchus de Nicée a amélioré ce système en déterminant la latitude à partir des mesures stellaires plutôt que l'altitude solaire et en déterminant la longitude par le minutage des éclipses lunaires, plutôt que par l'estime. Au 1er ou au 2e siècle, Marinus de Tyr a compilé un vaste répertoire géographique et une carte du monde mathématiquement tracée en utilisant des coordonnées mesurées à l'est à partir d'un méridien principal sur la terre connue la plus occidentale, désignée les îles Fortunées, au large des côtes de l'Afrique occidentale autour des Canaries ou du Cap. Iles Verde, et mesurée au nord ou au sud de l'île de Rhodes au large de l'Asie Mineure. Ptolémée lui a attribué l'adoption complète de la longitude et de la latitude, plutôt que de mesurer la latitude en fonction de la longueur du jour du milieu de l'été. [4]

IIe siècle de Ptolémée La géographie utilisé le même méridien principal mais mesuré la latitude à partir de l'équateur à la place. Après que leur travail ait été traduit en arabe au IXe siècle, l'ouvrage d'Al-Khwārizmī Livre de la Description de la Terre a corrigé les erreurs de Marinus et de Ptolémée concernant la longueur de la mer Méditerranée, [note 1] faisant en sorte que la cartographie arabe médiévale utilise un méridien principal autour de 10° à l'est de la ligne de Ptolémée. La cartographie mathématique reprit en Europe suite à la récupération par Maximus Planudes du texte de Ptolémée un peu avant 1300 le texte fut traduit en latin à Florence par Jacobus Angelus vers 1407.

En 1884, les États-Unis ont accueilli la Conférence internationale des méridiens, à laquelle ont participé des représentants de vingt-cinq nations. Vingt-deux d'entre eux ont accepté d'adopter la longitude de l'Observatoire royal de Greenwich, en Angleterre, comme ligne de référence zéro. La République dominicaine a voté contre la motion, tandis que la France et le Brésil se sont abstenus. [5] La France a adopté le temps moyen de Greenwich à la place des déterminations locales de l'Observatoire de Paris en 1911.

Afin d'être sans ambiguïté sur la direction de la surface "verticale" et "horizontale" au-dessus de laquelle ils mesurent, les cartographes choisissent un ellipsoïde de référence avec une origine et une orientation données qui correspondent le mieux à leur besoin de la zone à cartographier. Ils choisissent ensuite la cartographie la plus appropriée du système de coordonnées sphériques sur cet ellipsoïde, appelé système de référence terrestre ou système géodésique.

Les données peuvent être globales, ce qui signifie qu'elles représentent la Terre entière, ou elles peuvent être locales, ce qui signifie qu'elles représentent un ellipsoïde qui correspond le mieux à une partie seulement de la Terre. Les points à la surface de la Terre se déplacent les uns par rapport aux autres en raison du mouvement des plaques continentales, de l'affaissement et du mouvement diurne des marées terrestres causés par la Lune et le Soleil. Ce mouvement quotidien peut atteindre un mètre. Les déplacements continentaux peuvent atteindre 10 cm par an, ou 10 m par siècle. Une zone anticyclonique du système météorologique peut provoquer un enfoncement de 5 mm . La Scandinavie s'élève de 1 cm par an en raison de la fonte des calottes glaciaires de la dernière période glaciaire, mais l'Écosse voisine ne s'élève que de 0,2 cm . Ces changements sont insignifiants si une donnée locale est utilisée, mais sont statistiquement significatifs si une donnée globale est utilisée. [1]

Des exemples de systèmes géodésiques mondiaux comprennent le système géodésique mondial (WGS 84, également connu sous le nom d'EPSG:4326 [6] ), le système de référence par défaut utilisé pour le système de positionnement global, [note 2] et le système international de référence terrestre et le cadre (ITRF), utilisé pour estimer la dérive des continents et la déformation de la croûte. [7] La ​​distance au centre de la Terre peut être utilisée à la fois pour des positions très profondes et pour des positions dans l'espace. [1]

Les datums locaux choisis par une organisation cartographique nationale comprennent le datum nord-américain, l'ED50 européen et l'OSGB36 britannique. Étant donné un emplacement, la donnée fournit la latitude ϕ et la longitude λ . Au Royaume-Uni, trois systèmes communs de latitude, longitude et hauteur sont utilisés. Le WGS 84 diffère à Greenwich de celui utilisé sur les cartes publiées OSGB36 d'environ 112 m. Le système militaire ED50, utilisé par l'OTAN, diffère d'environ 120 m à 180 m. [1]

La latitude et la longitude sur une carte établie par rapport à un système géodésique local peuvent ne pas être les mêmes que celles obtenues à partir d'un récepteur GPS. La conversion de coordonnées d'une référence à une autre nécessite une transformation de référence telle qu'une transformation de Helmert, bien que dans certaines situations une simple traduction puisse être suffisante. [8]

Dans les logiciels SIG courants, les données projetées en latitude/longitude sont souvent représentées sous la forme d'un Système de coordonnées géographiques. Par exemple, les données en latitude/longitude si le datum est le datum nord-américain de 1983 sont désignées par 'GCS North American 1983'.


Lidar : Coordonnées et formats

Il y a eu récemment une question sur les avoirs en lidar de Digital Coast via Twitter. Les questions impliquaient pourquoi notre lidar avait des coordonnées horizontales de latitude et de longitude et pourquoi nous avions le format LAS 1.2 au lieu du format LAS 1.4. Les réponses sont trop longues pour Twitter, je vais donc essayer d'y répondre ici. Ou arrêtez simplement avec la réponse TLDR.

Réponse TLDR : nous gardons les choses dans leur format natif et tout le monde ne veut pas la même projection.

Commençons par le format. Essentiellement, nous ne les modifions pas par rapport à leur format lorsque nous les avons reçus. S'ils sont venus en tant que LAS 1.2, ils le restent. Nous ne les convertissons pas en LAS 1.4. Vous trouverez également des données LAS 1.1 et peut-être même 1.0. Pourquoi ne convertissons-nous pas ? Essentiellement, il n'y a pas de gain évident à le faire. Cela ne veut pas dire qu'il n'y a aucun avantage à utiliser LAS 1.4 lorsque vous travaillez avec les données, mais le simple fait de changer le format n'ajoute aucune information. Si vos processus doivent avoir des avantages avec LAS 1.4, utilisez certainement quelque chose comme las2las (partie open source de LAStools) pour convertir.

Il y avait une déclaration selon laquelle l'USGS servait toutes ses données dans LAS 1.4, mais je ne pense pas que ce soit vrai. Il est vrai que les données conformes à la spécification 3DEP seront en LAS 1.4, mais il y a beaucoup de données plus anciennes sur le site ftp rocheux qui est dans son format d'origine. Je n'ai pas pensé à une bonne raison pour justifier la conversion en gros de milliers de milliards de points de données. S'il y en a un, j'aimerais en entendre parler.

J'ai découvert que la question initiale tournait autour des données LAS 1.4 qui avaient été compressées en LAZ avant qu'un compresseur natif n'existe. Ces fichiers semblent être compressés en LAS 1.2, mais deviennent LAS 1.4 lorsque vous les décompressez avec laszip. Ils ont depuis été corrigés pour utiliser le compresseur natif.

Et ces systèmes de coordonnées ? La plainte qui est passée par mon bureau était que les coordonnées géographiques n'étaient pas utiles et pourquoi n'étaient-elles pas toutes en UTM comme l'USGS. Donc, j'admets que les coordonnées géographiques ont leurs inconvénients. Comme tout le reste. Les plus gros inconvénients que je connais sont : 1) vous n'obtiendrez qu'une précision d'environ 1 cm horizontalement et 2) les coordonnées horizontales et verticales ne sont pas les mêmes, donc certaines analyses 3D qui supposent les mêmes unités dans toutes les directions ne fonctionnera pas. Le premier n'est pas un gros problème car la précision horizontale du lidar aéroporté n'est nulle part proche d'un centimètre. La taille de la tache seule est beaucoup plus grande. Le second est un problème légitime.

Pourquoi ne pas tout faire en UTM ? Cela semble probablement très bien si vous utilisez UTM tout le temps. Beaucoup d'universitaires entrent dans cette catégorie. Cependant, beaucoup d'ingénieurs ne le font pas. Ils appartiennent à la catégorie des avions d'État et peuvent même être tenus d'utiliser l'avion d'État. Il n'y a pas une seule projection pour les gouverner tous. Tout ce que je choisirai rendra quelqu'un moins heureux. Pour les systèmes projetés, il y a aussi le problème des ensembles de données qui traversent les frontières. Cela ne se produit pas pour beaucoup de gens car leur domaine d'intérêt est assez limité, mais il existe des ensembles de données qui couvrent plusieurs états et plusieurs zones UTM. Cela rend les choses moche et vous devez faire beaucoup plus attention lors de la fusion de données. J'ai définitivement obtenu des données qui se trouvaient dans deux zones UTM différentes, puis converties comme si elles se trouvaient dans une zone, créant une séparation à six degrés des données.

Les coordonnées géographiques sont un système qui nous permet d'utiliser le même système partout. C'était plus important au début lorsque l'obtention de fichiers LAS correctement codés était un peu incertain. Tant que nous savions dans quoi se trouvaient les données, nous pouvions convertir en données géographiques et avoir confiance en ce que nous traitions au lieu de compter sur un logiciel médiocre pour le comprendre. Le logiciel est bien meilleur maintenant et nous pourrions probablement changer, mais cela apporterait toujours un changement de flux de traitement pour notre système personnalisé (Data Access Viewer DAV). À l'heure actuelle, nous pouvons fusionner les données nécessaires à une demande dans un seul fichier, puis reprojeter sur la projection demandée. Si les projections peuvent varier, nous devrions d'abord reprojeter chaque fichier, puis fusionner. Ce n'est pas grave, mais un peu plus compliqué avec plus de choses à suivre.

Depuis que j'ai évoqué le système DAV, je pense qu'il vaut la peine de parler des différentes intentions de la distribution des données que nous avons. À l'origine, nous n'avions que le système personnalisé et un utilisateur n'a jamais vu dans quelle projection se trouvaient les données de base. Il se trouve également qu'il s'agit de hauteurs ellipsoïdes. Le système DAV effectue de nombreux traitements pour vous, y compris la reprojection. Cependant, les limites de la taille d'un travail que nous pouvions traiter en toute sécurité sans causer de problèmes signifiaient que vous ne pouviez pas l'utiliser pour obtenir un comté de données. Nous avons commencé à mettre les données sur un site ftp pour les personnes qui avaient besoin d'un comté complet de données (après avoir appliqué le dernier modèle de géoïde). La théorie était que si vous aviez besoin d'un comté complet de données, vous aviez probablement les connaissances et la puissance nécessaires pour gérer les projections changeantes et bien d'autres choses. Il existe un tableau pratique de ces ensembles de données. Si nous n'avions pas tout transformé en coordonnées géographiques, nous aurions probablement besoin de faire la même chose que l'USGS et de l'afficher dans la projection reçue (pas tous UTM d'ailleurs).

Qu'en est-il du futur? Le plan de l'USGS et du 3DEP consiste à utiliser la projection à surface égale d'Albers comme norme. Il reste encore beaucoup de choses à régler et il est probable que nous verrons des données fournies dans plusieurs projections. Cependant, ils prennent essentiellement le même type de décision que nous avons pris avec la géographie et essaient de couvrir le terrain avec un seul système. Ils voulaient également une surface égale pour des tailles d'acquisition cohérentes. Il y a de bonnes raisons d'avoir des projections différentes et votre système préféré ne fera pas plaisir à tout le monde. Il est préférable de vous assurer que vous disposez des outils nécessaires pour passer d'un système à un autre. Soit dit en passant, toutes ces reprojections passent par la géographie pour passer d'un système à un autre. Vous ferez également des transformations géographiques pour effectuer des transformations de données.


Convertir les coordonnées géographiques en zone UTM 29N avec Javascript - Systèmes d'Information Géographique

Une collection de fonctions géodésiques qui résolvent une variété de problèmes de géodésie. Prend en charge une large gamme d'ellipsoïdes de référence communs et définis par l'utilisateur. La plupart des fonctions sont vectorisées. La version la plus récente est disponible sur <http://www.craymer.com/software/matlab/geodetic/>. Les fonctions incluent :
Conversions d'angles
deg2rad - Degrés en radians
dms2deg - Degrés,minutes,secondes en degrés
dms2rad - Degrés,minutes,secondes en radians
rad2deg - Radians en degrés
rad2dms - Radians en degrés, minutes, secondes
rad2sec - Radians en secondes
sec2rad - Secondes en radians
Conversions de coordonnées
ell2utm - Coordonnées ellipsoïdales (lat, long) à UTM (N, E)
ell2xyz - Coordonnées ellipsoïdales (lat,long) à cartésiennes (x,y,z)
sph2xyz - Coordonnées sphériques (az,va,dist) à cartésiennes (x,y,z)
xyz2sph - Coordonnées cartésiennes (x,y,z) à sphériques (az,va,dist)
xyz2ell - Coordonnées cartésiennes (x,y,z) à ellipsoïdales (lat,long,ht)
xyz2ell2 - xyz2ell avec formule de hauteur Bowring
xyz2ell3 - xyz2ell utilisant la version complète de Bowring
utm2ell - UTM (N, E) aux coordonnées ellipsoïdales (lat, long)
Coordinate Transformations
refell - Reference ellipsoid definition
ellradii - Various radii of curvature
ct2lg - Conventional terrestrial (ECEF) to local geodetic (NEU)
dg2lg - Differences in Geodetic (lat,lon) to local geodetic (NEU)
cct2clg - Conventional terrestrial to local geodetic cov. matrix
clg2cct - Local geodetic to conventional terrestrial cov. matrix
rotct2lg - Rotation matrix for conventional terrestrial to local geod.
rotlg2ct - Rotation matrix for local geod. to conventional terrestrial
lg2ct - Local geodetic (NEU) to conventional terrestrial (ECEF)
lg2dg - Local geodetic (NEU) to differences in geodetic (lat,lon)
direct - Direct geodetic problem (X1,Y1,Z1 + Az,VA,Dist to X2,Y2,Z2)
inverse - Inverse geodetic problem (X1,Y1,Z1 + X2,Y2,Z2 to Az,VA,Dist)
simil - Similarity transformation (translation,rotation,scale change)
Date Conversions
cal2jd - Calendar date to Julian date
dates - Converts between different date formats
doy2jd - Year and day of year to Julian date
gps2jd - GPS week & seconds of week to Julian date
jd2cal - Julian date to calenar date
jd2dow - Julian date to day of week
jd2doy - Julian date to year & day of year
jd2gps - Julian date to GPS week & seconds of week
jd2mjd - Julian date to Modified Julian date
jd2yr - Julian date to year & decimal year
mjd2jd - Modified Julian date to Julian date
yr2jd - Year & decimal year to Julian date

Error Ellipses
errell2 - Computes error ellipse semi-axes and azimuth
errell3 - Computes error ellipsoid semi-axes, azimuths, inclinations
plterrel - Plots error ellipse for covariance matrix

Divers
cart2euler- Converts Cartesian coordinate rotations to Euler pole rotation
euler2cart- Converts Euler pole rotation to Cartesian coordinate rotations
findfixed - Finds fixed station based on 3D covariance matrix
pltnet - Plots network of points with labels

DirInv - Simple partial GUI script for direct and inverse problems
DirProb - Example of direct problem
Dist3D - Example to compute incremental 3D distances between points.
InvProb - Example of inverse problem
PltNetEl - Example plot of network error ellipses
ToUTM - Example of conversion from latitude,longitude to UTM


Equations Used for Datum Transformations

The following formulae can be used to transform geographic coordinates between geodetic datums using three and seven parameter similarity transformations.

This conversion is a three-step process:

  • Convert geographic coordinates to their Cartesian equivalents
  • Apply similarity transformation to Cartesian coordinates
  • Convert Cartesian coordinates back to geographic values

Geographic coordinates to Cartesian coordinates

These formulae can be used to convert geographic coordinates, latitude ( Φ ), longitude ( λ ), and height ( h ), into Cartesian coordinates (X, Y, Z ):

Where a and f are obtained from the reference ellipsoid used for the respective geodetic datum and the h is the height of the computation point or approximated as zero

Three parameter transformation

The three parameter transformation is implemented using:

Seven parameter transformation

The Helmert seven parameter similarity transformation is implemented using:

Note: the rotation parameters ( R ) must be converted from arc-seconds to radians before being used in this equation.

Note: this is a simplified version of the Helmert formulae that applies for small rotation angles. This is the official formulae to use for the NZGD49-NZGD2000 seven parameter transformation.

Cartesian coordinates to geographic coordinates

These formulas can be used to convert Cartesian coordinates ( X, Y, Z ) into geographic coordinates latitude ( Φ ), longitude ( λ ), and height ( h ).

Where a and f are obtained from the reference ellipsoid used for the respective geodetic datum:

Note: because NZGD1949 is a horizontal datum the height resulting from the transformation will not be in terms of the output datum, the equation is shown here for completeness.


RNDOW

The rNDOW package is very early in development. I’ve added three functions that I use very frequently for exploratory analysis and visualization of animal movement data (post on that soon!). The package will be used by my colleagues to interface with our data management systems. Due to this, many functions will have hard-coded or default values that make data analysis easier for us. Most of the functions will encapsulate common data munging, exploration, and visualization procedures.

Les fonctions

There are 4 functions in the package, xyConv , moveParams , plotTraj and plot3DTraj . I’ve included a randomly sampled animal trajectory for examples muldDat . The packages can be found at the NDOW-ARG/rNDOW repository. The functions solve a common workflow I have in R, converting latlong coordinates to UTM, adding parameters for animal movement models, and plotting the animal trajectory for exploratory visualization. More information on the input and specifics of the functions can be found at the bottom of the post.

Here is an example of how I use the functions.

Line 1 calls the data from the package. Line 2 converts the coordinates to UTM Zone 11, the default options work for me. Line 3 adds the movement parameters to the dataframe. The timestamp isn’t class POSIXct so I convert it using the fasttime library. Line 4 and 5 plot the trajectory. Line 5 creates a “space time cube”, a visualization of the spatial and temporal distribution of the animals GPS locations.


Free Online Tools

Address Tools

BatchGeo
Generate a map from a list of addresses.

Noms géographiques
Look up coordinates for almost any place on Earth.

Texas A&M Geocoder
Free geocoding service for US addresses (max 2,500 addresses).

US Census Geocoder
Match addresses to geographic locations.

Coordinate System Tools

EPSG.io Coordinate Systems Worldwide
Look up information about over 6,000 coordinate systems.

FCC Degrees Minutes Seconds / Decimal Degrees Converter
Convert coordinates from DMS to/from decimal degrees.

State Plane Coordinate Systems
Look up the State Plane and UTM Zones for each county in the U.S.

Data Cleanup Tools

OpenRefine
Clean, transform, and explore large datasets.

Data Editing Tools

Mapshaper
Easy-to-use web interface for editing geospatial data (shapefile, GeoJSON, TopoJSON, CSV). Built-in console allows for command-line editing.

Georeferencing Tools

Map Warper
Search for georectified maps or upload and georeference your own.

David Rumsey Historical Map Collection Georeferencer
Georeference a historical map from the Rumsey collection.

Imagery Analysis Tools

Picterra
User-friendly platform for training an AI to detect objects in remotely sensed imagery.

Storytelling Tools

Leaflet Storymaps with Google Sheets
A do-it-yourself approach to storytelling using Github and Google Drive. Step-by-step instructions are available on the project’s Github page.

Mapbox Interactive Storytelling
For users with some knowledge of HTML and CSS, this is a configurable storytelling option that produces simple, chapter-based scrollable stories. Requires registering for a free Mapbox account. Generous free tier allows 50,000 map loads per month before charges begin accruing.

MapStory
Free story-telling app (alternative to ArcGIS StoryMaps) built by the nonprofit MapStory Foundation.

StoryMapJS
Free story-telling app (alternative to ArcGIS StoryMaps). One of several online tools by Northwestern University Knight Lab designed for telling compelling stories with multimedia, VR, numerical data, and timelines.

Web Mapping Tools

ColorBrewer 2.0
Get help picking out colors for your map.

Datawrapper
Easy-to-use interface for generating web-based maps, charts, and tables – no coding experience necessary. Generous free tier allows unlimited map creation.

Google Earth for Chrome
Only works in Chrome browser . Explore global satellite imagery from your browser.

Google Maps JavaScript API
Series of JavaScript tutorials for building a custom map in Google Maps. A demo API is provided for the tutorials, but an API must be purchased before sharing maps with the public.

Kartograph
Python and JavaScript libraries for designing interactive online vector maps.

Leaflet
JavaScript library for designing interactive mobile-friendly maps.

MaskMy.XYZ
Client-side JavaScript tool that masks sensitive location information.

Boîte à cartes
Suite of tools and libraries for designing eye-catching, interactive maps. Generous free tier allows 50,000 map loads per month before charges begin accruing.

Recogito
Annotate texts and images and link to digital gazetteers.


Voir la vidéo: Le geocoding et le geocoding inversé - Geocasts