Дорожинський О.Л., Тукай Р. Фотограмметрія

Подождите немного. Документ загружается.

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

внутрішнього та зовнішнього орієнтування та просторових координат точок місцевості,

взятих

з ЦМР.

Математичною основою слугує раніше наведене рівняння колінеарності,

яке для наочності подано ще раз:

де Χ,

Υ,

Ζ - координати центра піксела на ЦМР у вибраній просторовій системі

координат;

, у

, f - елементи внутрішнього орієнтування;

- лінійні елементи зовнішнього орієнтування;

,Ь

,с

- напрямні косинуси, обчислені за кутовими елементами зовнішнього

орієнтування α,ω,χ.

Зазначимо, що вихідним є цифровий знімок, де кожний піксел має коорди-

нати х

, у

; а у рівняннях (7.1) фігурують плоскі прямокутні координати точки

х,у. їхній взаємозв'язок подано у π. 4.3. Параметри такого переходу є найчастіше

коефіцієнтами поліноміальних функцій і є величинами відомими.Саме їх використо-

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

(х,

у) до циф-

рового зображення (х

, у

). У літературі рекомендується залежно від типу цифрового

знімка використовувати:

для

метричних знімків афінне, білінійне

або

проективне перетворення;

- для неметричних знімків подібне перетворення.

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

передати відповідну щільність (колір) від оригінального зображення. Про цю операцію

ми згадували у п. 4.6. (англ. "resampling") та вказали на кілька інтерполяційних

прийомів визначення оптичної щільності.

На четвертому етапі необхідно виконати "зшивання" окремих ортозображень

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

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

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

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

опрацювання. Найчастіше за остаточне значення приймається середнє арифметичне

із

100

- 200 пікселів, що потрапили в зону опрацювання. Складніші алгоритми дають

змогу плавно змінювати оптичну щільність, але їхнє застосування потребує більших

затрат машинного

часу.

Наведений загальний підхід у різних програмах має конкретну

реалізацію і не обов'язково повинен повторювати описану вище схему.

Нижче подаємо короткий опис створення цифрової ортофотокарти на ЦФС

фірми "Intergraph".

У~Уо =-/

-*Ч)

~ f

(Χ-Χ

)

^(Υ-Υ

)+€

(Ζ-Ζ

(X-X

)+b

(Y-Y

)+c

(Z-Z

)

(X-X

)+b

(Y-Y

)+c

(Z-Z

(7.1)

260

7.5.2. Γенерування ортозображень

Вихідні дані беруть з даних проекту (рис 7 ЗО)

ЯВІ

ФЇ^^ТІЇТ

Рис 7

ЗО

Вікно з даними для цілого проекту

New OrthoPro Project

идаШОПоІвСІ гІВІШіЕ

|ХЛТМСЕ0102mgeSproiek^ortoV0102mge op|

Шл

Рис 7 31 Вікно з базовими даними

261

Далі з вікна, в якому є ярлик OrthoPro Project Filename, вибирають файл, який

надалі буде базовим (рис. 7.31).

Окрім назви файла, необхідно вибрати систему для проекту загалом. Є три

варіанти:

а) Use Current Geomedia Workspace Coordinate System, із вибором системи

координат з поточного проекту Geomedia;

б) Coordinate System or Design File через вибір файла, з якого ці дані можна

отримати;

в) Define Explicitly, тобто вручну.

Далі визначаємо:

- рисунок з границями проекту (рис. 7.32);

-розмір піксела (рис. 7. 33);

- розміщення фотограмметричного проекту;

-розташування ЦМР (рис.7. 34);

- формат вихідних файлів: тип компресії, розмір файла (рис. 7.35).

ЩЯЩшЩЩ: щя

жіДхІ

ГіфіІ Імм ТшИпіифІІ ГЬшДшіГК Л fli|| Τ РюйосЗД

Cow. ЇХ (72636392 »

ЬтШЩш

C«»1Y|eiaW35e т

ШіШшшШют

Свм«?Х 177339799 «

ViRtfCMto» 1

Coew2Y|57042342 »

ш* <*$*«$« І

I^iiHBReWiort 1

IMSta іиооош» «

МММІІ 1

<* І

гятпівт· - "* Ш1——Г~ ЮВГ

Рис. 7.32. Вихідні дані для ортоперетворень: границі проекту

262

Рис. 7.33. Вихідні дані для ортоперетворень: розмір піксела

ттм*

_____

РГИЛІ , .<v.;·.

•MfW Щ\ - - - --.·; , V*: ^Ійї'і:: 4

І&М&М

jPraiKten *м*.*по*НРІ [763296.33.98782133

В» вч-. Jb f ,

^ /Ν

t* 1 ] Ргайк* Ага· - оіїсфгфй

t* riSebcMProtbcjs -atteerotect

Ifc

В Snfrct»d Product! l«twrt-ortcprt»d

tfc QS

АуШШ Products tjibft* - gctopro^ct

tfc : : Orthw - опоргоікл

k S)OiftoiUMi.ataNid

^ BMtcA fmiini »(MmiKi

t^ Hi

Dtwtton FcuQiiH*

UMt

ortni* ufrrt

Рис. 7.34. Вихідні дані для ортоперетворень: рисунок розміщення ЦМР

263

ІшІїШІІІІІІІІі

Intapototion jJ

Sp^bPjgate j32 Zl ?'

-VMBMteiC

Г U»£wtao*0M«fe№ |150

Γ UttAvataMref QotMtDTM

ΡΊ

Htfta^ajpctoti*·*)

f* fl^W '

(

Γ CfpMdilma9»£d9»

'MltiM^AdMlRirtt ·

NofewFMFlM** ЩГЗ

Ρ Adjust Xfяпврвмпсу

fntentfy

ДОoPA |U:\T«^\orto_p^<N|*<*»c<\

- і

1 Cancel І Нф

Рис. 7.35. Вихідні дані для ортоперетворень: відомості про вихідні файли

7.5.3. Генерування ортозображень

з виправленням оптичної щільності

Оскільки в межах одного зображення може проявлятись нерівномірна оптична

щільність, то, використовуючи функцію Dodge (рис.7.36), можна її поправити.

У групі Use Intensities

in Range

можна встановити нижню

верхню межі оптичної

щільності пікселів (відсікання гістограми).

Використовуємо:

1.1.1. Minimum 15

Minimum 240

Група Intensity Influence Limits дає змогу видозмінити темні та світлі області

з одночасним зазначенням розміру області, в межах якої виконана ця дія.

Використовуємо:

Kernel Size 33

Darken Maximum 10

Brighten Maximum 10

У групі Pixel Intensity Adjustments можна встановити значення пікселів, які

будуть прозорими.

У Dodge Destination, подібно, як і для ортовиправлення, окреслюється формат

і місце, де будуть записані результуючі файли.

264

['осі· jt' l'rup«:i' Jt_-4

Parameters

Use Intensities In Range —

Мдемгс І·

її

urn: 240 jJ

Target

Average Intensity—j

Auto Calculate

ifiteniity f"1

Dodged Destination

ШЯШ

. ι, ι ι .

Kernel Size:

33 J

Ра^Махітиг^ГИИіО jj

Ві^Матшк |ШЩЇо З

-Pixel

Intensity Adjustments

[TransparencyІЇШЩ

[jjo jr]

Dodged Path:

JO:\mSVstykftSTYKJJ^^

fiefine (

OK І I

Hefc

Рис. 7.36. Вікно для корекції оптичної щільності одного знімка

7.5.4. Визначення лінії "склеювання" мозаїки

Найкраще проводити цю лінію так:

а) лінію проводимо в смузі перекриття знімків як у маршруті, так і між

маршрутами;

б)

для плоскої місцевості бажано, щоб лінія була ототожнена на сусідніх знімках

і проходила по узбіччю доріг, по стежках, посередині смуги зелені на автострадах

чи інших двосмугових дорогах, на забудованих територіях вздовж вулиць (обходячи

автомобілі);

в) якщо

рисуємо лінію через велике

поле, то

лінія повинна

бути

зигзагоподібною

і оминати елементи озеленення, дерева тощо;

г) у лісових масивах чи плантаціях проводимо цю лінію посередині просік;

д) при переході через водні об'єкти лінію необхідно проводити у найвужчому

місці, а там, де є кущі, лінію проводять перпендикулярно до русла ріки.

Після виконання цієї дії необхідно впевнитись, що всі зображення будуть

використані при "склеюванні".

Також необхідно перевірити, чи полігони, утворені цими лініями, повністю

покриваються масивом знімків.

265

7.5.5. Вирівнювання кольорів для створених ортофотокарт

У межах усього масиву зображень радіометрію ортозображень виправляють

так, щоб не було істотних відхилень за тональністю і кольором.

Для цього використовують інструмент IRASC, вибираючи з головного меню

аплікацію Contrast Contrast|Brightness та Contrast Gamma.

Якщо ті самі звернення (contrast, brightness, gamma) можна застосовувати до

більшості файлів,

то

доцільно використовувати аплікацію, яка

дасть

змогу автоматично

вирівняти тональність.

Готують файл-зразок із розширенням (наприклад, гес);

CONVERT.REC - назва файла

IRASC OPEN

Де:

SOURCE DIRECTORY - каталог з вихідними файлами;

DESTINATION DIRECTORY - каталог, в якому будуть записані змінені

файли;

CONTRAST BRIGHTNESS - СІ5-12 - це відповідно контраст і яскравість;

IRASC GAMMA - r - с 0.83 - червоний;

IRASC GAMMA - g - c 0.72 - зелений;

IRASC GAMMA - в - c 0.83 - голубий.

Якщо файл-зразок підготовано, то з рівня MicroStation задіємо аплікацію mgs.

ma, вказуємо файл-зразок (рис. 7.37) та вибираємо файли, для яких хочемо виконати

радіометричну корекцію (рис. 7.38).

Wybierz plik гес

Г»»;

alenj.rec

alenJl.rec_outrec

convert. rec_out гес

con vert 1.rec

convert"!. rec_out re

convert2.rec

convert2.rec_out.re

eonvert3.rec

convert3.rec_out.re

dorotarec

Jst

Files

of Jype:

*rec

Diectories.'

сЛ

frc:\

Q cj№a

documents and settings

drawing files

geoworkspaees

ου

program files

rec

recycler

CD system volume information

- rl

Cancel 1

ΕΦ t

Рис. 7.37. Вікно, де вказуємо файл-зразок

266

Wybierz pliki tif

ШшШт

8.184.08.1.»

ЇЇЖ08 3~Й

1 85

08Tw~~ "

8185

08 3 bf

8186

081 hf

8 1

86 08 3 ttf

8 1 87

081 Μ

8.187.08.3. tif

ai 88.08.1. tif

8.188.08.3. tif

Ш Fie* of Ідо:

£?ortho

fe no«i<zasieg_ob1

W id

Dgte [

Смйі I

i.«rt.··.. .1

•50-

Ftename:

o:\adioVrKW^asiegLob148.184.083.tif

o:\orthoW^2asiesLob1 \8.185.08.! tif

o: SorthoSnowyzasiegLobl \8.185.08.3. tf

o:\wtho\ncwiy2adegLob1 \8.186.08.1 Xi

Remove |

Рис 7 38 Вікно з переліком файлів для радіометричної корекції

7.5.6. Створення ортофотокарт з обрізанням по рамках секцій

з одночасним тональним вирівнюванням

Під

час

"склеювання" отримаємо поліпшений тонально

масив, який

покриває всю

область проекту. Тепер ортофотокарти будуть "нарізані" на аркуші (по трапеціях).

Оскільки метою є отримання ортофотокарти високої

якості,

то на кожному етапі

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

має однакову тональність (рис. 7.39).

"Склеювання" можна розпочати після підготовки окремих ортозображень та

ліній склеювання.

До них входить база з рамками аркушів. Тоді з головного меню OrthoPro

переходимо на Project Planning та закладку Product(s) (рис. 7.40).

Щоб ввести лінію "склеювання" до бази проекту, необхідно у головному вікні

Geomedia з головного меню Tool Display Design Files (рис. 7.41) вибрати положення

файла, систему координат та шар, на якому розміщена лінія склеювання. Далі

імпортуємо

цю

лінію

автоматично через Assign

Mosaic

Polygons прив'язуємо ортообраз

до визначеного полігону з використанням лінії склеювання (рис. 7.42).

267

3203531 3203532

шл

Рис. 7.39. Загальний вигляд ортофотокарти з покриттям аркушів

карти

Projiut Ρΐ.ιικιιηη: О: LPIS vtyki SIYKI 0b.?_k0S2ahfV projekt

Project Area

ISPM Souroefc)

Elevation Data

Qrihofc)

ізяйшвіі

Ptodudtype-

Ґ Oteete one pfoduct ov8( the ectiw pcojcct жм.

і Γ Deate NDOP

DOQ products

rakfe

tfie

proiec* area.

Ρ ЦИзіІЇ

СМшеїdefinedptoducUbsadetheproject«e&

Product Deanibon Be: jO MP!S4081\d^eVdanej»oduk^btok«\ .J

J3K J Caned 1 Ha» (

Рис. 7.40. Вибір бази з рамками продуктів

268

GeoMedia - GeoWorkspacel

Fie Edit Vim J Took Angsts Warehouse

[Projection +east;+northfm1 MG5M Data Server • ρ

fl]V|0h| Ш

Display Design Files.

Measure Distance

WebBrowser

ψ A

Ο Φ #

Customize,,,

^•Ml^lrJH ш\

Рис 7 41 Вибір файла, системи координат, шару

з лінією "склеювання"

т Ш т* tocrt Toob Arwfrsfc Warehouse Upend Qrtfwftt» Шфт Нф

аіай^і t &НУ щмшвіям ol т

|Pro»ecte)n +ββ»».Ήχκ^ηι) "wj| |Τ623710.94 5

ІІШіУІММіП^

b Q Mosaic Polygons - new2000s7

Mosaic Polygons Latoete - new2000s7

Ifc

·/ ^artfres

wvJCO «s

fe I JPr

A Af

iWiMI ?

^ [X]

еИЧеі Ρ oci.ji-t

' Label-'

new2Q0f^·"

t^ ["» ^vato'eP оФ OoOs

artl ρ Pre luct

Umf* r-Vk^L 00

t* Orthoe new2000®7

£Sj

Orthoe Labels new2000s7

&eva*ton РГнД^іГї^

Ь Κ)

EtevaiionFjotpf nt

Ufceb <* GQ£4?

Ь • ^elerie·} !SPM Ojrce rv а *

)

000^"

Ife

[X] Selected ISPM Sources Ubei

neft2030s"*

Ik LI> чл c r^ У20П0 "

b (S

^'''Stiafcfe KPM "uurces usbe!*

π·*ν»?υ<\ι*

Рис 7 42 Визначення місцезнаходження ортозображення у полігоні

269