¿Cuál es el error aproximado del teorema de Pitágoras versus la fórmula de Haversine al medir distancias en la esfera a varias escalas?

30

Muchas personas cuando intentan calcular por primera vez las distancias entre dos pares de longitud / latitud preguntan si el teorema de Pitágoras funciona como una función de distancia apropiada.

La mayoría de las personas responden "no, el teorema de Pitágoras solo funciona en un plano euclidiano 2D". En raras ocasiones, sin embargo, las personas mencionan el efecto de la escala y la ubicación en la esfera sobre cuán inexacto es el teorema de Pitágoras.

La idea básica es a escalas muy pequeñas, la superficie de una esfera se parece mucho a un plano. A escalas muy grandes, las distancias a lo largo de la superficie son más curvas y, por lo tanto, la diferencia entre el teorema de Pitágoras incorrecto y la fórmula de Haversine correcta es mayor.

¿Alguien conoce una fórmula o regla general que le indique la diferencia entre las dos medidas de distancia en función de la escala de la distancia que está tratando de medir?

Creo que tener esto explícitamente ayudaría en:

  1. explicando por qué el teorema de Pitágoras no es perfecto; y
  2. para que las personas que buscan distancias más "aproximadas" sepan cuándo Pitágoras realmente cumplirá sus propósitos.
Amos Budde
fuente

Respuestas:

34

Usar la fórmula pitagórica en posiciones dadas en latitud y longitud tiene tan poco sentido como, por ejemplo, calcular el área de un círculo usando la fórmula para un cuadrado: aunque produce un número, no hay razón para suponer que debería funcionar.

Aunque a escalas pequeñas cualquier superficie lisa parece un plano, la precisión de la fórmula pitagórica depende de las coordenadas utilizadas. Cuando esas coordenadas son latitud y longitud en una esfera (o elipsoide), podemos esperar que

  1. Las distancias a lo largo de las líneas de longitud serán razonablemente precisas.

  2. Las distancias a lo largo del ecuador serán razonablemente precisas.

  3. Todas las demás distancias serán erróneas, en proporción aproximada a las diferencias en latitud y longitud.

El error depende del punto inicial y final de los cálculos de distancia. Sin embargo, debido a que tanto la esfera como el elipsoide tienen una simetría circular alrededor del eje, el error depende solo de la diferencia de las longitudes, por lo que para estudiar este error también podríamos tomar el punto de origen en el Primer Meridiano. Debido a que tanto la esfera como el elipsoide son simétricos bajo una reflexión norte-sur, solo necesitamos estudiar los puntos de origen en el hemisferio sur. Para cualquier punto de este tipo, podemos dibujar un mapa de contorno del error relativo, igual a [cálculo de Pitágoras] / [Distancia verdadera].

La fórmula pitagórica, que usa el radio medio de la tierra, es

Pythagorean distance =  6371000. * Sqrt[dx^2 + dy^2]] * pi / 180 meters

donde dx es la diferencia en longitudes y dy es la diferencia en latitudes, ambas en grados. (La diferencia en los valores de longitud se reduce en el módulo 360 para dar el valor correcto de dx cuando se cruza el antimeridiano; no hacerlo introduciría errores artificialmente grandes que no nos dicen nada sobre la fórmula pitagórica en sí).

Los siguientes gráficos muestran el error relativo en comparación con la distancia correcta en el elipsoide WGS 84 para latitudes de -70 a 0 en incrementos de 10 grados. La coordenada horizontal es la diferencia en longitudes y la coordenada vertical es la latitud del destino. Las regiones claras tienen un error relativamente pequeño: las líneas de contorno están en 1, 1.01, 1.02, 1.05, 1.1, 1.2, 1.5, 2, etc. (Las áreas blancas puras en las esquinas son lugares donde el error va más allá del rango de estos contornos .) Los puntos rojos muestran el punto de origen.

Parcelas

Las bandas blancas verticales atestiguan la exactitud de la expectativa (1): las distancias pitagóricas son precisas cuando hay una pequeña diferencia en las longitudes. Las bandas blancas horizontales en latitudes bajas confirman la expectativa (2): cerca del ecuador, las distancias horizontales son razonablemente precisas. De lo contrario, como lo atestiguan las extensas regiones más oscuras, en todas las demás distancias la fórmula pitagórica es mala.


Podemos hacer estimaciones cuantitativas del máximoerror alcanzado para pares de puntos cercanos (dentro, por ejemplo, a unos cientos de kilómetros entre sí). La escala, que usa un valor apropiado para el radio, es verdadera a lo largo del meridiano pero a lo largo de un círculo de latitud yerra aproximadamente por la secante de la latitud. Por ejemplo, a una latitud de 40 grados, la secante es 1.31, lo que implica que la fórmula de Pitágoras dará distancias aproximadamente 31% demasiado grandes en la dirección este-oeste. (Esto es evidente en el diagrama de contorno superior derecho, para un punto de origen a -40 grados de latitud, donde la región inmediatamente al este-oeste del punto rojo se encuentra entre los contornos 1.2 y 1.5). Las distancias cortas en todas las demás direcciones serán demasiado grande en alguna cantidad entre 0% y 31%; distancias más largas pueden errar aún más (como muestran los gráficos de contorno).

whuber
fuente
1
Realmente debe haber una funcionalidad de 'respuesta favorita' para respuestas como estas.
Devdatta Tengshe
2
@DevdattaTengshe: explícitamente requiere que esté cuerdo: "donde dx es la diferencia en longitudes (expresada entre -180 y 180) y dy es la diferencia en latitudes, tanto en grados".
lynxlynxlynx
1
¿Eso es 2, ya que 2 * 179 es mayor que 180?
lynxlynxlynx
1
@whuber: Sé esto y tú lo sabes, pero la mayoría de las personas que intentan ciegamente usar geometría pitagórica / euclidiana, ni siquiera piensan o saben sobre eso. Sería útil si este hecho (debe usar Mod 360) estuviera presente en su respuesta.
Devdatta Tengshe
1
@ToolmakerSteve Está bien, a menudo uso esta corrección, pero espero que su usuario entienda que es una aproximación y que para distancias grandes y algunas otras circunstancias puede estar muy lejos de la realidad.
whuber
8

Interpreté "distancia pitagórica" ​​como "distancia euclidiana". Entonces la respuesta es la misma que "¿cuál es la diferencia entre la longitud de un acorde de un círculo y el perímetro subtendido?" Deje que el radio sea R, el ángulo subtendido es A (radianes).

perimeter = L = A*R
chord = C = 2*sin(A/2)*R
diff = D = L - C
     = (A-2*sin(A/2))*R
     = A^3/24 * R  (for A small)
     = L^3/(24*R^2) (eliminating A)
relative error = D/L
               = (L/R)^2/24

Para la tierra, sustituya R = 6400 km. Por cierto, llámelo "distancia de gran círculo" (lo que es) no "distancia de haversine" (cómo se calcula). (Esto es similar a la distinción entre la distancia pitagórica y la distancia euclidiana).

cffk
fuente
Siguiendo su razonamiento, puede sustituir aún más L y hacer que la estimación solo requiera A.
lynxlynxlynx
¿Puedes elaborar la expresión que terminas? ¿Cómo salió A ^ 3/24 * R?
curioso
Expanda sin (A / 2) para A pequeño, usando sin (x) = x - x ^ 3/6, y obtendrá este resultado.
cffk
5

Para obtener una respuesta completa y rigurosa, consulte la respuesta de Whuber arriba. Voy a responder de una manera más visual y básica.

La razón por la cual los cálculos planos / pitagóricos son inapropiados es porque los cálculos se basan en el hecho de que moverse un paso en cualquier dirección es un cambio constante de magnitud, independientemente de dónde se encuentre en el gráfico.

gráfico simple

La longitud no se ajusta a este requisito. Las líneas de longitud convergen en los polos.

globo mostrando convergencia

Es por eso que cuando aplanamos la Tierra para reflejar las reglas de un gráfico plano, obtenemos distorsión.

Mapa de proyección de Mercator

Si observa ese mapa, parece que Groenlandia es aproximadamente del tamaño de África y la Antártida es aproximadamente del tamaño de Eurasia. Por supuesto que eso no es cierto. Groenlandia y la Antártida están extremadamente distorsionadas porque están cerca de los polos donde converge la longitud.

vista del globo del hemisferio norte

Como puede ver, Groenlandia es aproximadamente del tamaño de México.

vista del guante del hemisferio sur

Y la Antártida es aproximadamente del tamaño del sur de África (no Sudáfrica).

Como puede ver, los errores que obtendrá aplicando fórmulas pitagóricas dependen más de dónde están los puntos que de la distancia entre los puntos. Con la advertencia importante de que distancias más largas aumentarán cualquier error. Es por eso que las soluciones planas, aunque tentadoras, son una mala elección. Las distorsiones te morderán y no es tan simple como un desplazamiento. Los errores son el resultado de deformar la tierra para que se ajuste a reglas inapropiadas.

Erik
fuente
En realidad, lo que está mostrando es un tipo diferente de error. Utilizado correctamente, el teorema de Pitágoras calcula la distancia de longitud en función de la longitud a lo largo de la línea de latitud en la que se encuentra, multiplicada porcos(lat) . Utilizado de esta manera, los errores son pequeños para distancias pequeñas, en cualquier lugar de una esfera (excepto si pasa el polo N o S). Lo que está mostrando es una distorsión de una proyección de toda la tierra, donde inevitablemente algunas regiones están muy distorsionadas. "los errores que obtendrá ... dependen más de dónde ... que de la distancia" no es cierto si se usa * cos(lat).
ToolmakerSteve