Преобразование Гельмерта - Википедия - Helmert transformation

Преобразование из опорной рамы 1 к опорной раме 2 может быть описано с тремя переводами? X, Δy, DZ, три поворота Rx, Ry, Rz и масштабный параметром ц.

В Преобразование Гельмерта (названный в честь Фридрих Роберт Хельмерт, 1843–1917) - метод трансформации в трехмерном пространстве. Часто используется в геодезия производить преобразования без искажений из одного датум другому. Преобразование Гельмерта также называют семипараметрическое преобразование и является преобразование подобия.

Определение

Это можно выразить как:

{ Displaystyle X_ {T} = C + mu RX ,}

куда

$Икс Т$ преобразованный вектор
$Икс$ - начальный вектор

В параметры находятся:

$C$ – вектор перевода. Содержит три переводы по осям координат
$μ$ – масштаб, которая безразмерна; если это дано в промилле, его необходимо разделить на 1000000 и суммировать до 1.
$р$ – матрица вращения. Состоит из трех осей (малая вращения вокруг каждой из трех осей координат) $р Икс$ , $р у$ , $р z$ . Матрица вращения - это ортогональная матрица. Углы указаны либо в градусы или же радианы.

Вариации

Частный случай - двумерное преобразование Гельмерта. Здесь нужно всего четыре параметра (два перевода, одно масштабирование, одно вращение). Их можно определить по двум известным точкам; если доступно больше баллов, можно провести проверки.

Иногда достаточно использовать преобразование пяти параметров, состоящий из трех перемещений, только одного поворота вокруг оси Z и одного изменения масштаба.

Ограничения

Преобразование Хельмерта использует только один масштабный коэффициент, поэтому оно не подходит для:

Манипуляции с размерными чертежами и фотографии
Сравнение деформаций бумаги при сканирование старые планы и карты.

В этих случаях более общий аффинное преобразование предпочтительнее.

Заявление

Преобразование Гельмерта используется, среди прочего, в геодезия для преобразования координат точки из одной системы координат в другую. Используя его, становится возможным конвертировать региональные геодезия указывает на WGS84 локации, используемые GPS.

Например, начиная с Координата Гаусса – Крюгера, $Икс$ и $у$ , плюс высота, $час$ , преобразуются в 3D-значения поэтапно:

Отменить картографическая проекция: расчет эллипсоидальной широты, долготы и высоты ( $W$ , $L$ , $ЧАС$ )
Конвертировать из геодезические координаты к геоцентрические координаты: Расчет $Икс$ , $у$ и $z$ относительно опорный эллипсоид геодезии
7-параметрическое преобразование (где $Икс$ , $у$ и $z$ изменяются почти равномерно, максимум несколько сотен метров, а расстояния изменяются на несколько мм на км).
Из-за этого положения, измеренные на суше, можно сравнивать с данными GPS; затем они могут быть внесены в съемку как новые точки, преобразованные в обратном порядке.

Третий шаг состоит в применении матрица вращения, умножение на коэффициент масштабирования ${ displaystyle mu = 1 + s}$ (со значением около 1) и добавлением трех переводов, $c Икс$ , $c у$ , $c z$ .

Координаты системы отсчета B выводятся из системы отсчета A по следующей формуле:^[1]

{ Displaystyle { begin {bmatrix} X Y Z end {bmatrix}} ^ {B} = { begin {bmatrix} c_ {x} c_ {y} c_ {z} end {bmatrix}} + (1 + s times 10 ^ {- 6}) cdot { begin {bmatrix} 1 & -r_ {z} & r_ {y} r_ {z} & 1 & -r_ {x} -r_ {y} & r_ {x} & 1 end {bmatrix}} cdot { begin {bmatrix} X Y Z end {bmatrix}} ^ {A}}

или для каждого отдельного параметра координаты:

{ displaystyle { begin {align} X_ {B} & = c_ {x} + (1 + s times 10 ^ {- 6}) cdot (X_ {A} -r_ {z} cdot Y_ {A } + r_ {y} cdot Z_ {A}) Y_ {B} & = c_ {y} + (1 + s times 10 ^ {- 6}) cdot (r_ {z} cdot X_ { A} + Y_ {A} -r_ {x} cdot Z_ {A}) Z_ {B} & = c_ {z} + (1 + s times 10 ^ {- 6}) cdot (-r_ {y} cdot X_ {A} + r_ {x} cdot Y_ {A} + Z_ {A}). end {выравнивается}}}

Для обратного преобразования каждый элемент умножается на -1.

Семь параметров определены для каждого региона с тремя или более «идентичными точками» обеих систем. Чтобы привести их в соответствие, небольшие несоответствия (обычно всего несколько см) скорректированный используя метод наименьших квадратов - то есть исключены статистически достоверным образом.

Стандартные параметры

Примечание: углы поворота, указанные в таблице, указаны в угловые секунды и должен быть преобразован в радианы перед использованием в расчете.

Область, край	Начальная точка отсчета	Целевая точка отсчета	$c Икс$ (метр )	$c у$ (метр)	$c z$ (метр)	$s$ (промилле )	$р Икс$ (угловая секунда )	$р у$ (угловая секунда )	$р z$ (угловая секунда )
Словения ETRS89	D48	D96	409.545	72.164	486.872	17.919665	−3.085957	−5.469110	11.020289
Англия, Шотландия, Уэльс	WGS84	OSGB36^[2]	−446.448	125.157	−542.06	20.4894	−0.1502	−0.247	−0.8421
Ирландия	WGS84	Ирландия 1965	−482.53	130.596	−564.557	−8.15	1.042	0.214	0.631
Германия	WGS84	DHDN	−591.28	−81.35	−396.39	−9.82	1.4770	−0.0736	−1.4580
Германия	WGS84	Бессель 1841	−582	−105	−414	−8.3	−1.04	−0.35	3.08
Германия	WGS84	Красовский 1940	−24	123	94	−1.1	−0.02	0.26	0.13
Австрия (BEV)	WGS84	MGI	−577.326	−90.129	−463.920	−2.423	5.137	1.474	5.297
Соединенные Штаты	WGS84	Кларк 1866	8	−160	−176	0	0	0	0

Это стандартные наборы параметров для 7-параметрического преобразования (или преобразования данных) между двумя базами данных. Для преобразования в обратном направлении необходимо рассчитать параметры обратного преобразования или применить обратное преобразование (как описано в статье «О геодезических преобразованиях»^[3]). Переводы $c Икс$ , $c у$ , $c z$ иногда описываются как $т Икс$ , $т у$ , $т z$ , или же $dx$ , $dy$ , $дз$ . Вращения р_Икс, р_у, и р_z иногда также описываются как ${ displaystyle omega}$ , ${ displaystyle phi}$ и ${ displaystyle kappa}$ .^{[ВОЗ? ]} В Соединенном Королевстве основной интерес представляет преобразование между датумом OSGB36, используемым в обзоре артиллерийского оборудования для привязки к сетке на картах Landranger и Explorer, в реализацию WGS84, используемую технологией GPS. В Система координат Гаусса – Крюгера используемый в Германии обычно относится к Бессель эллипсоид. Еще одна интересная информация была ED50 (European Datum 1950) на основе Эллипсоид Хейфорда. ED50 был частью основ НАТО координаты до 1980-х годов, и многие национальные системы координат Гаусса – Крюгера определены ED50.

Земля не имеет идеальной эллипсоидальной формы, но описывается как геоид. Вместо этого геоид Земли описывается множеством эллипсоидов. В зависимости от фактического местоположения для съемки и картографии использовался «эллипсоид с наилучшим локальным выравниванием». Стандартный набор параметров дает точность около 7 мес. для преобразования OSGB36 / WGS84. Это недостаточно точно для съемки, и Ordnance Survey дополняет эти результаты с помощью таблицы поиска дальнейших переводов, чтобы достичь 1 см точность.

Оценка параметров

Если параметры преобразования неизвестны, их можно рассчитать с помощью опорных точек (то есть точек, координаты которых известны до и после преобразования. Поскольку необходимо определить в общей сложности семь параметров (три перемещения, один масштаб, три поворота), должны быть известны как минимум две точки и одна координата третьей точки (например, координата Z). Это дает систему с семью уравнениями и семью неизвестными, которые можно решить.

На практике лучше использовать больше очков. Благодаря такому соответствию достигается большая точность и возможна статистическая оценка результатов. В этом случае расчет корректируется с помощью гауссова наименьших квадратов метод.

Числовое значение точности параметров преобразования получается путем вычисления значений в контрольных точках и взвешивания результатов относительно центроид точек.

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

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

внешняя ссылка

http://www.w-volk.de/museum/mathex02.htm
https://www.webcitation.org/query?url=http://www.geocities.com/mapref/savpub/savpub-23.htm%23item40&date=2009-10-26+02:12:14 (Геометрия для обмена данными)
http://www.mapref.org/
TrafoStar гибкие преобразования 3D BestFit с: 3 переводами, 3 поворотами, 3 масштабами, 3 аффинными параметрами
Вычисление преобразований Гельмерта

[1] Уравнения преобразования базы http://www.linz.govt.nz/geodetic/conversion-coordinates/geodetic-datum-conversion/datum-transformation-equations/index.aspx

[OSguide-2] Руководство по системам координат в Великобритании, версия 1.7, октябрь 2007 г. D00659 Ordnance Survey

[3] О геодезических преобразованиях, Бо-Гуннар Рейт, 2009 г. https://www.lantmateriet.se/contentassets/4a728c7e9f0145569edd5eb81fececa7/rapport_reit_eng.pdf

[1]

[2]

[3]