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

Конспект до лекції 3

.pdf
Скачиваний:
7
Добавлен:
28.06.2022
Размер:
612.95 Кб
Скачать

-апріорі прийняте положення, що ЕВО невідомі;

-апріорі прийняте положення про тип функції, яка описує параметри деформації зображення.

Запишемо для точки лівого знімка рівняння (2.12) так:

x

а рівняння компланарності як

= x

0

f

 

 

y = y0

ytл ytп

r

 

x

,

r

 

z

 

fry , rz

=0

(3.18)

(3.19)

Тут ytл - трансформована ордината точки лівого знімка за кутами нахилу лівого знімка

стереопари,

yt

п

- трансформована ордината точки правого знімка за кутами нахилу правого

 

 

 

знімка стереопари.

Доповнимо рівняння (3.18) та (3.19) функціями параметрів деформації зображення ψx та ψy . Тоді

x

ytл

x

0

+

 

 

yy0

ytн

frx + x = rz

+f ry + yл rz

+ y

л

y

 

н

0

,

=0

=0 .

,

(3.20)

Лінеаризація перших двох рівнянь з (3.20) приводить до рівнянь поправок, аналогічних до (3.1), а останнього з (3.20) – до рівняння типу (2.5). В узагальненому матричному вигляді матимемо:

A S + B

x

Г + D

E + L

x

x

 

x

 

A

S + B

y

Г + D

E + L

y

y

 

 

y

 

A

 

S + B

 

Г + D E + L

 

q

 

q

q

q

=Vx

=Vy

=Vq

,

,

,

(3.21)

де δS – вектор поправок до елементів зовнішнього орієнтування знімків,

δΓ -вектор поправок до геодезичних (просторових) координат точок об’єкта, δE – вектор невідомих параметрів калібрування зображення,

Lx , Ly , Lq – вільні члени, обчислені з (3.20),

Vx , Vy , Vq поправки до виміряних величин (координат точки на лівому знімку і поперечного паралаксу),

Ax ,Aу , Bx ,Bу – матриці коефіцієнтів (3.2).

Якщо в задачі приймається, що елементи внутрішнього орієнтування треба уточняти, то матриці Dx , Dy ,Dq формують із частинних похідних, аналогічних до (3.2):

 

 

 

 

D

=

m

 

 

n

 

q

T

 

 

 

 

 

x

x

 

;

 

 

 

 

 

 

x

 

 

 

 

x

 

 

 

 

Dy

= m y

 

 

n y

 

T

;

 

 

 

 

 

q y

Dq = m y m y

 

 

n y n y

q y q y T ;

 

m

 

=

 

xл x0

 

 

; nx = 1; qx = 1;

 

x

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

=

y л y0

; ny = 1; qy = 1;

 

y

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

=

 

yп

y0

; ny΄ = 1; qy΄ = 1.

y

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тепер звернемось до функцій x

л

,

y

л

, y

п

, коефіцієнти яких теж входять в матриці

 

 

 

 

 

 

 

 

 

 

 

Dx , Dy ,Dq .

 

 

 

 

 

 

 

 

 

 

 

Однією з пропозицій (G.Schut, США) є застосування таких поліномів:

 

 

 

 

 

 

ψx = a1x + a2y + a5x3 + a6x2y + a7xy2 ,

 

 

 

 

 

 

ψx = a3x2 + a4xy +a5x2y + a6xy2 + a7xy3.

(3.22)

 

 

 

А. Grun (Німеччина) прoпонує поліноми такого типу:

 

 

 

 

 

 

ψx = a1x + a2y + a7x2y + a8xy2 ,

 

 

 

 

 

 

ψy = a3x2 + a4xy + a5x2y + a6xy2 .

(3.23)

 

 

 

Обидва типи поліномів (3.22) та (3.23) використовують при поперечному перекритті між

маршрутами 20-30%.

 

 

 

 

 

 

 

 

 

 

 

В. Дубіновський (Росія) пропонує такі поліноми:

 

 

 

 

 

 

ψx = c0 + c1x +c2y +c3x2 + c4xy + c5y2 + c6x3 + c7x2y + c8xy2 + c9x4 ,

(3.24)

 

 

 

 

 

ψx = c0΄ + c1΄x +c2΄y +c3΄x2 + c4΄xy + c5΄y2 + c6΄x3 + c7΄x2y + c8΄xy2 + c9΄x4 .

 

 

 

 

В формулах (3.22), (3.23), (3.24) опущено індекс знімка л, п. Зрозуміло, якщо йдеться про

x

л

,

y

л

,то в правій частині рівняння беруть координати хл , ул точки на лівому знімку, а

 

 

 

 

 

 

 

 

 

 

 

 

якщо йдеться про

yп

, то беруть координати хп , уп точки на правому знімку.

Далі знаходять розв’язок за класичною схемою: від рівнянь поправок (3.21) переходять

до нормальних рівнянь, знаходять поправки

S , Г , E

до невідомих параметрів. При

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

І. Антіпов (Росія) запропонував здійснювати вибір кількості членів поліномів з (3.24) на основі аналізу точності калібрування. Пояснимо ідею такого підходу. Спочатку візьмемо повний поліном (3.31) і розв’яжемо систему лінійних рівнянь n - го порядку за МНК. Отримана обернена матриця нормальних рівнянь дозволяє виконати оцінку точності по аналогії з формулами (4.19-4.21). На наступному кроці один елемент, наприклад c0 , виключається і розв’язок отримують для системи порядку (n-1). Почергово вилучають з розв’язку інші елементи, а отриманим вважається той розв’язок, в якому оцінка точності є найкращою.

Зовсім інший підхід запропонував В. Дубіновський (Росія), який назвав свій спосіб – принцип зон. Знімок розбивають на зони (рис.3.9).

 

 

 

y

 

 

 

 

 

 

 

 

1

2

3

. . . . . . . . . . . .

9

 

 

 

 

 

 

 

10

11

12

. . . . . . . . . . . .

1

 

 

 

 

 

 

 

.

.

.

.

.

 

.

.

.

.

. x

.

.

.

.

.

 

.

.

.

.

.

 

.

.

.

.

.

 

.

.

.

.

.

 

.

.

.

.

.

 

Рис.3.9. Розбиття знімка на зони для знаходження параметрів калібрування

Розмір зони для знімка 180х180 мм та їх кількість подані в таблиці 3.2.

Таблиця 3.2.Кількість зон на знімку при калібруванні

Параметр

Точність вимірювань, мм

 

 

 

 

 

0.005

0.010

0.05

 

 

 

 

Сторона квадрата (зони), мм

4.2

8.5

42.5

 

 

 

 

Число зон на знімку

1600

400

16

 

 

 

 

Аналіз показав, що для фототріангуляції достатньо знімок ділити на 9 зон. Для кожної зони визначають поправки х і Δy , вважаючи, що вони в межах однієї зони є незмінними

величинами. Тоді функції ψх , ψy мають вигляд:

 

 

Ψх

i

= хi

, ψy

i

= Δyi ,

(3.25)

 

 

 

 

 

де i – номер зони.

 

 

 

 

 

 

Точка вважається належною до i-ї зони, якщо

 

 

 

X – xi ≤ εx , y – yi ≤ εy ;

(3.26)

тут xi , yi - координати центра i-ї зони, εx , εy -допуск на відхилення, який дорівнює половині розміру зони.

Аналіз показав, що для фототріангуляції достатньо знімок поділити на 9 зон.

В подальшому рівняння (3.21) розв’язують під умовою мінімізації суми квадратів поправок, тобто складають нормальні рівняння та знаходять поправки до всіх шуканих невідомих. Як правило, організовується ітераційний процес, який триває до стабілізації шуканих невідомих за критеріями.

Ефективність фототріангуляції із самокалібруванням перевірялась у різних експериментальних умовах. За даними В. Дубиновського побудова фототріангуляції із самокалібруванням в 1,3 – 3,1 разів є точнішою від фототріангуляції без самокалібрування. Точність побудови не залежить від числа опорних точок в блоці фототріангуляції, а їх кількість може бути мінімальною (для розв’язання задачі геодезичного орієнтування блочної моделі та усунення залишкової деформації мережі).

3.2.4 Деформація фототріангуляційної мережі. Точність фототріангуляції

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

Сумарний вплив похибок на процес фототріангуляційних побудов є досить складним. Найлегше він простежується графічно на побудові висотної фототріангуляційної мережі способом моделей. Якщо кожна наступна ланка під’єднується до попередньої, побудованої з окремої стереопари, то елементарні похибки εі в висотах точок кожної моделі, викликані кутовими похибками α1 , α2 ,…αn , нагромаджуються за законом подвійного сумування:

= n 1 + (n 1) 2 + ... + n .

(3.27)

Графічно це представляється ламаною, подібною до параболи (рис.3.10). Похибки в азимутальному напрямку (кути k) спричиняють аналогічну до рис.3.10 деформацію мережі в плановому положенні. Нагромадження похибок в кутах ω1 спричиняють деформацію мережі у вигляді кручення.

Усунення деформації просторової мережі виконується, як правило, із використанням опорних точок, розміщених за певною схемою, та відповідних апроксимуючих функцій. Визначення координат опорних точок прийнято називати прив’язкою знімків.

Z

 

Z

 

3

 

 

3

Z

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Z

 

1

 

1

 

 

 

 

 

 

 

0

 

1

2

3

Рис.3.10. Деформація висот в фототріангуляційній мережі під впливом похибок поздовжнього кута нахилу моделей.

n

Для маршрутної фототріангуляції найбільш поширеними є такі схеми (рис.3.11).

Рис.3.11. Схеми розміщення опорних точок в маршруті знімків.

∆ - опорна точка.

Для блочної фототріангуляції опорні точки розташовують по периметру блока (рис.3.12).

Рис.3.12. Схема розміщення опорних точок в блоці знімків

Апроксимуючими функціями найчастіше виступають степеневі поліноми, в яких кількість членів залежить від кількості опорних точок, наприклад:

 

x

= a

0

+ a x + a

2

 

 

1

y= b0 + b1 x + b2 y +

z = c0 + c1 x + c2

y + a

3

xy

 

 

b xy + b

3

 

 

y + c

3

xy

 

 

+ a

 

 

x

2

+ a

 

 

x

3

+ a

 

 

y

2

,

4

 

5

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 x

2

 

+ b5 x

3

 

+ b6 y

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ c

 

 

x

2

+ c

 

 

x

3

+ c

 

 

y

2

.

 

4

 

 

5

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.28)

Якщо n ≥ 7 , то зберігаються всі члени в (3.28),

якщо 7 ≥ n ≥ 5 , то зберігають a0…..a4 , b0….b4 , c0…c4; якщо 5 ≥ n ≥ 4 ,тo зберігають a0….a3 , b0….b3 , c0…c3 .

Невідомими в рівняннях (3.28) є величини аі, bi, сі . В цих рівняннях х , у – просторові координати опорної точки, отримані з фототріангуляції після абсолютного (“геодезичного”) орієнтування маршрутної чи блочної моделі (для методу моделей) або ж після завершення ітераційного процесу для методу в’язок. В деяких випадках (переважно для гірських районів) в рівняннях (3.28) фігурують вирази типу схz, тобто використовують ще й висоти опорних точок.

Очікувану точність фототріангуляції можна заздалегідь підрахувати, використовуючи відповідні формули (їх вивід дається в деяких підручниках з фотограмметрії). Для вільної маршрутної мережі середні квадратичні похибки координат точки, що знаходиться в кінці маршруту, будуть такими:

mx

 

=1.2m mq

n

3

, m y

 

= 0.6m mq

n

3

, mz

 

= 0.9

f

m mq

n

3

.

 

 

 

 

n

 

n

 

n

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для схеми прив’язки з 4-х опорних точок (рис.3.10) найбільша похибка посередині маршруту, а середні квадратичні похибки обчислюють так:

(3.29)

виникає

m

x

 

n 2

= 0.3m m

 

n

3

q

 

 

 

 

,

m

 

= 0.15m m

 

n

3

y

q

 

 

n 2

 

 

 

 

 

 

 

, m

zn2

= 0.2

f

m m

 

n

3

 

 

b

q

 

 

 

 

 

 

 

 

 

 

. (3.30)

тут n - кількість знімків, f - фокусна віддаль, b - базис в масштабі знімка, квадратична похибка виміру паралаксу, m - знаменник масштабу знімка.

Для блочної фототріангуляції можна скористатись формулами Ф. Лисенка:

mx, y = 0.75m mq

3(2ns + kn rs) (2kns) ,

mz = (0.25H p)mq

3(2ns + kn rs) (2kns) .

mq - середня

(3.31)

(3.32)

тут n - кількість знімків в блоці, s - кількість знімків, на яких зобразилась точка (залежить від поздовжнього та поперечного перекрить), k - число точок мережі на кожному знімку, r - число опорних точок в блоці, H - висота фотографування, p - середнє значення поздовжнього паралаксу.

Аналіз формул (3.31) приводить до таких висновків:

збільшення числа опорних точок в блоці з 10 до 100 підвищує точність побудови лише на 6-7%;

якщо кількість точок на знімку збільшити з 5-ти до 12-ти, то точність підвищиться в

1.5раза;

збільшення поперечного перекриття з 20% до 60% підвищує точність побудов в 1,8

раза.

Детальні дослідження реальної точності фототріангуляції були досить масштабними в 1975-85р.р., а результати опубліковані в матеріалах Міжнародного фотограмметричного

товариства.