Преобразования из геоцентрических в геоцентрические координаты
Эти преобразования могут использоваться как посредник между преобразованием из географических в географические координаты по схеме:
географические в геоцентрические > геоцентрические в геоцентрические > геоцентрические в географические
Geocentric translations
EPSG code: 9603
Параметры: смещение по оси X, смещение по оси Y, смещение по оси Z
Если исходная и конечная система координат геоцентрические, оси эллипсоидов параллельны, главный меридиан — Гринвичский и нет разницы в масштабах, это преобразование позволяется вычислить координаты в конечной системе координат простым прибавлением смещения соответствующим координатам в исходной системе координат.
Xt = Xs + dX Yt = Ys + dY Zt = Zs + dZ
Position Vector
EPSG code: 9606
Параметры: смещение по оси X, смещение по оси Y, смещение по оси Z, поворот по оси X, поворот по оси Y, поворот по оси Z, масштабирование
Одно из 7-параметрических преобразований Гельмерта, использующее формулу Бурша-Вольфа.
где
Xs, Ys, Zs— координаты точки в исходной системе координат.
Xt, Yt, Zt— координаты точки в конечной системе координат.
dX, dY, dZ — вектор смещения, добавляемый к исходной точке, также является координатами начала координат исходной системы координат в конечной системе координат.
Rx, Ry, Rz — повороты, добавляемые к вектору смещения. Положительное значение означает поворот по часовой стрелке исходя из начала координат вдоль положительного хода соответствующей оси. Углы измеряются в радианах.
M — масштабирование вектора преобразования в исходной системе координат необходимое, чтобы получить правильный масштаб в конечной системе. M = 1+dS*10-6, где dS — масштабирование выражаемое в частях на миллион.
Coordinate frame rotation
EPSG code: 9607
Параметры: смещение по оси X, смещение по оси Y, смещение по оси Z, поворот по оси X, поворот по оси Y, поворот по оси Z, масштабирование.
Одно из 7-параметрических преобразований Гельмерта, использующее формулу Бурша-Вольфа.
Преобразование аналогичное Position Vector, но отличающееся инвертированными значениями поворотов Rx, Ry, Rz. Международная геодезическая ассоциация (IAG) и международный стандарт ISO 19111 (Geographic information — Spatial referencing by coordinates) рекомендуют для описания преобразования использовать Position Vector. В ArcGIS это преобразование эквивалентно преобразованию Бурша-Вольфа.
Это преобразование может использоваться как посредник между преобразованием из географических в географические координаты (см. Geocentric translations).
Две различные системы координат
Прежде чем погрузиться в код, представленный в этой статье, необходимо обсудить координатные системы, для поддержки которых этот код предназначен: известная система широты и долготы и универсальная поперечная проекция Меркатора (Universal Transverse Mercator — UTM). Также нужно коснуться военной системы координат (Military Grid Reference System — MGRS), которая основана на UTM.
Широта и долгота
Системы широты и долготы, вероятно, самый известный способ определения географических координат. В ней местоположение представляется двумя числами. Широта — это угол от центра земли к некоторой параллели на земной поверхности. Долгота — это угол от центра земли к некоторому меридиану на земной поверхности. Широта и долгота могут быть выражены в десятичных градусах (DD) или градусах, минутах и секунда (DMS); в последнем случае получаются числа в таком формате — 49°30’00» S 12°30’00» E. Этот формат обычно используется в GPS-навигаторах.
Земля разделена экватором (0° широты) на северное и южное полушария, и нулевым меридианом (0° долготы), воображаемой линией от северного к южному полюсу, которая проходитчерез город Гринвич в Великобритании и делит планету на восточное и западное полушарие. Диапазон широт в северном полушарии от 0 до 90 градусов, а в южном полушарии — от 0 до -90 градусов. Диапазон восточного полушария от 0 до 180 градусов, а западного полушария от 0 до -180 градусов.
Например, точка с координатами 61.44, 25.40 (в формате DD) или 61°26’24»N, 25°23’60»E (в формате DMS) находится в южной Финляндии. А точка с координатами -47.04, -73.48 (DD) или 47°02’24»S, 73°28’48»W (DMS) находится в южном Чили. На рисунке 1 приведено изображение Земли с перекрывающимися линиями параллелей и меридианов:
Рисунок 1. Земля с показанными линиями параллелей и меридианов
Дополнительную информацию можно найти в разделе .
Поперечная проекция Меркатора
Система координат UTM — это метод, использующий сетку для определения координат. Система UTM делит Землю на 60 зон, каждая из которых основана на поперечной проекции Меркатора. Проекция карты в картографии — это способ представить двумерную неровную поверхность на плоскости, как обычную карту. На приведена поперечная проекция Меркатора:
Рисунок 2. Поперечная проекция Меркатора
Зоны долготы в UTM пронумерованы от 1 до 60; все зоны кроме двух, о которых будет рассказано позже, имеют ширину 6° от востока до запада. Зоны долготы полностью покрывают поверхность Земли между широтами 80°S и 84°N.
Также есть 20 зон широты, каждая в 8° высотой; эти зоны пронумерованы от C до X, буквы I и O пропущены. Зоны A, B, Y и Z находятся за пределами этой системы, они покрывают Арктику и Антарктику. На приведены UTM зоны для Европы. На этом рисунке видны две нестандартные зоны долготы: зона 32V расширена для покрытия всей южной Норвегии, а зона 31V сокращена, чтобы покрывать только водное пространство.
Рисунок 3. Зоны UTM в Европе
Координаты в UTM представлены в формате зона долготы зона широты восточное склонение северное склонение, где восточное склонение — это проекционное расстояние от центрального меридиана зоны долготы, северное склонение — это проекционное расстояние от экватора. Значения восточного и северного склонений задаются в метрах. Например, координаты широты/долготы 61.44, 25.40 в UTM представлены как 35 V 414668 6812844; координаты широты/долготы -47.04, -73.48 соответствуют координатам 18 G 615471 4789269 в UTM.
В разделе приведена дополнительная информация по UTM и поперечной проекции Меркатора.
Военная система координат
Система координат MGRS — стандарт, используемый военными НАТО. MGRS основана на UTM и дальше делит каждую зону на квадраты 100х100 км. Эти квадраты идентифицируются двухбуквенным кодом, первая буква — восточно-западная позиция в зоне долготы, а вторая буква — северо-южная позиция.
Например, координата в UTM 35 V 414668 6812844 эквивалента координате MGRS 35VMJ1466812844. Точность координаты в MGRS задается с точностью в один метр и представлена с помощью 15 символов, где последние 10 символов представляют значения восточного и северного склонений в указанной сетке. В MGRS координаты могут быть представлены 15 символами, как в прошлом примере, или 13, 11, 9 или 7 символами; представленные таким образом значения будут соответственно иметь точность 1, 10, 100, 1000 или 10000 метров.
В этой статье подробно не разбирается MGRS, но скачиваемый код включает в себя преобразования между широтой/долготой и MGRS. В разделе приведена дополнительная информация.
Преобразования координат
Чтобы определить ширину и долготу — координаты местоположения на Земле, как минимум нужно иметь возможность видеть звезды или Солнце, иметь секстан и часы, показывающие время на меридиане GMT. Можно определить широту из угла между небесным телом и горизонтом, а долготу можно вычислить из вращения Земли. Эта статья не погружается в подобные подробности, но в разделе о них можно узнать больше. Вместо этого, предположим, что у нас уже есть координаты в формате DD, DMS или UTM.
Преобразование десятичных градусов в градусы/минуты/секунды и обратно
Крайне просто преобразовать координаты из формата DD в DMS. Ниже приведена формула для подобного преобразования:
DD: dd.ff DMS: dd mm ss dd=dd mm.gg=60*ff ss=60*gg
В этом примере — это дробная часть вычисления. Отрицательная широта означает местоположение в южном полушарии (S), а отрицательная долгота — в западном полушарии (W). Например, предположим, что имеются координаты в формате DD — 61.44, 25.40. Их можно преобразовать следующим образом:
lat dd=61 lat mm.gg=60*0.44=26.4 lat ss=60*0.4=24
Далее:
lon dd=25 lon mm.gg=60*0.40=24.0 lon ss=60*0.0=0
Таким образом, в формате DMS получаем следующие координаты — 61°26’24»N 25°24’00»E.
Ниже приведена формула для перехода от DMS к DD:
DD: dd.ff DMS: dd mm ss dd.ff=dd + mm/60 + ss/3600
Напомним, что места, расположенные в южном полушарии (S), имеют отрицательную широту, а места в западном полушарии (W) имеют отрицательную долготу.
Теперь выполним преобразование DMS координат 47°02’24»S, 73°28’48»W в формат DD:
lat dd.ff= - (47 + 2/60 + 24/3600 )=-47.04 lon dd.ff= - (73 + 28/60 + 48/3600)=-73.48
Таким образом, координаты в DD равны -47.04, -73.48.
Преобразование долготы/широты в UTM и обратно
В отличие от десятичных координат, которые можно определить с помощью хронометра и часов, координаты UTM невозможно определить без помощи вычислений. Хотя эти вычисления ничто иное, как простая тригонометрия и алгебра, формулы у них достаточно сложные. Если ознакомиться со статьей «The Universal Grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS)» (ссылка на неё приведена в разделе ), то станет понятно, что я имею ввиду.
Формулы для преобразования UTM здесь не приводятся, но исходный код в следующем разделе немного освещает эту проблему, а в разделе есть ссылки на дополнительную информацию.
Алгоритмы перевода географических координат в прямоугольные
Для быстрого пересчета географических координат в прямолинейные и обратно действуют особые алгоритмы, которые стали основой автоматических программ по такому сервису. Разработаны также онлайн конвертеры, пересчитывающие как координаты Гаусса — Крюгера, так и UTM, когда градус нахождения объекта, даже его минута и секунда превращаются в точные метры — и наоборот, когда метры трансформируются в градусы.
В программу либо конвертер вводятся параметры широты с долготой, на которых расположен наш объект, а на выходе имеем величины x (горизонтальный параметр) и y (вертикальный параметр). Аналогично делается обратный перевод.
Формула пересчета (ключ) учитывает:
- нумерацию зоны по Гауссу-Крюгеру (из имеющихся 60-ти);
- коэффициент масштаба (для Гаусса-Крюгера это единица, для UTM это 0,9996);
- тригонометрические функции;
- начальную параллель;
- осевой меридиан;
- большую и малую полуоси;
- условные смещения, присущие начальной параллели по северу, а также центральному меридиану по востоку;
- величину приплюснутости;
- эксцентриситет.
В спутниковой навигации ГЛОНАСС и GPS действует постоянное отслеживание координат любого заданного формата. Можно задать величины, чтобы показывалась широта и долгота, а одновременно отображались метры либо километры.
Общая характеристика картографических данных объектов недвижимости
Под картографическими данными объектов недвижимости следует понимать их координаты. На обычной географической карте мира данные координаты обозначаются широтой и долготой. Однако в таких масштабах точно определить местоположение таких относительно мелких объектов, как дом или участок в несколько соток невозможно.
В этих целях ещё в Советском Союзе была разработана система координат СК63 (система координат 1963 года). Однако она не определяет индивидуальные координаты объектов, а лишь является способом их обозначения. На основе неё вся территория РФ была поделена на зоны, условно обозначенные латинскими буквами.
В целом, система использует три показателя:
- ширина и длина, обозначенные условно;
- высота – согласно Балтийской системе высот.
Однако основные инструменты обозначения остались неизменными. За основу берётся масштаб, равный 1:100 000.
Масштаб может быть увеличен или уменьшен в зависимости от насыщенности местности объектами недвижимости.
Данные для перевода геодезических данных в географические
Так как для каждого региона определены свои МСК, то и перевести геодезические данные в географические можно лишь при наличии ключевых данных по соответствующей МСК. МСК представлена в виде плоскости, но с указанием высоты координат. Так, ключевыми данными по МСК являются:
- Масштаб топографии.
- Ширина карты относительно общемировой параллели.
- Длина карты относительно меридиана.
- Отклонение касательно эллипса.
- Ключ расчёта.
Все данные помимо ключа расчёта можно получить, сравнив топографическую карту и стандартную карту России. Необходимо сопоставить масштаб и выявить точную ширину и длину в градусах. Ключ расчёта до недавних пор имел статус государственной тайны, так как ещё в 1963 году правительством Советского Союза было принято такое решение.
Хотя на сегодняшний день информация официально считается общедоступной, в официальных источниках её не найти. Однако на таких ресурсах, как mapbasic.ru имеется более или менее свежая информация по ключам.
Ключ расчёта представляет собой набор точек отсчета и отклонений на топографической карте. Например, ключ МСК Республики Адыгея выглядит следующим образом: «8, 1001, 7, 37.98333333333, 0, 1, 1300000, -4511057.628». Чтобы понять, на какую точку указывает каждая из этих цифр нужно иметь на руках геодезическую карту местности.
Что представляют собой прямоугольные координаты
Основа проекций эллипса на плоскость — что по Гауссу-Крюгеру, что по системе UTM — это принцип прямолинейных исчислений Декарта.
Система плоских прямоугольных координат
- За горизонтальную ось X берется абсцисса (параллель), идущая на восток, за вертикальную Y — ордината (меридиан), идущая на север, за начало отсчета O — их пересечение.
- Точка, отмеченная на плоскости карты, измеряется вертикальным расстоянием до линии оси X (это будет величина y), плюс горизонтальным до линии оси Y (это будет величина x).
- Плоскость делится осями на 4 части — так называемых квадранта с нумерацией против часовой стрелки (I, II, III, IV): I квадрант верхний правый (северо-восток), II верхний левый (северо-запад), III нижний левый (юго-запад), IV нижний правый (юго-восток).
Величины имеют как плюсовое значение, так и минусовое, что зависит от положения относительно квадранта:
- I квадрант имеет обе положительные величины (x, y);
- II квадрант задает смешанные величины (-x, y);
- III квадранту присущи обе отрицательные величины (-x,-y);
- IV квадрант обладает также смешанными величинами (x,-y).
Далее системы имеют существенные различия.
Для проекции Гаусса-Крюгера отображаемая на карте территория разделена на 60 зон, где расстояние между меридианами приравнено к 6º. Отсчет идет от Гринвича к востоку и к экватору на север. За коэффициент масштаба взята единица. Точкой отсчета выступает пересечение выбранного меридиана с экватором.
Для разработанной американцами системы UTM характерны аналогичные деления на 60 зон, но расчетный меридиан иной — первая по нумерации зона ведет начало от меридиана 177º западной долготы. Также отличия касаются масштабного коэффициента — он равен 0,9996. В системе UTM отсутствуют отрицательные значения — для этого к западной абсциссе приплюсовывают 500 километров, а к южной ординате — 10 тысяч километров.
Способы перевода
Перевести геодезические координаты объекта недвижимости в географические на сегодняшний день представляется задачей сложной. Всё дело заключается в закрытости подробной информации и нелинейности картографических данных, из-за чего в итоге расчётов могут возникнуть сдвиги от нескольких метров до нескольких километров.
Однако разработано множество программ, которые облегчают процесс проведения пересчёта. Одной из них является GPSMapEdit. К сожалению, российских программ не существует, и для расчёта необходимо самостоятельно вводить данные ключевых точек (ключ расчёта), причём делать это нужно очень точно. Но даже в этом случае сдвигов в несколько метров не избежать, так как ключи большинства регионов до сих пор не доработаны.
К примеру, переложение данных из геодезической карты в географическую и наоборот по ключам МСК – 50 (Москва) зона 2 получается точным, чего не скажешь обо всех регионах.
Другим способом является осуществление расчётов на специальных сайтах, где установлены геокалькуляторы с учётом данный российских топографических карт. К примеру:
- latlong.ru – сайт довольно популярный и по основным регионам (Москва, Санкт-Петербург) выдаёт точные результаты. Но всё зависит от вводных данных. Нужно просто ввести имеющие координаты по ГСК-2011 (Основа МСК) или координаты по GPS. Однако возможно сдвиги в несколько километров.
- www.the-mostly.ru – простой и удобный калькулятор для перевода картографических координат в виде десятичных дробей в стандартные показатели широты и долготы в градусах, минутах, секундах. Точная ссылка на вкладку: http://the-mostly.ru/konverter_geograficheskikh_koordinat.html.
Поискав на просторах интернета, можно найти десятки подобных сайтов, но никто не гарантирует точность расчётов и уж тем более не обеспечивает правовой статус перевода. Ведь в основном перевод может пригодиться для представления в государственные органы или же исполнения указаний органов власти.
Под ответом органа ставиться печать государственного образца, которая подстрахует от ответственности в случае чего.
Для чего нужен перевод?
Необходимость перевода координат объекта недвижимости может понадобиться в разных жизненных ситуациях. Это может быть как определение границ участка, так и расположение будущих зданий и сооружений на участке земли.
Вот классический пример.: Строительная компания при оформлении разрешительных документов для строительства предоставила в государственный орган карту местности с указанием на границы будущей постройки и затрагиваемые смежные участки земли. При проведении проверки со стороны госоргана выяснилось, что в части участка имеется исторический памятник – древнее захоронение. При выдаче разрешения на строительство на карте, как правило – географической, указывается точки, которые нельзя затрагивать. Возникает необходимость перевода указанных координат на карту местности, по которой рассчитывались границы строительства.
Стоит отметить, что даже наличие специальных знаний и навыков не могут гарантировать точность расчётов, так как карты представляют собой плоскость, тогда как реальный объект недвижимости располагается на трёхмерной поверхности. Именно данная нелинейность реальных показателей являются причиной ошибок в несколько метров, порой даже километров при определении местоположения таких объектов на карте.