Расстояние большого круга - Great-circle distance

Эта статья про кратчайшее расстояние на сфере. Кратчайшее расстояние на эллипсоиде см. геодезические на эллипсоиде.

Диаграмма, показывающая расстояние по большому кругу (нарисованное красным) между двумя точками на сфере, P и Q. Также показаны две противоположные точки, u и v.

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

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

Через любые две точки на сфере, которые не прямо напротив друг друга, есть уникальный большой круг. Две точки разделяют большой круг на две дуги. Длина более короткой дуги - это расстояние по большому кругу между точками. Большой круг, наделенный таким расстоянием, называется Риманов круг в Риманова геометрия.

Между двумя точками, которые находятся прямо напротив друг друга, называются противоположные точки, существует бесконечно много больших окружностей, и все дуги больших окружностей между противоположными точками имеют длину, равную половине длины длина окружности круга, или ${\displaystyle \pi r}$, где р это радиус сферы.

В Земля почти сферический (см. Радиус Земли ), поэтому формулы расстояния по дуге большого круга дают расстояние между точками на поверхности Земли с точностью до 0,5%.^[1] (Увидеть Длина дуги § Дуги больших кругов на Земле.)

Формулы

Иллюстрация центрального угла Δσ между двумя точками P и Q. λ и φ - это продольный и широтный углы P соответственно.

Позволять ${\displaystyle \lambda _{1},\phi _{1}}$ и ${\displaystyle \lambda _{2},\phi _{2}}$ быть географическим долгота и широта в радианах двух точек 1 и 2, и ${\displaystyle \Delta \lambda ,\Delta \phi }$ быть их абсолютными различиями; тогда ${\displaystyle \Delta \sigma }$, то центральный угол между ними, дается сферический закон косинусов если один из полюсов используется как вспомогательная третья точка на сфере:^[2]

${\displaystyle \Delta \sigma =\arccos {\bigl (}\sin \phi _{1}\sin \phi _{2}+\cos \phi _{1}\cos \phi _{2}\cos(\Delta \lambda ){\bigr )}.}$

Проблема обычно выражается в определении центрального угла ${\displaystyle \Delta \sigma }$. Учитывая этот угол в радианах, фактическая длина дуги d на сфере радиуса р можно тривиально вычислить как

${\displaystyle d=r\,\Delta \sigma .}$

Расчетные формулы

В компьютерных системах с низким точность с плавающей запятой, формула сферического закона косинусов может иметь большие ошибки округления если расстояние небольшое (если две точки находятся на расстоянии километра друг от друга на поверхности Земли, косинус центрального угла составляет около 0,99999999). Для современных 64-битные числа с плавающей запятой формула сферического закона косинусов, приведенная выше, не имеет серьезных ошибок округления для расстояний, превышающих несколько метров на поверхности Земли.^[3] В формула гаверсина является численно лучше оснащенный на небольшие расстояния:^[4]

${\displaystyle {\begin{aligned}\Delta \sigma &=\operatorname {archav} \left(\operatorname {hav} \left(\Delta \phi \right)+\cos \phi _{1}\cos \phi _{2}\operatorname {hav} \left(\Delta \lambda \right)\right)\\\\\Delta \sigma &=2\arcsin {\sqrt {\sin ^{2}\left({\frac {\Delta \phi }{2}}\right)+\cos {\phi _{1}}\cos {\phi _{2}}\sin ^{2}\left({\frac {\Delta \lambda }{2}}\right)}}.\\\end{aligned}}}$

Исторически использование этой формулы было упрощено наличием таблиц для гаверсин функция: hav (θ) = грех²(θ/2).

Хотя эта формула точна для большинства расстояний на сфере, она также страдает ошибками округления для специального (и несколько необычного) случая противоположные точки (на противоположных концах сферы). Формула, точная для всех расстояний, является следующим частным случаем Формула Винсенти для эллипсоида с равными большой и малой осями:^[5]

${\displaystyle \Delta \sigma =\arctan {\frac {\sqrt {\left(\cos \phi _{2}\sin(\Delta \lambda )\right)^{2}+\left(\cos \phi _{1}\sin \phi _{2}-\sin \phi _{1}\cos \phi _{2}\cos(\Delta \lambda )\right)^{2}}}{\sin \phi _{1}\sin \phi _{2}+\cos \phi _{1}\cos \phi _{2}\cos(\Delta \lambda )}}.}$

Векторная версия

Другое представление аналогичных формул, но с использованием нормальные векторы вместо широты и долготы для описания позиций находится с помощью 3D векторная алгебра, с использованием скалярное произведение, перекрестное произведение, или комбинация:^[6]

${\displaystyle {\begin{aligned}\Delta \sigma &=\arccos \left(\mathbf {n} _{1}\cdot \mathbf {n} _{2}\right)\\&=\arcsin \left|\mathbf {n} _{1}\times \mathbf {n} _{2}\right|\\&=\arctan {\frac {\left|\mathbf {n} _{1}\times \mathbf {n} _{2}\right|}{\mathbf {n} _{1}\cdot \mathbf {n} _{2}}}\\\end{aligned}}}$

где ${\displaystyle \mathbf {n} _{1}}$ и ${\displaystyle \mathbf {n} _{2}}$ являются нормалями к эллипсоиду в двух положениях 1 и 2. Подобно приведенным выше уравнениям, основанным на широте и долготе, выражение, основанное на arctan, является единственным, которое хорошо обусловлено для всех углов. Выражение на основе arctan требует величины перекрестного произведения на скалярное произведение.

От длины хорды

Линия через трехмерное пространство между интересными точками на сферической Земле - это аккорд большого круга между точками. В центральный угол между двумя точками можно определить по длине хорды. Расстояние большого круга пропорционально центральному углу.

Длина хорды большого круга, ${\displaystyle C_{h}\,\!}$, может быть рассчитана следующим образом для соответствующей единичной сферы с помощью Декартово вычитание:

${\displaystyle {\begin{aligned}\Delta {X}&=\cos \phi _{2}\cos \lambda _{2}-\cos \phi _{1}\cos \lambda _{1};\\\Delta {Y}&=\cos \phi _{2}\sin \lambda _{2}-\cos \phi _{1}\sin \lambda _{1};\\\Delta {Z}&=\sin \phi _{2}-\sin \phi _{1};\\C&={\sqrt {(\Delta {X})^{2}+(\Delta {Y})^{2}+(\Delta {Z})^{2}}}\end{aligned}}}$

Центральный угол равен:

${\displaystyle \Delta \sigma =2\arcsin {\frac {C}{2}}.}$

Радиус сферической Земли

Экваториальный (а), полярный (б) и среднего радиуса Земли, как определено в 1984 г. Мировая геодезическая система доработка. (Не в масштабе.)

Основная статья: Радиус Земли

В форма Земли напоминает сплюснутую сферу ( сфероид ) с экваториальным радиусом ${\displaystyle a}$ 6378,137 км; расстояние ${\displaystyle b}$ от центра сфероида до каждого полюса 6356,7523142 км. При вычислении длины короткой линии север-юг на экваторе окружность, которая лучше всего приближается к этой линии, имеет радиус ${\displaystyle b^{2}/a}$ (что равно меридиану полу-латусная прямая кишка ), или 6335,439 км, а сфероид на полюсах лучше всего аппроксимируется сферой радиуса ${\displaystyle a^{2}/b}$, или 6399,594 км, разница в 1%. Пока предполагается сферическая Земля, любая формула для расстояния на Земле гарантированно верна только в пределах 0,5% (хотя возможна более высокая точность, если формула предназначена только для применения к ограниченной области). С использованием средний радиус Земли, ${\displaystyle R_{1}={\frac {1}{3}}(2a+b)\approx 6371.009\,\mathrm {km} }$ (для WGS84 эллипсоид) означает, что в пределе малого уплощения средний квадрат относительная ошибка в оценках расстояние сведено к минимуму.^[7]

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

• Аэронавигация
• Изоазимутал
• Центральный угол
• Кругосветное плавание
• Планирование полета
• Геодезия
• Геодезические на эллипсоиде
• Геодезическая система
• Географическое расстояние
• Навигация по большому кругу
• Формула гаверсина
• Дуга меридиана
• Линия румба
• Сферическая земля
• Сферическая геометрия
• Сферическая тригонометрия

Ссылки и примечания

1. Адмиралтейское руководство мореплавания, Том 1, Канцелярия, 1987, стр. 10, ISBN  9780117728806, Ошибки, вносимые принятием сферической Земли на основе международной морской мили, составляют не более 0,5% для широты и 0,2% для долготы.
2. Kells, Lyman M .; Керн, Уиллис Ф .; Блэнд, Джеймс Р. (1940). Плоская и сферическая тригонометрия. McGraw Hill Book Company, Inc., стр.323 -326. Получено 13 июля, 2018.
3. "Вычислить расстояние, азимут и многое другое между точками широты и долготы". Получено 10 августа 2013.
4. Синнотт, Роджер В. (август 1984). «Добродетели гаверсов». Небо и телескоп. 68 (2): 159.
5. Винсенти, Фаддей (1975-04-01). «Прямые и обратные решения геодезических на эллипсоиде с применением вложенных уравнений» (PDF ). Обзор обзора. Кингстон-роуд, Толворт, Суррей: Управление зарубежных исследований. 23 (176): 88–93. Дои:10.1179 / sre.1975.23.176.88. Получено 2008-07-21.
6. Гейд, Кеннет (2010). «Неособое представление горизонтального положения» (PDF). Журнал навигации. Издательство Кембриджского университета. 63 (3): 395–417. Дои:10.1017 / S0373463309990415.
7. Маккоу, Г. Т. (1932). «Длинные линии на Земле». Обзор обзора империи. 1 (6): 259–263. Дои:10.1179 / sre.1932.1.6.259.

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

• Большой круг в MathWorld