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

книги из ГПНТБ / Проворов К.Л. Радиогеодезия учеб. пособие

.pdf
Скачиваний:
49
Добавлен:
27.10.2023
Размер:
26.9 Mб
Скачать

величину). Скаляр и можно рассматривать как функцию

координат

х, у, определяющих положение точки, т. е.

 

и^и(х,у).

(368)

Если функция (368) однозначна и непрерывна в рассматриваемой области, то при некотором значении щ — const мы получим множе­ ство точек, составляющих линию положения, уравнение которой в прямоугольных координатах будет

и(х,у)

= щ.

(369)

Давая различные значения щ = и0 4- гАи (i = 0,

1, 2, . . .),

получим семейство линий положения, отображающих

скалярное

поле и.

 

 

Рассматривая щ как величину,

измеренную для определения

положения точки М, можем считать (369) уравнением линии положе­ ния, проходящей через эту точку. Так как уравнение имеет два неизвестных, то для определения положения точки необходимо не менее двух линий, принадлежащих к двум скалярным полям. Поля могут быть одинакового или различных свойств. Так, за первое можно взять поле и (например, направление, измеренное на исходной точке), а за второе — поле г- (например, расстояние от той же или от другой исходной точки). Тогда для определения неизвестных координат полу­ чим два уравнения:

и (х, у) = щ;

 

ѵ(х,у) = ѵк.

(370)

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

Расстояние AN между смежными линиями положения, считае­ мое по нормалям к ним, будет определять густоту линий, а.следо­ вательно, и точность графического получения координат точек. Величина AN зависит от разности скаляров Au, принятой при по­ строении сетки, и от градиента. Градиент скалярного поля и (grad и) является вектором, направленным по нормали к линии положения в сторону возрастания скаляра и. Градиент определяется по формуле

т

du

— du - ,

du -

g1 adu = 1

F

v = ^ - l +

^ 7 ,

т. е. модуль его gu будет

 

 

 

330

Здесь Î и / — единичные векторы координатных осей, a ѵ — еди­ ничный вектор нормали. Из формулы (371)-видно, что при малых приращениях

8и~ш.

(372)

т. е. модуль градиента можно рассматривать как коэффициент про­

порциональности между приращением Au функции

и расстоянием

AN между соответствующими линиями положения. Обычно

рассмат­

ривают линии положения величин их,.

и2, и3, . . .

лишь

вблизи

какой-либо точки. В пределах такой,

достаточно

малой,

области

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

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

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

В качестве величин и при построении линий положения можно принимать или измеренные величины, редуцированные на поверх­ ность эллипсоида, или проекции этих величин на поверхность шара или на плоскость. Наиболее просто строить сетку на плоскости в проекции Гаусса—Крюгера в случае, если измеренные величины редуцированы на эту же плоскость. Положение точки будет опре­ деляться также в системе координат Гаусса—Крюгера. Как указы­ валось в § 41, в этом случае можно одной координатной зоной ох­ ватить территорию до 1000 км с востока на запад и неограниченно с севера на юг.

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

331

Уравнение

линии положения

в этом случае в системе

плоских

координат

Гаусса—Крюгера

будет

 

 

 

T, = arctg

V—Ух

(373)

 

 

 

X — Х^

 

где Т[ — дирекционный угол линии положения, хх и ух — коорди­ наты исходной точки, а х и у — текущие координаты линии поло­ жения. Придавая различные значения дирекционному углу Tt в уравнении (373), например 0, 5, 10° и т. д., получим семейство линий положения, представляющее пучок полупрямых, построен­ ных в исходной точке под углами Г, к оси абсцисс (см. рис. 178). Азимутальную сетку на плоскости можно построить при помощи

любого углоначертательного

прибора

(на­

пример, транспортира), а также по коор­

динатам,

для

чего достаточно нанести на

основу

по

одной

точке

для

каждого

луча, если вся сетка размещается на од­

ной основе. В противном случае

для

построения каждого луча необходимо на

каждой основе иметь не менее двух точек.

Координаты

точек

вычисляются

по

фор­

мулам

 

 

 

 

 

 

 

x~x1JrdcosTi

 

 

y = yx + dsin

Т{

Рис . 180

где d — произвольное,

по возможности

 

наибольшее, расстояние.

Азимутальную сетку в системе географических координат, если величинами и являются направления (азимуты) на сфере или на эллипсоиде, наиболее удобно строить по координатам отдельных точек. В этом случае сетка будет представлять собой цучок геодези­ ческих линий. Координаты точек Мх, М2 и т. д. (рис. 180), необхо­ димых для построения данного луча, можно найти путем последо­ вательного решения прямой геодезической задачи, начиная от дан­ ной точки А. Расстояние s между смежными точками назначается из расчета необходимой частоты для построения линии. Прямой

азимут каждого

отрезка следует получать путем прибавления 180°

к вычисленному

обратному азимуту предыдущей линии.

В большинстве случаев точность вычисления координат, необ­ ходимых для построения сеток, может быть сравнительно невысокой. Ошибки определения широт и долгот 6L нетрудно рассчитать по формулам

ЬВ" =

mp"àx

ЮООМ

ЬЪ" =

 

mp"àx

1000ІѴ cos В

332

где Ьх — необходимая точность графического построения, которую примем равной 0,1 мм; m — знаменатель масштаба основы, на ко­ торой строится сетка; M vi N — радиусы кривизны меридиана и первого вертикала эллипсоида для данного района. Положив M — = N = 6,4 • 106 м, получим следующие простые соотношения:

 

 

 

05" = 2>,2т • Ю-6 ;

 

 

 

ÔL"=^8B" sec

В.

Так, для

масштаба

сетки 1 :

100 ООО под широтой 60° найдем

= 0,3"; ÔL = 0,6".

Для

сетки

в масштабе 1 : 2 000 000 допусти­

мые

ошибки

вычисления

координат

составляют соответственно

6,4 и 12,8". Вычисление географических координат с такой точ­

ностью можно выполнять

по приближенным

формулам. Так,

поль­

зуясь только главными

членами формул

со

средними аргументами

B2

= B1 +

 

(l)mscosAm

 

 

 

Л і

= А

±

180° + (2)m tg Bms

sin Am ,

(374)

L 2

= L 1 +

 

(2)msecBmssmAm

 

 

 

можно при расстояниях до 60 км получить географические коорди­ наты на эллипсоиде с ошибками 0,1". В формулах (374) Вт и Ат средние значения азимута и широты. Вычисление по этим формулам производится методом приближений. Географические координаты, отнесенные к поверхности шара, можно получить также по форму­

лам (374), заменив в них (1) и

(2) через

Расстояние s при вычи­

слении

координат

следует

рассчитать

так, чтобы на основе оно

было

в

пределах

5—10

см.

Для

масштаба 1 : 100 000

это соот­

ветствует

5—10 км; для

масштаба

1

: 1 000 000—50 —100

км.

Линии положения для расстояний представляют собой концент­ рические окружности (см. рис. 179), поэтому сетку таких линий называют круговой или стадиметрической. Графическое построение

круговой

сетки на

плоскости

выполняется при помощи циркуля

или

другого аналогичного прибора. Уравнение

линий

положения

в системе плоских

координат Гаусса—Крюгера

имеет

вид

 

 

 

 

 

(x-xxf

+ {y~yxf^d\

 

 

(375)

или

в

параметрической форме

 

 

 

 

 

 

 

 

 

X = хл

4- à, cos

T )

 

 

 

 

 

 

 

 

у = Уі + dt sin

T J

 

 

 

где

di

— расстояние от центра (исходной точки); хх, ух

— координаты

исходной

точки;

х,

у — текущие координаты

линии

положения;

Т — текущий параметр, который можно рассматривать как дирек-

ционный

угол

соответствующего

радиуса.

 

Наиболее просто координаты

точек для

построения круго­

вой сетки

как

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

и географических

333

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

Этот же способ можно применить для построения сетки гипербол, являющихся линиями положения разностей расстояний, измерен­

ных с

определяемой точки до

двух

известных

точек. По­

строив

описанным

способом (в

системе плоских или географи­ ческих координат) на каждой известной точке круговую сетку, соединим плавной линией точки

пересечения

окружностей

раз­

ных семейств, где значения

раз­

ности радиусов равны

величине

гг .

Каждое

значение

rt

даст

соответствующую

гиперболу.

На

рис. 181

изображена

ги­

перболическая сетка,

 

построен­

ная этим методом, для исход­ ных точек А и В. Гиперболы изображены на рисунке жирными линиями; окружности, по которым

строились гиперболы, показаны тонкими линиями.

Для построения по точкам, прямоугольные координаты х, у

точек

гиперболы с

разностью

расстояний г,- до

исходных точек

А (Яц Уі) и В (х2,

у2)

можно

найти по формуле

 

 

 

 

у = а ] / і - ( £ ) 2

(377)

или в

параметрической

форме

 

 

 

 

 

x = btgt; y = asect.

(378)

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

1

а = = 2 Г < >

Ь = \Ѵ {х2-х^

+ (у2-Уіу-гІ

(379)

334

По

найденным описанным способом условным координатам

можно

получить координаты X , Y в

системе Гаусса—Крюгера

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

 

 

 

X = X sin Ѳ + у cos Ѳ +

Х0

(380)

 

Y — у sin Ѳ — X cos Ѳ + Y0 f

 

 

где

X0 (xx -|- x2)

tgO: У2 — Vi

х2 — хг

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

 

В = В 0 - ^ Р % [ 1 1 - ^ ( 5

+ 3 ^ ) ]

 

L =

L 0

-

 

(381)

 

 

 

 

 

 

 

N0cos Во

 

где L Q

— долгота

осевого

меридиана; В0 — приближенная широта,

найденная по абсциссе

точки, считаемой

длиной дуги меридиана

от экватора, и t0

=

tg

В0.

 

Все

изложенное

относительно построения гиперболических се­

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

оси и к центру,

имеет вид

 

 

 

 

Я.2 ~

Ь2

'

(382)

а параметрические

уравнения

будут

 

 

X = a cos

t

(383)

 

y =

bsint

 

 

Азимутальную, круговую, гиперболическую или эллиптическую сетку в любой проекции эллипсоида можно построить, если предва­ рительно на соответствующую основу (или карту) нанести достаточно частую сетку координатных линий в системе прямоугольных ко­ ординат Гаусса—Крюгера. Вершины сетки можно наносить по гео­ графическим координатам, найденным с помощью таблиц, пред­ назначенных для расчета координат вершин съемочных трапеций. Линии положения наносятся в этом случае или по прямоугольным

335

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

Можно построить сетку в проекции Гаусса—Крюгера так, что местоположение точки можно определять величинами, не реду­ цированными на плоскость. Например, построив азимутальную сетку в системе прямоугольных координат, как описано выше, повернем ее на основе (или карте) на угол сближения меридианов у в начальной точке, взятый с обратным знаком. Оцифровка линий сетки при этом не изменяется. Поправка за кривизну изображения геодезической линии ô искривит линию положения незначительно, поэтому в большинстве случаев нет необходимости ее учитывать.

Для построения аналогичной круговой сетки в расстояния от исходной точки по линии положения с номером st (s, — расстояние нередуцированное на плоскость) вводится поправка As по формуле

As = slVm

2 Я 2

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

Для построения гиперболической сетки в системе Гаусса—Крю­ гера, рассчитанной на определение точек по разностям, не редуци­ рованным на плоскость, при построении линии сетки в каждую раз­ ность расстояний г вводится поправка Ar, которая находится по формуле

 

Аг =

d 2

( j h + Л У

dj

( Уі + у

y

 

2 Д 2

V

2

)

2Д2

V 2

/

где аг

ж d2 — расстояния

до

исходных

точек. Следовательно, для

каждой

точки одной

и

той же

линии

положения поправка будет

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

Вообще, для любой сетки, строящейся в проекции Гаусса— Крюгера, параметрические уравнения линий положения можно за­

писать в виде

 

 

 

 

 

 

 

(384)

где

и — измеренная

величина, редуцированная

на

эллипсоид, а

Au ее редукция на плоскость. Если координаты хну

вычислять

по

аргументу +

Au), присваивая ему целые

равноотстоящие

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

336

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

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

§ 47. О Ш И Б К А О П Р Е Д Е Л Е Н И Я П О Л О Ж Е Н И Я Т О Ч К И

При помощи линий положения нетрудно найти общие формулы ошибок уравненных координат определяемой точки. На основании формул (371), а также, согласно рис. 182, коэффициенты а и Ъ в урав­ нениях погрешностей (364) можно записать в виде

 

 

дщ

 

йщ

dN і

• gl sin Ѳ(

 

 

й і

дх

 

dNi

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

йщ

dNi

— gt

COS Ѳ,-

 

 

 

 

 

dNi

ду

 

 

 

 

 

 

 

 

где g{ — градиент величины щ;

 

 

 

Ѳг — дирекционный

угол

соот­

 

 

 

ветствующей линии положения.

 

 

 

Считая веса измеренных

вели-

 

 

 

чин

равными — - и

обозначая

 

 

 

 

 

 

 

(385)

 

 

 

напишем нормальные

уравне­

 

 

 

ния,

соответствующие

системе

 

 

 

(364)

 

 

 

 

 

 

 

 

[q sin2 Ѳ] I [g sin Ѳ cos Ѳ] г) 4-

 

Рис .

182

 

£ і = 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-[gsin0cos8] I 4

[gcos2 9] n + L 2 = 0.

(386)

Как известно, вес последнего неизвестного в системе нормальных уравнений равен коэффициенту при этом неизвестном после исклю­ чения из уравнений всех других неизвестных. Полагая рх = р% и ру = рц, из уравнений (386) найдем

рх

= [gsin2 9]-

ig sin Ѳ cos Ѳ ] 2 .

[q C0s2 Ѳ]

'

 

 

p „ = [ ç Cosw­

[q sin Ѳ cos

Ѳ ] 2

[q sinae]

;

 

 

22 Заказ 129

337

Или после необходимых

преобразований

 

_

\дт sin2 vife] .

 

 

 

[q COS2 Ѳ]

 

Ру'-

\дідк Sin Vife]

 

[?sin2 Ѳ]

 

В полученных формулах

 

= Ѳе — Bf — угол между

линиями

положения и,- и иА , причем индексы і и & берутся во всех

комбина­

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

пункта получим следующие

формулы:

 

 

 

2 _

1 _

[g cos2 Ѳ]

 

 

 

мі=

Рх

Ыт sin2

yik]

 

 

 

 

 

Ml

= ± ~ .

[ q s i n

l ö ]

, ;

(387)

y

Ру

[im sm2 yik]

>

V /

M* = M%+Ml

= -t

Щ—r.

 

 

 

,

[<?i<7fe S

l n 2 yik]

 

Найденные формулы являются

общими для ошибок

положения

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

величин длины MS), и направления Мі какой-либо

стороны sr

между исходным и определяемым пунктами нетрудно

найти по тем

же формулам (387). Для этого дирекционные углы линий положения

следует

отсчитывать

от

направления

этой стороны. Тогда

 

 

 

 

S r

{qm sin2

yik\

(388)

 

 

 

 

 

 

 

 

 

 

 

 

[ЯіЯк sin2 yik]

 

где ß„- = Qt — Tr

— угол,

образованный линией

положения изме­

ренной

величины

щ

со

стороной sr.

положения,

соответствующих

Для

получения

формул ошибок

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

1. Для азимута (направления) Т с исходного пункта на опреде­ ляемый [в системе координат с началом в исходном пункте и осью абсцисс, совпадающей с исходной стороной (см. рис. 175, а)] имеем:

Т = arctg У-

где X is. у — текущие

координаты точки на линии положения.

338

'

Отсюда

дТ

дх

дТ

 

ду

 

следовательно,

 

gr = Ç

(389)

Так как величина градиента не зависит

от направления осей

и начала координат, то формула (389) будет справедлива для любых

измеренных направлений, дирекционных углов или углов

между

направлениями на определяемый и твердый пункты.

 

 

 

2. Функция координат расстояния d между исходным и опре­

деляемым

пунктами в той же частной

системе

координат

(см. рис.

175, а) имеет

вид

 

 

 

 

 

 

 

откуда

 

 

 

d=\rrx*

+ y*,

 

 

 

 

 

 

dd

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх

d I

 

 

 

 

следовательно,

 

 

ду

d •

 

 

 

 

 

 

 

gd = l.

 

 

(390)

 

 

 

 

 

 

 

 

 

3. Функция

координат

вершины угла ß при определяемой

точке

на

две исходные

точки [в системе

координат

с началом

в точке О

и

осью ординат,

совпадающей

со стороной AB (см. рис. 175, а)]

будет

 

 

 

 

С

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У+~7Г

 

У—je

 

 

 

 

 

 

 

ß - a r c t g —

 

a r c t g — j - .

 

 

 

Откуда

найдем

 

 

 

 

 

 

 

 

 

 

 

_

У + 2

. Ѵ

2

 

 

 

 

 

 

 

~дх~

d\

 

d\ '

 

 

 

 

 

 

 

 

<?ß

X

X

 

 

 

 

 

Следовательно,

~ду ~~~dj ~ ~dj'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d\ + d\-2хЪ-2уг

+ ~-

 

 

 

 

 

 

 

^ß==

 

d*d\

 

 

 

 

Но из

треугольника

ABM:

 

 

 

 

 

 

 

 

 

dl + d%-2x*-2y*

= f ,

 

 

 

22*

339

Соседние файлы в папке книги из ГПНТБ