Geodesia |
---|
El estudio de las geodésicas en un elipsoide surgió en conexión con la geodesia específicamente con la solución de redes de triangulación . La figura de la Tierra está bien aproximada por un elipsoide achatado, una esfera ligeramente aplanada. Una geodésica es la ruta más corta entre dos puntos en una superficie curva, análoga a una línea recta en una superficie plana. La solución de una red de triangulación en un elipsoide es, por tanto, un conjunto de ejercicios de trigonometría esferoidal ( Euler 1755 ).
Si la Tierra se trata como una esfera , las geodésicas son grandes círculos (todos los cuales están cerrados) y los problemas se reducen a unos en trigonometría esférica . Sin embargo, Newton (1687) demostró que el efecto de la rotación de la Tierra hace que parezca un elipsoide ligeramente achatado: en este caso, el ecuador y los meridianos son las únicas geodésicas cerradas simples. Además, el camino más corto entre dos puntos del ecuador no corre necesariamente a lo largo del ecuador. Finalmente, si el elipsoide se perturba aún más para convertirse en un elipsoide triaxial (con tres semiejes distintos), solo se cierran tres geodésicas.
Hay varias formas de definir las geodésicas ( Hilbert y Cohn-Vossen 1952 , págs. 220-221). Una definición simple es como el camino más corto entre dos puntos en una superficie. Sin embargo, con frecuencia es más útil definirlos como caminos con curvatura geodésica cero , es decir, el análogo de las líneas rectas.sobre una superficie curva. Esta definición abarca geodésicas que viajan tan lejos a través de la superficie del elipsoide que comienzan a regresar hacia el punto de partida, de modo que otras rutas son más directas e incluye caminos que se cruzan o se vuelven a trazar. Los segmentos suficientemente cortos de una geodésica siguen siendo la ruta más corta entre sus puntos finales, pero las geodésicas no son necesariamente mínimas globalmente (es decir, las más cortas entre todas las rutas posibles). Cada camino globalmente más corto es una geodésica, pero no al revés.
A finales del siglo XVIII, un elipsoide de revolución ( también se utiliza el término esferoide ) era una aproximación bien aceptada a la figura de la Tierra . El ajuste de las redes de triangulación implicó reducir todas las mediciones a un elipsoide de referencia y resolver el problema bidimensional resultante como un ejercicio de trigonometría esferoidal ( Bomford 1952 , Cap. 3) ( Leick et al. 2015 , §4.5).
Es posible reducir los diversos problemas geodésicos en uno de dos tipos. Considere dos puntos: A en latitud φ 1 y longitud λ 1 y B en latitud φ 2 y longitud λ 2 (ver Fig. 1). La geodésica de conexión (de A a B ) es AB , de longitud s 12 , que tiene acimutes α 1 y α 2 en los dos extremos. [1] Los dos problemas geodésicos que se suelen considerar son:
Como puede verse en la Fig. 1, estos problemas implican resolver el triángulo NAB dado un ángulo, α 1 para el problema directo y λ 12 = λ 2 - λ 1 para el problema inverso, y sus dos lados adyacentes. Para una esfera, las soluciones a estos problemas son simples ejercicios de trigonometría esférica , cuya solución viene dada por fórmulas para resolver un triángulo esférico . (Consulte el artículo sobre la navegación en círculo máximo ).
Para un elipsoide de revolución, Clairaut (1735) encontró la constante característica que define la geodésica . Legendre (1806) y Oriani (1806) (y artículos posteriores en 1808 y 1810 ) dieron una solución sistemática para los caminos de las geodésicas . Bessel (1825) da la solución completa para el problema directo (con tablas computacionales y un ejemplo resuelto ) .
Durante el siglo XVIII, las geodésicas se denominaban típicamente "líneas más cortas". El término "línea geodésica" (en realidad, una curva ) fue acuñado por Laplace (1799b) :
Nous désignerons cette ligne sous le nom de ligne géodésique [Llamaremos a esta línea la línea geodésica ].
Esta terminología se introdujo en inglés como "línea geodésica" o como "línea geodésica", por ejemplo ( Hutton 1811 ),
Una línea trazada de la manera que hemos estado describiendo, o deducida de medidas trigonométricas, por los medios que hemos indicado, se llama línea geodésica o geodésica: tiene la propiedad de ser la más corta que se puede trazar entre sus dos extremos en la superficie de la Tierra; y, por tanto, es la medida de itinerario adecuada de la distancia entre esos dos puntos.
En su adopción por otros campos , se prefirió la línea geodésica , frecuentemente abreviada a geodésica .
Esta sección trata el problema en un elipsoide de revolución (tanto oblato como prolato). El problema de un elipsoide triaxial se trata en la siguiente sección.
Aquí se desarrollan las ecuaciones para una geodésica; la derivación sigue de cerca la de Bessel (1825) . Jordan y Eggert (1941) , Bagratuni (1962 , §15), Gan'shin (1967 , Cap.5 ), Krakiwsky y Thomson (1974 , §4), Rapp (1993 , §1.2), Jekeli (2012) , y Borre y Strang (2012) también proporcionan derivaciones de estas ecuaciones.
Considere un elipsoide de revolución con ecuatorial radio de una y polar semi-eje b . Defina el aplanamiento f = ( a - b ) / a , la excentricidad e = √ a 2 - b 2 / a = √ f (2 - f ) , y la segunda excentricidad e ′ = √ a 2 - b 2 / b = e / (1 - f ). (En la mayoría de las aplicaciones en geodesia, el elipsoide se considera achatado, a > b ; sin embargo, la teoría se aplica sin cambios a los elipsoides prolongados, a < b , en cuyo caso f , e 2 y e ′ 2 son negativos).
Supongamos que un segmento elemental de una trayectoria en el elipsoide tiene una longitud ds . De las Figs. 2 y 3, vemos que si su acimut es α , entonces ds está relacionado con d φ y d λ por
donde ρ es el radio de curvatura meridional , R = ν cosφ es el radio del círculo de latitud φ y ν es el radio de curvatura normal . Por tanto, el segmento elemental está dado por
o
donde φ ′ = d φ / d λ y la función lagrangiana L depende de φ a ρ (φ) y R (φ) . La longitud de una ruta arbitraria entre (φ 1 , λ 1 ) y (φ 2 , λ 2 ) está dada por
donde φ es una función de λ que satisface φ (λ 1 ) = φ 1 y φ (λ 2 ) = φ 2 . El camino más corto o geodésico implica encontrar esa función φ (λ) que minimiza s 12 . Este es un ejercicio de cálculo de variaciones y la condición de minimización viene dada por la identidad Beltrami ,
Sustituyendo L y usando las ecuaciones. (1) da
Clairaut (1735) encontró esta relación , utilizando una construcción geométrica; Lyusternik (1964 , § 10) presenta una derivación similar . [2] Al diferenciar esta relación se obtiene
Esto, junto con las Ecs. (1), conduce a un sistema de ecuaciones diferenciales ordinarias para una geodésica
Podemos expresar R en términos de latitud paramétrica , β , usando
y la relación de Clairaut se convierte en
Esta es la regla del seno de la trigonometría esférica que relaciona dos lados del triángulo NAB (ver Fig.4), NA = 1 ⁄ 2 π - β 1 , y NB = 1 ⁄ 2 π - β 2 y sus ángulos opuestos B = π - α 2 y A = α 1 .
Para encontrar la relación para el tercer lado AB = σ 12 , la longitud del arco esférico y el ángulo incluido N = ω 12 , la longitud esférica , es útil considerar el triángulo NEP que representa una geodésica que comienza en el ecuador; ver Fig. 5. En esta figura, las variables referidas a la esfera auxiliar se muestran con las cantidades correspondientes para el elipsoide mostradas entre paréntesis. Las cantidades sin subíndices se refieren al punto arbitrario P ; E , el punto en el que la geodésica cruza el ecuador en dirección norte, se utiliza como origen para σ , sy ω .
Si el lado EP se extiende moviendo P infinitesimalmente (ver Fig.6), obtenemos
Combinando Ecs. (1) y (2) da ecuaciones diferenciales para s y λ
La relación entre β y φ es
lo que da
de modo que las ecuaciones diferenciales de la geodésica se conviertan en
El último paso es usar σ como parámetro independiente en ambas de estas ecuaciones diferenciales y de este modo para expresar s y λ como integrales. Aplicando la regla del seno a los vértices E y G en el triángulo esférico EGP en la Fig.5 da
donde α 0 es el acimut en E . Sustituyendo esto en la ecuación para ds / d σ e integrando el resultado da
donde
y los límites de la integral se eligen de modo que s (σ = 0) = 0 . Legendre (1811 , p. 180) señaló que la ecuación para s es la misma que la ecuación para el arco en una elipse con semiejes b √ 1 + e ′ 2 cos 2 α 0 y b . Para expresar la ecuación de λ en términos de σ , escribimos
que se sigue de la Ec. (2) y la relación de Clairaut. Esto produce
y los límites de las integrales se eligen de modo que λ = λ 0 en el cruce del ecuador, σ = 0 .
Esto completa la solución de la trayectoria de una geodésica utilizando la esfera auxiliar. Mediante este dispositivo, un gran círculo se puede asignar exactamente a una geodésica en un elipsoide de revolución.
También hay varias formas de aproximar las geodésicas en un elipsoide terrestre (con un pequeño aplanamiento) ( Rapp 1991 , §6); algunos de ellos se describen en el artículo sobre la distancia geográfica . Sin embargo, estos son típicamente comparables en complejidad al método para la solución exacta ( Jekeli 2012 , §2.1.4).
La Fig. 7 muestra las geodésicas cerradas simples que consisten en los meridianos (verde) y el ecuador (rojo). (Aquí, la calificación "simple" significa que la geodésica se cierra sobre sí misma sin una auto-intersección intermedia). Esto se sigue de las ecuaciones para las geodésicas dadas en la sección anterior.
Todas las demás geodésicas están tipificadas por las Figs. 8 y 9 que muestran una geodésica que comienza en el ecuador con α 0 = 45 ° . La geodésica oscila alrededor del ecuador. Los cruces ecuatoriales se denominan nodos y los puntos de máxima o mínima latitud se denominan vértices ; las latitudes paramétricas de los vértices están dadas por β = ± ( 1 ⁄ 2 π - | α 0 |) . La geodésica completa una oscilación completa en latitud antes de que la longitud haya aumentado en360 ° . Por lo tanto, en cada cruce sucesivo del ecuador hacia el norte (ver Fig.8), λ no llega a un circuito completo del ecuador en aproximadamente 2π f sinα 0 (para un elipsoide alargado, esta cantidad es negativa y λ completa más que un circuito; consulte la Fig.10). Para casi todos los valores de α 0 , la geodésica llenará esa porción del elipsoide entre las dos latitudes de vértice (ver Fig. 9).
Si el elipsoide es suficientemente achatado, es decir, b ⁄ a < 1 ⁄ 2 , es posible otra clase de geodésicas cerradas simples ( Klingenberg 1982 , §3.5.19). Dos de tales geodésicas se ilustran en las Figs. 11 y 12. Aquí b ⁄ a = 2 ⁄ 7 y el acimut ecuatorial, α 0 , para la geodésica verde (resp. Azul) se elige para ser53,175 ° (resp.75.192 ° ), de modo que la geodésica completa 2 (resp. 3) oscilaciones completas alrededor del ecuador en un circuito del elipsoide.
La Fig.13 muestra geodésicas (en azul) que emanan A con α 1 un múltiplo de15 ° hasta el punto en el que dejan de ser caminos más cortos. (El aplanamiento se ha aumentado a 1 / 10 con el fin de acentuar los efectos elipsoidales.) También se muestra (en verde) son curvas de constante s 12 , que son los círculos geodésicas centrados A . Gauss (1828) mostró que, en cualquier superficie, la geodésica y el círculo geodésico se cruzan en ángulos rectos. La línea roja es el locus de corte , el lugar geométrico de puntos que tienen múltiples (dos en este caso) geodésicas más cortas desde A . En una esfera, el lugar de corte es un punto. En un elipsoide achatado (mostrado aquí), es un segmento del círculo de latitud centrado en el punto antípodaa A , φ = −φ 1 . La extensión longitudinal del lugar de corte es aproximadamente λ 12 ∈ [π - f π cosφ 1 , π + f π cosφ 1 ] . Si A se encuentra en el ecuador, φ 1 = 0 , esta relación es exacta y, como consecuencia, el ecuador es solo una geodésica más corta si | λ 12 | ≤ (1 - f ) π . Para un elipsoide alargado, el lugar de corte es un segmento del anti-meridiano centrado en el punto antípoda de A , λ 12 = π, y esto significa que las geodésicas meridionales dejan de ser caminos más cortos antes de que se alcance el punto antípoda.
Varios problemas relacionados con las geodésicas requieren conocer su comportamiento cuando están perturbados. Esto es útil en ajustes trigonométricos ( Ehlert 1993 ), determinando las propiedades físicas de las señales que siguen a las geodésicas, etc. Considere una geodésica de referencia, parametrizada por s , y una segunda geodésica a una pequeña distancia t ( s ) de ella. Gauss (1828) mostró que t ( s ) obedece a la ecuación de Gauss-Jacobi
donde K ( s ) es la curvatura gaussiana en s . Como ecuación diferencial homogénea, lineal y de segundo orden, su solución puede expresarse como la suma de dos soluciones independientes
donde
La cantidad m ( s 1 , s 2 ) = m 12 es la denominada longitud reducida , y M ( s 1 , s 2 ) = M 12 es la escala geodésica . [3] Sus definiciones básicas se ilustran en la Fig. 14.
La curvatura gaussiana para un elipsoide de revolución es
Helmert (1880 , ecuación (6.5.1.)) Resolvió la ecuación de Gauss-Jacobi para este caso permitiendo que m 12 y M 12 se expresen como integrales.
Como vemos en la Fig. 14 (subfigura superior), la separación de dos geodésicas que comienzan en el mismo punto con acimutes que difieren en d α 1 es m 12 d α 1 . En una superficie cerrada como un elipsoide, m 12 oscila alrededor de cero. El punto en el que m 12 se convierte en cero es el punto conjugado con el punto de partida. Para que una geodésica entre A y B , de longitud s 12 , sea un camino más corto, debe satisfacer la condición de Jacobi ( Jacobi 1837 ) ( Jacobi 1866 , §6) ( Forsyth 1927, §§26-27) ( dicha de 1916 ), que no hay ningún punto conjugado a A entre A y B . Si no se cumple esta condición, entonces hay un camino cercano (no necesariamente geodésico) que es más corto. Por lo tanto, la condición de Jacobi es una propiedad local de la geodésica y es solo una condición necesaria para que la geodésica sea un camino global más corto. Las condiciones necesarias y suficientes para que una geodésica sea el camino más corto son:
Las geodésicas de un punto particular A si continúan más allá del lugar de corte forman una envolvente ilustrada en la Fig.15. Aquí las geodésicas para las cuales α 1 es un múltiplo de3 ° se muestran en azul claro. (Las geodésicas solo se muestran para su primer paso cerca del punto antípoda, no para los posteriores). Algunos círculos geodésicos se muestran en verde; estos forman cúspides en el sobre. El lugar de corte se muestra en rojo. La envolvente es el lugar geométrico de los puntos que se conjugan con A ; los puntos de la envolvente pueden calcularse encontrando el punto en el que m 12 = 0 en una geodésica. Jacobi (1891) llama astroide a esta figura en forma de estrella producida por el sobre .
Fuera del astroide, dos geodésicas se cruzan en cada punto; por tanto, hay dos geodésicas (con una longitud de aproximadamente la mitad de la circunferencia del elipsoide) entre A y estos puntos. Esto corresponde a la situación en la esfera donde hay rutas "cortas" y "largas" en un gran círculo entre dos puntos. Dentro del astroide, cuatro geodésicas se cruzan en cada punto. En la Fig. 16 se muestran cuatro de estas geodésicas, donde las geodésicas están numeradas en orden de longitud creciente. (Esta figura usa la misma posición para A que la Fig. 13 y está dibujada en la misma proyección). Las dos geodésicas más cortas son estables , es decir, m 12 > 0, de modo que no haya un camino cercano que conecte los dos puntos que es más corto; los otros dos son inestables. Solo la línea más corta (la primera) tiene σ 12 ≤ π . Todas las geodésicas son tangentes a la envolvente que se muestra en verde en la figura.
El astroid es la (exterior) evolute de los círculos geodésicas centrado en A . Asimismo, los círculos geodésicos son involutas del astroide.
Un polígono geodésico es un polígono cuyos lados son geodésicos. Es análogo a un polígono esférico , cuyos lados son grandes círculos. El área de dicho polígono se puede encontrar calculando primero el área entre un segmento geodésico y el ecuador, es decir, el área del cuadrilátero AFHB en la Fig. 1 ( Danielsen 1989 ). Una vez que se conoce esta área, el área de un polígono se puede calcular sumando las contribuciones de todos los bordes del polígono.
Aquí se desarrolla una expresión para el área S 12 de AFHB siguiendo a Sjöberg (2006) . El área de cualquier región cerrada del elipsoide es
donde dT es un elemento de superficie y K es la curvatura gaussiana . Ahora, el teorema de Gauss-Bonnet aplicado a un polígono geodésico establece
donde
es el exceso geodésico y θ j es el ángulo exterior en el vértice j . Al multiplicar la ecuación para Γ por R 2 2 , donde R 2 es el radio autálico , y restar esto de la ecuación para T da
donde se ha sustituido el valor de K por un elipsoide . Aplicando esta fórmula al cuadrilátero AFHB , observando que Γ = α 2 - α 1 , y realizando la integral sobre φ da
donde la integral está sobre la línea geodésica (de modo que φ es implícitamente una función de λ ). La integral se puede expresar como una serie válida para f pequeña ( Danielsen 1989 ) ( Karney 2013 , §6 y apéndice).
El área de un polígono geodésico se obtiene sumando S 12 sobre sus bordes. Este resultado es válido siempre que el polígono no incluya un poste; si es así, se debe sumar 2π R 2 2 a la suma. Si los bordes se especifican por sus vértices, entonces una expresión conveniente para el exceso geodésico E 12 = α 2 - α 1 es
Resolver los problemas geodésicos implica mapear la geodésica en la esfera auxiliar y resolver el problema correspondiente en la navegación de círculo máximo . Al resolver el triángulo esférico "elemental" para NEP en la Fig. 5, se pueden emplear las reglas de Napier para triángulos cuadrantes ,
El mapeo de la geodésica implica evaluar las integrales para la distancia, s , y la longitud, λ , Ecs. (3) y (4) y estos dependen del parámetro α 0 .
Manejar el problema directo es sencillo, porque α 0 se puede determinar directamente a partir de las cantidades dadas φ 1 y α 1 .
En el caso del problema inverso, se da λ 12 ; esto no puede relacionarse fácilmente con el ángulo esférico equivalente ω 12 porque se desconoce α 0 . Por tanto, la solución del problema requiere que α 0 se encuentre iterativamente.
En aplicaciones geodésicas, donde f es pequeña, las integrales se evalúan típicamente como una serie ( Legendre 1806 ) ( Oriani 1806 ) ( Bessel 1825 ) ( Helmert 1880 ) ( Rainsford 1955 ) ( Rapp 1993 ). Para f arbitraria , las integrales (3) y (4) se pueden encontrar por cuadratura numérica o expresándolas en términos de integrales elípticas ( Legendre 1806 ) ( Cayley 1870 ).
Vincenty (1975) proporciona soluciones para los problemas directos e inversos; Estos se basan en una expansión en serie llevada a cabo a tercer orden en el aplanamiento y proporcionan una precisión de aproximadamente0,1 mm para el elipsoide WGS84 ; sin embargo, el método inverso no logra converger para los puntos casi antípodas. Karney (2013) continúa las expansiones hasta el sexto orden, lo que es suficiente para proporcionar una precisión de doble precisión completa para | f | ≤ 1 ⁄ 50 y mejora la solución del problema inverso para que converja en todos los casos. Karney (2013 , apéndice) amplía el método para utilizar integrales elípticas que se pueden aplicar a elipsoides con aplanamiento arbitrario.
Resolver el problema geodésico para un elipsoide de revolución es, desde el punto de vista matemático, relativamente simple: debido a la simetría, las geodésicas tienen una constante de movimiento , dada por la relación de Clairaut que permite reducir el problema a cuadratura . A principios del siglo XIX (con el trabajo de Legendre, Oriani , Bessel, et al.), Había una comprensión completa de las propiedades de las geodésicas en un elipsoide de revolución.
Por otro lado, las geodésicas en un elipsoide triaxial (con tres ejes desiguales) no tienen una constante obvia del movimiento y, por lo tanto, representaron un problema desafiante sin resolver en la primera mitad del siglo XIX. En un artículo notable, Jacobi (1839) descubrió una constante del movimiento que permitía reducir este problema a la cuadratura también ( Klingenberg 1982 , §3.5). [4]
Considere el elipsoide definido por
donde ( X , Y , Z ) son coordenadas cartesianas centradas en el elipsoide y, sin pérdida de generalidad, a ≥ b ≥ c > 0 . [5]Jacobi (1866 , §§26-27) empleó las coordenadas elipsoidales (triaxiales) (con latitud elipsoidal triaxial y longitud elipsoidal triaxial , β, ω ) definidas por
En el límite b → a , β se convierte en la latitud paramétrica de un elipsoide achatado, por lo que el uso del símbolo β es consistente con las secciones anteriores. Sin embargo, ω es diferente de la longitud esférica definida anteriormente. [6]
Las líneas de la cuadrícula de constante β (en azul) y ω (en verde) se dan en la Fig. 17. Éstas constituyen un sistema de coordenadas ortogonales : las líneas de la cuadrícula se cruzan en ángulos rectos. Las secciones principales del elipsoide, definidas por X = 0 y Z = 0, se muestran en rojo. La tercera sección principal, Y = 0 , está cubierta por las líneas β = ± 90 ° y ω = 0 ° o ± 180 ° . Estas líneas se encuentran en cuatro puntos umbilicales (dos de los cuales son visibles en esta figura) donde los principales radios de curvaturason iguales. Aquí y en las otras figuras de esta sección, los parámetros del elipsoide son a : b : c = 1.01: 1: 0.8 , y se ve en una proyección ortográfica desde un punto por encima de φ = 40 ° , λ = 30 ° .
Las líneas de cuadrícula de las coordenadas elipsoidales se pueden interpretar de tres formas diferentes:
Jacobi demostró que las ecuaciones geodésicas, expresadas en coordenadas elipsoidales, son separables. Así es como le contó su descubrimiento a su amigo y vecino Bessel ( Jacobi 1839 , Carta a Bessel):
Anteayer reduje a cuadratura el problema de las líneas geodésicas en un elipsoide con tres ejes desiguales . Son las fórmulas más simples del mundo, las integrales abelianas , que se convierten en las bien conocidas integrales elípticas si 2 ejes se igualan.
Königsberg , 28 de diciembre de '38.
La solución dada por Jacobi ( Jacobi 1839 ) ( Jacobi 1866 , §28) es
Como señala Jacobi, "una función del ángulo β es igual a una función del ángulo ω . Estas dos funciones son simplemente integrales abelianas ..." En la solución aparecen dos constantes δ y γ . Normalmente, δ es cero si los límites inferiores de las integrales se toman como el punto de partida de la geodésica y la dirección de las geodésicas está determinada por γ . Sin embargo, para las geodésicas que comienzan en un punto umbilical, tenemos γ = 0 y δ determina la dirección en el punto umbilical. La constante γ puede expresarse como
donde α es el ángulo que forma la geodésica con líneas de constante ω . En el límite b → a , esto se reduce a sinα cosβ = const. , la familiar relación de Clairaut. Darboux (1894 , §§583–584) da una derivación del resultado de Jacobi ; da la solución encontrada por Liouville (1846) para superficies cuadráticas generales.
En un elipsoide triaxial, solo hay tres geodésicas cerradas simples, las tres secciones principales del elipsoide dadas por X = 0 , Y = 0 y Z = 0 . [7] Para estudiar las otras geodésicas, es conveniente considerar las geodésicas que intersecan la sección principal central, Y = 0 , en ángulos rectos. Tales geodésicas se muestran en las Figs. 18-22, que utilizan los mismos parámetros de elipsoide y la misma dirección de visión que la Fig. 17. Además, las tres elipses principales se muestran en rojo en cada una de estas figuras.
Si el punto de partida es β 1 ∈ (−90 °, 90 °) , ω 1 = 0 y α 1 = 90 ° , entonces γ> 0 y la geodésica rodea el elipsoide en un sentido "circumpolar". La geodésica oscila al norte y al sur del ecuador; en cada oscilación completa un poco menos de un circuito completo alrededor del elipsoide, lo que da como resultado, en el caso típico, que la geodésica llene el área delimitada por las dos líneas de latitud β = ± β 1 . Se dan dos ejemplos en las Figs. 18 y 19. La Figura 18 muestra prácticamente el mismo comportamiento que para un elipsoide achatado de revolución (porque a ≈ b); compárese con la Fig. 9. Sin embargo, si el punto de partida está en una latitud más alta (Fig. 18), las distorsiones resultantes de a ≠ b son evidentes. Todas las tangentes a una geodésica circumpolar tocan el hiperboloide confocal de una sola hoja que cruza el elipsoide en β = β 1 ( Chasles 1846 ) ( Hilbert y Cohn-Vossen 1952 , págs. 223-224).
Si el punto de partida es β 1 = 90 ° , ω 1 ∈ (0 °, 180 °) y α 1 = 180 ° , entonces γ <0 y la geodésica rodea el elipsoide en un sentido "transpolar". La geodésica oscila al este y oeste de la elipse X = 0 ; en cada oscilación completa un poco más de un circuito completo alrededor del elipsoide. En el caso típico, esto da como resultado que la geodésica llene el área delimitada por las dos líneas de longitud ω = ω 1 y ω = 180 ° - ω 1 . Si a = b , todos los meridianos son geodésicos; el efecto de un≠ b hace que dichas geodésicas oscilen de este a oeste. Se dan dos ejemplos en las Figs. 20 y 21. La constricción de la geodésica cerca del polo desaparece en el límite b → c ; en este caso, el elipsoide se convierte en un elipsoide alargado y la figura 20 se parecería a la figura 10 (girada de lado). Todas las tangentes a una geodésica transpolar tocan el hiperboloide confocal de doble hoja que corta al elipsoide en ω = ω 1 .
Si el punto de partida es β 1 = 90 ° , ω 1 = 0 ° (un punto umbilical) y α 1 = 135 ° (la geodésica deja la elipse Y = 0 en ángulo recto), entonces γ = 0 y la geodésica repetidamente intersecta el punto umbilical opuesto y vuelve a su punto de partida. Sin embargo, en cada circuito el ángulo en el que se cruza Y = 0 se acerca más a0 ° o180 ° de modo que asintóticamente la geodésica se encuentra en la elipse Y = 0 ( Hart 1849 ) ( Arnold 1989 , p. 265), como se muestra en la Fig. 22. Una sola geodésica no llena un área en el elipsoide. Todas las tangentes a las geodésicas umbilicales tocan la hipérbola confocal que cruza el elipsoide en los puntos umbilicales.
Las geodésicas umbilicales gozan de varias propiedades interesantes.
Si el punto inicial A de una geodésica no es un punto umbilical, su envoltura es un astroide con dos cúspides en β = −β 1 y las otras dos en ω = ω 1 + π . El lugar de corte de A es la parte de la línea β = −β 1 entre las cúspides.
Los problemas geodésicos directos e inversos ya no juegan el papel central en la geodesia que alguna vez tuvieron. En lugar de resolver el ajuste de las redes geodésicas como un problema bidimensional en la trigonometría esferoidal, estos problemas ahora se resuelven mediante métodos tridimensionales ( Vincenty y Bowring 1978 ). Sin embargo, las geodésicas terrestres todavía juegan un papel importante en varias áreas:
Por el principio de acción mínima , muchos problemas de la física pueden formularse como un problema variacional similar al de las geodésicas. De hecho, el problema geodésico es equivalente al movimiento de una partícula obligada a moverse en la superficie, pero por lo demás no sujeta a fuerzas ( Laplace 1799a ) ( Hilbert y Cohn-Vossen 1952 , p. 222). Por esta razón, las geodésicas en superficies simples como elipsoides de revolución o elipsoides triaxiales se utilizan con frecuencia como "casos de prueba" para explorar nuevos métodos. Ejemplos incluyen: