Меридианная дуга
В геодезии и навигации дуга меридиана — это кривая между двумя точками на поверхности Земли, имеющими одинаковую долготу . Этот термин может относиться либо к сегменту меридиана . либо к его длине ,
Целью измерения дуг меридианов является определение фигуры Земли .Одно или несколько измерений дуг меридианов можно использовать для определения формы опорного эллипсоида , который лучше всего аппроксимирует геоид в области измерений. Измерения дуг меридианов на нескольких широтах вдоль многих меридианов по всему миру можно объединить, чтобы приблизить геоцентрический эллипсоид, предназначенный для всего мира.
Самые ранние определения размера сферической Земли требовали одной дуги. Точные геодезические работы, начавшиеся в 19 веке, потребовали нескольких дуговых измерений в регионе, где должна была проводиться съемка, что привело к распространению опорных эллипсоидов по всему миру. В последних определениях используются астрогеодезические измерения и методы спутниковой геодезии для определения опорных эллипсоидов, особенно геоцентрических эллипсоидов, которые сейчас используются для глобальных систем координат, таких как WGS 84 (см. числовые выражения ).
История измерений
[ редактировать ]Ранние оценки размеров Земли зафиксированы в Греции в 4 веке до нашей эры и у ученых халифа Дома мудрости в Багдаде в 9 веке. Первое реалистическое значение было рассчитано александрийским ученым Эратосфеном около 240 г. до н. э. Он подсчитал, что длина меридиана составляет 252 000 стадий с погрешностью реального значения от -2,4% до +0,8% (при условии, что значение стадиона составляет от 155 до 160 метров). [1] Эратосфен описал свою технику в книге под названием «О мере Земли» , которая не сохранилась. Подобный метод использовал Посидоний примерно 150 лет спустя, а несколько лучшие результаты были рассчитаны в 827 г. методом дугового измерения , [2] приписывают халифу Аль-Мамуну . [ нужна ссылка ]
Эллипсоидальная Земля
[ редактировать ]В ранней литературе термин «сплюснутый сфероид» используется для описания сферы, «сжатой на полюсах». используется термин «эллипсоид вращения» В современной литературе вместо термина «сфероид» , хотя уточняющее слово «революция» обычно опускается. Эллипсоид , не являющийся эллипсоидом вращения, называется трехосным эллипсоидом. В этой статье сфероид и эллипсоид используются как взаимозаменяемые, при этом подразумевается сплюснутость, если не указано иное.
17 и 18 веков
[ редактировать ]было известно Хотя с классической античности , что Земля имеет сферическую форму , к 17 веку накопились доказательства того, что она не является идеальной сферой. В 1672 году Жан Рише нашел первое свидетельство того, что гравитация над Землей не является постоянной (как это было бы, если бы Земля была сферой); он взял маятниковые часы в Кайенну , Французская Гвиана , и обнаружил, что они потеряли 2 + 1 / 2 минуты в день по сравнению с ценой в Париже . [3] [4] Это указывало на то, что ускорение силы тяжести в Кайенне было меньше, чем в Париже. Маятниковые гравиметры начали брать с собой в путешествия в отдаленные части мира, и постепенно было обнаружено, что гравитация плавно увеличивается с увеличением широты , причем гравитационное ускорение примерно на 0,5% больше, на географических полюсах чем на экваторе .
В 1687 году Исаак Ньютон опубликовал в « Началах» доказательство того, что Земля представляет собой сплюснутый сфероид , , уплощенный равный 1 / 230 . [5] Это оспаривалось некоторыми, но не всеми французскими учёными. Дуга меридиана Жана Пикара была продлена до более длинной дуги Джованни Доменико Кассини и его сыном Жаком Кассини в период 1684–1718 годов. [6] Дуга была измерена как минимум с тремя определениями широты, поэтому они смогли определить среднюю кривизну северной и южной половин дуги, что позволило определить общую форму. Результаты показали, что Земля представляет собой вытянутый сфероид (с экваториальным радиусом меньше полярного радиуса). Для решения вопроса Французская академия наук (1735) предприняла экспедиции в Перу ( Бугер , Луи Годен , де ла Кондамин , Антонио де Ульоа , Хорхе Хуан ) и в Лапландию ( Мопертюи , Клеро , Камю , Ле Монье , аббат Утье , Андерс Цельсий ). Результаты измерений в экваториальных и полярных широтах подтвердили, что Землю лучше всего моделировать сплюснутым сфероидом, поддерживающим Ньютона. [6] Однако к 1743 году теорема Клеро полностью вытеснила подход Ньютона.
К концу века Жан Батист Жозеф Деламбр переизмерил и расширил французскую дугу от Дюнкерка до Средиземного моря ( дуга меридиана Деламбра и Мешена ). Он был разделен на пять частей четырьмя промежуточными определениями широты. Объединив измерения с измерениями дуги Перу,Были определены параметры формы эллипсоида, а расстояние между экватором и полюсом по парижскому меридиану рассчитано как 5 130 762 туаза в соответствии со стандартным туазом Парижа. Определение этого расстояния ровно в 10 000 000 м привело к построению новой стандартной метровой планки, равной 0,513 0762 туазов. [6] : 22
19 век
[ редактировать ]В XIX веке многие астрономы и геодезисты занимались детальными исследованиями кривизны Земли по различным дугам меридианов. В результате анализа было получено очень много моделей эллипсоидов, таких как Плесси 1817, Эйри 1830, Бессель 1841 , Эверест 1830 и Кларк 1866 . [7] Полный список эллипсоидов приведен в разделе « Земной эллипсоид» .
Морская миля
[ редактировать ]Исторически морская миля определялась как длина дуги в одну минуту вдоль меридиана сферической Земли. Модель эллипсоида приводит к изменению морской мили в зависимости от широты. Эту проблему удалось решить, определив морскую милю равной ровно 1852 метра. Однако для всех практических целей расстояния измеряются по широтной шкале карт. Как говорит Королевская яхтенная ассоциация в своем руководстве для шкиперов : «1 (минута) широты = 1 морской миле», за которой следует «Для большинства практических целей расстояние измеряется по шкале широты, предполагая, что одна минута широты равна одной морской миле». миля». [8]
Расчет
[ редактировать ]На сфере длина дуги меридиана — это просто длина дуги окружности .На эллипсоиде вращения для коротких дуг меридианов их длину можно аппроксимировать, используя меридиональный радиус кривизны Земли и формулировку дуги окружности.Для более длинных дуг длина получается из вычитания двух меридиональных расстояний , расстояния от экватора до точки на широте φ .Это важная проблема в теории картографических проекций, в частности поперечной проекции Меркатора .
Основными параметрами эллипсоида являются , b , f , но в теоретической работе полезно определить дополнительные параметры, в частности эксцентриситет e n и третье сплющивание a . Только два из этих параметров являются независимыми и между ними существует множество связей:
Определение
[ редактировать ]Можно показать, что радиус кривизны меридиана равен: [9] [10]
Длина дуги бесконечно малого элемента меридиана равна dm = M ( φ ) dφ (где φ в радианах). Следовательно, меридиональное расстояние от экватора до широты φ равно
Формула расстояния становится проще, если ее записать в терминах параметрическая широта ,
где tan β = (1 − f )tan φ и e ′ 2 = e 2 / 1 - е 2 .
Хотя широта обычно ограничивается диапазоном [− π / 2 , π / 2 ] , все приведенные здесь формулы применимы для измерения расстояния вокруг полного меридианного эллипса (включая антимеридиан). Таким образом, диапазоны φ , β и выпрямляющей широты µ не ограничены.
Связь с эллиптическими интегралами
[ редактировать ]Приведенный выше интеграл относится к частному случаю неполного эллиптического интеграла третьего рода . В обозначениях онлайн- NIST. справочника [11] ( Раздел 19.2(ii) ),
Его также можно записать в терминах неполных эллиптических интегралов второго рода справочника NIST (см. раздел 19.6(iv) ),
Вычисление (с произвольной точностью) эллиптических интегралов и аппроксимаций также обсуждается в справочнике NIST. Эти функции также реализованы в программах компьютерной алгебры, таких как Mathematica. [12] и Максима. [13]
Расширения серии
[ редактировать ]Вышеупомянутый интеграл можно выразить как бесконечный усеченный ряд, разложив подынтегральную функцию в ряд Тейлора, выполнив полученные интегралы почленно и выразив результат в виде тригонометрического ряда. В 1755 году Леонард Эйлер вывел разложение по квадрату третьего эксцентриситета . [14]
Разложения по эксцентриситету ( e )
[ редактировать ]Полуразрушенный в 1799 году. [15] вывел широко используемое разложение по e 2 ,
где
Ричард Рапп дает подробный вывод этого результата. [16]
Расширения в третьем уплощении ( n )
[ редактировать ]Ряды со значительно более быстрой сходимостью можно получить, разложив их по третьему уплощению n вместо эксцентриситета. Они связаны
В 1837 году Фридрих Бессель получил одну такую серию: [17] который был облечен в более простую форму Гельмертом , [18] [19]
с
Поскольку n меняет знак, когда a и b меняются местами, и поскольку начальный множитель 1 / 2 ( a + b ) при таком обмене постоянен, половина членов в разложениях H 2 k обращается в нуль.
Ряд можно выразить с помощью a или b в качестве начального множителя, записав, например,:
и разложим результат в ряд по n . Несмотря на то, что это приводит к более медленно сходящимся рядам, такие ряды используются в спецификации поперечной проекции Меркатора Национальным агентством геопространственной разведки. [20] и Служба боеприпасов Великобритании . [21]
Ряд по параметрической широте
[ редактировать ]В 1825 году Бессель [22] вывел разложение меридионального расстояния по параметрической широте β в связи со своими работами по геодезии ,
с
Поскольку этот ряд обеспечивает разложение эллиптического интеграла второго рода, его можно использовать для записи длины дуги через геодезическую широту как
Обобщенная серия
[ редактировать ]Вышеуказанные серии до восьмого порядка по эксцентриситету или четвертого порядка по третьему уплощению обеспечивают точность до миллиметра. С помощью систем символьной алгебры их можно легко расширить до шестого порядка при третьем уплощении, что обеспечивает полную точность двойной точности для наземных приложений.
ветхий [15] и Бессель [22] оба написали свои ряды в форме, позволяющей их обобщать в произвольном порядке. Коэффициенты ряда Бесселя можно выразить особенно просто:
где
и к !! — это двойной факториал , расширенный до отрицательных значений с помощью рекурсивного соотношения: (−1)!! = 1 и (−3)!! = −1 .
Коэффициенты в ряду Гельмерта аналогичным образом могут быть выражены в общем виде:
Этот результат был предположен Фридрихом Гельмертом. [23] и доказал Кадзусигэ Кавасэ. [24]
Дополнительный множитель (1 − 2 k )(1 + 2 k ) возникает в результате дополнительного расширения появляется в приведенной выше формуле и приводит к худшей сходимости ряда по φ по сравнению с рядом по β .
Числовые выражения
[ редактировать ]Приведенный выше тригонометрический ряд удобно оценить с помощью суммирования Кленшоу . Этот метод позволяет избежать вычисления большинства тригонометрических функций и позволяет быстро и точно суммировать ряды. Этот метод также можно использовать для оценки разницы m ( φ 1 ) − m ( φ 2 ) при сохранении высокой относительной точности.
Замена значений большой полуоси и эксцентриситета эллипсоида WGS84 дает
где φ ( ° ) = φ / 1° — это φ, выраженная в градусах (и аналогично для β ( ° ) ).
На эллипсоиде точное расстояние между параллелями в точках φ 1 и φ 2 равно m ( φ 1 ) − m ( φ 2 ) . Для WGS84 приблизительное выражение для расстояния Δ m между двумя параллелями на расстоянии ±0,5° от окружности на широте φ определяется выражением
Четверть меридиана
[ редактировать ]Вы можете помочь расширить этот раздел текстом, переведенным из соответствующей статьи на немецком языке . (Май 2021 г.) Нажмите [показать], чтобы просмотреть важные инструкции по переводу. |
Расстояние от экватора до полюса, четверти меридиана (аналог четверти круга ), также известного как земной квадрант , равно
Это было частью исторического определения метра и морской мили , а также использовалось в определении гебдомометра .
Четверть меридиана можно выразить через полный эллиптический интеграл второго рода :
где – первый и второй эксцентриситеты .
Четвертьмеридиан также задается следующим обобщенным рядом:
(Формулу c 0 см. в разделе #Обобщенный ряд выше.)Этот результат был впервые получен Джеймсом Айвори . [25]
Числовое выражение для четверти меридиана на эллипсоиде WGS84:
Окружность полярной Земли просто равна четырем четвертям меридиана:
Периметр 2π меридионального эллипса также можно переписать в виде периметра спрямляющего круга C p = M r . Следовательно, выпрямляющий радиус Земли равен:
Его можно оценить как 6 367 449,146 м .
Обратная задача меридиана для эллипсоида
[ редактировать ]В некоторых задачах нам нужно уметь решить обратную задачу: по заданному m определить φ . Эту проблему можно решить методом Ньютона , повторяя
до сближения. Подходящее начальное предположение задается формулой φ 0 = µ , где
– это выпрямляющая широта . Заметим, что при этом нет необходимости дифференцировать ряды по m ( φ ) формулу для радиуса кривизны меридиана M ( φ ) , поскольку вместо этого можно использовать .
Альтернативно, ряд Гельмерта для меридианного расстояния можно вернуть, чтобы получить [26] [27]
где
Точно так же ряд Бесселя для m через β можно преобразовать, чтобы получить [28]
где
Адриен-Мари Лежандр показал, что расстояние по геодезической на сфероиде такое же, как расстояние по периметру эллипса. [29] По этой причине выражение для m через β и обратное ему, приведенное выше, играют ключевую роль в решении геодезической задачи с m заменой на s , расстояние вдоль геодезической, и β, замененного на σ , длину дуги на вспомогательная сфера. [22] [30] Необходимые серии, расширенные до шестого порядка, даны Чарльзом Карни, [31] уравнения (17) и (21), где ε играет роль n , а τ играет роль µ .
См. также
[ редактировать ]Ссылки
[ редактировать ]- ^ Руссо, Лусио (2004). Забытая революция . Берлин: Шпрингер. п. 273 -277.
- ^ Торге, В.; Мюллер, Дж. (2012). Геодезия . Учебник Грютера. Грютер. п. 5. ISBN 978-3-11-025000-8 . Проверено 2 мая 2021 г.
- ^ Пойнтинг, Джон Генри; Джозеф Джон Томпсон (1907). Учебник физики, 4-е изд . Лондон: Charles Griffin & Co., с. 20 .
- ^ Виктор Ф., Ленцен; Роберт П. Мультауф (1964). «Документ 44: Развитие гравитационных маятников в 19 веке» . Бюллетень 240 Национального музея США: Материалы Музея истории и технологий перепечатаны в Бюллетене Смитсоновского института . Вашингтон: Издательство Смитсоновского института . п. 307 . Проверено 28 января 2009 г.
- ^ Исаак Ньютон: Начала , Книга III, Предложение XIX, Проблема III , перевод на английский Эндрю Мотта. Современный перевод с возможностью поиска доступен на сайте 17 Centurymaths . Найдите в следующем PDF-файле слово «сфероид».
- ^ Перейти обратно: а б с Кларк, Александр Росс (1880). Геодезия . Оксфорд: Кларендон Пресс. ОСЛК 2484948 . . Доступно бесплатно на сайтах Archive.org и Forgotten Books ( ISBN 9781440088650 ). Кроме того, книга была переиздана издательством Nabu Press ( ISBN 978-1286804131 ), первая глава посвящена истории ранних опросов.
- ^ Кларк, Александр Росс; Джеймс, Генри (1866). Сравнения эталонов длины Англии, Франции, Бельгии, Пруссии, России, Индии, Австралии, произведенные в Артиллерийском бюро Саутгемптона . Лондон: GE Eyre и W. Spottiswoode для канцелярского офиса HM. стр. 281–87. OCLC 906501 . Приложение к рисунку Земли.
- ^ Хопкинсон, Сара (2012). Справочник шкипера RYA - парус . Хэмбл: Королевская яхтенная ассоциация. п. 76. ИСБН 9781-9051-04949 .
- ^ Рапп, Р., (1991): Геометрическая геодезия, Часть I , §3.5.1, стр. 28–32.
- ^ Осборн, Питер (2013), Проекции Меркатора , doi : 10.5281/zenodo.35392 Раздел 5.6. Этот справочник включает вывод формул кривизны из основных принципов и доказательство теоремы Менье. (Дополнения: файлы Maxima , код и рисунки Latex )
- ^ FWJ Olver, DW Lozier, RF Boisvert и CW Clark, редакторы,2010, Справочник NIST по математическим функциям (Кембридж)Университетское издательство).
- ^ Руководство Mathematica: Эллиптические интегралы
- ^ Maxima , 2009, Система компьютерной алгебры, версия 5.20.1.
- ^ Эйлер, Л. (1755). «Элементы сфероидальной тригонометрии, взятые из метода максимумов и минимумов» . Мемуары Королевской академии наук Берлина 1753 г. (на французском языке). 9 : 258–293. Цифры .
- ^ Перейти обратно: а б Деламбре, JBJ (1799): Аналитические методы определения дуги меридиана ; им предшествуют мемуары А. М. Лежандра на ту же тему , De L'Imprimerie de Crapelet, Париж, 72–73.
- ^ Рэпп, Р., (1991), §3.6, стр. 36–40.
- ^ Бессель, Ф.В. (1837). «Оценка осей эллипсоида посредством измерений дуги меридиана» . Астрономические новости (на немецком языке). 14 (333): 333–346. Бибкод : 1837AN.....14..333B . дои : 10.1002/asna.18370142301 .
- ^ Гельмерт, Франция (1880): Математические и физические теории высшей геодезии , введение и 1 часть , напечатано и опубликовано Б. Г. Тойбнером, Лейпциг, § 1.7, стр. 44–48. Английский перевод (Аэронавигационного картографического и информационного центра, Сент-Луис) доступен по адресу: дои : 10.5281/zenodo.32050
- ^ Крюгер, Л. (1912): Равноугольное изображение земного эллипсоида на плоскости . Королевский прусский геодезический институт, новая серия 52, стр. 12
- ^ Дж. В. Хагер, Дж. Ф. Бехенский и Б. В. Дрю, 1989. Технический отчет TM 8358.2 Агентства оборонных карт. Универсальные сетки: Универсальная поперечная Меркатора (UTM) и Универсальная полярная стереографическая (UPS).
- ^ Руководство по системам координат в Великобритании , Артиллерийское управление Великобритании.
- ^ Перейти обратно: а б с Бессель, ФР (2010). «Расчет долготы и широты по геодезическим измерениям (1825 г.)». Астрон. Нахр . 331 (8): 852–861. arXiv : 0908.1824 . Бибкод : 2010AN....331..852K . дои : 10.1002/asna.201011352 . S2CID 118760590 . Английский перевод Астрона. Нахр. 4 , 241–254 (1825), §5.
- ^ Гельмерт (1880), §1.11
- ^ Кавасе, К. (2011): Общая формула для расчета длины дуги меридиана и ее применение для преобразования координат в проекции Гаусса-Крюгера , Бюллетень Управления геопространственной информации Японии , 59 , 1–13
- ^ Айвори, Дж. (1798). «Новая серия по исправлению многоточия» . Труды Королевского общества Эдинбурга . 4 (2): 177–190. дои : 10.1017/s0080456800030817 . S2CID 251572677 .
- ^ Гельмерт (1880), §1.10
- ^ Адамс, Оскар С. (1921). Разработки Latitude, связанные с геодезией и картографией . Специальная публикация береговой и геодезической службы США № 67. с. 127.
- ^ Гельмерт (1880), §5.6
- ^ Лежандр, AM (1811 г.). Упражнения по интегральному исчислению с различными порядками трансцендентов и квадратурами [ Упражнения по интегральному исчислению ] (на французском языке). Париж: Курьер. п. 180 . OCLC 312469983 .
- ^ Гельмерт (1880), гл. 5
- ^ Карни, CFF (2013). «Алгоритмы геодезических». Журнал геодезии . 87 (1): 43–55. arXiv : 1109.4448 . Бибкод : 2013JGeod..87...43K . дои : 10.1007/s00190-012-0578-z . S2CID 119310141 . Дополнения .