Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Бугаевский Л.М. - Математическая картография - 1998.pdf
Скачиваний:
633
Добавлен:
10.03.2016
Размер:
12.22 Mб
Скачать

РАЗДЕЛ 1. ОБЩАЯ ТЕОРИЯ МАТЕМАТИЧЕСКОЙ КАРТОГРАФИИ

В отличие от других учебников, в которых рассматрива­ ется только общая теория картографических проекций, в настоящей книге излагается общая теория математической картографии, которая включает общую теорию картографи­ ческих проекций (п.п. 1-4) и основные положения теории других элементов математической основы карт (п. 5).

1.1. ЗАКОНОМ ЕРНОСТИ И ОБЩИЕ ПОЛОЖ ЕНИЯ ОТОБРАЖ ЕНИЯ ПОВЕРХНОСТЕЙ НЕБЕСНЫХ ТЕЛ НА ПЛОСКОСТИ

1.1.1.ПОНЯТИЕ О ФИЗИЧЕСКОЙ ПОВЕРХНОСТИ ЗЕМЛИ

ИПОВЕРХНОСТЯХ ОТНОСИМОСТИ

Физическая поверхность Земли, как и других небесных тел, имеет сложную форму. Её изучение является первоосно­ вой для многих наук, в том числе математической картогра­ фии.

При этом используется понятие о поверхности геоида, кото­ рое ввел в 1873 году немецкий физик Листинг. В настоящее время под поверхностью геоида понимают уровенную поверх­ ность, проходящую через точку начала отсчета высот.

Уровенной поверхностью называется поверхность, ортого­ нальная к отвесным линиям, по которым в каждой точке по­ верхности данного небесного тела направлен вектор силы т я ­ жести. Строгое определение геоида связано со знанием строе­ ния земной коры.

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

икоторая совпадает с поверхностью геоида на морях и океанах

иотступает от неё до 2 метров на континентальной части Земли.

Представления ученых о фигуре Земли менялись с тече­ нием времени. С VI века до н. э. до конца XVII века Землю принимали за шар, с конца XVII века до второй половины

XIX века — фигуру, близкую к эллипсоиду вращения; со второй половины XIX века до сороковых годов XX столетия

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

В геодезии измерения, выполненные на физической по­ верхности, переносят на математическую, наиболее близкую к ф и зи ческой , которая может быть описана с о о тв ет с т ­ вующими уравнениями. В этой связи изучают и используют общеземной эллипсоид и референц-эллипсоиды.

Эллипсоид вращения, плоскость экватора и центр которого совпадает с плоскостью экватора и центром масс Земли и наи­ лучшим образом аппроксимирует поверхность геоида (квази­ геоида) в планетарном масштабе, называется общеземным эллипсоидом.

Эллипсоид, на поверхность которого отображаются мате­ риалы астрономо-геодезических работ и топографических съемок, и который наиболее полно соответствует поверхности геоида на соответствующие территории Земли, называется референц-эллипсоидом. Эти поверхности называются так же поверхностями относимости. В разных странах приняты свои референц-эллипсоиды, различающиеся своими параметрами (значения некоторых из них приведены в приложении 1).

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

Втеоретических исследованиях рассматриваются вопросы отображения и более сложных, регулярных поверхностей.

Вместе с тем, в настоящее время разрабатываются и ис­ пользуются картографические проекции, дающие отображе­ ние реальных поверхностей.

Применение этих проекций не даёт возможности в ы ­ полнять измерительные работы по картам, но они позволяют получить представление о формах картографируемых по­ верхностей небесных тел.

При изучении картографических проекций в математичес­ кой картографии исходят из того, что параметры исполь­ зуемого референц-эллипсоида, исходные геодезические даты

(и их аналогичные величины), а при необходимости и данные о фигурах и параметрах реальных небесных тел, полученные

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

1.1.2.СИСТЕМЫ КООРДИНАТ, ПРИМЕНЯЕМЫЕ В МАТЕМАТИЧЕСКОЙ КАРТОГРАФИИ

В математической картографии применяются пространст­ венные прямоугольные, криволинейные, плоские п рям о ­ угольные и полярные системы координат.

1.1.2.1. ГЕОГРАФИЧЕСКАЯ (ГЕОДЕЗИЧЕСКАЯ) И ГЕОЦЕНТРИЧЕСКАЯ СИСТЕМЫ КООРДИНАТ

Из бесчисленного множества параметрических линий, которые можно установить на поверхности эллипсоида (сферы), выбирем семейство географических параллелей и меридианов, составляющих географическую систему коорди­ нат: = const и Л = const.

В произвольной точке A (<ptA) эллипсоида проведем нор­ маль АО' к этой поверхности (рис.1), через которую можно провести бесчисленное множество нормальных сечений. Из них мы выбирем два главных: сечение, совпадаю щ ее с

Рис.1 Геодезическая и геоцентрическая пространственная система координат

плоскостью меридиана РАР\ называемое меридианным, и

сечение,

ортогональное первому в точке А,

н азы ваем ое

сечением

первого вертикала.

 

Радиусы кривизны этих нормальных сечений будут равны:

 

М

=

а ( 1 - е 2) / ( 1 - е 2 sin2 (p)3/2;

(1)

 

N

=

а / ( 1 - е 2 sin2 <р)|/2 ;

(2)

1/2

- первый эксцентриситет,

а, b — большая и малая полуоси эллипсоида вращения.

Рассмотрим пространственную геоцентрическую ситему координат ОХг Y Z , в которой начало совмещено с центром масс Земли (с центром эллипсоида вращ ения), ось Zr - направлена на средний северный полюс Земли, ось Хг - в точку пересечения Гринвичского меридиана с экватором, ось Yг - на восток.

Тогда связь геоцентрической и географической систем

координат может быть записана в виде:

 

X

= N cos^cosA ,

(3)

Y

= N cos^sinA ,

 

Z p = N ( l - e 2)sin^.

 

1.1.2.2. ТОПОЦЕНТРИЧЕСКАЯ (ГОРИЗОНТНАЯ)

И ПОЛЯРНАЯ

СФЕРОИДИЧЕСКАЯ (СФЕРИЧЕСКАЯ) СИСТЕМЫ КООРДИНАТ

Т о п о ц е н тр и ч е с к о й горизонтной системой к о о р д и н а т (рис. 2) будем называть систему, в которой начало совмещено

с произвольной точкой пространства

Q0 (<р0,Л0, Н 0), ось X лежит

в плоскости

меридиана

точки Q и

направлена на северный

(средн и й )

полюс, ось

Z с о в п а д а е т с н о р м ал ь ю 0'Q {) к

поверхности эллипсоида в точке Q, ось Y - дополняет систему

до левой.

 

 

 

Для установления связи этой системы QqXYZ и геоцент­

рической системы координат OXrYrZ r перенесем параллельно

,

i

l l

каждую из них в точку О . Тогда

:

Zr = Zr + e2 N0sm%; Z = Z+N{) + H0.

F*

Рис. 2. Топоцентрическая горизонтная и сфероидическая (сферическая) полярные системы координат

Повернув систему координат Q'XvYrZ[ на углы и

9 0 ° - ^ , получим формулы связи указанных систем координат в следующем виде:

( Х Л

X

'

0

Уг = А

Y

-

0

\ Z fj

KZ + N 0 + H QJ

<e2N0sin(po,

и

 

 

 

 

 

 

o

Y = A'

yr

~

0

 

 

(5)

KZ + e ’N t* n 9,j

где А, А' — соответственно матрица и транспонированная к А матрицы вращения:

- cos А,0 sin ф0

- sin XQ

cos X0 cos cp0

(6)

А = - sin X0 sin ф0

cos XQ

sin X0 cos ф0

cos ф0

0

sin ф0

 

С учетом (3 ) - (6) ф ормулы связи топоцентрической и географической систем координат принимают вид:

X = (N + А^Шфсозфо - cosф sin ф0 cos(k - А.0)] + + e2(Nо sin фо - N 8Шф)со8ф0;

Y = (N + h) cos q>sin(A - Я0) ;

(7)

Z = (N + /i)[sin ф sin фо + cosф cosф0 cos(k -

X0)\ +

+e2( N 0 sin фо - N sin ф) sin ф0 - (jV0 + # 0) ;

где h - превышения точек (при отображении поверхности Земли все И=0)

Т е п е р ь

введем

с ф е р о и д и ч е с к у ю

п о л я р н у ю

с и с те м у

координат

z = const,

а = const, где а - углы

между

нормаль­

ными плоскостями в точке полюса

Q,

Z

- углы м еж ду

нормалью O'Q0и направлениями в точке

О 'на

текущие точки

С. , леж ащ ие в соответствующих нормальных плоскостях

( рис.2 ). Обозначим СО'- N'Qи из рисунка

2 следует:

X = N о sin z cos а

 

Y = Nq sin z sin a

(8)

Z = Nq cos z - N 0

'

 

Если в точке С провести нормаль СО"к эллипсоиду, которая пересечется с осью вращения эллипсоида в точке О, то образу­ ется треугольник О О С .

Учитывая широту данной точки, значения сторон О С= N Q; 0"C=N', 0 0* = e2(Nsin<p-A^0sin^0) и значение N из формулы

(2), получаем по теореме

косинусов:

2

е2

N(> = N 0{1 - — (siny - sinф0)2[1 + — (5sin2 ф + sin2 ф0 +

(9)

+ 2sin<psin<p0 - 4]+...}.

Аналогично

N Q = N 1 + e 2(sin (p- sin<p0)jsin<p + ^-(sin9 ~ sin 9 0)(l - 2sin 2 ф|

Приравняв выражения (7) и (8) и учитывая (9) и (10), по­ лучаем следующие формулы связи полярных сфероидических координат z, Д и географических (геодезических) ср и А:

(H)

где:

/j

= sin фсо8ф0 - c o sф sin ф0 cos(>, -

A,0);

t2

= sin<psin<p0(sin<p + sinq>o) + (мпф - sin<p0)(3sin2 ф - 1);

r3

= sin2 ф - l/2 sin фо (sin ф - sinф0);

 

t4 = cos<psin(A. - A,0);

( 12)

t5

= sin ф sinф0 + созфсо5ф0 cos(A, -

Х0);

т= sin ф - sinф0.

Вбольшинстве случаев при выполнении вычислений дос­

т а т о ч н о

у д е р ж а т ь

в ф о р м у л а х (9)- (11) ч л е н ы до е2

включительно.

 

В этом

случае:

 

sin z cos а = tj

+ e2r{t] sin q>- cos <p0)+- • • »

sin^sinfl = r4(l+ e 2rsin^)+... ;

cosz = /5 +

sin cp - sin ^0)+*• •

При отображении поверхности шара получим: sin z cos а = г,; sin г sin а = /4; cos z = t5-

1.1.2.3. ПОЛЯРНАЯ ГЕОДЕЗИЧЕСКАЯ СИСТЕМА КООРДИНАТ

Полярными геодезическими координатами точки С ( (р,А ) называется длина геодезической линии s от полюса полярной системы координат Q0 ( ^ 0>>^0) до данной точки и азимут а линии Q0C в точке Q0 (рис 3). В этой системе семействами ко­ ординатных линий являются:

а * const — пучок геодезических линий из полюса ,

s const — геодезические окружности, ортогональные первому семейству, не являющиеся геодезическими линиями

ипредставляющие собой сложные кривые двойной кривизны.

1.1.2.4.ЭЛЛИПТИЧЕСКИЕ КООРДИНАТЫ

Любую картографируемую поверхность можно определить при помощи уравнений вида:

Ф ( х , у , г ) = 0 ;

где х, у, z , — прямоугольные пространственные координаты:

x = F l(u,v) ;

y = F2(u,v) ;

Z= Fi(u,v).

Рис. 3. Полярная геодезическая сис-

Рис. 4. Система эллиптических

тема координат

координат

Независимые переменные и, v являются криволинейными координатами, определяющими положение точки на карто­ графируемой поверхности.

В основу получения эллиптических координат положены две системы софокусных сферических эллипсов (Н.А.Урмаев, 1962 г., Л.А. Вахрамеева, 1986 г.).

Сферическим эллипсом называется геометрическое место точек на поверхности шара, сумма расстояний которых от двух неподвижных точек (фокусов) постоянна.

На рис. 4 фокус F является общим для сферического эл­ липса МС, который имеет второй фокус в точке F , и для сферического эллипса M B , второй фокус которого находится в точке Fv

Положение произвольной точки М определяется у д а ­ лением ее от ближайших фокусов: FM= а и F М = Ь. Если Я долгота точки Л/, отнесенная к плоскости начального мериди­

ана

РСАРХ перпендикулярной плоскости

чертежа FPPr то

по

формулам сферической тригонометрии:

 

cos я = sin sin (ро - cos(pcos<pQsin Я ;

 

 

cosb = sin <ps\n<pQ+ cos^cos^0 sin Я .

 

 

где q>Q- широта точки фокуса.

 

 

Примем дугу АС = и и дугу A B = v за эллиптические коор­

динаты. Зная а и Ь} можно найти эллиптические координаты

и и v по формулам:

 

 

а + b

 

 

sin и sin (pQ= cos—-— ;

 

 

а - b

 

 

sin v COS (pn = sin------- .

 

 

2

 

 

Контрольная формула : cos и cos v = cos cpcos Я.

 

Приведенные формулы показывают,

что эллиптические

координаты зависят от положения на поверхности шара ф о ­ кусов сферических эллипсов (F\F ,FXyFx ). По этому признаку эллиптические координаты подразделяются на различные системы, в частных случаях определяют координаты Гюйу, Пирса, Адамса.

Э ллиптические координаты Гюйу, в которых широта

л

полюса <р0 = , вычисляются по формулам:

4

 

cosa = —

(sin 40 - cospsinA) ;

 

 

2

 

 

 

 

v 2 /.

v

 

cos b = —

(sin <p+ cos <pcos A) ;

 

sin и =

-

 

a + b

 

V 2 cos

;

 

 

 

2

 

 

sinu =

r-

a - b

 

 

V2 sin------- .

Если

фокусы сферического 2эллипса расположить на эк­

ваторе

симметрично относительно

среднего или начального

 

 

 

з

_ л

меридиана АР и положить, что *о - —, то получают систему

4

эллиптических координат Пирса, вычисляемые по формулам:

cosa = cos<pcos^ + Х\;

a + b

sin w = V2 cos-

b

nr

a -

cos b = COS ф cos

 

sint; = V2 cos

-------

.

В случае, когда фокусы сферического эллипса располага­ ю тся в полю се и на э к в а т о р е , п о л у ч а ю т у р а в н е н и я поверхности шара в эллиптических координатах Адамса:

X= cos ^ cos А= coswcost? ;

У=cos^sin А = sinuy/l + cos2v + sinWl + cos2и j ;

Z = sin(p - —(sinuyl1 + cos21; - sinWl + cos2w J .

1.1.2.5. ЭЛЕМЕНТЫ СФЕРОИДИЧЕСКОЙ ТРАПЕЦИИ, ИЗОМЕТРИЧЕСКАЯ СИСТЕМА КООРДИНАТ

Рассмотрим на поверхности эллипсоида вращения вы ­ пуклую трапецию A B C D , ограниченную двумя бесконечно близкими меридианами и параллелями с разностью широт d(p

идолгот dA соответственно (рис.5).

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

Элементами этой трапеции являются:

бесконечно малый отрезок меридиана -

dsx - АВ = Mdcp,

(15)

бесконечно малый отрезок параллели -

ds2 = ВС = rdA,

(16)

линейный

элемент эллипсоида

-

 

 

 

(17)

азимут

линейного элемента -

 

,

( rdA \

(18)

а = arctg

-------

 

Mdcp’

 

площадь бесконечно малой трапеции -

dS = dslds2 = MrdcpdA,

(19)

где г - Ncosq? - радиус кривизны параллели с широтой (р. Учитывая (16), длина дуги параллели равна

 

 

(20)

Из

этих формул видно, что при

равенстве д и ф ф е р е н ­

циалов

d<p=dA бесконечно малые дуги

ds] и ds2 не равны,

т.к.

г * М.

Это обстоятельство в ряде

случаев, особенно

при

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

кой, в которой при равенстве дифференциалов аргументов равны между собой соответствующие бесконечно малые дуги меридианов и параллелей.

Запишем формулу квадрата линейного элемента (17) в виде:

Введем обозначение:

М

(21)

dq = — d(p,

г

 

тогда ф о р м у л а

(20)

принимает вид:

 

ds2 = r 2[dq2 + d%]

Теперь, при равенст­ ве дифференциалов dq = dA д л и н ы б е с ­

 

к о н е ч н о

м а л ы х

 

отрезков

меридиана

 

и п а р а л л е л и будут

Рис. 5. Элементы сфероидической трапе­

равны.

q, Л - изо ­

Здесь

ции

метрические коорди­

 

 

наты. При

этом Я —

одновременно изометрическая и геодезическая долгота. Изо­ метрическую широту q можно найти, проинтегрировав выра ­ жение (21) :

q,-=Jj—* d(p+ C .

J г

Учтя значения М и г; получим

J

а (1 - е 2)

(l —е2 sin2 <е>)'П d<p

 

\3/2

+ С =

 

( l - e 2sin2ср}'

acoscp

(1- e 2)d<p

л( l - e 2sin2<p)cos<p + С .

Умножим в числителе е 2 на тригонометрическую единицу и составим два интеграла:

dcp [ е cos (р d<p

+ С

cos " ( l - e 2sin2^)

Введем обозначение:

sin у = е sin

(22)

cos if/dip = е cos (pdcp