Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

FEATURES OF CALCULATING TRANSFORMATION PARAMETERS AND REGRESSION EQUATION FROM WGS84 TO CLARKE 1880 IN THE SYRIAN LAMBERT COORDINATE SYSTEM

Abaas Gafaar 1
1 Moscow State University of Geodesy and Cartography
The Lambert equiangular conic projection has been used in the Syrian Arab Republic since the 1920s and was adopted until 1958 as the basis for the creation of maps of Syrian lands by the Military Geodetic Directorate. The maps place in this projection, as well as the coordinates of their Control Points, are still used for research and large-scale civil engineering projects. This article presents the current state of the Lambert coordinate system in the Syrian Arab Republic and the calculation of the transformation parameters and the definition of the regression equation in excel for converting geodetic coordinates measured using GPS into geodetic coordinates on the Clark 1880 ellipse in Syria and the possibility of creating a spatial satellite geodetic network. Based on the results of the comparison it was found that the differences between the calculated and known local coordinates do not exceed 0.53 meters, taking into account that the GPS device used (Ashtic) and the monitoring method are designed to determine one point. And, too, the differences between the geodetic coordinates calculated using the program Leica_GEO_Office and the coordinates obtained from the regression equations do not exceed 0.004 seconds, and therefore the regression equations can be used to convert geodetic coordinates measured using GPS into geodetic coordinates on the Clark 1880 ellipse, and we convert geodetic coordinates to rectangular coordinates in the Lambert system using transformation equations.
equiangular conic Lambert projection
geodetic coordinate system
points of the geodetic network
ellipsoid
transformation parameters
regression equation
Scale factor

Картографическая проекция Ламберта в Сирии использует равноугольную проекцию, которая сохраняет углы при переходе от эллипса к поверхности проекции. Эта проекция была разработана французом (Дураффурад), который использовал базовые точки (тригонометрическую сетку) первого и второго ранга в расчетах, он использовал географические координаты, из которых он вывел прямоугольные координаты, используя формулу преобразования [1].

В рамках этой системы все карты были созданы до 1958 г. Военно-геодезическим управлением (в настоящее время Главное управление геодезических работ) в масштабе 1:25000 и 1:50 000, из которых были получены карты меньшего масштаба. Контрольные точки системы (которая простиралась от северного Алеппо до равнин Хорана на юге и от побережья на западе до окраин пустыни на востоке) послужили отправной точкой для их уплотнения российскими геодезистами, создавшими карты всех сирийских земель в конце 1950-х гг. на масштаб 1:200 000 по контракту с Министерством промышленности и торговли. Был использован эллипсоид Кларка 1880 г. [2].

В работе ставится цель совершенствования координатной основы Республики Сирии. На основе спутниковых наблюдений показать, как можно определить состояние сирийской системы координат Ламберта; параметры трансформации и уравнение регрессии для преобразования геодезических координат, измеренных с помощью GPS, в геодезические координаты на эллипсе Кларка 1880.

Материалы и методы исследования

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

Анализ современного состояния системы координат Ламберта в Сирии

Проекционная поверхность. Прямая коническая касательная к широте 38,5°, а поскольку сирийская территория расположена в пределах шести широт, трех северных и южных широт, то есть в пределах 35,5° и 41,5°, и был использован K0 = 0,99962560 для снижения линейных искажений, и деформации нет на двух широтах рядом с городами (Алеппо и Аль-Санамайн) [1, 3]:

40.25 = 1.75 + 38.5 = 1φгр.

36,75 = 1,75 – 38 = 2φгр.

Точка перспективы расположена в центре Земли, и, поскольку эта проекция является равноугольной, она зависит от математических формул в процессе проектирования. Географические координаты центра проекции:

Abaas001.wmf

Прямоугольная система. В сирийской конической проекции Ламберта ось y применяется к долготе, спроецированной графикой, проходящей через центр системы, а направление x с запада на восток перпендикулярно ей, что относительно широты, проходящей через центр системы, И для различения прямоугольной системы в сирийской конической проекции Ламберта из прямоугольной стереографической проекционной системы, так что все территории находятся в первой четверти и все координаты точек становятся положительными, центр системы координат (x, y) был перемещен на 300 км к западу и югу [3].

X = x + 300KM; Y = y + 300Km;

x = R sinγ; y = R0-Rcosγ.

Координаты с помощью этой системы с учетом коэффициента искажения следующим образом [4, 5]:

Abaas002.wmf

Abaas003.wmf

Abaas004.wmf = 0.99962560.

Abaas005.wmf.

Экваториальный радиус кривизны первого вертикала Abaas006.wmf

Abaas007.wmf,

Abaas008.wmf

Abaas009.wmf,

Abaas010.wmf.

И тоже можно вычислить с помощью двух формул

Abaas011.wmf,

Abaas012.wmf.

Масштабный коэффициент для точки или расстояния в несколько километров составляет

Abaas013.wmf,

Abaas014.wmf,

Abaas015.wmf

Abaas016.wmf,

где a – большая полуось; b – малая полуось; е2 – первый эксцентриситет; N – экваториальный радиус; f – геометрическое (полярное) сжатие; k0 – масштабный коэффициент в начале координат; (R, γ) – элементы полярных координат точки; γ – полярный угол, зажатый между разверткой конуса, проходящей через точку, и развертки конуса, проходящей от центра предложения; R – длина развертки конуса, проходящей через точку; R0 – длина развертки конуса, проходящей от центра предложения; ε – продольное искажение в км.

И масштабный коэффициент, равный единице на расстоянии 174728 м от φ0 севера и юга. А это примерно приводит к тому, что масштабный коэффициент равен единице. То есть нет искажений при ширине 36,75° и 40,25°. А линейное искажение на широте, проходящей через центр системы, равно -37,44 см/км.

Формирование уравнения регрессии

Регрессионный анализ – это статистический метод анализа различных факторов и понимания того, какие из них могут повлиять на достижение цели, а какие можно игнорировать. В нашем случае его можно использовать для формирования уравнений регрессии геодезических координат от WGS84 до эллипса Кларка 1880, которые можно использовать вместо параметров преобразования [6].

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

R-квадрат – это коэффициент корреляции, который измеряет силу линейной зависимости между двумя переменными. Чем больше абсолютное значение, тем сильнее взаимосвязь.

− 1 означает прочные позитивные отношения;

− -1 означает сильную отрицательную связь;

− 0 означает отсутствие каких-либо отношений вообще.

Вычисление параметров преобразования (Leica_GEO_Office)

 

У нас есть шесть точек, распределенных вблизи района Пальмиры, который находится близко к центру сирийских координат, как видно на рис. 1. В табл. 1 показаны их местные координаты и с помощью GPS.

На рис. 2 показано геометрическое сжатие и начальные координаты системы координат Ламберта.

Используя четыре точки и с помощью метода Бурсы – Вольфа (семь параметров), мы вычисляем параметры перехода от WGS84 на Кларка 1880.

Локальные координаты и местные геодезические координаты (Кларк 1880) оставшихся двух точек, рассчитанные с использованием параметров преобразования, показанных в табл. 2.

Таблица 1

Геодезические и декартовы координаты точек

Геодезические координаты Кларк 1880

Декартовы координаты

по Ламберту

Геодезические координаты WGS84

φ֯

λ֯

X (m)

Y (m)

H (m)

φ֯

λ֯

H (m)

1

34,43971069

37,20772322

286928,26

276690,96

914,879

34,43948567

37,20887906

952,27

2

34,57962625

37,2064295

286831,47

292205,73

696,523

34,57940711

37,20760986

733,12

3

34,53606069

37,31068669

296392,19

287365,64

759,118

34,53582172

37,31187142

796,85

4

34,33346022

36,94926556

263135,11

264973,23

933,472

34,33325933

36,95036692

969,22

5

34,42030367

37,09905006

276939,07

274558,71

779,585

34,42009283

37,10019569

816,13

6

34,41447153

37,30001428

295406,26

273884,84

725,067

34,41423958

37,30117783

763,36

Abaas1.pdf

Рис. 1. Пункты геодезической сети

Abaas2.pdf

Рис. 2. Параметры системы координат Ламберта

Параметры перехода от WGS 84 на Кларка 1880 следующие:

− смещение по оси X TX = 758.2309 m;

− смещение по оси Y TY = -176.8499 m;

− смещение по оси Z TZ = 323.1733 m;

− поворот вокруг оси X RX'' = -6.82987'';

− поворот вокруг оси Y RY'' = -16.88788'';

− поворот вокруг оси Z RZ'' = -29.98073'';

− масштабный фактор S (PPM)= -98.6265 PPM.

Таблица 2

Декартовы координаты и местные геодезические координаты

φ֯

λ֯

X (m)

Y (m)

5

34,42030762

37,09905859

276939,6057

274558,9367

6

34,4144807

37,30001663

295406,3671

273885,246

Формирование уравнения регрессии геодезических координат от WGS84 до эллипса Кларка 1880 в Excel

Для формирования уравнения регрессии геодезических координат от WGS84 до эллипса Кларка 1880 были использованы четыре точки (табл. 3) и вычислены коэффициенты регрессии долготы и широты (табл. 4), и с использованием этих формул будут рассчитаны геодезические координаты оставшихся точек в локальном эллипсоиде Кларка 1880. Эти координаты указаны ниже в табл. 5.

Таблица 3

Геодезические координаты точек с географическими разностями координат

WGS84

WGS84 – CLARK1880

φ֯

λ֯

Δφ֯

Δλ֯

1

34,43948567

37,20887906

-0,000225028

0,001155833

2

34,57940711

37,20760986

-0,000219139

0,001180361

3

34,53582172

37,31187142

-0,000238969

0,001184722

4

34,33325933

36,95036692

-0,000200883

0,001101361

Таблица 4

Коэффициенты регрессии долготы и широты

Регрессионная статистика

Δφ֯

Δλ֯

Множественный R

0,987686401

0,999732071

R-квадрат

0,975524426

0,999464213

Нормированный R-квадрат

0,926573277

0,998392639

Стандартная ошибка

4,27666E-06

1,53595E-06

Наблюдения

4

4

Коэффициенты

Y-пересечение

0,002792441

-0,009947338

φ֯

5,35117E-05

0,000181012

λ֯

-0,000130701

0,000130834

И уравнение регрессии для геодезических координат от 84-й рабочей группы до эллипса Кларка 1880 г. будет выглядеть следующим образом:

Abaas017.wmf

φКларк 1880Abaas018.wmf,

Abaas019.wmf,

λКларк 1880 Abaas020.wmf

Таблица 5

Рассчитанные геодезические координаты

φ֯ WGS84

λ֯ WGS84

Δφ֯

Δλ֯

φ֯ CLARK1880

λ֯ CLARK1880

5

34,42009283

37,10019569

-0,0002147

0,001137084

34,42030753

37,09905861

6

34,41423958

37,30117783

-0,000241282

0,001162319

34,41448087

37,30001551

Таблица 6

Локальные различия в координатах

X (m) выч

Y (m) выч

X (m) Исх

Y (m) Исх

ΔX (m)

ΔY (m)

5

276939,6057

274558,9367

276939,07

274558,71

0,5357

0,2267

6

295406,3671

273885,246

295406,26

273884,84

0,1071

0,406

Таблица 7

Различия геодезических координат

Δφ֯

Δλ֯

Δφ''

Δλ''

5

8,58E-08

-1,90E-08

0,000309

-0,0000688

6

-1,70E-07

1,12E-06

-0,00061

0,004029809

Результаты исследования и их обсуждение

Сличая локальные координаты, рассчитанные с использованием программы Leica_GEO_Office, с известными координатами, мы обнажаем следующие разницы (табл. 6).

Сличая местные геодезические координаты, рассчитанные с использованием уравнений регрессии, и координаты, рассчитанные с использованием программы Leica_GEO_Office, мы обнажаем следующие разницы (табл. 7).

Из результатов сличения мы замечаем следующее:

– Разницы между вычисленными локальными координатами и известными не превосходят 0,53 м, мы замечаем, что результаты находятся в допустимых пределах, и поэтому параметры преобразования, рассчитанные с помощью программы Leica_GEO_Office, могут быть использованы для преобразования измеренных координат в локальные координаты.

− Различия между геодезическими координатами, рассчитанными с помощью программы Leica_GEO_Office, и координатами, полученными из уравнений регрессии, не превышают 0,004 с, и поэтому уравнения регрессии можно использовать для преобразования геодезических координат, измеренных с помощью GPS, в геодезические координаты на эллипсе Кларка 1880, и мы преобразуем геодезические координаты в прямоугольные координаты в системе Ламберта с использованием уравнений преобразования.

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

Заключение

На основе результатов анализа были рассмотрены современное состояние системы координат Ламберта в Сирийской Арабской Республике и вычисление параметров преобразования с использованием программного обеспечения Leica_GEO_Office и определение уравнения регрессии в Excel для преобразования геодезических координат, измеренных с помощью GPS, в геодезические координаты на эллипсе Кларка 1880 в Сирии.

Эта система использовалась в Сирийской Арабской Республике с 1920-х гг. и была принята до 1958 г. в качестве основы для создания карт сирийских земель Военно-геодезическим управлением. Основываясь на результатах предыдущего анализа, мы можем использовать эту проекцию для обновления сирийской геодезической сети, чтобы после войны сделать реконструкцию страны.