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

Фотограмметрія Дорожинський

.pdf
Скачиваний:
61
Добавлен:
28.06.2022
Размер:
15.32 Mб
Скачать

Професор І. Антіпов довів, що цього негативного явища можна позбутись, якщо до рівнянь взаємного орієнтування додається умова рівності координат зв'язкових точок.

2.7.5.Сумісне визначення елементів взаємного орієнтування та передавання масштабу для двох суміжних моделей

Нехай визначені елементи взаємного орієнтування j-ϊ моделі за формулою (2.83) та розв'язана пряма фотограмметрична засічка, тобто в у-й моделі визначено координати точки А (початок координат у точці S2). Необхідно одночасно знайти елементи взаємного орієнтування (/+1)-ї моделі та виконати приведення масштабу (/+1)-ї моделі до масштабу попередньоїу-ї моделі (рис. 2.29). Як відомо, визначення елементів взаємного орієнтування виконується на підставі виміряних поперечних паралаксів, а передавання масштабу виконується за зв'язковими точками. Математичне формулювання задачі полягає у використанні рівняння (2.83) взаємного орієнтування та рівності векторів

S2A — S2A!,

тобто

(2.85)

S3

Si

А'

Рис. 2.29. З'єднання двох моделей за зв'язковими точками

60

Строго кажучи, рівняння (2.85) є нелінійним відносно шуканих елементів г, ν, Δα, БА, ΔΑ:, Δ6 І тому треба було б застосувати описаний вище прийом лінеаризації (розклад функції у ряд Тейлора).

З практичного погляду, враховуючи, що \R\»|z|

та те, що розв'язання задачі

відбувається ітеративно, перепишемо (2.85) так:

 

^Sb

+ Zj - Zj+i = vb,

(2.86)

або

a b - Sb + lb=Vb.

(2.87)

Отже, для всіх точок (/+1)-ї стереопари складається система рівнянь (2.83), до неї долучаються рівняння (2.87), записані для зв'язкових точок, і одночасно розв'язується система нормальних рівнянь з 6 невідомими

Χ = [δτ δν δΑα δΑω δΑκ 0 b f .

(2.88)

Уразі застосування МНК послідовність розв'язання така:

1.Для кожної зв'язкової точки розв'язують пряму фотограмметричну засічку

у7-й моделі.

2.Задаються наближеними значеннями елементів взаємного орієнтування та

базису в (/+1)-й моделі

ν°, Δα 0 , Αω°, А к b ° .

3.Для кожної пари виміряних точок обчислюють коефіцієнти та вільні члени (2.74), (2.75).

4.Для кожної зв'язкової точки розв'язують пряму фотограмметричну засічку

у(/+1)-й моделі, використовуючи спочатку наближені (а потім уточнені) параметри (елементи взаємного орієнтування та базис).

5.Для кожної зв'язкової точки складають рівняння (2.87).

6.Складають нормальні рівняння, отримують розв'язок (2.88).

7.Обчислюють уточнені значення невідомих

Γ = + (<fr)( 0 , V = v° + (δν)Μ,

 

 

Δα = Δα 0 + (<5Δα)(,), Αω = Αω° + (δΑω)(ί),

(2.89)

ΔΑ: = ΔΑ:0 + (δΑκ)\ b = b° +

.

 

8.Перевірка ітераційного процесу на збіжність аналогічна до (2.81), тільки додається перевірка

9.Оцінка точності розв'язання задачі аналогічна до (2.82), тільки в знаменнику стоятиме величина (п-6).

61

Зазначимо, що геометричну сутність рівнянь (2.85) вперше використав проф. А. Скірідов, назвавши їх "умовою бокових сторін". Він також запропонував використати цю умову для уточнення кутового орієнтування моделі на універсальних стереоприладах.

Як вже зазначалось, проф. І. Антіпов застосував цю саму умову у побудові мережі маршрутної фототріангуляції методом залежних моделей. Таке одночасне розв'язання задачі визначення елементів взаємного орієнтування та передавання масштабу принципово змінило проблему "невизначеності елементів взаємного орієнтування". Оскільки детермінант сукупної системи нормальних рівнянь відрізняється від нуля, то розв'язок задачі є завжди стабільним.

2.7.6. Геодезичне орієнтування фотограмметричної моделі

У результаті розв'язання прямої фотограмметричної засічки, описаного в п. 2.7.1, отримують просторові координати точки А моделі об'єкта (рис. 2.25). Якщо система SXYZ є абсолютною (геодезичною), і кути нахилу знімків та базису визначені в цій системі, то зрозуміло, що модель будується в абсолютній системі координат.

Проте найчастіше система SXYZ є відносною, і тому треба знайти елементи

геодезичного орієнтування моделі:

 

 

 

 

 

υ = [х

γ

Ζ

а0

ω0 к0

t]T.

(2.90)

Згеометричного погляду ця задача складається з трьох дій (рис. 2.24):

-перенесення початку координат з центра ОΦ у центр О;

-"поворот простору";

-зміна масштабу моделі, тобто приведення його до величини t. Формули перетворення координат є такими:

Yr =

+ A

Y

 

 

ζ

Оскільки функції (2.91) є нелінійними відносно вектора (2.90), то приведемо їх до лінійного вигляду:

ЯА

+ bx5Yo + сх&0

+ dxSa0

+ βχδω0 + fxStc0

+ gxSt

+ / , =

ν,,

aySX0

+ bySY0 + cy5Z0

+ dySa0

+ eySa>0 + /у0

+ gy8t

+ ly=vy9

(2.92)

azSXQ + bzSY0 +czSZ0+ dzSa0

+ βζδω0 + /ζδκ0

+gz5t + lz = vz.

62

Тут:

 

 

ах=1

ьх=0>

сх=°>

а у = 0,

Ьу = 1,

су =0,

а2= 0,

bz= 0,

c z =l,

dx = -ΔΖ,

 

= -Δ7sin a0,

 

fx=

(a2X - ax Y)t,

= AX,

dy =0,

e

= -AYtgco0

—Ζ

· U

fy=

Φ2Χ

- b, Y)t,

gy = AY,

 

7

COS^o

 

 

 

 

 

dz=AX,

 

ez=A7cosa0,

/z

=(c2X-q7V,

 

gz=AZ,

 

AX = (axX + a2Y + a3Z)t,

 

 

 

lx=X0+AX-Xr,

 

AY = (blX + b2Y + b3Z)t,

 

 

 

ly=Y0+AY-Yr,

 

AZ = (cxX + c2Y + czZ)t,

 

 

/z

= Ζ 0 +

Ζ - Ζ Γ .

 

Урівняння (2.92) входять сім невідомих. Для однозначного розв'язання задачі треба мати дві планово-висотні опорні точки (вони дають шість рівнянь) та одну висотну (дає одне рівняння). Тоді задача розв'яжеться без контролю.

Найчастіше беруть чотири планово-висотні точки, розташовані по кутах стереопари, і розв'язують задачу, використовуючи МНК.

Знаходження вектора (2.90) є ітераційним і принципово не відрізняється від описаного у п. 2.7.2 чи 2.7.4.

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

2.8.Геометричні перетворення, вживані у фотограмметри

Уфотограмметрії доводиться оперувати такою математичною категорією, як перетворення простору. Найуживанішим є проективне перетворення двоабо тривимірного простору. Подамо в стислій формі найпоширеніші у фотограмметрії типи перетворень.

1.Проективне; має вигляд:

у' = Ату,

(2.93)

де А - матриця перетворень, (detv4*0); у - вектор початковий; у - вектор перетворений.

2. Афінне; має вигляд:

х' = Ах + Ь,

(2.94)

де А- матриця перетворення; х,х\Ь- вектори-стовпчики.

Ї

і

63

Властивості цього перетворення: площина, лінія і відрізок переходять відповідно у площину, лінію і відрізок; не порушується паралельність прямих і площин; зберігається метрика векторів; трикутник переходить в трикутник.

Особливий випадок цього типу перетворення:

x' = RAx + b,

(2.95)

за умови, що

А~1Т,

тобто А - ортогональна матриця, а вектор R є діагональною матрицею, тобто

Rn

R =

R22

 

Rmr?

3. Афінне у разі збігу осей координат; має вигляд:

χ*= Ах,

дейсїА

4. Афінне зі збереженням подібності; має вигляд:

х' = кАх + Ь,

д,ек>0,А~1т.

Властивості цього перетворення:

-кути зберігаються;

-зберігається подібність фігур;

-зберігається пропорційність відрізків;

-не порушується перпендикулярність векторів.

5.Афінне подібне середнє:

(2.96)

(2.97)

(2.98)

х' = кАх,

(2.99)

тобто отримане з (2.98) за умови Ь=0. У цьому разі зберігається початок системи координат.

6. Гомотетія:

xf = Ax,

(2.100)

тобто з (2.99) к=1. 7. Ізометрія:

х' = Ах + Ь.

(2.101)

Тут зберігаються довжини векторів.

64

8. Ізометрія середня:

 

х' = Ах.

(2.102)

Тут збігаються початки систем координат.

 

9. Поворот простору:

 

х' = Ах,

(2.103)

де A~ l =A T , det = 1.

 

10. Перенесення:

 

х' = х + Ь,

(2.104)

А = Е,к = 1.

 

Властивість: Змінюється початок системи координат.

 

11. Симетричне відображення:

 

=

(2.105)

А - Е , к = 1, випливає з (2.99).

 

Властивість: Змінюється орієнтування у просторі.

 

12. Ідентичне перетворення:

 

X' = JC,

(2.106)

тут зберігаються всі елементи.

Звернемо тепер увагу на деякі аспекти геометричних перетворень. Спочатку розглянемо перетворення площини Ε у площину Ρ (рис. 2.30).

Для визначення положення точки а на площині Ρ прийнята довільна (не прямокутна) система координат 0'ххх2. Для так званих невласних точок (перетин двох паралельних прямих) ці координати не дають змогу визначити їхнє розміщення. Тому у проективній геометрії положення точки визначається його напрямком. Напрямок вектора визначається співвідношенням його складових к : т . Для будь-якої точки а (рис. 2.30) визначальним є відношення:

Уі'У2-Уз=хі-х2:1-

(2-107)

Έ

65

Координати yl9 у2, у3 називають однорідними координатами точки а у системі О'Х\Х2 · Співвідношення (2.107) означає пропорційність складових к.

Для кожної точки площини можна записати

 

кУі у2 у^ = хх2Л.

(2.108)

Звідси можна зробити висновок, що системі однорідних координат у (л+1)- вимірному просторі може відповідати довільна система координат и-вимірного простору.

З (2.107) маємо:

 

jc 1 =A 5

j c 2 = ^ - .

(2.109)

 

Уз

Уз

 

 

Для тривимірного простору необхідно взяти чотири числа, тобто (2.107)

запишеться так:

 

 

 

 

 

Уі'>У2:Уз:У4=хі2з:1>

(2.110)

 

кух '· ку2 : ку3 : ky4 = х\ '·х2 '· хз : 1 ·

 

Якщо врахувати (2.93) - лінійне перетворення, то для

четвірки чисел

запишемо:

 

 

 

 

 

= а\\Х1 + а\2х2 +

α 13*3 + а\4х4>

 

ку2

= а2\х\ + α22χ2 + а23*3 + *24*4>

(2-111)

ку3

= ЛзЛ + а32х2

+ а33х3 +

а34х4,

 

кул

= ахх + а42х2

+ а43л;3 +

а^х4.

 

З (2.111) та (2.110) легко отримуємо формули для обчислення координат при геометричних проективних перетвореннях.

Для тривимірного простору:

кУ\ =

а\\х\

+ **12*2 +

+ ^14*4

 

ку4

#41*1 + ^42^2 + а43х3 + а44Х4

 

кУ2 =

а2\х\

+ д 22*2 + а23х3 + а24х4

(2.112)

ку4

а4\х\ + а42х2 + а43х3 + а44х4

 

кУЗ =

а3\х\

+ а32х2 + **33*3 + **34*4

 

ку4

а4\х\ + а42х2 + а43х3 + а44х4

 

Ставиться вимога, щоб для матриці перетворень det А * 0.

66

Рівняння (2.112) можна спростити, розділивши чисельник і знаменник у (2.112) на 044*4.

Ввівши позначення для величин

 

 

 

α44

Л44 Х4

 

Q44χ4

 

 

- Μ -

= Αηζ

іт.д.,

(2.113)

отримаємо:

 

 

 

 

 

 

^Axxx

+ Al2y + Ax3z + Ai4

 

Уі

A4Xx + A42y + A43z + l

 

у

= A l l X +

Л г 2 У +

+ А і 4 ,

(2.114)

 

А41х + А42у + Α43Ζ +1

 

у3 зхх + А32у + Α33ζ + А34 .

 

 

А41х + А42у + Α43Ζ +1

 

Для двовимірного простору

(перетворення площини

на площину) маємо

у3 = 0, ζ = 0 і відповідно

 

 

 

 

 

 

 

Аххх + Ах2у + Ахз ·

 

 

 

Ах + А42у

+1

 

 

у2==А21х

+ А22у + А23

( 2 п 5 )

 

 

Ах + А42у + \

 

Для одновимірного простору (перетворення прямої на пряму) відповідно запишемо:

( 2 Л 1 6 )

АХ + 1

Серед описаних вище перетворень назвемо декілька найвживаніщих.

Перетворення тривимірного простору на тривимірний (3D * 3D) у випадку евклідового простору (прямокутної системи координат). У цьому разі матриця А з (2.103) є ортогональною, а величини у4 та х4 з (2.112) відсутні. Тоді (2.103) запишеться так:

ух хххх

Х2х2 +aX3x3,

 

Уг =а2\х\

+ а22х2 + ^23*3 >

(2л 17)

у331хх32х233х3.

 

67

Матриця

аи

°12

ап

 

А = а21

а22

°23

(2.118)

_°31

°32

а 3 3

 

є матрицею напрямних косинусів і описує поворот однієї системи координат відносно іншої.

Про цю матрицю ми вже подали відомості у π. 2.3.2. За формулами (2.7) обчислюють компоненти матриці напрямних косинусів, а позначення літерами

може бути таке, як в (2.8) або в (2.118).

 

 

 

 

 

 

 

 

Перетворення тривимірного простору на двовимірний (3D

* 2D) (англ.

Direct Linear Transformation — DLT).

 

 

 

 

 

 

 

 

Тоді уз з (2.114)

відсутнє і маємо

 

 

 

 

 

 

 

 

 

 

 

_ А п Х + AnY + Al3Z + Al4

 

 

 

 

 

 

Уі = А41Х + A42yY + A43Z +1

 

 

 

 

 

 

 

 

_ А2\Х + a22Y + A23Z + А24

 

 

(2.119)

 

 

 

 

" А41Х + А42у Υ + Α43ζΖ +1

 

 

 

 

 

 

 

 

 

 

Це перетворення з 11 невідомими параметрами.

 

 

 

 

Перетворення площини на площину

(2D *

2D)

як афінне перетворення зі

збереженням подібності (2.98).

 

 

 

 

 

 

 

 

 

Тоді маємо, якщо к=1

у{ = апХ + al2Y + l\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—і

 

τ

 

у2 = a2lX

+ α22Υ + Ь2 .

 

 

(2.120)

 

 

 

 

 

 

 

 

 

 

 

 

Оскільки А

,

то

матриця

А є

матрицею

напрямних

косинусів

і є

функцією тільки

одного

кута

κ повороту

системи координат. У

такому

разі

ап = а22, а2\ = -аі2

і

бачимо, що це

є

перетворенням з чотирма

невідомими

параметрами.

 

 

 

 

 

 

 

 

 

 

 

 

 

Перетворення площини на площину (2D *2D) з незбіжністю масштабів за

осями координат. Це перетворення описується рівняннями (2.95) - (2.96). Тоді

 

 

 

 

 

 

y ^ R ^ X

+ a ^

+ by

 

 

 

 

 

 

 

 

 

y2=R2(a2lX

 

+ a22Y) + b2

 

 

 

або ж, якщо R\=R2=l

 

 

 

 

 

 

 

 

 

 

(2.121)

 

 

 

 

 

ух пХ + aX2Y +

 

 

 

 

 

Уг =a2lX + a22Y + b2

і це є перетворення з шістьма параметрами.

68

Перетворення Гельмерта описується рівняннями

 

0

1

2

(2.122)

у = Ь0+ а2Х + α{Υ

 

і є чотирипараметричним двовимірним перетворенням, виконує зміну масштабу, "поворот простору (площини)" та паралельне перенесення.

Нелінійні перетворення у двовимірному просторі (зокрема такі, що використовуються у цифровій фотограмметрії при допасуванні оптичних щільностей для пікселів - див. т. зв. операція Resampling).

Одним з них є білінійне перетворення типу

 

х = аь+ОуХ + a2Y + α3ΧΥ,

і

y = b0+blX + b2Y + b3XY.

 

Це перетворення з 8 невідомими параметрами.

Інші типи нелінійних перетворень (кубічна інтерполяція, поліном Лагранжа, сплайн-інтерполяція) подано у п. 4.6.

2.9. Математичні моделі похибок фотограмметричних вимірів. Вплив кривини Землі

2.9.1. Моделі похибок фотограмметричних вимірів

Фотограмметричні вимірювання супроводжуються значним впливом різних за своєю природою чинників, які викликають похибки у вимірах. У загальному вигляді виділимо такі основні групи чинників:

1)що спричиняють відхилення реальної моделі центральної проекції від ідеальної;

2)що супроводжують вимірювання фотокоординат;

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

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

похибок, запропоновану Ю.В. Кемніцом [8] (адитивна -

дослівно "отримана

підсумовуванням"):

 

£ = Θ + Δ,

(2.124)

де Є - повна, Θ - систематична, - випадкова похибки.

Кожна з вищеназваних трьох груп наявна як у векторі похибок©, так і у векторі похибок .

Проаналізуємо кожну з цих груп окремо.

69