Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
генике_побединский_ГНСС.pdf
Скачиваний:
1035
Добавлен:
27.03.2016
Размер:
7.22 Mб
Скачать

3.2.2.Геодезические системы координат и их преобразования

Вгеодезической системе координат положение точки определя-

ется высотой Н над принятым референцтэллипсоидом, широтой В и долготой L.

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

Такая эллипсоидальная система координат применяется при обработке наземных геодезических измерений. В космической геодезии при создании спутниковых геодезических сетей, которые являются пространственными и физически не связаны с какой-либо отсчетной поверхностью, более удобна система пространственных прямоугольных координат X, Y, Z. Преобразование эллипсоидальных геодезических координат в прямоугольные осуществляется по формулам:

 

X = (N + Н) cos В cos L;

(3.15)

 

Y = (N + H)cosBsm L;

 

 

 

h2

 

где N= .

; а и b — соответственно большая и малая

Va2 cos2 В + b2 sin2

В

 

полуоси референц-эллипсоида (см. рис. 3.2).

Зависимость между геодезической широтой В и геоцентрической широтой определяется формулой:

tg<p = (l-e2)tgB.

Переход от X, УиНкВ, ЬиН неизбежно связан с итерациями при вычислении широты В и высоты Н. Удобные для вычисления на ЭВМ формулы можно получить из выражения (3.15):

6

X

 

 

tgB = . Z

+

Ne2

tgB = с + dtgB.

(3.16)

Формулы (3.16) позволяют организовать цикл приближений. Для

первого приближения принимается

 

 

 

 

tgB = с.

 

Для второго и последующих приближений по значению В из пре-

дыдущего приближения вычисляют значения

 

н = fx

2 +

v2 sec B - N ;

 

-

 

 

 

122

N= I

л/1-е2 sin2 В

после чего повторяют вычисления по формуле (3.16). Процесс приближений заканчивается при пренебрежимо малом расхождении АВ из двух последних приближений. Практически для достижения расхождения АВ = 0,03" во всех случаях достаточно трех приближений.

И эллипсоидальная, и пространственная прямоугольная системы координат являются геодезическими по определению, несмотря на различия между ними.

При обработке континентальных и региональных геодезических сетей используются различные геодезические координатные системы. Сведения о некоторых геодезических системах приведены в табл. 3.1.

Референц-

Год

Большая

Сжатие а

Таблица 3.1

Страны, ис-

эллипсоид

вычис-

полуось а,

 

пользующие

 

ления

м

 

референц-

Бесселя

1841

6 377 397

1:299,2

эллипсоид

Европы, Азии

Эйри

1849

6 377 563

1:299,3

Великобритания,

Деламбра

1810

6 376 428

1:311,5

Ирландия

Бельгия

Датский

-

6 377 104

1:300,0

Дания, Исландия

Плессиса

-

6 376 523

1:308,6

Франция

Струве

-

6 378 298

1:294,7

Испания

Хейфорда

1909

6 378 388

1:297,0

Европы, Азии,

 

 

 

 

Южной Америки,

Эвереста

1830

6 377 276

1:300,8

Антарктида

Индия, Пакистан,

Кларка

1858

6 378 293

1:294,3

Непал, Шри Ланка

Австралия, Ирлан-

Кларка

1866

6 378 206

1:295,0

дия

Северной и Цент-

Кларка

1880

6 378 249

1:293,5

ральной Америки

Африки, Барбадос,

 

 

 

 

Ямайка, Израиль,

Красовского

1940

6 378 245

 

Иордания, Иран

1:293,8

Страны Восточной

(СК-42)

 

 

 

Европы и Средней

Австралий-

1984

 

 

Азии, Антарктида

6 378 160

1:298,2

Австралия, Папуа-

ский

 

 

 

Новая Гвинея

GR80

1980

США

123

 

 

 

Продолжение табл. 3.1

Референц-

Год

Большая

Сжатие а

Страны, ис-

эллипсоид

вычис-

полуось я,

 

пользующие

 

ления

м

 

референц-

 

 

 

1:298,3

эллипсоид

WGS-72

1972

6 378 135

США

WGS-84

1984

6 378 137

1:298,3

США

ПЗ-90

1990

6 378 136

1:298,3

СНГ

Каждая из референцных систем имеет свой исходный пункт, на котором осуществлено «внешнее» ориентирование «своего» рефе- ренц-эллипсоида. Сведения о исходных пунктах основных эллипсоидов приведены в табл. 3.2 [6, 74].

 

 

 

Таблица 3.2

Страна

Эллипсоид,

Широта

Долгота

 

начальный пункт

 

+ 133°12'30,1"

Австралия

Австралийский,

-25°56'54,6"

 

Джонстон Оригон

 

+ 10°19'20,6"

Алжир, Тунис

Кларка 1880,

+36°5Г06,5"

 

Карфаген

+51°28'39,7"

0°00'00,0"

Великобритания

Эйри, Гринвич

Египет

Хейфорда, Гельмерта,

+30°01'42,9"

+31 1637,0"

 

Трэнзит-оф-Винес

+24°07'11,3"

+77°39'17,6"

Индия, Бирма,

Эвереста, Калианпур

Пакистан

Хейфорда, Рим,

+41°55'25,5"

+ 12°27'08,4"

Италия

 

Монте-Марио

+39°13'26,7"

 

Канада, США,

Кларка 1866, Мидс-

-98°32'30,5"

Мексика

Рэнч, штат Канзас

+32°03'26,4"

+ 118°46'56,0"

КНР

Бесселя, Нанкин

Польша

Бесселя, Борова Гура

+52°28'32,8"

+21°02'12,1"

Бывший СССР

Красовского,

+59°46'18,5"

+30°19'38,6"

 

Пулково

 

 

Франция

Кларка 1880, Париж

+48°50'46,5"

+2°20'48,6"

Чехословакия,

Бесселя,

+48°16'15,3"

+ 16°17'56,0"

Югославия,

Германскогль

 

 

Австрия

Бесселя, Стокгольм

+59°20'32,7"

+ 18°03'29,6"

Швеция

Страны Европы

Хейфорда, Потсдам

+52°22'51,4"

+ 13°03,58,9"

При таком ориентировании совмещается отвесная линия на исходном пункте с нормалью к референц-эллипсоиду, а плоскость мери-

124

диана исходного пункта устанавливается параллельно оси вращения Земли по астрономическому азимуту. При этом отвесную линию и азимут на исходном пункте определяют по астрономическим наблюдениям. Неучет уклонений отвесных линий на исходном пункте, так же как и отличие формы и размеров принятого референц-эллипсоида от общеземного, приводит к сдвигу начала системы координат (центра ре- ференц-эллипсоида) относительно центра масс Земли, не нарушая параллельности малой оси референц-эллипсоида с осью вращения Земли [6]. Поэтому для каждого исходного пункта будет своя геодезическая система координат. Если на исходном пункте азимут был определен с ошибкой dA и имеются уклонения отвеса, отличные от нуля

"-4

' рл

\

 

...

А-Ь)соз(р

,

(3.17)

. НА

)

 

то будет иметь место не только сдвиг начала системы координат на г н о и перекос координатных осей. Направляющие косинусы новых осей относительно старых в этом случае имеют вид [6]:

 

 

-(d/*cos2?sinL+>

1

(d/*sini?-7CosZ?)

+ £cos L +

 

 

 

+ 7] sin В sin L)

 

 

(dA cos В sin

L-

sin Я - 7 7 cos 2?)

1

- £ s i n 1 +

(3.18)

(<£4cosi?sinZ,+

-(flWcos^sinL-

+ /7 sin В cos L)

 

 

+ £cosL +

- £ cos L +

1

 

+ 77 sin В sin L)

+ 77 sin BcosL)

 

 

Общеземная система координат

Как и всякая другая геодезическая система, общеземная система координат определяется параметрами земного эллипсоида, гравитационным полем Земли и координатами пунктов на физической поверхности Земли. Начало системы координат располагается в центре масс Земли. Направления ее осей были зафиксированы на XIV Генеральной ассамблее Международной ассоциации геодезии (МАГ) в 1967 г. [18, 48]. Малая ось общего земного эллипсоида совмещена со средней осью вращения Земли. Пространственное направление средней оси вращения Земли закреплено относительно звезд координатами среднего полюса на среднюю эпоху 1900 - 1905 гг. - Международное условное начало (МУН) — Convensional International Origin (СЮ). После введения поправок за нутацию эта точка в соответствии с рекомендациями Международной службы вращения Земли IERS (International Earth Rotation Service) известна как условный земной по-

125

люс СТР (Conventional Terrestrial Pole). Его перемещение относительно мгновенного полюса вначале определялось Международной службой полюса по наблюдениям на шести широтных станциях, расположенных на разных долготах примерно на параллели 39° северной широты. Координаты широтных станций приведены в табл. 3.3 [7].

Таблица 3.3

Страна

Название

Широта

Долгота

 

станции

 

 

Сардиния

Карлсфорте

+39°08'08,9"

-Ю8°18'43,5"

Узбекистан

Китаб

+39°08'01,7"

+66°52'55,5"

Япония

Мудзусава

+39°08'03,4"

+ 141°07'51,9"

США

Юкайа

+39°08'12,0"

-123°12'34,5"

США

Цинциннати

+39°08'19,8"

-84°25'21,0"

США

Гейтерсберг

+39°08'13,2"

-77°1Г57,0"

В настоящее время положение полюса относительно МУН (координаты хр, ур) определяют с ошибкой в несколько сантиметров из регулярных наблюдений выполняемых методом радиоинтерферометрии со сверхдлинной базой (РСДБ) - Very Long Baseline Interferometry (VLBI).

Плоскость начального меридиана устанавливается определением начала отсчета долгот как результат обработки долготных наблюдений национальных служб времени, сотрудничающих в рамках Международного бюро времени (МБВ) - Bureau International de ГНеиге (BIH). Начало отсчета долгот задается принятыми долготами национальных служб после учета поправок за движение полюса и соответствует точке на среднем экваторе в период 1900-1905 гг. вблизи Гринвичского меридиана [18, 69].

По форме общеземная система координат может быть эллиптической ( В L % , #*) и пространственной прямоугольной (X, У, Z). Переход между ними осуществляется по формулам (3.15) и (3.16). В космической геодезии более удобно использовать пространственную прямоугольную систему.

Составной частью общеземной системы координат являются координаты пунктов, закрепленных на физической поверхности Земли и объединенных в геодезические сети. Различия общеземных координатных систем связаны с особенностями построения и обработки геодезических сетей.

Наиболее известной из современных является общеземная геоцентрическая координатная система ITRS (IERS Terrestrial Reference System), которая поддерживается Международной службой вращения

126

Земли IERS. Ежегодно, начиная с 1989 г., новейшими методами космической геодезии и измерениями с использованием спутниковых приемников формируется сеть пунктов ITRF (IERS Terrestrial Reference Frame). Сеть ITRF с высокой точностью закрепляет начало координат в центре масс Земли и ориентирует координатные оси. Точность положения пунктов оценивается погрешностью до 10 см. Со временем геоцентрические координаты пунктов сети ITRF вследствие непрерывного совершенствования сети и геодинамических процессов изменяются. Поэтому каталоги координат обновляют и указывают их эпоху, например, ITRF-89, ITRF-94 и т.д. Сеть 1994 г. содержала около 150 пунктов, расположенных на всех материках и на островах всех океанов.

Топоцентрическая система координат

Иногда наряду с геоцентрической геодезической системой координат удобно использовать топоцентрическую систему. В этой системе (рис. 3.7) определяется взаимное положение пунктов. Прямоугольные координаты АX, AY, AZB этой системе будут равны разности геоцентрических координат пункта и начала координат топоцентрической системы:

(АХ^

г-х?

 

AY

Yi-Yt

(3.19)

A Z

Z 2 - Z ,

системе могут

Сферические координаты

в топоцентрической

быть определены длиной стороны D12 и углами Л и Ф аналогичными углам у и 5 в звездной системе координат. Если за основную плоскость X0Y принята плоскость, параллельная экватору, то в этой плоскости угол Л12 отсчитывается от меридиана Гринвича до проекции линии а угол Ф12 характеризует наклон линии /;/2 к экватору. Углы Л и Ф в эк-

ваториальной топоцентрической системе иногда называют ориентирующими углами хорды D12.

Из рис. 3.7 очевидно:

ДГ

 

Л = a

rctg—\

 

 

Ф = arctg-

АЛ

 

 

(3.20)

AZ

 

 

 

Y1'

 

 

'УДЛГ'+А

 

 

D = \1АХ2 + ДУ2 + AZ2.

 

 

Обратный переход выполняется по формулам:

 

 

 

'cos Ф cos А^

 

(3.21)

 

AY

 

sin

 

 

= D cos Ф n Л

 

 

Если же за основную

плоскость взять плоскость горизонта пунк-

 

,

°

,

 

та РО, то сферические координаты будут определены углом h наклона

127

линии iji2 к плоскости горизонта (дополнением до 90° зенитного расстояния z) и азимутом А, отсчитанным в плоскости горизонта от направления на север до проекции линии /;/2. Переход от углов А и h к углам А и Ф осуществляется по формулам:

sin А

 

 

tgA = -tg/7Cos£ + sini?cos/4

г

(3.22)

sin Ф = sin /isin h - cos Вcosh cos A,J

 

где 5 - широта пункта i2.

Рис. 3.7. Топоцентрическая система координат

В топоцентрической горизонтальной системе, ось Z которой направлена по нормали к эллипсоиду, длина и направление хорды (азимут и зенитное расстояние) могут быть выражены через геодезические координаты пунктов - концов хорды [39]:

Z)2

= (N + Я)2

+ (N + Н)]

- 2(N + Я), (N + Я)2

cos ц/ -

-

sin В2 -

Nx sin В,)2

- 2e2(N2

sin Я2 - ЛГ, sin Вх

х (Я2 sin

- Я, sin Вх),

 

 

 

где cos у/ = sin Вх sin В2

+ cos Вх cos В2 cos(I2 - Z,,);

 

 

я

 

 

 

 

 

 

 

 

I2

2 (N, sin Д, -

N. sin В.) cos

L

 

 

- e ^

?

!

 

= crga

 

 

 

 

 

 

 

 

(N + Я)2 cosЯ2 sin(Z2 -Lx)'

128