Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Огородова Л.В. - Высшая геодезия. Часть III. Теоретическая геодезия - 2006.pdf
Скачиваний:
1017
Добавлен:
10.03.2016
Размер:
9.96 Mб
Скачать

§35. РЕДУКЦИЯ ЛИНЕЙНЫХ ИЗМЕРЕНИЙ

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

Редукция пространственного отрезка. Пусть измерено расстоя­ ние PQ = D (см. рис. 6.3), нужно получить длину pq геодезической линии. Редуцирование выполняют в два этапа: сначала находят хор­ ду p q - D эллипсоида и затем переходят от хорды к геодезической линии.

Редуцирование длины отрезка к хорде эллипсоида можно вы­ полнить так. Длина D связана с прямоугольными координатами X , У , Ζ , X , Υ , Z точек Р и Q известным соотношением

Р9 Р Р я9 я4 я

^

 

D2 = (Xq - X p)2+(Yq - Y p)2 +(Zq - Z qf .

(6.15)

Используя зависимости (2.6) между прямоугольными X, У, Z и геодезическими В, L, Н координатами, можно выразить длину D отрезка через геодезические координаты точек

D2 = (N p + H p)2 +(Nq +Hq)2- 2(Np + Н р)(Nq + Hq)cos^ +

(6.16)

e

где ψ - угол между нормалями к эллипсоиду в точках Р и Q, опреде­ ляемый формулой (2.39); определена формулой (6.5).

Поскольку плановые координаты точек и их проекций на эл­ липсоид одинаковы, длину D хорды можно получить, положив в уравнении (6.16) Нр = Hq = О

Вычитая выражение (6.17) из (6.16) и используя тождество

cos ψ= 1 - 2sin2 ψ!2,

(6.18)

177

получим зависимость между длиной D пространственного отрезка и хордой D эллипсоида

D 2 - D2 - {Hq - Н р)2 - 2{Hq - H p)(Nq - N p)~

(6.19)

- 4kN pNqsin2 у + 2A(Hqsin Bq - H p sin Bp),

H P

, H o

, H PH 4

(6.20)

N p

Nq

N pNq

 

Исключив из выражения (6.16) cosy/ с помощью равенства (6.17), находим иную запись формулы для перехода от простран­ ственного отрезка к хорде эллипсоида

(1 + k)D 2 =D2 - (H q - H p)2 - 2 (Hq - H p){Nq - N p) +

+2A(HqsinBq - H psinBp) +k[(Nq - N p)2 + ^ ( e 2 -2)]. (6'21)

Согласно уравнению (6.21) для нахождения длины хорды нуж­ но знать геодезические высоты и широты точек Р и Q.

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

К

(6.22)

Р = ^(Рр + Рч1

Np

Nq

^ р 1 + e 2cos2 Врcos2 А ^ 4 l +e'2cos2 Bqcos2 А

где е' - первый эксцентриситет.

Для сферы уравнения (6.20) и (6.21) примут вид

(l +k0)D 2

=D2 - ( H q - H p)2,

Я .

Я

я . я

к0 = -^ - +^ - + ^ ±

Р

Р

Р2 ’

(6.24)

(6.25)

178

откуда

P 2 -{ H q - H p)2

D:

(6.26)

(1 + ^ ) ( 1

+ ^ )

P

P

Оценим точность формулы (6.26). Введем обозначение

D\ =D2 ~(H q - H p)2,

(6.27)

(1+ k)D 2 = D2.

(6.28)

Логарифмируя и дифференцируя выражение (6.28), получим

dP_ dP{

D А

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

Оценим погрешность dk. Так как

Η « N ,

Н .

Н 0

к =——+ —- и

Н в

н а

N P

N<,

 

 

* dN B-----TdNa4 .

 

Ν 1 ρ

Ν 2

4

 

Ошибки dN и dN возникли из-за замены эллипсоида сферой радиуса р, поэтому d& = p - N p, dNq= p - Nq. Используя выражения (2.7) для радиуса кривизны первого вертикала и (6.22) - (6.23) для радиуса сферы, с точностью до членов порядка е2 получим

p - N = Ne2[^s\n{Bm-B )sin 2 5 -co s2 Втcos2 А],

где Вт- средняя широта отрезка.

Таким образом,

 

 

 

Н ве1

1

,

,

dk = —

bsi n(Bm- B p)sin2Bm- cos2 Втcos2 А] -

1Ур

1

 

 

- ^ — [jsin(Bm - Bq)sin2Bm- cos2 Bmcos2 А].

179

Положим Ν ρ в N q, s ' m ( B - B p ) = -sin(Bm- B ^ ) = sin(fi? - B p ) / 2, тогда

Λ

dk - Н« - Нр e2 sin sin 2В. (6.29)

2N

V /

Согласно этой формуле, ошибка dk зависит от превышения и длины линии (разности широт). Максимальной она будет на ши­ роте 45°, поскольку в этом случае радиусы кривизны эллипсоида изменяются наиболее быстро. На этой широте

dk = 0,5 · 10‘6 (Я - Я )sin(£y- Вр)/2,

где превышение Я - Я выражено в километрах.

Для линии длиной 2000 км {{Bq - Вр)/2= 9°) dk = 0,078 · 10'6 (Hq - Яр и даже при Hq-H p = 10 км ошибка редуцирования будет меньше 10Л Оценим разности длин хорды эллипсоида и окружности. Эта разность равна нулю, если точки Р и Q лежат на одной параллели и максимальна, когда они лежат на одном меридиане. Вычитая из длины D0 хорды окружности радиуса Мсрдлину хорды эллипсоида

(6.17), при ψ= Bq- Вр, найдем максимальное отличие

Ъ

аг - D 2 - е 2 (2 -

e 2 ){Nqsin Bq - N p sin Bp) 2 -

 

- 4(NpNq - M 20)sin2

- (Na - N p)2,

 

<■P)

 

 

где M - средний радиус кривизны меридиана.

 

Оценка по этой формуле дает такой результат:

 

D0, км

200

1000

1700

2000

D0 - D , м

0,0031

0,425

1,668

2,692

D0- D

1,6 · ю-8

4,2 · ΙΟ'8

1,0 · 106

1,3 · ю-6

D0

 

 

 

 

Таким образом, для расстояний до 1700 км ошибка из-за заме­ ны хорды эллипсоида хордой окружности меньше Ю-6 и ею можно пренебречь.

Для редуцирования с точностью 10-6 формулу (6.26) можно ис­ пользовать, если расстояние не превышает 1500 км, а превышение меньше 10 км.

180

Переход от хорды D эллипсоида к длине σ дуги нормального сечения или длине s геодезической линии. Будем считать, что σ и 5 равны между собой и равны длине дуги окружности радиуса р, оп­ ределенного выражениями (6.22), (6.23). Тогда

s - о - 2Rarcsin — = D

(6.30)

 

Эта формула имеет точность порядка---- , поэтому, если ошиб- р 4

ка редуцирования не должна превышать ЮЛ ее можно использо­ вать для расстояний до 200 км.

Найдем также зенитное расстояние z и азимут Apq хорды D. Обратимся вновь к рис. 6.3. Спроектируем ломаную PQOqOpHa за­ мыкающую РОр.

DcoszpQ = (N + H )q cosy/ - (N + Н )р -Asin BpJ

(6.31)

где zPQ- зенитное расстояние пространственного отрезка Д Нр и H(j - геодезические высоты точек Р и Q; cos допределен формулой (6.14), а функция А - формулой (6.5).

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

Dcoszpq = Nq cosy/ —N р - Asin В .

Исключая отсюда с помощью выражения (6.31) произведение A sin#p, найдем выражение, связывающее зенитное расстояние про­ странственного отрезка и хорды эллипсоида

_ Dcoszpg - Hq cosy/ + Η р

(6.32)

COSZpq -

D

 

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

Получим азимут А хорды, соединяющей проекции точек Р и Q на эллипсоид. Геодезический азимут Арс} хорды D отличается от азимута ApQизмеренного отрезка D вследствие отличия плоскости PQOq, в которой лежит D , от нормальной плоскости PQOp, про­

181

ходящей через нормаль к эллипсоиду в точке Р и точку Q. Поэтому для редуцирования азимута можно воспользоваться поправкой δ2

(6.4) за высоту наблюдаемого предмета. Согласно рис. 6.3

qq' = H

sine = Η

e (NasinBq - N p sinBD)

 

A cos В

QOp

— cos B„ - H,

 

 

 

4

4

Q ° p

QOp = Np·,

pq = D,

для поправки в азимут получим

 

 

 

 

qq’ sin Aq

HqA sin Aqcos Ba

 

 

Apq ~ ApQ ~ δ 2 -

(6.33)

PQ NqD

Точная формула редукции азимута A pQ пространственного от­ резка к азимуту Apq хорды эллипсоида имеет вид [20]

sin(APQ - Apq) = -sin δ2 =

= - e 2

cos Bp sin APQ(cos Bpcos APQ+sin Bpctgzpq).

(6,34)

Если исключить из этой формулы зенитное расстояние, то

 

 

Hqe2Acos Bp

 

 

sin δ, =

 

-sin Apn sin^n/J.

(6.35)

2

N q(N + H )qcosBqsm(Lq - L pf ......PQ

pr

 

Редуцирование по этим формулам требует последовательных приближений, поскольку в правую часть формулы (6.35) входит редуцированный азимут А . Поскольку поправка <5, мала, выраже­ ние (6.35) можно упростить. Положим Np = Nq = (N + H)q= N, Apq =

= ApQ= A, COSBp= cosBq, тогда

 

sin52 = —M sintf - s in £

)— sin2 A

 

N

'

sin(Lq - L p)

Ho

sinBq-sinBp=(Bq-B p)cosBm,

 

CQCJ

Bq-B p=—^ D , sin(Lq - L p) =

_ r

r _ DsinA . e e

 

 

= Lq - L p = N w s B ’ sin °1 ~ °2’ поэтомУ ПРИ сделанных предполо­ жениях формула (6.32) совпадает с (6.7).

Формулы (6.19), (6.26) и (6.32) - (6.35) решают задачу редуци­ рования пространственного отрезка к хорде эллипсоида.

182

Редукция базиса к поверхности эллипсоида. Базис в триангуля­ ции является ломаной линией, звенья которой - пролеты D бази­ са - пространственные отрезки малой длины (при измерениях ба­ зиса обычно использовали инварные проволоки длиной 24 м). На рис. 6.4. изображен один пролет базиса.

Рис. 6.4. К редуцированию базиса на поверхность эллипсоида

Базисные измерения сопровождали геометрическим нивелирова­ нием и получали проекции ΔΖ измеренных отрезков D сначала на уровенную поверхность, а затем находили длину А1оотрезка РР2, па­ раллельного эллипсоиду. После этого проектировали базис на эллип­ соид и находили хорду Δ&. Получим формулу для вычисления АЬ.

Так как точки Р и Pj близки между собой, можно не принимать во внимание непараллельность нормалей к эллипсоиду, проходящих через эти точки, и считать треугольник PPjP2 прямоугольным. Тогда

Alo = M)cos(a - ϋ 41) = AD(cosa COST}41' - sina sin^O»

где a - угол наклона; ϋΑΓ - составляющая астрономо-геодезическо­ го уклонения отвеса в азимуте отрезка РРГ

183

Согласно рис. 6.4 произведение ADcoscc равно А1, а D sina - это превышение Ah. Положив cosifu= 1, s in ta s ϋ 4Γ, получим

AIо = A I - iVrAh.

Хорду Ab найдем из пропорции

А1о/ (Я + R) = АЫ R,

откуда

Ab = Alo - Alo(H /R - (Я /Я)2),

где Я -геодезическая высота пролета; R - радиус кривизны эллип­ соида.

Введя в последнее выражение значение Δ/ο, получаем формулу для редуцирования пролета базиса к хорде эллипсоида

Ab = AI - bArAh - A/(Я /R - (Я /Я)2).

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

Я

\2

Ъ= ^ А / - ^ $ ^ гА /г - ^ ^ А / + ^Г

Δ/,

V R

/

или, если использовать средние вдоль базиса значения уклоне­ ния отвеса и Нтгеодезической высоты,

b ^ A l - ^ A h - ^ l ^ l - ^

(6.36)

Согласно формуле (6.36) для определения длины b базиса на эллипсоиде нужно знать геодезическую высоту и уклонение отвес­ ной линии.

Внастоящее время базисы не измеряют. Формулы (6.26), (6.30)

и(6.36) используют для сравнения измеренного ранее базиса и на­ клонной дальности.

184