Задачи на сфере: угловая засечка

Эта страница является черновиком статьи.

Линейная засечка — это нахождение положения точки по координатам двух исходных пунктов и значениям азимутов направлений с этих пунктов на определяемую точку.

Общие положения

В качестве модели Земли принимается сфера с радиусом R, равным среднему радиусу земного эллипсоида. Аналогом прямой линии на плоскости является геодезическая линия на поверхности. На сфере геодезическая линия — дуга большого круга.

Введём следующие обозначения:

• φ — географическая широта,
• λ — географическая долгота,
• α — азимут дуги большого круга,
• σ — сферическое расстояние (длина дуги большого круга, выраженная в долях радиуса шара).

Линейное расстояние по дуге большого круга s связано со сферическим расстоянием σ формулой s = R σ.

Постановка задачи

Исходные данные

координаты пунктов Q₁, Q₂ — φ₁, λ₁, φ₂, λ₂,

начальные направления с пунктов Q₁, Q₂ на точку Q₃ — α₁₃, α₂₃.

Определяемые величины

координаты точки Q₃ — φ₃, λ₃.

Алгоритм

Файл:Sph ang.png

Угловая засечка

Решение любого вида засечек сводится к нахождению полярных координат искомой точки, т.е. начального направления и расстояния на неё с одного из исходных пунктов. На конечном этапе координаты находятся из решения прямой геодезической задачи. Поскольку в угловой засечке направления α₁₃ и α₂₃ уже заданы, остаётся определить расстояние σ₁₃ или σ₂₃.

На рисунке синим цветом выделены заданные элементы сферического треугольника, красным цветом неизвестные, зелёным — вспомогательные элементы. Очевидно, в треугольнике Q₁Q₂Q₃ нет ни одного известного элемента. Однако из решения обратной геодезической задачи для пунктов Q₁, Q₂ могут быть получены расстояние σ₁₂, а также азимуты α₁₂ и α₂₁, после чего углы β₁ и β₂ вычисляются как разности азимутов при соответствующих пунктах. Далее из решения треугольника Q₁Q₂Q₃ найдём сторону σ₁₃.

Последовательность действий:

1. решить обратную геодезическую задачу для Q₁, Q₂: по φ₁, λ₁, φ₂, λ₂ получить α₁₂, α₂₁, σ₁₂;
2. вычислить углы β₁, β₂;
3. в треугольнике Q₁Q₂Q₃ по σ₁₂, β₁, β₂ вычислить σ₁₃;
4. решить прямую геодезическую задачу для Q₁, Q₃: по φ₁, λ₁, α₁₃, σ₁₃ вычислить φ₃, λ₃.

Действия по первому и последнему пунктам рассмотрены в статьях Задачи на сфере: обратная геодезическая задача и Задачи на сфере: прямая геодезическая задача.

Углы β₁, β₂ и длина σ₁₃ вычисляются по формулам:

${\begin{array}{rcl}\beta _{1}&=&\alpha _{12}-\alpha _{13}\\\beta _{2}&=&\alpha _{23}-\alpha _{21}\\\sigma _{13}&=&\operatorname {arctg\,} {\dfrac {\sin \beta _{2}\sin \sigma _{12}}{\cos \beta _{2}\sin \beta _{1}+\sin \beta _{2}\cos \beta _{1}\cos \sigma _{12}}}\end{array}}$

Правда, до вычисления длины σ₁₃ необходимо проанализировать полученные значения углов β₁ и β₂. Ниже в коде функции можно увидеть пример такого анализа:

• если линии Q₁Q₃ и Q₂Q₃ совпадают с Q₁Q₂, решение не определено, т.к. решением может быть любая точка геодезической линии Q₁Q₂;
• если одна из линий Q₁Q₃ и Q₂Q₃ совпадает с Q₁Q₂, а другая нет, решением является пункт, из которого выходит другая;
• если линии Q₁Q₃ и Q₂Q₃ уходят в разные полушария от Q₁Q₂, функция находит ближайшее к Q₁Q₂ «ложное пересечение» этих линий.

Здесь необходимо пояснить, что на сфере две несовпадающие геодезические линии всегда пересекаются в двух точках, которые являются антиподами. В традиционной постановке направление на нужное пересечение задаётся явно. Если же задача ставится так, что прямое и обратное направления равнозначны, возникает вопрос выбора одного из антиподов: φ₃, λ₃ или φ₃′ = −φ₃, λ₃′ = λ₃ ± 180°.

Ссылки

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