Jede Software, die Koordinaten genau projizieren kann, kann genaue Tissot-Indikatoren berechnen .
Eine gute Quelle für die Formeln ist Snyder, John, Map Projections - A Working Manual , hauptsächlich auf den Seiten 20-26. (Ich werde sie hier nicht reproduzieren, da diese Site keine geeigneten Werkzeuge zur Übermittlung mathematischer Formeln enthält.) Sie erfordern alle vier ersten Ableitungen der projizierten Koordinaten (x, y) in Bezug auf die sphärischen Koordinaten (lat, lon) = (Phi, Lambda):
dx / d(phi), dx / d(lambda);
dy / d(phi), dy / d(lambda).
Alles andere an den TIs wird in diesen Begriffen berechnet (unter Verwendung einiger arithmetischer und trigonometrischer Funktionen: dem Cosinus, dem inversen Hauptsinus und dem inversen Haupttangens). Die Berechnungen erfordern eine Beschreibung der Erdform. Verwenden Sie für die höchste Genauigkeit ein Ellipsoid-Datum mit der Halbschwerachse a und der Exzentrizität e. (Diese werden der Software bekannt sein.)
Snyders Buch enthält Anweisungen, wie man alles außer diesen Derivaten berechnet. Mach es numerisch. Ich habe hervorragende Ergebnisse erzielt, wenn ich Schätzungen für zentrale finite Differenzen erster Ordnung in einer Entfernung von h = 10 ^ (- 5,2) Radiant (typischerweise um die 50 Meter) verwendet habe: Dies ist ein guter Kompromiss zwischen dem Versuch, unendlich nahe zu kommen und dem Verlust von zu viel Präzision Gleitkomma-Rundung (unter der Annahme einer doppelten Genauigkeit), da der gemachte Fehler proportional zu (10 ^ (- 5.2)) ist. ^ 2 = 10 ^ (- 10.4) und 10 ^ (- 5.2) entspricht dem 10 ^ 10.4-fachen der Genauigkeit nach IEEE mit doppelter Genauigkeit von 10 ^ (- 15,6) und es ist immer noch viel größer als die typische Präzision in Projektionen, die normalerweise von 10 ^ (- 10) bis etwa 10 ^ (- 14) laufen.
Wie berechnen Sie also die Schätzungen für endliche Differenzen? Dieser Teil ist überraschend einfach. Bitten Sie Ihr GIS, die Punkte zu projizieren, um dx / d (phi) an einem Punkt (phi, lambda) zu erhalten
(phi - h/2, lambda) --> (x0,y0),
(phi + h/2, lambda) --> (x1,y1).
Verwenden Sie die Schätzungen
dx / d(phi) = (x1 - x0)/h,
dy / d(phi) = (y1 - y0)/h.
In ähnlicher Weise projizieren Sie die Punkte
(phi, lambda - h/2) --> (x2,y2),
(phi, lambda + h/2) --> (x3,y3)
und verwenden Sie die Schätzungen
dx / d(lambda) = (x3 - x2)/h,
dy / d(lambda) = (y3 - y2)/h.
Das braucht vier Projektionen und ein bisschen Arithmetik. (Sie können es auf drei reduzieren, indem Sie nicht-zentrale Unterschiede verwenden, aber die Genauigkeit nimmt ein wenig ab. Es ist ratsam, eine hohe Genauigkeit anzustreben, ohne dass sie zu klein wird, es sei denn, Sie sind sicher, dass Ihr GIS Vermessungsgrad (Millimeter) verwendet.) Genauigkeit in seinen Projektionsformeln.)
Aus diesen Ableitungen können Sie zusammen mit Snyders Formeln (unter Berücksichtigung der unter 4-19 und 4-21 beschriebenen Modifikationen) die Längen der Achsen der Tissot-Indikatrix bei (phi, lambda) und ihre Ausrichtung erhalten. Auf Karten im Weltmaßstab ist der TI so klein, dass er unsichtbar ist. Als letztes müssen Sie entscheiden, wie viel Sie für jeden TI neu skalieren möchten. Ich bestimme den Skalierungsfaktor, indem ich herausfinde, wie groß die Karte sein wird, die Größen typischer TIs auf der Karte ermittle und skaliere, sodass diese TIs ungefähr 6% so breit sind wie die Karte. Es ist jedenfalls ein guter Anfang; Ich lasse den Benutzer die Größe des TI von dort aus anpassen. Natürlich werden Sie alle TIs um den gleichen Betrag neu skalieren, damit sie verglichen werden können, und jeder wird um sein eigenes Zentrum neu skaliert (was durch eine fünfte Projektion (phi, lambda) -> (x, y) erhalten wird. ).
Eine schöne Ergänzung zur elliptischen Darstellung des TI ist die Darstellung der Richtungen des lokalen Meridians und der Parallelen: Auf einen Blick können Sie dann die Gitterkonvergenz beurteilen . Ich zeige auch einen Standardkreis (der keine Verzerrung darstellt), der konzentrisch zu jedem TI ist, da dies die Fähigkeit des Lesers verbessert, das Ausmaß der von jeder Ellipse dargestellten Verzerrung zu messen.
Bemerkenswert bei dieser Mollweide-Projektion ist der extreme TI in der Nähe des Südpols. Es ist immer noch eine perfekte Ellipse und beschreibt genau die Kartenverzerrung dort.
