Преобразование_геодезических_систем_координат

Преобразование геодезических систем координат Материал из Википедии — свободной энциклопедии В геодезии задача перехода между различными системами координат возникает из-за существования множества систем координат, возникающих во всем мире на протяжении долгого времени. Применение различных систем координат при решении практических задач геодезии, картографии, навигации и в геоинформационных системах неизбежно. Различают несколько типов преобразования координат: переход между различными форматами координат, переход между различными системами координат и картографическими проекциями, а также преобразование датумов. Все перечисленные виды преобразования будут рассмотрены в данной статье.[1] Содержание Изменение формата и единиц Переход между различными системами координат Из геодезической системы координат в прямоугольную Из прямоугольной системы координат в геодезическую Метод Ньютона-Рафсона Решение Ferrari Применение решения Феррари Степенной ряд Переход из геодезической системы координат в ENU и обратно Преобразование координат из прямоугольной в топоцентрическую систему координат ENU Преобразование координат из топоцентрической системы координат ENU в прямоугольную. Переход к другой картографической проекции Преобразования датумов Проекционные преобразования Преобразование Молоденского Уравнения множественной регрессии Преобразование Гельмерта Преобразование Молоденского — Бадекаса См. также Ссылки на литературу Примечания Изменение формата и единиц Обозначить географическое место обычно значит передать широту и долготу места. Числовые значения для широты и долготы могут быть представлены в нескольких различных видах единиц и форматах:[2] шестидесятеричная система: градусы, минуты и секунды: 40° 26′ 46″ N 79° 58′ 56″ W градусы и десятичные минуты: 40° 26.767′ N 79° 58.933′ W десятичные градусы: 40,446° N 79,982° W

В градусе 60 минут, а в минуте 60 секунд. Поэтому для перевода из формата градусы/минуты/секунды в формат десятичных градусов можно использовать формулу: десятичные градусы=градусы+минуты/60+секунды/3600. Для обратного перевода из формата десятичных градусов в формат градусы/минуты/секунды можно использовать формулы: градусы=[десятичные градусы] минуты=[60*(десятичные градусы-градусы)] секунды=3600*(десятичные градусы-градусы)-60*минуты где обозначение [x] значит, что надо взять целую часть x, и обратиться к «полочной функции». Переход между различными системами координат Преобразование системы координат — это переход из одной системы координат в другую, причём обе системы координат основаны на одном и том же геодезическом датуме. Часто задача преобразования заключается в переходе от геодезической системы координат к прямоугольным координатам либо переход от одной картографической проекции к другой. Из геодезической системы координат в прямоугольную Прямоугольные координаты точек в пространстве можно вычислить по известным геодезическим координатам этих точек (широта B, долгота L, высота H) по формулам:[3] где где и — экваториальный (большая полуось) и полярный радиусы (малая полуось), соответственно. — квадрат первого эксцентриситета эллипсоида. радиус кривизны первого вертикала — расстояние по нормали к эллипсоиду от точки пересечения поверхности эллипсоида нормалью до оси oZ (Рис. 1). Из прямоугольной системы координат в геодезическую При переходе от прямоугольных пространственных координат к геодезической системе координат (такой как WGS84) геодезические широты B и высоты H зачастую требуется вычислять итеративными методами, то есть выполнением последовательных приближений. Что касается долгот L, то они вычисляются обычным путём. Существует несколько методов для вычислений геодезических широт и высот, рассмотрим два из них. Метод Ньютона-Рафсона

Следующее иррациональное уравнение Боуринга[4] для геодезической широты решается итерационным методом Ньютона — Рафсона:[5][6] где , . Широту B можно найти из уравнения Высота H рассчитывается как: Итерация может быть преобразована к следующему виду: где Постоянная является хорошим начальным значением для итерации, когда . Боуринг показал, что в таком случаи уже первая итерация дает достаточно точное решение. Он использовал дополнительные тригонометрические функции в своей первоначальной формулировке. Решение Ferrari Вышеизложенное уравнение для можно решить методом Феррари:[7][8] Применение решения Феррари Существует ряд методов и алгоритмов, но наиболее точным, согласно Чжу[9], является следующая последовательность, установленная Хейккиненом[10]. Предполагается, что геодезические параметры известны.

Примечание: arctan2 [Y, X] — обратная касательная к четырем квадрантам. Степенной ряд начинается с Для малых е2 степенной ряд Переход из геодезической системы координат в ENU и обратно Преобразование из геодезических координат в топоцентрические координаты ENU состоит их двух этапов: 1. Преобразование координат из геодезической системы в прямоугольную. 2. Преобразование координат из прямоугольной в топоцентрическую систему координат ENU. Преобразование координат из прямоугольной в топоцентрическую систему координат ENU Чтобы преобразовать прямоугольные координаты в топоцентрические необходимо знать начальную точку топоцентрической системы координат, обычно она располагается в некоторой точке наблюдений. Если наблюдение производится в точке , а наблюдаемый объект в тогда радиус-вектор этого направления в системе координат ENU имеет вид:

Преобразование координат из топоцентрической системы координат ENU в прямоугольную. Путём обратного преобразования координат из прямоугольной системы получим топоцентрическую систему координат: Переход к другой картографической проекции Преобразование координат и положений на карте между различными картографическими проекциями карты, привязанные к одной и тоже геодезической поверхности, может быть выполнено либо с помощью формул прямого перехода из одной проекции в другую, либо сначала проекция преобразуется в промежуточную систему координат, такую как прямоугольная, а уже из неё в проекцию . Используемые формулы могут быть сложными, в некоторых случаях преобразование не имеет решения в замкнутой форме, и необходимо использовать приближенные методы. Обычно для выполнения задач преобразования координат используются компьютерные программы, например, с программой GEOTRANS, поддерживаемой DoD и NGA.[11] Преобразования датумов Преобразования между датумами могут быть выполнены различными способами. Существуют преобразования, которые позволяют совершить прямой переход от геодезических координат одного датума к геодезическим координатам другого датума. Существуют менее прямые переходы, которые преобразуют геодезические координаты в геоцентрические (ECEF), преобразуют геоцентрические координаты из одного датума в другой, затем преобразуют геоцентрические координаты другого датума обратно в геодезические. Также существуют проекционные преобразования, которые позволяют совершить прямой переход из одной (датум, проекция) пары к другой (датум, проекция) паре. Проекционные преобразования Проекционные преобразования позволяют совершить прямой переход от координат на карте для одной (картографическая проекция, датум) пары к координатам на карте для другой (картографическая проекция, датум) пары. Как пример можно привести метод NADCON для преобразования от Северо — Американского датума (North American Datum) (NAD) 1927 года к датуму NAD 1983 года[12]. Высокоточная эталонная сеть (The High Accuracy Reference Network) (HARN), высокоточная версия преобразований NADCON, имеет точность приблизительно 5 сантиметров. Национальное преобразование версия 2(The National Transformation version 2) (NTv2 (https://z5h64q92x9.net/proxy_u/en-ru.ru/https/en.wikipedia.org/wiki/NTv2)) является канадской версией NADCON для перехода между NAD 1927 и NAD 1983. Методы HARN также известны, как NAD 83/91 и Высокоточные проекционные сети (High Precision Grid Networks) (HPGN)[13]. Впоследствии, Австралия и

Новая Зеландия адаптировали для себя формат NTv2 для того, чтобы создать методы проекционных преобразований для переходов между их собственными местными датумами. Как и преобразования с помощью уравнений Базовое смещение между NAD27 и NAD83 — На карте показана множественной регрессии, проекционные приблизительная величина (в метрах) общего сдвига между методы используют интерполяцию низких одним и тем же набором значений координат на север и восток в порядков для преобразования координат карты, NAD27 и NAD83 но в двух пространствах вместо трех. NOAA предоставляет программное обеспечение (как часть NGS Geodetic Toolkit) для производства преобразований NADCON.[14][15] Преобразование Молоденского Преобразование Молоденского позволяет совершить прямой переход между геодезическими координатами разных датумов без необходимости промежуточного перехода к геоцентрическим координатам.[16] Для него требуются три смещения между центрами систем координат и разницы между большими полуосями и параметрами сжатия референц- эллипсоидов. Преобразование Молоденского используется Национальным агентством геопространственной разведки (National Geospatial-Intelligence Agency) (NGA) в их техническом сообщении TR8350.2, а также в программе GEOTRANS, поддерживаемой агентством.[17] Преобразование Молоденского было популярно до прихода эры современных компьютеров, и метод является частью многих геодезических программ. Уравнения множественной регрессии Преобразования датумов с использованием эмпирических методов множественной регрессии были созданы для достижения большей точности для маленьких географических регионов чем стандартные преобразования Молоденского. Данные преобразования используются для преобразования местных датумов, которые создаются для континентов или меньших регионов, в глобальные датумы, такие как WGS 84.[18] В стандарте NIMA TM 8350.2, Приложение D[19] перечислены преобразования с помощью уравнений множественной регрессии от нескольких местных датумов к WGS 84, с точностью около 2 метров.[20] Метод уравнений множественной регрессии позволяет совершить прямое преобразование геодезических координат без промежуточного перевода в геоцентрические координаты. Геодезические координаты в новом датуме B моделируются как полиномы вплоть до девятой степени в геодезических координатах изначального датума A. Например, приращение может быть разложено как (показано разложение только до квадратичных членов): где параметры подобранные множественной регрессией фактор масштаба начало датума для и выстраиваются схожие уравнения. При достаточном количестве пар координат (A, B) для точек в обоих датумах для хорошей статистики используются множественные методы регрессии для подбора параметров этих полиномов. Полиномы, вместе с подобранными коэффициентами, формируют уравнения

множественной регрессии. Преобразование Гельмерта Использование преобразования Гельмерта при переходе от геодезических координат датума к геодезическим координатам датума происходит в три шага: 1 Преобразование геодезических координат датума в геоцентрические; , для 2 Применение преобразования Гельмерта, с соответствующими параметрами преобразования для перехода от геоцентрических координат датума к геоцентрическим координатам датума ; 3 Преобразование геоцентрических координат в геодезические координаты для датума . Для геоцентрических координат XYZ преобразование Гельмерта имеет форму:[21] Преобразование Гельмерта это преобразование с семью элементами, с тремя параметрами смещения , тремя параметрами разворота и одним масштабным параметром . Преобразование Гельмерта это приближенный метод, который можно считать точным только, когда параметры преобразования малы по сравнению с величинами векторов геоцентрической системы координат. С такими условиями преобразование можно считать обратимым.[22] Преобразование Гельмерта с четырнадцатью параметрами, с линейной зависимостью от времени для каждого параметра, может быть использовано для наблюдений за изменением во времени географических координат из- за геоморфологических процессов, таких как континентальный дрейф[23] и землетрясения.[24] Оно было преобразовано в программное обеспечение, такое как инструмент Horizontal Time Dependent Positioning (HTDP) в программном обеспечении U.S. NGS.[25] Преобразование Молоденского — Бадекаса Для исключения связи между смещениями и разворотами преобразования Гельмерта могут быть использованы три дополнительных параметра, чтобы получить новый XYZ центр разворота ближе к преобразуемым координатам. Это преобразование с десятью параметрами называется преобразованием Молоденского — Бадекаса, его не стоит путать с более простым преобразованием Молоденского. Как при использовании преобразования Гельмерта, использование преобразования Молоденского — Бадекаса состоит из трех шагов: 1. Преобразование геодезических координат датума в геоцентрические. 2. Применение преобразования Молоденского — Бадекаса, с соответствующими параметрами преобразования для , для перехода от геоцентрических координат датума к геоцентрическим координатам датума . 3. Преобразование геоцентрических координат в геодезические координаты для датума . Преобразование имеет форму[26]:

где — начало отсчёта для разворота и масштабного преобразования, а — масштабный фактор. Преобразование Молоденского — Бадекаса используется для преобразования местных геодезических датумов в глобальные датумы, такие как WGS 84. Применяется при расчёте посадок воздушных судов на палубу морских судов в достаточно статичном состоянии (в случае присутствия сильной качки или ветра формула не применяется). В отличие от преобразования Гельмерта, преобразование Молоденского — Бадекаса необратимо из-за того, что начало отсчёта для разворота относится к изначальному датуму. См. также Геодезическая система координат Прямоугольная система координат Список картографических проекций Проекция Гаусса — Крюгера Универсальная поперечная проекция Меркатора Переход между различными системами координат Ссылки на литературу 1. Roger Foster Dan Mullaney Basic Geodesy Article 018: Conversions and Transformations (http://earth-inf o.nga.mil/GandG/coordsys/geoarticles/pdfs/Article018_Conversions_and_Transformations.pdf) (4 марта 2014). 2. Ordnance Survey Great Britain. Coordinate transformer (http://www.ordnancesurvey.co.uk/gps/transformati on). 3. B. Hofmann-Wellenhof H. Lichtenegger J. Collins. GPS - theory and practice. — 282 с. — ISBN 3-211- 82839-7. 4. Bowring B. R. Transformation from Spatial to Geographical Coordinates // Surv. Rev.. — 1976. — Т. 23, № 181. — С. 323—327. — doi:10.1179/003962676791280626 (https://dx.doi.org/10.1179%2F003962676 791280626). 5. Fukushima, T. Fast Transform from Geocentric to Geodetic Coordinates (англ.) // J. Geod. : journal. — 1999. — Vol. 73, no. 11. — P. 603—610. — doi:10.1007/s001900050271 (https://dx.doi.org/10.1007%2Fs0 01900050271). (Appendix B) 6. Sudano, J. J. (1997). «An exact conversion from an earth-centered coordinate system to latitude, longitude and altitude». doi:10.1109/NAECON.1997.622711 7. Direct Transformation from Geocentric to Geodetic Coordinates // Vermeille, H.H J. Geod.. — 2002. — Т. 76. — С. 451—454. — doi:10.1007/s00190-002-0273-6 (https://dx.doi.org/10.1007%2Fs00190-002-027 3-6). 8. Irene PoloBlanco Gonzalez-Vega. A symbolic analysis of Vermeille and Borkowski polynomials for transforming 3D Cartesian to geodetic coordinates // J. Geod.. — 2009. — Т. 83. — С. 1071—1081. — doi:10.1007/s00190-009-0325-2 (https://dx.doi.org/10.1007%2Fs00190-009-0325-2). 9. J.Zhu. Conversion of Earth-centered Earth-fixed coordinates to geodetic coordinates // IEEE Transactions on Aerospace and Electronic Systems. — 1994. — Т. 30. — С. 957—961. — doi:10.1109/7.303772 (http s://dx.doi.org/10.1109%2F7.303772). 10. M.Heikkinen. Geschlossene formeln zur berechnung räumlicher geodätischer koordinaten aus rechtwinkligen koordinaten // Z. Vermess.. — 1982. — Т. 107. — С. 207—211. 11. MSP GEOTRANS 3.3 (Geographic Translator) (http://earth-info.nga.mil/GandG/geotrans/). NGA: Coordinate Systems Analysis Branch. 12. ArcGIS Help 10.1: Grid-based methods (http://resources.arcgis.com/en/help/main/10.1/index.html#//003r00 000013000000). ESRI. 13. NADCON/HARN Datum ShiftMethod (http://www.bluemarblegeo.com/knowledgebase/geocalc/classdef/d atumshift/datumshifts/nadcon.html). bluemarblegeo.com. 14. NADCON - Version 4.2 (http://www.ngs.noaa.gov/PC_PROD/NADCON/). NOAA. 15. Donald M. Mulcare. NGS Toolkit, Part 8: The National Geodetic Survey NADCON Tool (http://www.profsur v.com/magazine/article.aspx?i=1193) (недоступная ссылка). Professional Surveyor Magazine.

Архивировано (https://web.archive.org/web/20140306001134/http://www.profsurv.com/magazine/article.aspx?i=1193) 6 марта 2014 года. 16. ArcGIS Help 10.1: Equation-based methods (http://resources.arcgis.com/en/help/main/10.1/index.html#//0 03r00000012000000). ESRI. 17. Datum Transformations (http://earth-info.nga.mil/GandG/coordsys/datums/index.html). National Geospatial-Intelligence Agency. 18. (https://web.archive.org/web/20160412230130/http://www.iho.int/iho_pubs/standard/S60_Ed3Eng.pdf) U Handbook on Datum Transformations Involving WGS 84 (https://web.archive.org/web/20160412230130/ht tp://www.iho.int/iho_pubs/standard/S60_Ed3Eng.pdf) (3rd ed.), Special Publication No. 60, Monaco: International Hydrographic Bureau, August 2008, <https://web.archive.org/web/20160412230130/http://www.iho.int/iho_pubs/standard/S60_Ed3Eng.pdf >. Проверено 10 января 2017. 19. DEPARTMENT OF DEFENSE WORLD GEODETIC SYSTEM 1984 Its Definition and Relationships with Local Geodetic Systems (http://earth-info.nga.mil/GandG/publications/tr8350.2/wgs84fin.pdf) . National Imagery and Mapping Agency (NIMA). 20. Taylor Chuck. High-Accuracy Datum Transformations (http://home.hiwaay.net/~taylorc/bookshelf/math-scie nce/geodesy/datum/transform/high-accuracy/). 21. Equations Used for Datum Transformations (http://www.linz.govt.nz/geodetic/conversion-coordinates/geod etic-datum-conversion/datum-transformation-equations/index.aspx). Land Information New Zealand (LINZ). 22. Geomatics Guidance Note Number 7, part 2 Coordinate Conversions and Transformations including Formulas (http://info.ogp.org.uk/geodesy/guides/docs/G7-2.pdf) (недоступная ссылка). International Association of Oil and Gas Producers (OGP). Архивировано (https://web.archive.org/web/20140306005736/htt p://info.ogp.org.uk/geodesy/guides/docs/G7-2.pdf) 6 марта 2014 года. 23. Paul Bolstad. GIS Fundamentals, 4th Edition (https://web.archive.org/web/20160202201558/http://www.pa ulbolstad.net/4thedition/samplechaps/GISFundChap3.pdf) . — Atlas books. — 93 с. — ISBN 978-0- 9717647-3-6. 24. Addendum to NIMA TR 8350.2: Implementation of the World Geodetic System 1984 (WGS 84) Reference Frame G1150 (http://gis-lab.info/docs/nima-tr8350.2-addendum.pdf) . National Geospatial-Intelligence Agency. 25. HTDP - Horizontal Time-Dependent Positioning (https://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml). U.S. National Geodetic Survey (NGS). 26. Molodensky-Badekas (7+3) Transformations (http://earth-info.nga.mil/GandG/coordsys/datums/molodensk y.html). National Geospatial Intelligence Agency (NGA). Примечания Источник — https://ru.wikipedia.org/w/index.php?title=Преобразование_геодезических_систем_координат&oldid=112244056 Эта страница в последний раз была отредактирована 8 февраля 2021 в 04:40. Текст доступен по лицензии Creative Commons Attribution-ShareAlike; в отдельных случаях могут действовать дополнительные условия. Wikipedia® — зарегистрированный товарный знак некоммерческой организации Wikimedia Foundation, Inc.