stringtranslate.com

Геодезические на эллипсоиде

Геодезическая на сплющенном эллипсоиде

Изучение геодезических на эллипсоиде возникло в связи с геодезией, в частности, с решением триангуляционных сетей . Фигура Земли хорошо аппроксимируется сплющенным эллипсоидом , слегка сплющенной сферой. Геодезическая — это кратчайший путь между двумя точками на искривленной поверхности, аналогичный прямой на плоской поверхности. Решение триангуляционной сети на эллипсоиде представляет собой, таким образом, набор упражнений по сфероидальной тригонометрии (Эйлер 1755).

Если рассматривать Землю как сферу , то геодезические линии представляют собой большие окружности (все из которых замкнуты), и проблемы сводятся к проблемам сферической тригонометрии . Однако Ньютон (1687) показал, что эффект вращения Земли приводит к тому, что она становится похожей на слегка сплющенный эллипсоид: в этом случае экватор и меридианы являются единственными простыми замкнутыми геодезическими. Более того, кратчайший путь между двумя точками на экваторе не обязательно проходит по экватору. Наконец, если эллипсоид еще больше возмущен, чтобы стать трехосным эллипсоидом (с тремя различными полуосями), замкнутыми будут только три геодезические.

Геодезические на эллипсоиде вращения

Существует несколько способов определения геодезических (Hilbert & Cohn-Vossen 1952, стр. 220–221). Простое определение — как кратчайший путь между двумя точками на поверхности. Однако часто бывает полезнее определять их как пути с нулевой геодезической кривизной — т. е. аналог прямых линий на искривленной поверхности. Это определение охватывает геодезические, проходящие так далеко по поверхности эллипсоида, что они начинают возвращаться к начальной точке, так что другие маршруты являются более прямыми, и включает пути, которые пересекаются или повторяют сами себя. Достаточно короткие сегменты геодезических по-прежнему являются кратчайшим путем между своими конечными точками, но геодезические не обязательно являются глобально минимальными (т. е. кратчайшими среди всех возможных путей). Каждый глобально кратчайший путь является геодезическим, но не наоборот.

К концу XVIII века эллипсоид вращения (также используется термин сфероид ) был общепринятым приближением к фигуре Земли . Корректировка сетей триангуляции повлекла за собой сведение всех измерений к опорному эллипсоиду и решение полученной двумерной задачи в качестве упражнения по сфероидальной тригонометрии (Бомфорд 1952, гл. 3) (Лейк и др. 2015, §4.5).

Рис. 1. Геодезическая AB на эллипсоиде вращения. N — северный полюс, а EFH лежит на экваторе.

Различные геодезические задачи можно свести к одному из двух типов. Рассмотрим две точки: A на широте φ 1 и долготе λ 1 и B на широте φ 2 и долготе λ 2 (см. рис. 1). Соединительная геодезическая (от A до B ) — это AB , длиной s 12 , которая имеет азимуты α 1 и α 2 в двух конечных точках. [1] Обычно рассматриваются две геодезические задачи:

  1. прямая геодезическая задача или первая геодезическая задача , заданные A , α 1 и s 12 , определяют B и α 2 ;
  2. обратная геодезическая задача или вторая геодезическая задача , при заданных A и B , определяет s 12 , α 1 и α 2 .

Как видно из рис. 1, эти задачи включают решение треугольника NAB по одному углу, α 1 для прямой задачи и λ 12 = λ 2λ 1 для обратной задачи, и его двум смежным сторонам. Для сферы решения этих задач являются простыми упражнениями по сферической тригонометрии , решение которых дается формулами для решения сферического треугольника . (См. статью о навигации по большому кругу .)

Для эллипсоида вращения характеристическая постоянная, определяющая геодезическую, была найдена Клеро (1735). Систематическое решение для путей геодезических было дано Лежандром (1806) и Ориани (1806) (и последующими работами 1808 и 1810 годов). Полное решение прямой задачи (с вычислительными таблицами и разработанным примером) дано Бесселем (1825).

В XVIII веке геодезические линии обычно называли «кратчайшими линиями». Термин «геодезическая линия» (на самом деле, кривая ) был придуман Лапласом (1799b):

Nous désignerons cette ligne sous le nom de ligne géodesique [Мы назовем эту линию геодезической линией ].

Этот термин был введен в английский язык как «геодезическая линия» или как «геодезическая линия», например (Хаттон 1811, стр. 115),

Линия, проведенная описанным нами способом или выведенная из тригонометрических измерений указанными нами способами, называется геодезической или геодезической линией: она обладает тем свойством, что является кратчайшей линией, которую можно провести между ее двумя концами на поверхности Земли; и поэтому она является надлежащей мерой расстояния между этими двумя точками.

При принятии его в других областях науки предпочтение отдавалось геодезической линии , часто сокращаемой до геодезической .

В этом разделе рассматривается задача на эллипсоиде вращения (как сплющенном, так и вытянутом). Задача на трехосном эллипсоиде рассматривается в следующем разделе.

Уравнения геодезической

Здесь разрабатываются уравнения для геодезической; вывод близко следует выводу Бесселя (1825). Джордан и Эггерт (1941), Багратуни (1962, §15), Ганшин (1967, гл. 5), Кракивский и Томсон (1974, §4), Рапп (1993, §1.2), Джекели (2012) и Борре и Стрэнг (2012) также приводят выводы этих уравнений.

Рассмотрим эллипсоид вращения с экваториальным радиусом a и полярной полуосью b . Определим сплющивание f , эксцентриситет e и второй эксцентриситет e :

(В большинстве приложений в геодезии эллипсоид считается сплюснутым, a > b ; однако эта теория применима без изменений к вытянутым эллипсоидам, a < b , и в этом случае f , e 2 и e2 отрицательны.)

Пусть элементарный отрезок пути на эллипсоиде имеет длину ds . Из рис. 2 и 3 видно, что если его азимут равен α , то ds связан с и соотношением

(1)

где ρмеридиональный радиус кривизны , R = ν cos φ — радиус круга широты φ , а νнормальный радиус кривизны . Элементарный сегмент, таким образом, задается как

или

где φ ′ = / и функция Лагранжа L зависит от φ через ρ( φ ) и R ( φ ) . Длина произвольного пути между ( φ 1 , λ 1 ) и ( φ 2 , λ 2 ) определяется как

где φ — функция λ, удовлетворяющая φ ( λ 1 ) = φ 1 и φ ( λ 2 ) = φ 2 . Кратчайший путь или геодезическая влечет за собой нахождение той функции φ ( λ ), которая минимизирует s 12 . Это упражнение по вариационному исчислению , а условие минимизации задается тождеством Бельтрами ,

Подставляя L и используя уравнения (1) , получаем

Клеро (1735) нашел это соотношение , используя геометрическое построение; аналогичный вывод представлен Люстерником (1964, §10). [2] Дифференцирование этого соотношения дает

Это, вместе с уравнениями (1) , приводит к системе обыкновенных дифференциальных уравнений для геодезической

Мы можем выразить R через параметрическую широту β , используя

и отношение Клеро тогда становится

Это правило синусов сферической тригонометрии, связывающее две стороны треугольника NAB (см. рис. 4), NA = 12 πβ 1 , и NB = 12 πβ 2 , и их противолежащие углы B = π − α 2 и A = α 1 .

Чтобы найти соотношение для третьей стороны AB = σ 12 , длины сферической дуги и прилежащего угла N = ω 12 , сферической долготы , полезно рассмотреть треугольник NEP, представляющий геодезическую, начинающуюся на экваторе; см. рис. 5. На этом рисунке переменные, относящиеся к вспомогательной сфере, показаны с соответствующими величинами для эллипсоида, показанными в скобках. Величины без нижних индексов относятся к произвольной точке P ; E , точка, в которой геодезическая пересекает экватор в северном направлении, используется в качестве начала координат для σ , s и ω .

Рис. 6. Дифференциальный элемент геодезической на сфере.

Если сторону EP расширить, перемещая P бесконечно мало (см. рис. 6), то получим

(2)

Объединение уравнений (1) и (2) дает дифференциальные уравнения для s и λ

Соотношение между β и φ следующее:

что дает

так что дифференциальные уравнения для геодезической становятся

Последний шаг — использовать σ как независимый параметр в обоих этих дифференциальных уравнениях и тем самым выразить s и λ как интегралы. Применение правила синусов к вершинам E и G в сферическом треугольнике EGP на рис. 5 дает

где α 0 — азимут в точке E. Подставляя это в уравнение для ds / d σ и интегрируя результат, получаем

(3)

где

и пределы интеграла выбираются так, что s ( σ = 0) = 0 . Лежандр (1811, стр. 180) указал, что уравнение для s такое же, как уравнение для дуги эллипса с полуосями b 1 + e2 cos 2 α 0 и b . Чтобы выразить уравнение для λ через σ , запишем

что следует из уравнения 2 и соотношения Клеро. Это дает

(4)

и пределы интегралов выбираются так, что λ = λ 0 при пересечении экватора, σ = 0 .

Это завершает решение пути геодезической с использованием вспомогательной сферы. С помощью этого устройства большой круг может быть точно отображен на геодезическую на эллипсоиде вращения.

Существует также несколько способов аппроксимации геодезических на земном эллипсоиде (с небольшим сплющиванием) (Rapp 1991, §6); некоторые из них описаны в статье о географическом расстоянии . Однако они, как правило, сопоставимы по сложности с методом точного решения (Jekeli 2012, §2.1.4).

Поведение геодезических линий

Рис. 7. Меридианы и экватор являются единственными замкнутыми геодезическими. (Для очень сплющенных эллипсоидов существуют другие замкнутые геодезические; см. рис. 11 и 12).
Геодезическая на сплющенном эллипсоиде ( f = 150 ) с α 0 =45° .
Рис. 10. Геодезическая на вытянутом эллипсоиде ( f = − 150 ) с α 0 =45° . Сравните с рис. 8.

На рис. 7 показаны простые замкнутые геодезические, которые состоят из меридианов (зеленые) и экватора (красные). (Здесь определение «простые» означает, что геодезические замыкаются на себя без промежуточного самопересечения.) Это следует из уравнений для геодезических, приведенных в предыдущем разделе.

Все остальные геодезические проиллюстрированы на рис. 8 и 9, на которых показана геодезическая, начинающаяся на экваторе с α 0 = 45° . Геодезическая колеблется вокруг экватора. Экваториальные пересечения называются узлами , а точки максимальной или минимальной широты называются вершинами ; параметрические широты вершин задаются как β = ±( 12 π − |α 0 |) . Геодезическая завершает одно полное колебание по широте до того, как долгота увеличится на360° . Таким образом, при каждом последующем пересечении экватора на север (см. рис. 8) λ не дотягивает до полного оборота экватора примерно на 2 π f sinα 0 (для вытянутого эллипсоида эта величина отрицательна, и λ совершает более одного полного оборота; см. рис. 10). Почти для всех значений α 0 геодезическая заполнит ту часть эллипсоида, которая находится между двумя широтами вершин (см. рис. 9).

Две дополнительные замкнутые геодезические для сплюснутого эллипсоида, ba = 27 .

Если эллипсоид достаточно сплющен, т. е. ba < 12 , возможен другой класс простых замкнутых геодезических (Клингенберг 1982, §3.5.19). Две такие геодезические показаны на рис. 11 и 12. Здесь ba = 27 и экваториальный азимут, α 0 , для зеленой (соответственно синей) геодезической выбран равным53.175° (соотв.75,192° ), так что геодезическая совершает 2 (соответственно 3) полных колебания вокруг экватора за один оборот эллипсоида.

Рис. 13. Геодезические (синие) из одной точки для f = 110 , φ 1 =−30° ; геодезические окружности показаны зеленым цветом, а место сечения — красным.

На рис. 13 показаны геодезические (синего цвета), исходящие из точки А с α 1 кратным15° до точки, в которой они перестают быть кратчайшими путями. (Сглаживание было увеличено до 110 , чтобы подчеркнуть эллипсоидальные эффекты.) Также показаны (зеленым цветом) кривые постоянной s 12 , которые являются геодезическими окружностями с центром в точке A . Гаусс (1828) показал, что на любой поверхности геодезические и геодезические окружности пересекаются под прямым углом.

Красная линия — это разрезное место , место точек, которые имеют несколько (в данном случае две) кратчайших геодезических из A. На сфере разрезное место — это точка. На сплющенном эллипсоиде (показанном здесь) это сегмент круга широты с центром в точке, антиподальной к A , φ = − φ 1 . Продольная протяженность разрезного места приблизительно равна λ 12 ∈ [ π (1 − f cos φ 1 ), π (1 + f cos φ 1 )] . Если A лежит на экваторе, φ 1 = 0 , это соотношение является точным, и, как следствие, экватор является кратчайшей геодезической, только если | λ 12 | ≤ π (1 − f ) . Для вытянутого эллипсоида точкой разреза является отрезок антимеридиана с центром в точке, антиподальной к A , λ 12 = π , и это означает, что меридиональные геодезические перестают быть кратчайшими путями до того, как будет достигнута антиподальная точка.

Дифференциальные свойства геодезических

Различные проблемы, связанные с геодезическими, требуют знания их поведения при возмущении. Это полезно в тригонометрических корректировках (Ehlert 1993), определении физических свойств сигналов, которые следуют за геодезическими и т. д. Рассмотрим опорную геодезическую, параметризованную s , и вторую геодезическую на небольшом расстоянии t ( s ) от нее. Гаусс (1828) показал, что t ( s ) подчиняется уравнению Гаусса-Якоби

Рис. 14. Определение приведенной длины и геодезического масштаба.

где K ( s )гауссова кривизна в точке s . Как линейное однородное дифференциальное уравнение второго порядка, его решение может быть выражено как сумма двух независимых решений

где

Величина m ( s 1 , s 2 ) = m 12 является так называемой приведенной длиной , а M ( s 1 , s 2 ) = M 12 является геодезическим масштабом . [3] Их основные определения проиллюстрированы на рис. 14.

Гауссова кривизна для эллипсоида вращения равна

Гельмерт (1880, уравнение (6.5.1.)) решил уравнение Гаусса-Якоби для этого случая, что позволило выразить m 12 и M 12 в виде интегралов.

Как мы видим из рис. 14 (верхний подрисунок), разделение двух геодезических, начинающихся в одной точке с азимутами, отличающимися на d α 1 , равно m 12 d α 1 . На замкнутой поверхности, такой как эллипсоид, m 12 колеблется около нуля. Точка, в которой m 12 становится равным нулю, является точкой, сопряженной начальной точке. Для того чтобы геодезическая между A и B , длиной s 12 , была кратчайшим путем, она должна удовлетворять условию Якоби (Якоби 1837) (Якоби 1866, §6) (Форсайт 1927, §§26–27) (Блисс 1916), что между A и B нет точки, сопряженной с A . Если это условие не выполняется, то есть соседний путь (не обязательно геодезический), который короче. Таким образом, условие Якоби является локальным свойством геодезической и является лишь необходимым условием для того, чтобы геодезическая была глобальным кратчайшим путем. Необходимые и достаточные условия для того, чтобы геодезическая была кратчайшим путем:

Оболочка геодезических линий

Геодезические из одной точки ( f = 110 , φ 1 =−30° )

Геодезические из определенной точки A, если их продолжить за пределы сечения, образуют огибающую, показанную на рис. 15. Здесь геодезические, для которых α 1 кратно показаны светло-голубым цветом. (Геодезические показаны только для их первого прохождения вблизи антиподальной точки, а не для последующих.) Некоторые геодезические окружности показаны зеленым цветом; они образуют выступы на оболочке. Геометрическое место точек разреза показано красным цветом. Оболочка является геометрическим местом точек, которые сопряжены с A ; точки на оболочке могут быть вычислены путем нахождения точки, в которой m 12 = 0 на геодезической. Якоби (1891) называет эту звездообразную фигуру, созданную оболочкой, астроидой .

Снаружи астроиды две геодезические пересекаются в каждой точке; таким образом, между A и этими точками есть две геодезические (длиной примерно в половину окружности эллипсоида) . Это соответствует ситуации на сфере, где есть «короткие» и «длинные» пути на большом круге между двумя точками. Внутри астроиды четыре геодезические пересекаются в каждой точке. Четыре такие геодезические показаны на рис. 16, где геодезические пронумерованы в порядке возрастания длины. (Этот рисунок использует то же положение для A , что и рис. 13, и нарисован в той же проекции.) Две более короткие геодезические устойчивы , т. е. m 12 > 0 , так что нет соседнего пути, соединяющего две точки, который короче; другие два неустойчивы. Только самая короткая линия (первая) имеет σ 12π . Все геодезические касаются оболочки, которая показана зеленым на рисунке.

Астроида является (внешней) эволютой геодезических окружностей с центром в точке A. Аналогично, геодезические окружности являются эволютами астроиды.

Площадь геодезического многоугольника

Геодезический многоугольник — это многоугольник, стороны которого являются геодезическими. Он аналогичен сферическому многоугольнику , стороны которого являются большими окружностями. Площадь такого многоугольника можно найти, сначала вычислив площадь между геодезическим сегментом и экватором, т. е. площадь четырехугольника AFHB на рис. 1 (Danielsen 1989). Как только эта площадь известна, площадь многоугольника можно вычислить, суммируя вклады всех ребер многоугольника.

Здесь выражение для площади S 12 AFHB разработано в соответствии с Sjöberg (2006). Площадь любой замкнутой области эллипсоида равна

где dT — элемент площади поверхности, а Kгауссова кривизна . Теперь теорема Гаусса–Бонне, примененная к геодезическому многоугольнику, утверждает

где

— геодезический избыток, а θ j — внешний угол при вершине j . Умножая уравнение для Γ на R 2 2 , где R 2аутентичный радиус , и вычитая это из уравнения для T, получаем

где было подставлено значение K для эллипсоида . Применяя эту формулу к четырехугольнику AFHB , отмечая, что Γ = α 2 − α 1 , и выполняя интеграл по φ, получаем

где интеграл берется по геодезической линии (так что φ неявно является функцией λ ). Интеграл можно выразить в виде ряда, справедливого для малых f (Danielsen 1989) (Karney 2013, §6 и приложение).

Площадь геодезического многоугольника определяется суммированием S 12 по его ребрам. Этот результат справедлив при условии, что многоугольник не включает полюс; если он включает, то к сумме необходимо добавить 2 π R 2 2. Если ребра заданы их вершинами, то удобное выражение для геодезического избытка E 12 = α 2 − α 1 имеет вид

Решение прямой и обратной задачи

Решение геодезических задач влечет за собой отображение геодезической на вспомогательную сферу и решение соответствующей задачи в навигации по большому кругу . При решении «элементарного» сферического треугольника для NEP на рис. 5 можно использовать правила Непера для квадрантных треугольников ,

Картографирование геодезической включает оценку интегралов для расстояния s и долготы λ , уравнения (3) и (4) , которые зависят от параметра α 0 .

Решение прямой задачи не вызывает затруднений, поскольку α 0 можно определить непосредственно из заданных величин φ 1 и α 1 ; пример расчета см. в работе Карни (2013).

В случае обратной задачи λ 12 задано; это не может быть легко связано с эквивалентным сферическим углом ω 12, поскольку α 0 неизвестен. Таким образом, решение задачи требует, чтобы α 0 находили итеративно ( нахождение корня ); подробности см. в Karney (2013).

В геодезических приложениях, где f мало, интегралы обычно вычисляются как ряды (Legendre 1806) (Oriani 1806) (Bessel 1825) (Helmert 1880) (Rainsford 1955) (Rapp 1993). Для произвольного f интегралы (3) и (4) можно найти с помощью числовой квадратуры или выразив их через эллиптические интегралы (Legendre 1806) (Cayley 1870) (Karney 2024).

Винсенти (1975) предлагает решения для прямой и обратной задач; они основаны на разложении в ряд, выполненном до третьего порядка при выравнивании, и обеспечивают точность около0,1 мм для эллипсоида WGS84 ; однако обратный метод не сходится для почти антиподных точек.

Карни (2013) продолжает разложения до шестого порядка, что достаточно для обеспечения полной двойной точности для | f | ≤ 150 и улучшает решение обратной задачи, так что оно сходится во всех случаях. Карни (2013, приложение) расширяет метод для использования эллиптических интегралов, которые можно применять к эллипсоидам с произвольным сплющиванием.

Геодезические на трехосном эллипсоиде

Решение геодезической задачи для эллипсоида вращения математически просто: из-за симметрии геодезические имеют постоянную движения , заданную соотношением Клеро, что позволяет свести задачу к квадратуре . К началу 19 века (благодаря работам Лежандра , Ориани , Бесселя и др.) было достигнуто полное понимание свойств геодезических на эллипсоиде вращения.

С другой стороны, геодезические на трехосном эллипсоиде (с тремя неравными осями) не имеют очевидной константы движения и, таким образом, представляли собой сложную нерешенную проблему в первой половине 19 века. В замечательной статье Якоби (1839) открыл константу движения, позволяющую свести эту задачу к квадратуре (Клингенберг 1982, §3.5). [4]

Система координат трехосного эллипсоида

Рис. 17. Трехосные эллипсоидальные координаты.

Рассмотрим эллипсоид, определяемый формулой

где ( X , Y , Z ) — декартовы координаты с центром на эллипсоиде и, без потери общности, abc > 0. [ 5]Якоби (1866, §§26–27) использовал (триаксиальные) эллипсоидальные координатытриаксиальной эллипсоидальной широтой и триаксиальной эллипсоидальной долготой , β , ω ), определяемые как

В пределе ba , β становится параметрической широтой для сплющенного эллипсоида, поэтому использование символа β согласуется с предыдущими разделами. Однако ω отличается от сферической долготы, определенной выше. [ 6]

Линии сетки постоянной β (синим цветом) и ω (зеленым цветом) показаны на рис. 17. Они образуют ортогональную систему координат: линии сетки пересекаются под прямым углом. Главные сечения эллипсоида, определяемые X = 0 и Z = 0, показаны красным цветом. Третье главное сечение, Y = 0 , покрыто линиями β = ±90° и ω = 0° или ±180° . Эти линии встречаются в четырех точках умбилики (две из которых видны на этом рисунке), где главные радиусы кривизны равны. Здесь и на других рисунках в этом разделе параметры эллипсоида равны a : b : c = 1,01:1:0,8 , и он рассматривается в ортогональной проекции из точки выше φ = 40° , λ = 30° .

Линии сетки эллипсоидальных координат можно интерпретировать тремя различными способами:

  1. Они являются «линиями кривизны» на эллипсоиде: они параллельны направлениям главной кривизны (Монж, 1796).
  2. Они также являются пересечениями эллипсоида с софокусными системами однополостных и двуполостных гиперболоидов (Дюпен 1813, часть 5).
  3. Наконец, это геодезические эллипсы и гиперболы, определяемые с помощью двух соседних точек пупка (Hilbert & Cohn-Vossen 1952, стр. 188). Например, линии постоянной β на рис. 17 могут быть получены с помощью знакомой конструкции струны для эллипсов с концами струны, прикрепленными к двум точкам пупка.

Решение Якоби

Якоби показал, что геодезические уравнения, выраженные в эллипсоидальных координатах, являются разделимыми. Вот как он рассказал о своем открытии своему другу и соседу Бесселю (Якоби 1839, Письмо Бесселю),

Позавчера я свел к квадратуре задачу геодезических линий на эллипсоиде с тремя неравными осями . Это простейшие формулы в мире, абелевы интегралы , которые становятся известными эллиптическими интегралами, если две оси сделать равными.

Кёнигсберг , 28 декабря 1938 года.

Решение, данное Якоби (Якоби 1839) (Якоби 1866, §28), следующее:

Как отмечает Якоби, «функция угла β равна функции угла ω . Эти две функции являются просто абелевыми интегралами...» В решении появляются две константы δ и γ . Обычно δ равно нулю, если нижние пределы интегралов принимаются за начальную точку геодезической, а направление геодезических определяется γ . Однако для геодезических, начинающихся в точках пупочной области, мы имеем γ = 0 , а δ определяет направление в точке пупочной области. Константа γ может быть выражена как

где α — угол, который геодезическая образует с линиями постоянной ω . В пределе ba это сводится к sinα cos β = const. , известному соотношению Клеро. Вывод результата Якоби дан Дарбу (1894, §§583–584); он приводит решение, найденное Лиувиллем (1846) для общих квадратичных поверхностей.

Обследование трехосных геодезических

Циркумполярные геодезические, ω 1 = , α 1 =90° .

На трехосном эллипсоиде есть только три простые замкнутые геодезические, три главных сечения эллипсоида, заданные как X = 0 , Y = 0 и Z = 0 . [7] Для обзора других геодезических удобно рассматривать геодезические, которые пересекают среднее главное сечение, Y = 0 , под прямым углом. Такие геодезические показаны на рис. 18–22, которые используют те же параметры эллипсоида и то же направление взгляда, что и рис. 17. Кроме того, три главных эллипса показаны красным на каждом из этих рисунков.

Если начальная точка β 1 ∈ (−90°, 90°) , ω 1 = 0 и α 1 = 90° , то γ > 0 и геодезическая охватывает эллипсоид в «циркумполярном» смысле. Геодезическая колеблется к северу и югу от экватора; при каждом колебании она совершает немного меньший полный оборот вокруг эллипсоида, в результате чего в типичном случае геодезическая заполняет область, ограниченную двумя линиями широты β = ± β 1 . На рис. 18 и 19 приведены два примера. Рисунок 18 показывает практически то же поведение, что и для сплющенного эллипсоида вращения (потому что ab ); сравните с рис. 9. Однако, если начальная точка находится на более высокой широте (рис. 18), искажения, возникающие из-за ab, очевидны. Все касательные к циркумполярной геодезической касаются конфокального однополостного гиперболоида, который пересекает эллипсоид в точке β = β 1 (Chasles 1846) (Hilbert & Cohn-Vossen 1952, стр. 223–224).

Трансполярная геодезическая, β 1 =90° , α 1 =180° .

Если начальная точка β 1 = 90° , ω 1 ∈ (0°, 180°) и α 1 = 180° , то γ < 0 и геодезическая охватывает эллипсоид в «трансполярном» смысле. Геодезическая колеблется на восток и запад от эллипса X = 0 ; при каждом колебании она совершает немного больше полного оборота вокруг эллипсоида. В типичном случае это приводит к тому, что геодезическая заполняет область, ограниченную двумя линиями долготы ω = ω 1 и ω = 180° − ω 1 . Если a = b , все меридианы являются геодезическими; эффект ab заставляет такие геодезические колебаться на восток и запад. Два примера приведены на рис. 20 и 21. Сужение геодезической вблизи полюса исчезает в пределе bc ; в этом случае эллипсоид становится вытянутым эллипсоидом и рис. 20 будет напоминать рис. 10 (повернутый на бок). Все касательные к трансполярной геодезической касаются конфокального двухполостного гиперболоида, который пересекает эллипсоид в точке ω = ω 1 .

На рис. 18–21 геодезические (почти) замкнуты. Как отмечено выше, в типичном случае геодезические не замкнуты, а заполняют область, ограниченную предельными линиями широты (в случае рис. 18–19) или долготы (в случае рис. 20–21).

Рис. 22. Геодезическая омбиликальная линия, β 1 =90° , ω 1 = , α 1 =135° .

Если начальная точка β 1 = 90° , ω 1 = 0° (точка пупка) и α 1 = 135° (геодезическая линия покидает эллипс Y = 0 под прямым углом), то γ = 0 и геодезическая линия многократно пересекает противоположную точку пупка и возвращается в свою начальную точку. Однако на каждом витке угол, под которым она пересекает Y = 0, становится ближе к или180° , так что асимптотически геодезическая лежит на эллипсе Y = 0 (Hart 1849) (Arnold 1989, стр. 265), как показано на рис. 22. Одна геодезическая не заполняет область на эллипсоиде. Все касательные к омбилическим геодезическим касаются конфокальной гиперболы, которая пересекает эллипсоид в омбилических точках.

Геодезические пуповины обладают несколькими интересными свойствами.

Если начальная точка A геодезической не является точкой пупка, то ее огибающая является астроидой с двумя точками возврата, лежащими на β = − β 1 , а двумя другими на ω = ω 1 + π . Точкой разреза для A является часть линии β = − β 1 между точками возврата.

Приложения

Прямые и обратные геодезические задачи больше не играют той центральной роли в геодезии, как когда-то. Вместо решения уравнивания геодезических сетей как двумерной задачи сфероидальной тригонометрии, эти задачи теперь решаются трехмерными методами (Vincenty & Bowring 1978). Тем не менее, наземная геодезическая система по-прежнему играет важную роль в нескольких областях:

По принципу наименьшего действия многие проблемы в физике можно сформулировать как вариационную задачу, похожую на задачу для геодезических. Действительно, геодезическая задача эквивалентна движению частицы, вынужденной двигаться по поверхности, но в остальном не подверженной никаким силам (Лаплас 1799a) (Гильберт и Кон-Фоссен 1952, стр. 222). По этой причине геодезические на простых поверхностях, таких как эллипсоиды вращения или трехосные эллипсоиды, часто используются в качестве «тестовых случаев» для исследования новых методов. Вот некоторые примеры:

Смотрите также

Примечания

  1. ^ Здесь α 2прямой азимут в точке B. Некоторые авторы вместо этого вычисляют обратный азимут; он задается как α 2 ± π .
  2. ^ Лаплас (1799a) показал, что частица, вынужденная двигаться по поверхности, но не подверженная никаким другим силам, движется вдоль геодезической линии этой поверхности. Таким образом, соотношение Клеро является просто следствием сохранения момента импульса для частицы на поверхности вращения.
  3. ^ Багратуни (1962, §17) использует термин «коэффициент сходимости ординат» для геодезической шкалы.
  4. ^ Этот раздел адаптирован из документации GeographicLib (Karney 2015, Geodesics on a triaxial ellipsoid)
  5. ^ Это обозначение для полуосей несовместимо с тем, которое использовалось в предыдущем разделе об эллипсоидах вращения, в котором a и b обозначали экваториальный радиус и полярную полуось. Таким образом, соответствующие неравенства имеют вид a = ab > 0 для сплюснутого эллипсоида и ba = a > 0 для вытянутого эллипсоида.
  6. ^ Предел bc дает вытянутый эллипсоид, где ω играет роль параметрической широты.
  7. ^ Если ca < 12 , существуют другие простые замкнутые геодезические, подобные показанным на рис. 11 и 12 (Клингенберг 1982, §3.5.19).

Ссылки

Внешние ссылки