Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Geodeziya_1_chastina_26-09-2011

.pdf
Скачиваний:
98
Добавлен:
23.02.2016
Размер:
10.12 Mб
Скачать

501

Рис. 185. Координати точки А в різних зонах

Протягом кількох років різними авторами (А.М. Віровець, Б.Н. Рабінович, А.А. Тімофєєв, Є.Є. Бірюков, В.Л. Коган та ін.) складено багато таблиць. Але всі вони дозволяють виконувати перетворення координат із зони в зону для пунктів, віддалення яких від осьового меридіана не перевищує 3°30, тобто у смузі шириною біля 30на схід і на захід від роздільного меридіана.

При підготовці вихідних даних для розв’язання задач на великі відстані необхідність у перетворенні координат пунктів виникає у більш широкий смузі на стиках зон. У цьому випадку застосування вищезазначених таблиць або обмежене, або взагалі неможливе.

Найбільш зручні для цієї цілі таблиці В.П. Морозова. Вони дозволяють виконувати переобчислення координат пунктів із однієї шестиградусної зони проекції Гаусса в суміжну шестиградусну зону в смузі перекриття, що дорівнює ширині зони.

При користуванні зазначеними таблицями для перетворення прямокутних координат із однієї шестиградусної зони в іншу використовують такі формули:

502

xII xI co c x c (k)x k (c2 )x c2 (ck)x ck (k 2 )x k 2(c2k)x c2k (ck2 )x ck2 (k3 )x 0.1k3

yII yI ko (c)y c (k)y k (c2 )y c2 (ck)y ck (k 2 )y k 2

(c2k)y c2k (ck2 )y ck2 (k3 )y 0.1k3.

де хІ, уІ – відомі координати пункту, які необхідно перетворити; хІІ, уІІ – невідомі координати того ж пункту в системі сумі-

жної зони;

с = (хІ – х0) · 10 -5

k = (уІ – у0) · 10 -5 – інтерполяційні множники; х0 – найближча до хІ таблична абсциса;

у0 – таблична ордината, зазначена під значенням х0. Величини с0, k0, (с)х, (с)у, (k)х, (k)у, …, (k3)у є членами

рядів розкладу функції координат за ступенями різниці довгот осьових меридіанів. Вони дозволяють увести в координати пунктів, які необхідно перетворити в систему суміжної зони, поправки на поворот координатних осей і на зміну масштабу зображення в проекції Гаусса.

Значення величин с0, k0, (с)х, (с)у, (k)х, (k)у, …, (k3)у вибирають із таблиць без інтерполяції за аргументами х0, у0.

У коефіцієнтів, поданих у таблиці з подвійними знаками, верхні знаки використовуються при переході з західної зони в суміжну східну.

Для контролю обчислень перетворення координат доцільно здійснювати двічі, тобто після перетворення із західної зони в східну виконують зворотний перехід – із східної в західну з метою отримання вже відомих координат.

Надійний контроль надає також використання при обчисленнях двох систем коефіцієнтів, які відповідають двом аргументам, один з яких х0І, а другий х0І.

Точність переобчислення координат за вказаними вище таблицями залежить від величини k (принаймні при переході в суміжну зону вони дозволяють перетворити коорди-

503

нати з точністю до 1 м). Існують і більш точні формули та відповідні їм таблиці.

3.2. Перетворення геодезичних координат у плоскі прямокутні координати Гаусса

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

іL деякої точки А еліпсоїда і невідомими плоскими прямокутними координатами Гаусса х і у тієї ж точки. Скористаємося для цього системою прямокутних сферичних координат як перехідною, тобто спочатку встановимо зв’язок відомих і невідомих величин із прямокутними сферичними координатами Х і Y.

Залежність між плоскими прямокутними і сферичними координатами однієї і тієї ж точки, як уже зазначалося, виражається формулою (7.13). Якщо тепер прямокутні сфероїдичні координати виразити через геодезичні (географічні)

іпотім прийняти сфероїдичні прямокутні координати точки на еліпсоїді такими, що збігаються зі сферичними прямокутними координатами, та підставити ці вирази у рівняння (7.13), можна вважати, що задача буде розв’язана. Такий шлях виведення формул, що визначають залежність між географічними і плоскими прямокутними координатами, був наданий Б.С. Кузьміним.

Поверхня земного еліпсоїда близька до сферичної, тому для спрощення виведення приймемо її за поверхню кулі з радіусом, що дорівнює довжині нормалі N в даній точці. Тоді, вважаючи, що сферична ордината Y точки А (рис.

186) приблизно дорівнює довжині дуги паралелі АА0, можна записати

l"" Y .

r

504

Оскільки із трикутника АСС1 r N cos B, то

Y l"" N cos B.

Наведемо без виведення більш строгий вираз сферичної ординати через координати В і l:

Y

l"

N cos B

l"3

Nsin

2

B cos B,

(7.37)

"

6 3

 

 

 

 

 

 

 

де l′′ – довгота точки А відносно осьового меридіана ОР, вираженого в секундах; В – геодезична широта точки А;

N – довжина нормалі під широтою В;

ρ′′ – величина радіана, виражена в секундах градусної міри.

Рис. 186. Залежність між

Рис. 187. Залежність між

прямокутними сферичними

прямокутними сферичними

і плоскими координатами

і плоскими абсцисами

Підставивши тепер значення Y із виразу (7.37) в другу формулу (7.13) і зберігаючи при цьому члени третього порядку і менш відносно l, отримаємо значення плоскої ординати у, виражене через геодезичні координати, тобто

505

 

l

"

 

l

" 3

N

 

3

 

2

 

 

 

 

 

 

 

 

 

 

 

 

(7.38)

y

 

 

 

N cos B

 

 

 

 

 

cos

 

B(1 tg

 

B).

 

"

 

"

6

 

 

 

 

 

 

 

 

 

 

 

 

 

Сферичну абсцису Х точки А можна зобразити у вигляді суми двох дуг: дуги ОА0=S (рис. 187) і дуги

A A'

X S. Тоді будемо мати тотожність

 

0

0

Х = S + (Х – S).

(7.39)

 

 

 

Значення довжини дуги S, виражене в метрах, можна

вибрати за аргументом В даної точки зі спеціальних таблиць.

Отже, щоби найти Х, достатньо відшукати різницю Х –

S.

Рис. 188. Визначення різниці Х – S

Для отримання Х – S скористаємося фігурою АА0М, де А0М і АМ – дотичні в точках А0 і А до відповідних меридіанів, які створюють при точці М кут γ, що збігається з кутом зближення меридіанів даних точок.

З’єднавши точки А і А0 прямою, розглянемо отриманий трикутник АА0М на окремому рисунку (рис. 188). Прове-

демо AA0' перпендикулярно А0М. Перпендикуляр AA0"

506

приймемо рівним сферичній ординаті Y, а відрізок A0 A0'

різниці Х – S.

У рівнобедреному трикутнику АА0М кут

A 180o

90o .

 

0

2

 

 

 

2

 

 

 

 

 

 

 

 

У прямокутному трикутнику AA A'

кут

 

 

 

 

0

0

 

 

 

A 90o (90o

)

,

 

а

 

 

2

2

 

 

X S Y tg .

 

 

 

 

 

 

 

 

 

2

 

 

 

Замінивши через малість кута tg

 

 

на

, будемо ма-

2

2

 

 

 

 

 

ти

 

 

 

 

 

 

 

X S Y

.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

(7.40)

 

 

 

 

 

 

 

Підставивши у вираз (7.40) значення Y із формули (7.37), а γ із формули (7.32), отримаємо

X S Nl2 sin B cos B, 2 2

звідки

 

N

l 2

 

X S

 

 

 

 

sin B cos B.

(7.41)

2

 

 

 

 

 

 

Але за першою формулою (7.13) х = Х, тому

 

N

 

l

"

2

 

 

x S

 

 

 

 

sin B cos B.

(7.42)

 

 

 

 

2

 

 

"

 

 

 

 

 

 

При виведенні формул (7.38) і (7.42) були зроблені спрощення, тому ці формули дають приблизне вираження

507

плоских прямокутних координат через геодезичні координати. Вони забезпечують перетворення координат з точні-

стю ±5 м.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Якщо ввести позначення

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

8

 

 

a

 

 

 

sin B cos B 10 ,

 

 

 

2 2

 

 

 

 

c

 

 

N

cos B 104 ,

 

(7.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

l

"

 

3

 

3

 

2

 

y

 

 

 

 

 

 

cos

B(1 tg

B),

 

 

 

 

 

 

"

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l l" 10 4 ,

 

l

" L L , Х = S,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

то рівності (7.38) і (7.42) матимуть вигляд

 

 

 

x X al2 ,

 

(7.44)

 

 

 

yo

cl y.

 

 

 

 

 

 

Уформулах (7.44) х і уо – плоскі прямокутні координати Гаусса, при цьому уо – ордината, що показує відстань від даної точки до осьового меридіана зони, виражена в

метрах. Знак уо визначається знаком l.

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

Удеяких спеціальних посібниках розміщені таблиці, складені за формулами (7.44) з позначенням (7.43). Величини Х, а, с обирається з них за аргументом В (широта точки) із застосуванням лінійного інтерполювання. Поправка

δу (в метрах) обирається за аргументами В і lо (довгота точки від осьового меридіана). Вона алгебраїчно додається до сl. При від’ємних значеннях l табличні знаки δу необхідно змінити на зворотні, оскільки таблиця складена для додатних значень l.

Щоб уникнути зайвих помилок при перетворенні геодезичних координат у плоскі табличні інтерполяційні коефіцієнти можна ввести безпосередньо в робочі формули

508

(7.44). Посібником з обчислення азимута передбачені для цього такі позначення:

a a m

 

N

 

 

sin B cos 108 ,

 

 

 

 

 

 

 

1

1

 

2

2

 

 

 

 

 

 

 

 

 

 

 

b b m

N

cos B 104 ,

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

N

cos3

B(1 tg2 B) 1012 ,

(7.45)

6 3

2

 

 

 

 

 

 

 

 

 

n l" 10 4 ,

 

 

l L L ,

 

 

 

 

 

 

 

 

 

 

o

 

S X X1 Xm.

Тоді формули (8.38) і (8.42) набудуть вигляду

x (a a m)n2

X

1

Xm,

 

 

 

1

1

 

 

 

y

o

(b

b m b n2 )n,

(7.46)

 

1

1

2

 

 

де т = (В – Втабл.)′′ – інтервал інтерполяції;

а1, ∆b1 і ∆Х – зміни табличних величин а1, b1 і Х при зміні широти В на 1′′;

х і уо – шукані абсциса і умовна ордината.

Значення Х1, а1, b1, b2 та інтерполяційні коефіцієнти ∆Х, а1 і ∆b1 обираються без інтерполяції із таблиць за аргументом Втабл (найближчим меншим значенням широти В, закругленим до 10).

Формули (7.44) і (7.46) – ідентичні. Їх зовнішня відмінність викликана лише позначеннями, що відрізняються одне від одного.

У таблиці 35 наведений приклад перетворення координат за формулами (7.46).

Для перетворення геодезичних координат у плоскі з більш високою точністю можуть бути застосовані інші таблиці, складені за більш строгими формулами.

509

Таблиця 35 Перетворення геодезичних координат у плоскі прямокутні

координати

Елементи

Пункт Пул-

Пояснення

дій

формул

кове

 

1

В

52°2007′′,6

 

2

L

27 31 35.0

 

3

Lо

27°

 

4

l = L – Lо

+0°3135′′,0

Переведення в секунди здій-

5

п =l′′· 10-4

0,18950

 

 

 

снюється за допомогою таб-

6

п2

0.0359

лиць

 

7

т

7′′,6

Інтервал інтерполяції

 

 

 

т=52°2007′′,6–52°2000=7,6

8

а1

3633

За аргументом Втабл= 52°20

 

 

 

за таблицями. Там же ∆а1 =–

9

а1т

0

0.009

 

10

1+а1т)п2

130

 

11

Х1

5 800 536

30,910×7.6=235

12

∆Хт

235

13

х

5 800 901

До цілих метрів

14

b1

189 355

За аргументом В = 52°20

 

 

 

Там же ∆b1 = - 1.88

15

b1т

- 9

b2 = - 19

 

16

b2п2

- 1

 

17

b1+b1т+

189 345

 

 

+b2п2

 

 

18

(b1+b1т+

+35 881

Значення ординати від осьо-

19

+b2п2)п

5 535 881

вого меридіана

у

 

При переобчисленні координат із точністю до 0.001 м використовують таблиці для обчислення координат Гаусса

– Крюгера в межах широт від 30 до 80°, складені під керів-

510

ництвом Д.А. Ларіна. При складанні цих таблиць для обчислення абсциси й ординати використані рівняння, які, крім вказаних у формулах (7.38) і (7.42) двох членів, містять додаткові члени п’ятого і шостого порядків малої ве-

личини l"" .

Наведемо більш точні вирази для координат х і у, виведення яких надається в курсах сфероїдної геодезії:

x S

N

(

l"

)2 sin B cos B

N

(

l"

)4 sin B cos3 B(5

2

"

24

"

 

 

 

 

 

 

 

 

 

9 2 4 4 )

N

(

l"

)6 sin B cos5 B(61 58t2 t4 ),

720

"

 

 

 

 

 

 

 

 

 

 

 

yN ( l"" )cos B N ( l"" )3 cos3 B(1 t2 2 )

6

 

 

N

(

l"

)

5

cos

5

B(5 18t

2

t

4

14

2

58

2

t

2

),

120

"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

(7.47)

де N – радіус кривизни перерізу першого вертикала під широтою В;

t tgB ;

e' cos B ;

e' a2 2 b2 – другий ексцентриситет меридіанного еліпса; b

аі b – велика і мала півосі еліпсоїда.

Узгаданих таблицях Д.А. Ларіна довжина дуги меридіана позначена через Х, коефіцієнти при l2 і l4 у виразі для х – відповідно через а1 і а2, а коефіцієнти при l і l3 у виразі для у – через b1 і b2; члени п’ятого і шостого порядків від-

носно

l

позначені через x і

y . З цими позначеннями

 

 

 

 

формули (7.47) матимуть вигляд

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]