Пропустить навигацию.
Главная

Пересчет координат

        Спутниковые системы определяют координаты в общеземных системах: GPS - в WGS-84, ГЛОНАСС - в ПЗ-90. Используют же координаты в референцной, например в СК-42, или местной системе. Возникает необходимость пересчета координат из одной системы в другую. Общеземные геоцентрические координаты WGS-84 пересчитывают в геодезические широты B, долготы L и высоты H. Геодезические высоты трансформируют в ортометрические высоты Hg. По геодезическим широтам и долготам вычисляют плоские прямоугольные координаты, например, в проекции UTM. От вычисленных таким образом координат необходимо перейти к соответствующим координатам, применяемым в России. Переход может быть выполнен на любом этапе. Естественно в самом начале перейти от WGS-84 к ПЗ-90, затем на эллипсоид Красовского, вычислить геодезические координаты, нормальные высоты и плоские координаты Гаусса-Крюгера.
        Пересчет геоцентрических координат Взаимосвязь между WGS-84 и ПЗ-90 устанавливают при помощи формул Гельмерта:

RWGS = DRO + (1+m)W.Rпз,

RWGS = (XWGS, YWGS, ZWGS)T, DRO = (DXO, DYO, DZO)T,

Rпз = (Xпз, Yпз, Zпз)T,

где RWGS и RПЗ - векторы в соответствующих координатных системах, DRO - вектор
начала координатной системы ПЗ-90 в системе WGS-84; m - различие линейных масштабов в этих системах; W - матрица поворота координатных осей, зависит от трех малых углов
wX, wY, wZ. Обычно угловые параметры не превышают 1". Поэтому матрицу W, где ее элементы в радианах, представляют в виде:

 

Таким образом, для пересчета координат надо знать семь параметров трансформирования: DXO, DYO, DZO, m, wX, wY, wZ. Их появление обусловлено точностью установки общеземных геодезических систем отсчета. Эти параметры многократно уточнялись. На момент подготовки учебного пособия параметры трансформирования официального статуса не имели. Для примера несколько их числовых значений приведено в таблице 11 (см. также: Параметры Земли 1990).

Параметры ПЗ-90
WGS-84
ПЗ-90
WGS-84
СК-42
ПЗ-90
СК-42
ПЗ-90
DXO, м -0,3 0 +26,6 +25,00
DYO, м +2,2 0 -134,8 -141,00
DZO, м +1,0 +1 -77,3 -80,00
m.106 -0,06 0 -0,06 0,00
wX'' -0,05 0 -0,17 0,00
wY'' -0,01 0 -0,39 -0,35
wZ'' -0,07 -0,2 -0,83 -0,66

Формулу Гельмерта, учитывая, что W-1 = WT, перепишем так, чтобы она была удобной для перехода от WGS-84 к ПЗ-90:

Rпз = (1- m)WT(RWGS - DRO) » (1- m)WTRWGS - DRO .

       Если координаты не трансформировать, то возникнет погрешность:

|d p| < (D XO2 + D YO2 + D ZO2)1/2 + Rз(3m2 + 2(w X2 + w Y2 + w Z2))1/2,

где Rз - радиус Земли, а w X, w Y, w Z выражены в радианах. В настоящее время, как следует из данных табл. 11, эта погрешность может достигать 10 м.

        Пересчет геоцентрических координат в квазигеоцентрические. Связь координат референцных и общеземных систем устанавливается формулой Гельмерта (Пеллинен, 1978):

R = DR + (1+m)W.r,

r = (Xr, Yr, Zr)T, R = (XR, YR, ZR)T, DR = (DXO, DYO, DZO)T,

где r и R - векторы соответственно в референцной и в общеземной системах координат, DR - вектор начала национальной системы в общеземной системе, m - разница в линейных масштабах систем; W - матрица поворота координатных осей от референцной системы к общеземной.
Ориентировочные значения всех семи параметров взаимосвязи референцной и общеземной систем приведены в таблице 11 (Базлов и Бовшин и др., 1995; Бойков и др., 1993; Пеллинен, 1992). Несколько другие значения даны в (Параметры Земли 1990).
С целью пересчета координат из общеземной в референцную систему уравнения связи
должны быть записаны в виде:

r = (1- m)WT(R - DR) .

        Предположим, определены координаты двух пунктов: RA пункта A и RВ пункта B. Тогда
уравнение связи для разности координат принимает следующий вид:

rB - rA = (1-m) WT(RB - RA) .

        В этой формуле нет линейных параметров DR . Вместе с ними исчезли и их погрешности.
Поэтому в относительных определениях имеет смысл трансформировать не координаты, а измеренные с базовой станции приращения координат.
        Пересчет координат в местную систему. Если с помощью спутниковой системы позиционирования определены координаты трех - четырех опорных пунктов, то появляется возможность пересчитать координаты текущих точек в координатную систему окружающих пунктов. Пересчет тем точнее, чем меньше площадь объекта работ.
         Для каждой опорной точки с известными GPS прямоугольными геоцентрическими координатами R в общеземной системе и вычисленными по геодезическим широтам, долготам и высотам координатам r в местной системе, используя преобразованные формулы трансформирования, можно составить следующие уравнения поправок
vi:

vX = DXO - Zr wY + Yr wZ + Xr m - (XR - Xr),

vY = DYO + Zr wX - Xr wZ + Yr m - (YR - Yr),

vZ = DZO - Yr wX + Xr wY + Zr m - (ZR - Zr).

В уравнениях 7 неизвестных параметров, образующих вектор:

t = (DXO, DYO, DZO, wX, wY, wZ, m)T.

        Чтобы его найти, нужны координаты в 2 системах минимум лоя 3 пунктов. Тогда будем иметь матричное уравнение поправок:

v = A t - L,

где в матрицу A собраны коэффициенты, стоящие в уравнениях поправок vi перед искомыми параметрами трансформирования, а в вектор L - разности координат в общеземной и местной системах. Пусть P весовая матрица поправок v. Решая задачу по MНК находим параметры:

t = (ATPA)-1. ATPL.

        Задачу можно упростить, исключив вектор DR. Для этого координаты некоторого опорного пункта А следует вычесть из соответствующих координат других опорных пунктов и составить следующего вида уравнения поправок:

vXi = - (Zri - ZrA)wY + (Yri- YrA)wZ + (Xri- XrA)m - [(XRi - Xr i) - (XRA - XrA)],

vYi = +(Zri - ZrA)wX - (Xri - XrA)wZ + (Yri - YrA)m - [(YRi - Yri) - (YRA- YrA)],

vZi = - (Yri - YrA)wX + (Xri- XrA)wY + (Zri - ZrA)m - [(ZRi - Zri) - (ZRA - ZrA)].

       Решив эти уравнения по МНК, определим вектор искомых параметров
e = (wX, wY, wZ, m)T

       В России создана новая Государственная геодезическая сеть. На смену координатам СК-42 могут быть введены новые координаты на эллипсоиде Красовского СК-95. Чтобы перейти от общеземной к новой референцной системе достаточно будет к геоцентрическим координатам прибавить поправки (Базлов и др., 1996):

   D X = -25,90 м,     D Y = +130,94 м,      D Z = +81,76 м.

        Пересчет геодезических координат. Геодезические широту B, долготу L и высоту H легко вычислить по трансформированным пространственным прямоугольным координатам X, Y, Z и по формулам раздела 2. Однако, чтобы найти X, Y, Z, надо иметь надежные параметры трансформирования. Кроме того, для отдельных регионов модель Гельмерта с 7 параметрами трансформирования может оказаться недостаточно точной. Поэтому применяют и более сложные модели трансформирования, например, полиномиальные:

 

Br, Lr и BW, LW - широты и долготы, отнесенные соответственно к референц-эллипсоиду и общеземному эллипсоиду, BWO, LWO - средние для данного региона значения координат BW, LW, au, bv - коэффициенты пропорциональности, выбирают произвольно, но так, чтобы U и V стали небольшими величинами, ast, bst - эмпирические коэффициенты. Их общее число равно (n+1)(n+2). Для их определения, надо иметь не менее (n+1)(n+2)/2 точек, для которых известны и Br, Lr и BW, LW. Полиномиальные модели эффективно исключают разного рода систематические, свойственные только для данного региона, искажения, которые не возможно учесть иными, более простыми методами ( Featherstone, 1997).

        Пересчет геодезических координат по дифференциальным формулам В географической практике иногда GPS измеряют геодезические координаты B, L, H. Их следует пересчитать из WGS-84 в СК-42. Это можно выполнить по дифференциальным формулам, учитывающим различия в размерах эллипсоидов, а также смещение начал и разворот осей координатных систем (Герасимов и др., 1993). В формулах геодезические координаты WGS-84 даны без индексов, а в СК-42 - с индексами r; в формулах размерности и знаки параметров трансформирования изменены так, чтобы они соответствовали данным табл. 11. Формулы имеют вид:

Br = B + [e2N sinB cosB da/a + (1+N2/a2)N sinB cosB de2/2 +

+ (DXo cosL + DYo sinL) sinB - DZo cosB]/(M+H) +

+ (1 + e2 cos2B)(wX sinL - wY cosL)/r" +

+ me2 sinB cosB;

Lr = L + (DXo sinL - DYo cosL)/(N+H) cosB -

- (1 - e2) tgB (wX cosL + wY sinL)/r" + wZ/r";

Hr = H - a.da/N + N sin2B de2/2 -

- (DXo cosL + DYo sinL) cosB - DZo sinB +

+ e2N sinB cosB (wX sinL - wY cosL)/r" -

- (N + H - e2N sin2B)m;

da = ar - a; de2 = er2 - e2, r" - число угловых секунд в радиане.

        Пересчет плоских прямоугольных координат. Поясним суть решения задачи на примере пересчета плоских прямоугольных UTM координат, найденных при помощи GPS, в произвольную локальную систему.
Между двумя плоскими прямоугольными координатными системами существует следующая связь:

 

        В формулах локальные координаты даны без штрихов, а координаты UTM - со штрихами.
Необходимо найти 6 неизвестных параметров трансформирования:
xO, yO, xO', yO' -
координаты точки, относительно которой координатные оси развернуты на угол g, m различие линейных масштабов. Пусть для n > 1 пунктов координаты известны в обеих системах. Примем

xI’= [x']/n , yI’= [y']/n .

        Квадратными скобками обозначены по Гауссу суммы заключенных в них величин. Иными словами, в качестве xO' принимается среднее из абсцисс, а в качестве yO' - среднее из ординат опорных точек. Применяя МНК, получаем:

xI = [x]/n , yI = [y]/n ,

C = m cosg = [Dx Dx' + Dy Dy']/[( Dx'2 + Dy'2)],

S = m sing = [Dx Dy' - Dx' Dy]/[( Dx'2 + Dy'2)],

m = (C2 + S2)1/2, tgg = S/C ,

где для i-ой опорной точки

Dxi' = xi- xI’, Dyi' = yi'- yI’, Dxi = xi - xI, Dyi = yi - yI.

        Этот метод можно использовать для пересчета прямоугольных координат UTM в координаты Гаусса-Крюгера. При этом для точек в пределах листа топографической карты, например масштаба 1:100 000, погрешности трансформирования останутся в пределах графической точности масштаба.
        Пересчет высот. По координатам X,Y,Z определяют геодезические высоты H. В повседневной практике в России пользуются не геодезическими, а нормальными
Hg . При этом

Hg = H - z ,

где z - высота квазигеоида. Таким образом, для вычисления нормальных высот по геодезическим высотам нужны карты высот квазигеоида или гравиметрические определения. Другим путем решения задачи являются измерения спутниковой системой на ряде пунктов, нормальные высоты которых уже известны. Тогда можно образовать разности из высот для двух пунктов:

D Hg = D H - D z .

         Как видно, разность нормальных высот отличается от разности геодезических высот на
величину Dz. Неучет величин Dz приведет к погрешностям в передаче нормальных высот. Погрешности тем больше, чем больше расстояния между пунктами. Приближенно можно полагать, что с каждым километром разности Dz изменяется на 1- 2 см. При расстояниях между пунктами 50 км эти погрешности могут превысить 1 м. Для учета изменений высот квазигеоида применяют методы интерполяции, например, по формуле (Непоклонов и др., 1996):

D z = Ax + By + Cxy + D,

где x, y - плоские прямоугольные координаты, A, B, C, D - коэффициенты. Чтобы их вычислить, надо знать высоты H и Hg для 5 или большего числа пунктов. Достаточно прост способ среднего весового, восходящий к Гауссу. Пусть для ряда пунктов известны высоты HGPS в местной системе и геодезические высоты HGPS. Алгоритм пересчета следующий
1. Определяем веса измерений Pi = (1/Si)q, где Si   - расстояние от текущей точки до i-го
пункта с известными высотами в двух системах, q > 0 (подбирается экспериментально).
2. Находим преобразованную высоту текущего пункта по формуле
HЛОК = HGPS + a Pi.( HЛОК - HGPS)i / aPi.