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

788

.pdf
Скачиваний:
1
Добавлен:
07.01.2021
Размер:
684.5 Кб
Скачать

В2 = 45˚02'34", 8726;

В1 =33˚01'18",2632.

По широтам B1, B2, B3 и B4 найти дуги меридиана X от экватора до точек с заданными широтами по формулам

 

 

 

 

 

 

X = a0 B

a2

 

sin 2B +

a4

sin B

a6

sin B +...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

= m

m2

 

+

 

 

3 m +

5

 

m , m = a(1e2 );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

2

 

 

 

 

 

 

 

8

 

 

4

16

 

6

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

=

m2

 

+

 

m4

 

+

15 m ,

m =

3 e2m

 

;

 

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

2

 

 

 

 

 

 

32

6

 

2

 

 

2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

=

m4

 

+

 

 

3

m , m =

 

5 e2m

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

2

 

 

 

 

 

 

 

6

 

 

 

4

 

 

4

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

6

 

 

=

m6

,

m = 7 e2m

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32

 

 

 

2

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

n2

 

n4

 

 

3

 

 

 

 

 

 

 

3

 

 

 

 

 

 

3

 

 

 

 

 

 

15

 

 

2

 

15

 

 

4

 

35

 

3

 

X =

 

1

+

 

+

 

 

 

B

n

 

 

 

 

 

 

n

 

sin 2B +

 

 

n

 

 

 

 

n

 

sin 4B

 

n

 

sin 6B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ n

 

 

4

 

64

 

 

2

 

 

 

 

 

 

 

16

А

 

16

 

 

 

 

 

64

 

 

 

 

 

 

48

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =

 

α

 

 

 

α2

α3

 

 

α4

 

 

 

 

 

1

,

 

 

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

+

4 + 8 +

16 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 f 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a – ольшая полуось эллипсоида; α

 

 

 

 

 

 

С

 

 

 

где

 

 

 

 

 

 

 

 

бсжатие; f – знаменатель сжатия.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 3 дуги меридианов обозначе-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ны как ХII и ХI. Разность Х равна дуге в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

метрах между соответствующими точками.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эта же дуга вычислена, кроме того, как BI"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в секундах. Таким же образом найдём дугу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SII

в метрах и BII" в секундах.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При разности широт в пределах не-

 

скольких градусов длину дуги

 

 

 

 

X найти по способу Симпсона (9).

 

 

 

 

 

X =

Bm B0 (M0 + M m + 4(M1 + M3 +... + M m1 )+

 

 

 

 

 

 

 

 

 

 

 

 

3m

 

 

 

+ 2(M 2 + M 4 +... + M m2 )),

 

 

 

 

 

 

(9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где B0 и Bm – широты начальной и конечной точек дуги меридиана;

11

Mi – радиус кривизны меридиана в заданной точке вычисляют по формуле (1);

m – число равных интервалов, на которое делят дугу меридиана (m – чётное число).

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

(9) выражают в радианах.

По формулам (7) или (8) вычислить дуги меридиана от экват ора до заданных точек для эллипсоида Ф. Н. Красовского и ГСК–2011 г.

Вычислить разности дуг.

По формулам Симпсона (9) с учётом (1) вычислить дуги второй

раз.

К сдаче представить:

И

 

 

 

– два варианта вычисления дуг, расхождения между вариантами для эллипсоида Ф. Н. Красовского, Д

– два варианта вычисления дуг, расхождения между вариантами для эллипсоида ГСК–2011.

Контрольные вопросы и задания

1.Как вычислить длиныбдуг меридианов?

2.Как рассчитатьидлины дуг от экватора до каждой из четырех

точек?

3.Что такоеСвысшая геодезия?

4.Напишите формулу С мпсона.АРешение геодезических задач

на поверхности земного эллипсоида

Цель: научиться решать геодезические задачи на поверхности земного эллипсоида.

Главные геодезические задачи, решаемые методами сфероидической геодезии, позволяют определить геодезические координаты B и L некоторой точки или длину геодезической линии и ее азимут на поверхности эллипсоида (рис. 4).

12

Решение ГГЗ для пунктов P и Q сводится к решению сфероидического треугольнка PDQ. При решении ПГЗ известными являются

 

 

 

 

 

P

сторона DR = 900 B1, длина геодезиче-

 

 

 

 

 

 

 

ской линии s и азимут этой линии A1. При

 

 

 

 

L

решении ОГЗ известными являются сто-

 

 

 

 

 

 

 

роны DR = 900 B1, DQ = 900 B2 и угол

 

 

 

 

 

 

 

L = L2 L1 . Следовательно,

решение

 

 

 

 

 

 

 

ГГЗ в свете сфероидической геодезии за-

 

 

 

 

 

A2-1

 

ключается в вычислении двух углов и

 

 

 

 

 

 

 

стороны сфероидического треугольника,

 

 

 

 

 

 

 

A1-2

 

 

 

 

 

 

 

Q2 (B2, L2)

по известным двум его сторонам и углу

 

 

Q1 (B1, L1)

 

 

 

 

 

 

 

 

 

между ними.

 

 

 

 

Рис. 4. Полярный

 

 

 

 

Формулы сферической тригономет-

сфероидический треугольник

И

решения

рии используются только для

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

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

Прием 1. Прямой путь решения, при котором используется ка-

кая-либо проекция поверхности эллипсоида на шар. По трем элемен-

там сфероидического треугольника находят элементы соответствую-

шара на эллипсоид. Одной из лучных реализаций этого приема счита-

ется способ Бесселя.

 

А

Прием 2. Косвенный путь решения предполагает нахождение

 

б

приращений геодез ческ х координат и азимута, от которых осуще-

и

 

С

 

 

ствляется переход к скомым величинам. Этот прием был распространен для решения ГГЗ для рядов триангуляции при построении АГС в XX столетии. В связи с развитием принципиально новых методов создания ГГЗ и вычислительной техники этот прием потерял актуальность.

Прием 3. Сводится к редуцированию исходных величин на плоскость в некоторой проекции, например в проекции Гаусса– Крюгера. ГГЗ решаются по формулам плоской тригонометрии, после чего выполняется обратное редуцирование.

Прием 4. Решение ГГЗ в пространстве. Для этого от известных величин на поверхности эллипсоида переходят к их значениям в пространственной прямоугольной системе координат. Дальнейшие вычисления ведутся по формулам аналитической геометрии с обратным преобразованием найденных элементов на поверхность эллипсоида.

13

Прямая геодезическая задача

Исходные данные: геодезические координаты B1 и L1 точки Q1 ,а также азимут A1-2 с точки Q1 на точку Q2 и расстояние S между пунк-

тами Q1 и Q2 (рис.4).

Требуется определить широту B2 и долготу L2 точки Q2 и обратный азимут A2-1 с точки Q2 на точку Q1.

Прямая геодезическая задача решается главным образом при вычислении геодезических координат пунктов триангуляции 1 класса. Зная координаты начального пункта, а также расстояния и азимуты на соседние пункты, вычисляют геодезические координаты и обратные азимуты для всех других пунктов, непосредственно связанных с на-

чальным пунктом.

 

 

 

 

 

 

 

 

 

И

 

 

Исходные данные: по своему порядковому номеру в списке

журнала преподавателя из табл. 4 взятьВ1

; L1;S и А1 2 ..

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задание

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

Вычислить геодезические координаты L1

и L2

и обратный ази-

мут А21

с пункта Q2

на пункт Q1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

Порядок выполнения задания

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

Порядок решен я прямой геодезической задачи по способу

Бесселя на эллипсо де Ф. Н. Красовского

 

 

 

 

 

Исходными данными являются широта B1, долгота L1, длина

геодезической линии s и азимут A1

этой линии в начальной точке. Не-

обходимо найти широту B2

и долготу L2

конечной точки и обратный

азимут A2 в этойСточке. При вычислениях угловые величины необхо-

димо выражать в радианах.

 

 

 

 

 

 

 

 

 

 

 

 

1.

Вычисляем приведённую широту начальной точки

 

 

 

 

tgU1

 

 

 

 

 

 

 

11

 

 

 

 

sinU

1

=

 

 

 

; cosU

1

=

 

 

 

; tgU

1

=

1e2tgB

. (10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+tg2U1

 

1+tg2U1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2. Для перехода от s к σ найдём вспомогательные величины A0 и σ1 по формулам

14

tgα =

tgU1

, tgA = sin A1 cos A1 .

 

 

 

 

 

(11)

 

 

 

 

 

 

1

 

 

 

 

cosA1

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

tgU1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введём обозначение k 2 = +e2 cos A0

 

и вычислим коэффициенты

α, β, γ по формулам (12) и α1, β1

 

по формулам (13)

 

 

 

 

 

 

 

 

 

 

A =1+ k 2

 

 

 

3

 

k

4 +

 

 

 

5

 

k6 ...;

 

 

 

 

 

 

 

 

64

 

256

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B = k 2

 

k 4

 

+

 

 

15

 

 

k6

...;

 

 

 

 

 

 

 

 

 

 

 

 

 

512

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

16

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 4 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C =

 

+512k6

+

...;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

128

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

7

 

 

 

4

 

 

 

15

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

k

 

 

;

 

α =

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bA

b

1

 

4

 

64

 

 

 

256

 

 

+...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

4 4

 

 

 

 

37

 

 

 

 

 

 

 

 

 

 

 

 

 

β = B = k

 

 

 

k

 

 

 

 

 

+

 

 

k6

...

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

4

 

 

 

 

8

 

 

 

 

 

 

512

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

k 4

 

4

 

 

 

 

k6

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

+...

 

;

 

 

 

 

 

(12)

 

 

 

 

 

γ =

A

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

128

 

 

 

 

128

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

б6

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

e

 

 

 

 

 

 

 

 

2

A0

 

 

 

 

 

 

 

 

 

 

 

4

 

;

(13)

β1 =

16

16

... cos

 

 

 

32

+... cos

 

A0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ = αS + βsin σcos(2σ1 +σ)+ γsin 2σcos(4σ1 + 2σ +...) .

(14)

За первое приближение примем σ0 = αS . Примем обозначения

Iσ = βsin σ0 cos(2σ1 +σ0 ); IIσ = γsin 2σ0 cos(4σ1 +σ0 ); IIIσ = Iσβcos(2σ1 +σ0 ).

Для вычисления σ" с погрешностью менее 0,001" получили формулу

15

σ = σ0 + Iσ + IIσ + IIIσ.

(15)

Находим сферическое расстояние σ по формуле и затем σ2 = σ1 + σ. Решаем прямую задачу на шаре по формулам

tgλ1 = sin A0tgσ1 ;

 

tgλ2 = sin A0tgσ2 ;

tgA =

tgA0

 

;

tgU

2

= −cos A tgσ

2

;

 

 

2

 

cos σ2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ = λ2 λ1.

 

 

(16)

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

 

И

 

 

 

 

 

 

 

 

 

 

 

 

l = λsin A0 (α1σ + β1 sin σcos 2σ1

+σ).

 

(17)

 

2

 

 

Д2

 

 

 

Вычисляем геодезическую широту конечной точки

 

 

 

А

 

 

 

 

 

 

б

2

 

 

 

 

 

 

 

tgB = 1+e tgU

 

 

 

и геодезическую долготуиL2=L1+ l.

Способ Бесселя по простоте и строгости решения превосходит любые другие,Си его можно пр менять при любых расстояниях между точками на эллипсо де. Погрешности по способу Бесселя зависят от количества верных значащих цифр, удерживаемых при вычислениях, так как ряды в формулах ( 14) и (17) быстро сходятся и могут быть вычислены с любой необходимой точностью.

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

Решив прямую геодезическую задачу, получим координаты второго пункта Q2 геодезической линии Q1Q2, то есть В2 и L2 , которые

будут исходными данными для решения обратной геодезической задачи.

16

Обратная геодезическая задача

К решению обратной геодезической задачи прибегали в некоторых методах уравнивания полигонов триангуляции 1-го класса. В прошлом различали:

1.Малые расстояния до 30–45 км.

2.Средние расстояния до 600 км.

3.Большие – до 5 000 км.

4.Очень большие – до 19 000–20 000 км.

В настоящее время для решения прямой и обратной геодезической задачи на эллипсоиде целесообразнее применять способ Бесселя. Точность получения прямого и обратного азимута и расстояния между пунктами Q1 и Q2 геодезической линии будут зависеть от погрешностей применяемого математического аппарата, исходных данных и

погрешностей вычислений.

 

 

 

Д

 

Рассмотрим решение обратной геодезической задачи.

 

Исходные данные: В1 ; L1

и В2 ; L2

; требуется получить S1 2 ; А21

и

А1 2 .

 

 

 

 

А

И

 

 

 

 

Порядок выполнения задания

 

 

 

 

 

б

 

 

 

 

Решение обратной геодезической задачи на эллипсоиде по спо-

 

 

 

и

 

 

 

 

 

 

собу Бесселя сложнее прямой задачи, так как по исходным данным

можно вычислить на шаре только приведённые широты U1 и U2

по

С

 

tgU

 

 

1e2 tgB .

 

формулам tgU

1

=

1e2 tgB ;

2

=

 

 

 

 

1

 

 

 

2

 

Предварительные вычисления

Разность долгот l выражается в радианах.

Для определения долготы λ на шаре необходимо знать A0, σ1 и δ, которые зависят от s и A1. Следовательно, применим последовательные приближения.

W =

 

1e2 sin2 B

 

;

W =

 

1e2 sin2 B

;

 

1

1

 

 

2

 

 

2

 

 

 

sin B

1e2

 

 

 

 

 

 

sin B

1e2

 

sin u1 =

1

 

 

;

sin u2

=

2

 

;

 

W1

W2

 

 

 

 

 

 

 

 

 

17

cosu

= cos B1 ; cosu

2

= cos B2 ; l = L

L

;

 

 

1

W1

2

1

 

 

 

 

W2

 

 

 

 

 

 

a1 = sin u1 sin u2 ;

 

a2 = cosu1 cosu2 ;

 

 

 

 

b1 = cosu1 sin u2 ;

b1 = sin u1 cosu2 .

 

 

 

 

λ= l+ δ.

 

 

 

 

 

В первом приближении принимаем δ = 0 радиан.

 

 

 

 

 

p = cosu2 sin λ; q = b1 b2 cos λ; A1′ = arccos

 

 

q

 

;

 

 

 

 

 

 

p2

+ q2

 

 

 

 

 

 

 

 

(18)

(19)

(20)

(21)

 

 

если p > 0 или A1 = 2

 

 

 

 

= 360

0

 

 

 

если p < 0.

где A1 = A1

π− A1

 

A1

Вычисляем дугу на шаре

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arccosσ = α1 2 cos λ

,

 

 

 

 

(22)

далее A0

и коэффициенты α,

 

β и поправку δ

И

 

 

sin A = cosu

sin A

;

cos

2

A =1

 

 

 

2

 

 

 

 

 

 

 

2

A cos σ

 

sinДA ; x = 2 α cos

 

 

0

 

 

1

 

 

 

1

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

1

 

0

 

 

 

e2

 

e4

 

 

e6

 

e4

 

 

e6

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

cos2

 

 

 

 

 

 

e6 cos4 A ;

 

 

α =

 

+

 

+

 

 

 

 

 

+

 

 

 

A +

 

 

 

 

 

 

2

 

4

 

 

 

 

 

 

 

 

 

А

0

 

 

128

 

 

0

 

 

 

 

 

 

 

16

 

 

16

 

16

 

 

 

 

 

 

 

 

 

 

 

e4

С

 

б

 

 

 

 

 

 

 

 

 

 

 

 

β

 

 

+

 

 

 

 

 

 

cos

2

A0

;

 

 

δ = (α σ β

x sin σ) sin A0 . (23)

 

=

16

16

32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

β и поправку δ получаем в радианах.

Вычисленное значение δ подставляем в (20) и повторяем все вычисления по формулам (21), (22) и (2 3) и т.д. Приближения заканчиваем, когда δi δi 1 < 48·10-11.

Значения λ, A1, σ, x и sinA0, полученные в последнем приближении, принимаем за окончательные.

Вычисляем коэффициенты A, B', C' и геодезическую линию на эллипсоиде s.

Для нахождения коэффициентов A, B', C' возьмем формулы

k

2

2

cos

2

A0

;

(24)

 

= e

 

18

B′ = e2

4

A =1+

k

2

3k 4

+

5k6

;

4

64

256

e4

 

 

 

 

 

cos2

 

A0

+

 

15

 

e6 cos4 A0 ;

 

512

 

16

 

 

 

 

 

 

 

 

 

C′ =

e4

 

3e6

 

 

 

64

128 .

 

(25)

(26)

(27)

Длину геодезической линии на эллипсоиде вычисляем по формулам

 

 

 

 

y = (cos4 A0 2x2 ) cos σ ;

 

 

 

 

 

 

s = (A σ +(Bx +Cy) sin σ)

 

1

 

 

 

 

 

α 1

 

,

 

(28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

где a – большая полуось эллипсоида, f – знаменатель сжатия.

 

 

Обратный азимут найдем из формул

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

p′ = cosu1 sin λ;

q′ = b1

cos

λb2 ;

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

= arccos

 

 

 

 

 

 

 

 

 

;

 

 

 

 

(29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

q

+ p

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где A

=180

0

, если p'>0 или

 

 

 

 

0

= π+ A

 

+ A

 

A = π− A =180

 

A , если

2

2

 

 

2

А

 

2

2

 

 

 

2

p'<0.

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Погрешности выч сленбй прямой и обратной задач не превы-

шают 2 мм в коорд натах

в длинах линий и 0,001” в азимутах.

Контрольные вопросы и задания

1. Какие приемы используются для решения сфероидического треугольника?

2. Как решается прямая геодезическая задача?

3. Опишите порядок решения прямой геодезической задачи по способу Бесселя на эллипсоиде Ф. Н. Красовского.

4. Что такое обратная геодезическая задача?

19

Лабораторная работа 4

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

Гауссового сближения меридианов

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

Необходимо вычислить прямоугольные координаты x и y в проекции Гаусса–Крюгера по геодезическим координатам B и L в системе СК95 и ГСК2011.

Каждый студент пользуется данными своего индивидуального задания, т. е. геодезическими координатами пункта. Индивидуальные задания составлены так, что все пункты находятся в полосах перекрытия шестиградусных зон. Шестиградусная зона проекции Гаусса– Крюгера представляет собой сфероидический двуугольник, построенный от одного полюса до друго-

го и ограниченный меридиана-

И

ми.

Д

В каждой зоне изображе-

ние осевого меридиана прини-

мается за ось абсцисс,

а изо ра-

жение экватора за ось ординат.

 

А

В проекции Гаусса–Крюгера

осевой меридиан

экваторб

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

Точка а (рис. 5) − изображение точки А (эллипсоида) на

плоскости. Положение a определяется прямоугольными координатами x и y. Каждый студент вычисляет координаты точек от осевых меридианов двух смежных зон, например от осевого меридиана

75˚ и 81˚.

Исходные данные: Для точки с геодезическими координатами В1 и L1, которые выбирают из табл. 4, по номеру студента в журнале

20

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