Геодезические на эллипсоиде - Geodesics on an ellipsoid

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

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

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

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

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

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

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

Можно свести различные геодезические задачи к одному из двух типов. Учтите два момента: А в широтаφ1 и долгота λ1 иB на широте φ2 и долготаλ2 (см. рис. 1). Соединительная геодезическая (от А к B) является AB, длиныs12, у которого есть азимуты α1 иα2 на двух конечных точках.[1] Обычно рассматриваются две геодезические задачи:

  1. то прямая геодезическая задача или же первая геодезическая задача, данный А, α1, и s12, определять B и α2;
  2. то обратная геодезическая задача или же вторая геодезическая задача, данный А и B, определять s12, α1, и α2.

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

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

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

Nous désignerons cette ligne sous le nom de Ligne Géodésique [Мы будем называть эту линию геодезическая линия].

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

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

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

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

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

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

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

Рассмотрим эллипсоид вращения с экваториальным радиусома и полярная полуось б. Определить сглаживание ж = (аб)/а, эксцентриситете = а2б2/а = ж(2 − ж), а вторая децентриситет е′ = а2б2/б = е/(1 − ж). (В большинстве случаев в геодезии эллипсоид считается сжатым,а > б; однако теория без изменений применима к вытянутым эллипсоидам, а < б, в таком случае ж, е2, и е2 отрицательны.)

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

(1)

куда ρ этомеридиональный радиус кривизны,р = ν cosφ это радиус круга широтыφ, и ν этонормальный радиус кривизны Таким образом, элементарный сегмент имеет вид

или же

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

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

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

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

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

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

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

Рис. 4. Геодезическая задача с привязкой к вспомогательной сфере.
Рис. 5. Элементарная геодезическая задача на вспомогательной сфере.

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

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

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

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

(2)

Комбинируя уравнения. (1) и (2) дают дифференциальные уравнения для s и λ

Связь между β и φ является

который дает

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

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

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

(3)

куда

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

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

(4)

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

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

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

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

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

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

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

Две дополнительные замкнутые геодезические для сжатого эллипсоида, ба = ​27.
Рис. 11. Вид сбоку.
Рис. 12. Вид сверху.

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

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

На Рис.13 показаны геодезические (синим цветом), исходящиеА с α1 кратный15° вплоть до того момента, когда они перестают быть кратчайшими путями. (Сглаживание увеличено до110 чтобы подчеркнуть эллипсоидальные эффекты.) Также показаны (зеленым цветом) кривые постоянной s12, которые представляют собой геодезические окружности с центром А.Гаусс (1828) показал, что на любой поверхности геодезические и геодезическая окружность пересекаются под прямым углом. Красная линия - этовырезать место, геометрическое место точек, которые имеют кратчайшие (в данном случае две) кратчайшие геодезические из А. На сфере разрез - это точка. На сплющенном эллипсоиде (показанном здесь) это сегмент круга широты с центром в точке противоположный к А, φ = −φ1. Продольная протяженность локуса разреза составляет приблизительноλ12 ∈ [π - ж π cosφ1, π + ж π cosφ1]. ЕслиА лежит на экваторе, φ1 = 0, это соотношение является точным, и, как следствие, экватор является кратчайшей геодезической, только если| λ12| ≤ (1 − ж) π. Для вытянутого эллипсоида место разреза - это сегмент антимеридиана с центром в точке, противоположной точке. А, λ12 = π, а это означает, что меридиональные геодезические перестают быть кратчайшими путями до того, как будет достигнута точка противоположности.

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

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

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

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

куда

Количество м(s1, s2) = м12 так называемыйуменьшенная длина, и M(s1, s2) = M12 этогеодезическая шкала.[3]Их основные определения проиллюстрированы на рис. 14.

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

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

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

  • для сплющенного эллипсоида, | σ12| ≤ π;
  • для вытянутого эллипсоида, | λ12| ≤ π, если α0 ≠ 0; если α0 = 0, дополнительное условие м12 ≥ 0 требуется, если | λ12| = π.

Конверт геодезических

Геодезические из одной точки (ж = ​110, φ1 = −30°)
Рис. 15. Огибающая геодезических от точки. А в φ1 = −30°.
Рис. 16. Четыре геодезические, соединяющие А и точка B, φ2 = 26°, λ12 = 175°.

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

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

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

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

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

Здесь выражение для площади S12 из AFHBразработан после Сьёберг (2006). Площадь любой замкнутой области эллипсоида равна

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

куда

- геодезический эксцесс и θj внешний угол при вершине j. Умножая уравнение на Γк р22, куда р2 этоаутентичный радиус, и вычитая это из уравнения для Т дает

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

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

Площадь геодезического многоугольника определяется суммированием S12по его краям. Этот результат верен при условии, что многоугольник не включает полюс; если это так, р22 необходимо добавить к сумме. Если ребра задаются своими вершинами, тоудобное выражение за геодезический избыток E12 = α2 - α1 является

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

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

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

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

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

В геодезических приложениях, где ж мала, интегралы обычно вычисляются как ряд (Legendre 1806 г. )(Ориани 1806 ) (Бессель 1825 ) (Гельмерт 1880 г. )(Рейнсфорд 1955 ) (Рэпп 1993 ). Для произвольныхж, интегралы (3) и (4) могут быть найдены с помощью числовой квадратуры или выражением их черезэллиптические интегралы (Legendre 1806 г. ) (Кэли 1870 ).

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

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

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

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

Трехосная система координат

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

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

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

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

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

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

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

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

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

Кенигсберг, 28 декабря '38.

Решение, данное Якоби (Якоби 1839 )(Якоби 1866, §28) является

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

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

Съемка трехосных геодезических

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

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

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

Заполярные геодезические, β1 = 90°, α1 = 180°.
Рис.20. ω1 = 39.9°.
Рис.21. ω1 = 9.966°.

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

Рис. 22. Пупочная геодезическая, β1 = 90°, ω1 = , α1 = 135°.

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

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

  • Через любую точку эллипсоида проходят две омбилические геодезические.
  • Геодезическое расстояние между противоположными точками шлангокабеля одинаково независимо от начального направления геодезической.
  • Тогда как замкнутые геодезические на эллипсах Икс = 0 и Z = 0 устойчивы (геодезическая, изначально близкая к эллипсу и почти параллельная ему, остается близкой к эллипсу), замкнутая геодезическая на эллипсе Y = 0, который проходит через все 4 точки пуповины, экспоненциально нестабильный. Если его потревожить, он вылетит из плоскости Y = 0 и перевернитесь, прежде чем вернуться, чтобы приблизиться к самолету. (Это поведение может повторяться в зависимости от природы начального возмущения.)

Если отправная точка А геодезической не является точкой пупка, ее оболочка - астроида с двумя вершинами, лежащими наβ = −β1 и два других наω = ω1 + π. Место разреза А это часть строки β = −β1 между бугорками.

Приложения

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

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

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

Примечания

  1. ^ Здесь α2 это вперед азимут на BНекоторые авторы рассчитывают назад азимут вместо; это даетсяα2 ± π.
  2. ^ Лаплас (1799a) showed that a particle constrained to move ona surface but otherwise subject to no forces moves along a geodesic forthat surface. Thus, Clairaut's relation is just a consequence ofсохранение углового момента for a particle on a surface ofrevolution.
  3. ^ Bagratuni (1962, §17) uses the term "coefficient ofconvergence of ordinates" for the geodesic scale.
  4. ^ This section is adapted from the documentation for GeographicLib(Karney 2015, Геодезические на трехосном эллипсоиде )
  5. ^ This notation for the semi-axes is incompatible with that used in theprevious section on ellipsoids of revolution in which а иб stood for the equatorial radius and polar semi-axis.Thus the corresponding inequalities are а = аб > 0 foran oblate ellipsoid and ба = а > 0 for a prolateellipsoid.
  6. ^ Лимит бc gives a prolate ellipsoid withω playing the role of the parametric latitude.
  7. ^ Если cа < ​12, there are other simple closed geodesicssimilar to those shown in Figs. 11 and 12(Klingenberg 1982, §3.5.19).

Рекомендации

внешняя ссылка