Comment créer une Tissot Indicatrix précise?


31

Un Tissot Indicatrix est une méthode utile pour communiquer en un coup d'œil les types de distorsion auxquels une projection donnée est sujette (dans la figure ci-dessous, chacun des cercles rouges occupe la même zone). On m'a dit que les méthodes populaires pour générer des TI ont leurs propres problèmes, au point d'être parfois terriblement inexactes.

Quel est le problème avec les méthodes populaires et quelle est la manière la plus correcte de générer un TI accessible à votre type SIG moyen (ette)?

Mercator et globes avec Tissot


2
grande question. J'aimerais aussi le savoir.
George Silva

Réponses:


23

Tout logiciel capable de projeter des coordonnées avec précision peut calculer des indicatrices Tissot précises .

Une bonne source pour les formules est Snyder, John, Map Projections - A Working Manual , principalement aux pages 20-26. (Je ne les reproduirai pas ici car ce site ne dispose pas d'outils appropriés pour communiquer des formules mathématiques.) Ils nécessitent les quatre premières dérivées des coordonnées projetées (x, y) par rapport aux coordonnées sphériques (lat, lon) = (phi, lambda):

dx / d(phi), dx / d(lambda);
dy / d(phi), dy / d(lambda).

Tout le reste sur les TI est calculé en fonction de ceux-ci (en utilisant certaines fonctions arithmétiques et trigonométriques: le cosinus, le sinus inverse principal et la tangente inverse principale). Les calculs nécessitent une description de la forme de la terre. Pour la plus grande précision, utilisez une donnée ellipsoïdale d'axe semi-principal a et d'excentricité e. (Ceux-ci seront connus du logiciel.)

Le livre de Snyder contient des instructions sur la façon de tout calculer sauf ces dérivés. Faites-le numériquement. J'ai obtenu d'excellents résultats en utilisant des estimations de différences finies centrales de premier ordre à une distance de h = 10 ^ (- 5,2) radians (généralement autour de 50 mètres): c'est un bon compromis entre essayer de se rapprocher infinitésimalement et perdre trop de précision de arrondi en virgule flottante (en supposant une double précision), car l'erreur commise est proportionnelle à (10 ^ (- 5.2)) ^ 2 = 10 ^ (- 10.4) et 10 ^ (- 5.2) est égal à 10 ^ 10.4 fois la précision de la double précision IEEE de 10 ^ (- 15,6) et il est encore beaucoup plus grand que la précision typique dans les projections, qui vont généralement de 10 ^ (- 10) à environ 10 ^ (- 14).

Alors, comment calculez-vous les estimations des différences finies? Cette partie est étonnamment facile. Pour obtenir dx / d (phi) en un point (phi, lambda), demandez à votre SIG de projeter les points

(phi - h/2, lambda) --> (x0,y0),
(phi + h/2, lambda) --> (x1,y1).

Utilisez les estimations

dx / d(phi) = (x1 - x0)/h,
dy / d(phi) = (y1 - y0)/h.

De même, projetez les points

(phi, lambda - h/2) --> (x2,y2),
(phi, lambda + h/2) --> (x3,y3)

et utiliser les estimations

dx / d(lambda) = (x3 - x2)/h,
dy / d(lambda) = (y3 - y2)/h.

Cela prend quatre projections et un tout petit peu d'arithmétique. (Vous pouvez le réduire à trois en utilisant des différences non centrales, mais la précision diminue un peu. Il est sage de viser une grande précision, sans laisser h trop petite, à moins que vous ne soyez sûr que votre SIG utilise une qualité d'enquête (millimètre) précision dans ses formules de projection.)

A partir de ces dérivés, ainsi que des formules de Snyder (en faisant attention aux modifications décrites en 4-19 et 4-21), vous pouvez obtenir les longueurs des axes du Tissot Indicatrix en (phi, lambda) et son orientation. Sur les cartes à l'échelle mondiale, le TI sera si petit qu'il sera invisible, donc la dernière chose à faire est de décider combien vous voulez redimensionner chaque TI. Je détermine le facteur d'échelle en déterminant la taille de la carte, en trouvant les tailles des TI typiques sur la carte et en les mettant à l'échelle de sorte que ces TI soient environ 6% aussi larges que la carte. C'est un bon début, de toute façon; Je laisse l'utilisateur ajuster la taille de la TI à partir de là. Bien sûr, vous redimensionnerez tous les TI du même montant, afin qu'ils puissent être comparés, et chacun sera redimensionné autour de son propre centre (qui est obtenu par une cinquième projection, (phi, lambda) -> (x, y) ).

Un bon ajout à la représentation elliptique du TI est de montrer les directions du méridien local et parallèle: puis, en un coup d'œil, vous pouvez évaluer la convergence de la grille . Je montre également un cercle standard (ne représentant aucune distorsion) concentrique avec chaque TI, car il améliore la capacité du lecteur à évaluer la quantité de distorsion représentée par chaque ellipse.

texte alternatif

Dans cette projection de Mollweide, il convient de noter l'extrême TI près du pôle sud. C'est toujours une ellipse parfaite et décrit avec précision la distorsion de la carte là-bas.


2
whuber est le meilleur: P.
George Silva

1
J'ai remarqué qu'ESRI a publié un article sur la création d'une indicatrice tissot avec des tampons, la méthode du tampon est-elle correcte car une indicatrice et un tampon «ne sont pas du tout les mêmes»? blogs.esri.com/Support/blogs/mappingcenter/archive/2011/03/28/…
SaultDon

1
@Sault Les tampons remplacent les TI corrects. Ils fonctionneront assez bien en général, à condition qu'ils ne s'étendent pas sur plus de quelques degrés, sauf lorsque les distorsions deviennent extrêmes, où elles ne fonctionnent pas du tout (car le tampon se répand sur une région de distorsion infinie, comme indiqué sur cette carte ESRI ) ou ils donnent des formes non elliptiques. Un léger changement rendrait cette approche beaucoup plus efficace: calculez de petits tampons, tels que des tampons de 50 m, et agrandissez-les uniformément (sur la carte) autour de son centre pour le rendre visible.
whuber

1
Comment avez-vous calculé l'orientation des ellipses?
Jason Davies

3
@Jason L'indicatrice Tissot n'est rien d'autre qu'une représentation picturale de la dérivée de la projection appliquée à un cercle dans le plan tangent en un point. Le moyen le plus simple de le dessiner est alors de calculer cette dérivée (c'est une matrice A 2 x 2 ), de paramétrer le cercle et d'appliquer la dérivée aux sommets paramétrés. Voila, vous avez une ellipse. (C'est exactement ainsi que les ellipses de ma figure ont été dessinées.) Ses axes peuvent être calculés comme les vecteurs propres de A'A ; il s'agit d'un calcul simple n'impliquant rien de plus compliqué que la formule quadratique.
whuber
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.