¿Cómo crear una Tissot Indicatrix precisa?

31

Una Tissot Indicatrix es un método útil para comunicar de un vistazo los tipos de distorsión a los que es propensa una proyección dada (en la figura siguiente, cada uno de los círculos rojos ocupa la misma área). Me han dicho que los métodos populares para generar TI tienen sus propios problemas, hasta el punto de ser a veces lamentablemente inexactos.

¿Cuál es el problema con los métodos populares y cuál es la forma más correcta de generar un TI que sea accesible para su tipo promedio de SIG (ette)?

Mercator y globos con telas

wilkie mate
fuente
2
gran pregunta Me gustaría saber eso también.
George Silva

Respuestas:

23

Cualquier software que pueda proyectar coordenadas con precisión puede calcular indicadores Tissot precisos .

Una buena fuente para las fórmulas es Snyder, John, Map Projections - A Working Manual , principalmente en las páginas 20-26. (No los reproduciré aquí porque este sitio no tiene herramientas apropiadas para comunicar fórmulas matemáticas). Requieren las cuatro primeras derivadas de las coordenadas proyectadas (x, y) con respecto a las coordenadas esféricas (lat, lon) = (phi, lambda):

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

Todo lo demás sobre los TI se calcula en términos de estos (utilizando algunas funciones aritméticas y trigonométricas: el coseno, el seno inverso principal y la tangente inversa principal). Los cálculos requieren una descripción de la forma de la tierra. Para obtener la mayor precisión, utilice un dato elipsoidal con semieje mayor a y excentricidad e. (Estos serán conocidos por el software).

El libro de Snyder tiene instrucciones sobre cómo calcular todo excepto estos derivados. Hazlo numéricamente. He tenido excelentes resultados usando estimaciones de diferencias finitas centrales de primer orden a una distancia de h = 10 ^ (- 5.2) radianes (típicamente alrededor de 50 metros): este es un buen compromiso entre tratar de acercarse infinitamente y perder demasiada precisión de redondeo de punto flotante (suponiendo doble precisión), porque el error cometido es proporcional a (10 ^ (- 5.2)) ^ 2 = 10 ^ (- 10.4) y 10 ^ (- 5.2) es igual a 10 ^ 10.4 veces la precisión de doble precisión IEEE de 10 ^ (- 15.6) y todavía es mucho más grande que la precisión típica en proyecciones, que generalmente van de 10 ^ (- 10) a aproximadamente 10 ^ (- 14).

Entonces, ¿cómo se calculan las estimaciones de diferencias finitas? Esta parte es sorprendentemente fácil. Para obtener dx / d (phi) en un punto (phi, lambda), solicite a su SIG que proyecte los puntos

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

Usa las estimaciones

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

Del mismo modo, proyecte los puntos

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

y usar las estimaciones

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

Eso requiere cuatro proyecciones y un poco de aritmética. (Puede reducirlo a tres mediante el uso de diferencias no centrales, pero la precisión disminuye un poco. Es sabio apuntar a una alta precisión, sin dejar que h sea demasiado pequeña, a menos que esté seguro de que su SIG está usando grado de encuesta (milímetro) precisión en sus fórmulas de proyección.)

A partir de estos derivados, junto con las fórmulas de Snyder (prestando atención a las modificaciones descritas en 4-19 y 4-21), puede obtener las longitudes de los ejes de la Tissot Indicatrix en (phi, lambda) y su orientación. En los mapas a escala mundial, el TI será tan pequeño como invisible, por lo que lo último que debe hacer es decidir cuánto desea reescalar cada TI. Determino el factor de escala averiguando qué tan grande será el mapa, encontrando los tamaños de TI típicos en todo el mapa y escalando para que esas TI tengan aproximadamente un 6% del ancho del mapa. Es un buen comienzo, de todos modos; Dejo que el usuario ajuste el tamaño de la TI desde allí. Por supuesto, volverá a escalar todas las TI en la misma cantidad, para que puedan ser comparadas, y cada una será reescalada alrededor de su propio centro (que se obtiene mediante una quinta proyección, (phi, lambda) -> (x, y) )

Una buena adición a la representación elíptica de la TI es mostrar las direcciones del meridiano local y el paralelo: luego, de un vistazo, puede evaluar la convergencia de la cuadrícula . También muestro un círculo estándar (que no representa distorsión) concéntrico con cada TI porque mejora la capacidad del lector para medir la cantidad de distorsión representada por cada elipse.

texto alternativo

En esta proyección de Mollweide, cabe destacar la extrema TI cerca del polo sur. Sigue siendo una elipse perfecta y describe con precisión la distorsión del mapa allí.

whuber
fuente
2
Whuber es el mejor: P.
George Silva el
1
Me di cuenta de que ESRI publicó un artículo sobre la creación de una matriz de tejido con buffers, ¿es correcto el método del buffer ya que una matriz de indicadores y un buffer "no son lo mismo"? blogs.esri.com/Support/blogs/mappingcenter/archive/2011/03/28/…
SaultDon
1
@Sault Las memorias intermedias son un sustituto de las TI correctas. Funcionarán bastante bien en general, siempre que no se extiendan más de unos pocos grados, excepto donde las distorsiones se vuelven extremas, donde tampoco funcionarán en absoluto (porque el búfer se extiende a través de una región de distorsión infinita, como se muestra en este mapa ESRI ) o dan como resultado formas no elípticas. Un ligero cambio haría que este enfoque funcionara mucho mejor: calcule pequeños buffers, como los buffers de 50 m, y expanda cada uno uniformemente (en el mapa) alrededor de su centro para hacerlo visible.
whuber
1
¿Cómo calculó la orientación de las elipses?
Jason Davies
3
@Jason La matriz de Tissot no es otra cosa que una representación gráfica de la derivada de la proyección aplicada a un círculo en el plano tangente en un punto. La forma más simple de dibujarlo, entonces, es calcular esta derivada (es una matriz A de 2 por 2 ), parametrizar el círculo y aplicar la derivada a los vértices parametrizados. Voila, tienes una elipse. (Así es exactamente como se dibujaron las elipses en mi figura). Sus ejes se pueden calcular como los vectores propios de A'A ; Este es un cálculo simple que no implica nada más complicado que la fórmula cuadrática.
Whuber