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

120104.62 Конспект лекций

.pdf
Скачиваний:
225
Добавлен:
27.03.2016
Размер:
1.7 Mб
Скачать

 

X a ( X X

 

) b (Y Y ) c (Z Z

 

)

 

 

1

S

1

S

1

S

 

 

 

Y a2 ( X X S ) b2 (Y YS ) c2 (Z ZS )

x f

X * Z a ( X X ) b (Y Y ) c (Z Z )

Z *

S

3

S

3

S

 

 

 

3

 

 

 

 

 

 

 

U

 

U V UV

 

 

 

 

 

 

 

V 2

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

a Z * a X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

Z

*

 

 

 

fa

 

f

 

X

*

 

a

 

fa

 

 

a x

 

fa a x

 

 

x

 

 

 

 

 

 

 

 

f

 

1

Z * 2

3

 

 

1

 

 

 

 

 

 

1

 

 

3

 

 

1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

X S

 

 

 

X S

 

 

 

 

 

 

 

 

 

 

 

 

Z *

 

 

 

Z * 2

3

 

 

Z *

 

 

 

Z *

 

 

 

Z *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

b Z * b X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

Z *

 

 

fb

f

 

X *

 

b

 

fb

 

 

b x

 

 

fb b x

 

x

 

 

 

 

 

 

 

f

1

 

3

 

 

1

 

 

 

 

 

1

 

 

3

 

 

 

1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z * 2

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

x

YS

 

 

 

YS

 

 

 

 

 

 

 

 

 

 

Z *

 

 

Z * 2

Z *

 

 

 

Z *

 

 

 

 

Z *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

c Z * c X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

Z *

 

 

 

fc

 

с

 

 

 

 

 

 

 

 

 

f

1

Z * 2

3

 

 

1

f

x

ZS

 

ZS

 

 

 

Z *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим пример дифференцирования по

 

 

 

 

 

 

 

X *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

*

 

 

 

Z

*

 

*

 

 

 

 

 

f

 

 

 

 

 

Z

*

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

Z

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

Z * 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1 cos cos sin sin sin a2 cos sin sin sin cos a3 sin cos

b1 cos sin b2 cos cos b3 sin

c1 sin cos cos sin sin c2 sin sin cos sin cos c3 cos cos

X *

 

 

fc

 

c x

 

fc

c x

 

c

 

1

3

1

3

 

 

 

 

 

 

Z * 2

3

 

Z *

 

Z *

 

 

Z *

sin cos cos sin

Продифференцируем

X *

и

Z *

отдельно

 

 

 

 

 

X *

a1

c1

X *

 

a1 X X

S

c1 Z Z

S

 

 

 

 

 

 

 

 

 

cos ( sin ) cos sin sin c1

cos cos sin sin sin a1

с1 X X S a1 Z ZS

41

abc deg

 

 

 

Z *

a3 X X

S

c3

Z Z

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a3

cos cos c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

sin cos a

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z *

c3 X X S

a3 Z ZS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X *

 

 

 

Z *

 

 

 

 

 

 

 

 

x

 

 

Z *

 

X *

 

 

 

d

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

x

 

 

Z *

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z * c X X

a Z Z X * c X X

 

a Z Z

 

 

 

f

1

S

1

S

3

S

3

S

 

 

 

 

 

 

 

 

 

 

 

 

 

Z * 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f c1 X X S a1 Z ZS x c3 X X S a3 Z ZS

 

 

 

 

 

 

 

 

Z *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотренные

ранее

теоретические

операции были выполнены с

целью преобразования строгого уравнения для определения элементов внешнего ориентирования и приведения уравнений к линейному виду путем разложения в ряд Тейлора.

Решение же этих уравнений не рассматривалось.

Выше были получены формулы для вычисления коэффициентов

и lxly .

После этого формируют матрицу коэффициентов

 

 

 

a1

 

b1

c1

d1

e1

g1

 

An6

 

 

a2

 

b2

c2

d2

e2

g2

 

 

 

 

 

 

 

 

 

 

 

 

an

 

bn

cn

dn

en

gn

 

 

 

 

X S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

YS

 

 

 

 

 

 

 

X 61

 

Z S

 

матрица неизвестных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

l2

 

матрица свободных членов уравнений

n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

 

42

 

 

 

v1

 

 

 

V

 

v2

матрица поправок

 

n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vn

 

 

 

n – число уравнений

 

AX L V уравнение поправок

ax X S

bx YS

cx ZS

dx ex gx xвыч xизм vx

ay X S

by YS

cy ZS

d y ey g y yвыч yизм vy

 

Чтобы

решить эти уравнения нужно перейти к нормальным

уравнениям. Для этого формируют транспонированную матрицу AT (столбцы

меняем на строки).

AT AX AT L 0

B66 X C61 0 нормальное уравнение

 

aa

ab

ac

ad

ae

ag

B

ba

bb

bc

cd

ce

cg

66

 

 

 

 

 

 

 

 

ga

gb

gc

gd

ge

gg

 

 

 

 

 

 

 

al

L bl

gl

Решать нормальные уравнения можно любым предложенным в математике способом, но в фотограмметрии предпочитают решать по методу

Гаусса.

X Q 1L

Q 1 – обратная матрица диагональных коэффициентов нормальных уравнений.

В МНК предусмотрена оценка точности решенных уравнений и определенных уравнений.

Сначала вычисляют СКО единицы веса

 

vx vx vy vy

n k

 

n – число уравнений (число опорных точек умноженное на 2)

43

k – число неизвестных

Далее вычисляют СКО неизвестных

mX S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

QX S X S

mY

 

 

QY Y

 

 

S

 

 

 

 

S S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mZS

 

 

QZS ZS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

Q

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

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

чтобы добиться точности используют итерационный метод Ньютона.

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

X S X S0 X S

0

Y Y 0

Y

0

S S

S

 

ZS ZS0 ZS

0

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

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

44

2.7 Влияния рельефа, угла наклона снимка и кривизны Земли на положение точек на снимке

Вывод формул смещения точек из-за влияния рельефа и угла наклона снимка.

n (o) a0 a

снимок

S

 

 

 

H

 

A

 

N’

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

h

 

 

 

a0

 

 

N

средняя плоскость

r0

A0

 

 

 

r

 

 

 

n

 

 

 

 

 

 

 

местность

 

 

 

 

 

 

Рисунок 31 - Смещения точек из-за влияния рельефа

So=f

SN=H

N’N=h=AA0

Введем обозначения AN’=A0N=R na=r, na0=r0

aa0=δrh – смещение точки за влияние рельефа. r и r0 – радиусы векторы.

Из подобия треугольников

na0S и

NA0S

naS и

N'AS

na0 NA0

na N A

r0 H f

 

So

 

r0

 

f

 

R

r0 H

 

 

SN

R

H

f

 

 

 

 

 

So

 

r

 

 

 

f

R

r(H h)

SN

R

 

 

H h

 

 

f

r(H h)

f

r0 H rH rh

r0 r rh

45

rH rh H rH rh

rh rhH

Эта формула позволяет вычислить смещение точки за рельеф, если на снимке известно её положение, n, H и h.

Пример r=100мм

H=1000м h=20м (столб)

r

 

rh

 

100мм * 20м

2мм

 

 

h

 

H

 

1000м

 

 

 

Известно, что если местность рельефная, то точки за влияние рельефа будут смещаться, и изображение не будет подобно местности. Следовательно,

чтобы получить на снимке подобное изображение надо в точки ввести поправки за рельеф.

v

 

0

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

c’

o’

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

0

0

 

 

ot

 

 

0

a

 

 

 

 

 

P

 

 

0

 

 

k

c

 

 

 

 

 

 

 

 

P

v

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 32 – Смещение точек из-за влияния угла наклона снимка Пусть из одной точки фотографирования получены наклонный и горизонтальный снимки.

Р и Р0 – наклонный и горизонтальный снимки соответственно

So=Sot=f

o и ot – главные точки наклонного и горизонтального снимков hchc – линия пересечения Р и Р0

46

а и а0 – изображение точки А на Р и Р0

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

ca=r, ca0=r0 – радиусы векторы определения положения точки на снимках относительно т. с.

Известно, что точки на наклонном снимке смещены за влияние угла наклона снимка либо к т. с, либо от неё.

Из

 

 

 

 

Sa0оt и

 

Sаo’

 

 

 

ac

 

 

So

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0c

 

So

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

Обозначим o ot

f

ak

 

Тогда с учетом обозначений

 

r

 

 

 

f f

 

 

 

 

ac= ac’=r , так как сaс’ равнобедренный

 

r 0

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из

 

 

 

 

акс

 

 

 

 

 

 

 

 

 

f

 

sin

0

 

 

 

 

f

r sin

0

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

f r sin 0

 

 

 

 

 

r 0

 

 

 

 

f

 

 

 

 

 

 

r 0

 

 

 

 

 

 

fr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

f r sin

 

 

 

Из формулы видно, что если 0 0 , то r=r0. Следовательно, r -r0 =δrα, будет зависеть от величины угла наклона и радиуса вектора.

Здесь рассмотрен случай когда радиусы векторы лежат на линии главного вертикала, но это не всегда так, поэтому учитывается угол φ – угол отсчитываемый от линии нулевых искажений против часовой стрелки до r

или r0.

 

 

v

a

 

 

 

 

 

hc

 

 

 

c

 

hc

v

 

 

 

Тогда формула примет следующий вид:

r 0

fr

 

– строгая формула

 

 

f r sin 0

sin

 

 

47

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

Разделим на f

r 0

 

 

 

 

 

 

r

 

 

 

 

если принять

 

 

 

 

 

за а, то зная, что при разложении

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r sin 0 sin

 

 

 

r sin 0 sin

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в ряд

 

 

 

1

1 a , т.о. можем записать следующее

 

 

 

 

1 a

r 0 r(1

 

r sin

0

sin

) r

r 2 sin

0

sin

 

 

 

 

 

f

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r- r0=δrα

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r 2

sin

0

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

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

Влияния кривизны Земли на положение точек на снимке

a0

a

n

P

 

 

δr

 

f

 

 

 

S

 

H

N

A0

E

 

 

 

A

 

R

 

 

O

Рисунок 33 - Влияния кривизны Земли на положение точек на снимке

R – радиус Земли; n – точка надира;

N – точка на местности, соответствующая точке надира;

48

А0 – ортогональная проекция точки А на горизонтальную плоскость E на;

а0 – изображение точки А0 на горизонтальном снимке; r – радиус-вектор от точки надира до точки а;

r0– радиус-вектор от точки надира до точки а0;

δr - смещение точки вызванное влиянием кривизны Земли.

r r r0 ;

rh

 

r 3 H

ф

 

 

 

 

2R f

2

 

 

 

 

 

 

 

.

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

3. Трансформирование аэроснимков. Фотосхемы. Фотопланы.

3.1 Назначение и методы трансформирования снимков. Цифровое трансформирование снимков.

Назначение и методы трансформирования снимков. Оптико-

механическое трансформирование снимков. Цифровое трансформирование снимков. Вывод формул связи координат плоского и наклонного снимков.

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

Существуют два основных способа трансформирования:

1.снимки исправляются только за угол наклона и приводятся к заданному масштабу

2.снимки исправляются за угол наклона, приводятся к заданному масштабу и исправляются за влияние рельефа

49

Теоретически первый способ применим, когда местность плоская и горизонтальная. В действительности такой местности не бывает и практически первый способ применим, когда смещение точек за рельеф не превышает заданного допуска. Например, rh 0.3мм в масштабе карты. Так как масштаб аэроснимков как правило в 2-3 раза мельче масштаба карты,

следовательно допуск на снимке должен быть меньше в k Mm раз. Если

смещение за рельеф превышает допуск, то выполняют ортотрансформирование.

Способы трансформирования:

аналитическое;

фотомеханическое;

цифровое трансформирование;

ортофототрансформирование.

Аналитическое трансформирование.

Как известно,

X X S

(Z Z S )

a (x x

0

) a

2

( y y

0

) a

3

f

1

 

 

 

 

 

c1 (x x0 ) c2 ( y y0 ) c3

f

 

 

 

 

 

 

b1 (x x0 ) b2 ( y y0 ) b3 f

 

 

(Y YS ) (Z Z S )

 

 

c1 (x x0 ) c2 ( y y0 ) c3 f

 

 

 

 

 

 

 

 

 

Выведем формулу вычисления координат точек местности для горизонтального снимка.

Если α = = æ = 0, то a1 = b2 = c3 = 1, a2 = a3 = b1 = b3 = c1 = c2 = 0.

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

X X S

(Z Z S )

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

y 0

 

 

(Y YS ) (Z Z S )

 

 

f

 

 

 

 

 

 

Так как левые части уравнений равны, приравняем их правые части и решим относительно x0 и y0:

50

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