Преобразование некоторых систем координат
Преобразование геодезических (эллипсоидальных) координат B, L, H в прямоугольные (пространственные) координаты X, Y, Z.
(1)
где B и L – соответственно геодезические широта и долгота пункта, H – высота пункта над принятым референц-эллипсоидом;
N – радиус кривизны первого вертикала:
(2)
e – эксцентриситет меридианного эллипса:
(3)
a и b – большая и малая полуоси эллипсоида.
Преобразование прямоугольных (пространственных) координат X, Y, Z в геодезические (эллипсоидальные) координаты B, L, H.
Для обратного перехода от пространственных координат X, Y, Z к геодезическим B, L, H выполняются итерации при вычислении широты B и высоты H с использованием следующего алгоритма.
1) Вычисляется вспомогательная величина D по формуле:
. (4)
2) Анализируют значение D:
если D = 0, то
, , (5)
если D > 0, то
; (6)
при этом
если , то
если , то
если , то
если , то
3) Анализируют значение Z:
если Z = 0, то
Во всех других случаях используется следующая схема вычислений:
1) Находят вспомогательные величины r, c, p по формулам:
, (7)
, (8)
. (9)
2) реализуют итеративный процесс:
,
, (10)
;
если модуль разности меньше заданного ε, то
,
. (11)
В противном случае приравнивают s1 = s2 и вычисления повторяют, начиная с вычисления b. Значение ε во все случаях принимается равным 0,0001², что обеспечивает погрешность вычисления H не более 0,003 м.
Преобразование геодезических координат B и L в плоские прямоугольные координаты x и y
Плоские прямоугольные координаты x и y с погрешностью не более 0,001 м в принятой на территории Узбекистана проекции Гаусса-Крюгера по геодезическим координатам на эллипсоиде Красовского вычисляются по формулам:
(12)
(13)
В формуле (12) значение B выражают в радианной мере.
В формуле (13) n – номер шестиградусной зоны в проекции Гаусса-Крюгера. Величина l вычисляется в радианной мере по формуле
(14)
Преобразование пространственных прямоугольных координат из системы А в систему Б
Пространственные прямоугольные координаты двух систем (А и Б), у которых отличаются начало координат, ориентирование осей и линейный масштаб, связаны следующими соотношениями:
(15)
Обратное перевычисление, из системы Б в систему А, выполняется по формулам:
(16)
где m – масштабный коэффициент при переходе от системы А к системе Б;
wX , wY , wZ – угловые элементы трансформирования, углы вращения осей системы А, обеспечивающие их параллельность осям системы Б;
DX, DY, DZ – линейные элементы трансформирования, смещение начала системы координат А относительно системы Б.
Приведенные выше формулы называются формулами Гельмерта, а собственно преобразование – преобразованием Гельмерта по 7-ми параметрам.
Преобразование высот
Различие систем высот, применяемых в геодезии, связано с отсчетной поверхностью (эллипсоид, квазигеоид, геоид). Высоты отсчитываются от точки на земной поверхности до выбранной поверхности по отвесной линии (нормали) к ней.
Геодезическая (эллипсоидальная) высота – расстояние от точки на земной поверхности до поверхности референц-эллипсоида.
Нормальная высота – расстояние от точки на земной поверхности до поверхности квазигеоида.
Ортометрическая высота – расстояние от точки на земной поверхности до поверхности геоида.
Геодезическая высота H и нормальная высота Hg связаны соотношением:
(17)
где z – высота квазигеоида (геоида) над отсчетным эллипсоидом.
Высоты квазигеоида (геоида) над отсчетным эллипсоидом вычисляются по моделям гравитационного поля Земли, входящими в состав параметров общеземных систем координат.
Для вычисления нормальных высот, обычно используемых в практике геодезических работ, необходимо иметь карты высот квазигеоида или гравиметрические определения, точность которых определяет точность перевычисления.