¿Calculando la distorsión del área fuera de la zona UTM?

26

Uno de mis colegas está trabajando con datos que se distribuyen en dos zonas UTM. La mayoría de los datos se encuentran en una zona, con algunos valores atípicos en otra zona. Le gustaría saber cuál sería la distorsión del área de esos valores extremos si estuvieran en la zona UTM principal.

¿Existe una fórmula para calcular la distorsión del área sabiendo qué tan lejos estaban las características en la otra zona UTM?

    
pregunta kenbuja 17.08.2012 - 21:50

2 respuestas

30

UTM utiliza una proyección de Mercator transversal con un factor de escala de 0.9996 en el meridiano central. En el Mercator, el factor de escala distancia es la secante de la latitud (una fuente: enlace ), por lo que el factor de escala del área es el cuadrado de este factor de escala (porque se aplica en todas las direcciones, siendo el Mercator conforme). Entendiendo la latitud como la distancia esférica al ecuador , y aproximando el elipsoide con una esfera, podemos aplicar esta fórmula a cualquier aspecto de la proyección de Mercator. Así:

  

El factor de escala es 0.9996 veces la secante de la distancia (angular) al meridiano central. El factor de escala del área es el cuadrado de esta cantidad.

Para encontrar esta distancia, considere el triángulo esférico que se forma al viajar a lo largo de una geodésica desde un punto arbitrario en (lon, lat) = (lambda, phi) directamente hacia el meridiano central en la longitud mu, a lo largo de ese meridiano hasta el polo más cercano , y luego de vuelta a lo largo del meridiano lambda hasta el punto original. El primer giro es un ángulo recto y el segundo es un ángulo de lambda-mu. La cantidad recorrida a lo largo de la última porción es de 90 phi. La Spherical Law of Sines se aplica a los estados de este triángulo

  

sin (lambda-mu) / sin (distancia) = sin (90 grados) / sin (90-phi)

con solución

  

distancia = ArcSin (sin (lambda-mu) * cos (phi)).

Esta distancia se da como un ángulo, lo cual es conveniente para calcular la secante.

Ejemplo

Considere la zona UTM 17, con el meridiano central en -183 + 17 * 6 = -81 grados. Deje que la ubicación periférica sea de longitud -90 grados, latitud 50 grados. Entonces

Paso 1: la distancia esférica desde (-90, 50) al meridiano de -81 grados es igual a ArcSin (sin (9 grados) * cos (50 grados)) = 0.1007244 radianes.

Paso 2: la distorsión del área es igual a (0.9996 * s (0.1007244 radianes)) ^ 2 = 1.009406.

(Los cálculos numéricos con el elipsoide GRS 80 dan el valor como 1.009435, mostrando que la respuesta que calculamos es 0.3% demasiado baja: es el mismo orden de magnitud que el aplanamiento del elipsoide, lo que indica que el error se debe a la esférica aproximación.)

Aproximaciones

Para tener una idea de cómo cambia el área, podemos usar algunas identidades trigonométricas para simplificar la expresión general y expandirla como una serie de Taylor en lambda-mu (el desplazamiento entre la longitud del punto y la longitud del meridiano central UTM ). Funciona para

  

Factor de escala de área ~ 0.9992 * (1 + cos (phi) ^ 2 * (lambda-mu) ^ 2).

Como con todas estas expansiones, el ángulo lambda-mu debe medirse en radianes. El error es menor que 0.9992 * cos (phi) ^ 4 * (lambda-mu) ^ 4, que está cerca del cuadrado de la diferencia entre la aproximación y 1, es decir, el cuadrado del valor después del punto decimal .

En el ejemplo con phi = 50 grados (con un coseno de 0.642788) y lambda-mu = -9 grados = -0.15708 radianes, la aproximación da 0.9992 * (1 + 0.642788 ^ 2 * (-0.15708) ^ 2) = 1.009387. Mirando más allá del punto decimal y la cuadratura, deducimos (incluso sin saber el valor correcto) que su error no puede ser mayor que (0.009387) ^ 2 = menor que 0.0001 (y de hecho el error es solo una quinta parte de ese tamaño).

A partir de este análisis, es evidente que en latitudes altas (donde cos (phi) es pequeño), los errores de escala siempre serán pequeños; y en latitudes más bajas, los errores de escala de área se comportarán como el cuadrado de la diferencia en longitudes.

    
respondido por el whuber 18.08.2012 - 00:17
4

La herramienta GeographicLib de GeoConvert

enlace

permite una superposición generosa entre las zonas UTM (específicamente, la conversión a se permite una zona vecina siempre que el este resultante sea en el rango [0km, 1000km]). GeoConvert también puede reportar el meridiano. convergencia y escala y, como señala Whuber, la distorsión del área es la cuadrado de la escala.

Por ejemplo, su zona "principal" es 42 y se le da un punto dado

41N 755778 3503488

(Universidad de Kandahar) que está a unos 29 km al oeste de la zona 42. Para convertir esto a la zona 42, use

echo 41N 755778 3503488 | GeoConvert -u -z 42   == > 42N 186710 3505069

Para determinar la convergencia y la escala de los meridianos en la zona 42, agregue -c bandera

echo 41N 755778 3503488 | GeoConvert -u -z 42 -c   == > -1.73405 1.0008107

Por lo tanto, la distorsión del área es 1.0008107 ^ 2 = 1.0016221.

    
respondido por el cffk 01.09.2012 - 21:41

Lea otras preguntas en las etiquetas