Matematické základy fotogrammetrie, souřadnicové soustavy, transformace Pořízení snímků + změření hodnot • Fotogrammetrické zpracování • - transformace • - vyrovnání • - korelace ■ ■ ■ • Fotogrammetrické výstupy Snímkové orientace Fotogrammetrie řeší přepočet polohy bodu ze snímkových souřadnic do reálných souřadnic požadovaného systému. Před vlastní transformací souřadnic a následným vyhodnocováním LMS je nutné provést rekonstrukci polohy snímacího zařízení v době pořízení snímku - tzv. orientaci snímků. Negativ jS* Pro měřické úkoly a pro následnou tvorbu map je nezbytné přesně znát: > polohu středu promítání vzhledem k rovině snímku (vnitřní orientaci) > polohu středu promítání k vnějším souřadnicím a orientaci osy záběru v prostoru (vnější orientaci). Kontaktní pozitiv m. ] obni skutá letecká (pozemní) FM • počátek souřadnicové soustavy = střed snímku M', kde M'je průsečík spojnic rámových značek • souřadnicové osy • osa x'spojnice horizontálních RZ; + vpravo • osa y' (z) _L na osu x' v rovině snímku ;orientace v matematickém smyslu (česká konvence os = letecký x pozemní snímek) Snímkové souřadnice-schéma 0' osa záběru • Charakteristiky: 2D souřadnice v rovině snímku vysoká přesnost měření => řád um H' = M' ideál H ^ M' obecný případ => PVniO - dx', dy' (dz) Orientace os - pozemní snímek tři modelové souřadnice x, y, z - v měřítku, pravoúhlé, prostorové počátek - střed promítání O souřadnicové osy: - osa x ve směru letu - osa y kolmá v pravotočivém směru - osa z + orientace k zenitu (směr nadm. výšek tg negativy ✓> ks s ---^4-*----✓ • rovina xy - horizontální rovina (rovnoběžná se snímkovou rovinou, nebo stejná rovina jako geodetický systém) OK • rotace co kolem osy x (primárni) cp kolem osy y (sekundárni) K kolem osy z (terciální) kladný směr => pravotočivý Pozn.: stereodvojice - počátek leží ve středu promítání levého snímku negativ^ [A s ------ě JF ^ \ \y pozitiv x,y,z „ modelové souř. x'fx",z',z" .. snímkové souř. A B „stanoviska FTG základny pozemní FTG • tři geodetické souřadniceX, Y, Z nejčastěji prostorové, pravoúhlé • výsledná souřadnicová soustava objektu, nejčastěji S-JTSK (pro katastr, základní mapy apod.), WGS-84, ETRS89 aj. (GPS), u pozemní FM možné i vlastní místní souřadnicové soustavy Prvky vnitřní orientace Prvky vnitřní orientace přesně definují polohu středu promítání vzhledem k rovině snímku. Umožňují rekonstruovat svazek paprsků, který v okamžiku expozice vytvořil měřický snímek. Střed promítání je střed výstupní pupily O' > konstanta komory /- určuje se s přesností na setiny mm. > poloha hlavního snímkového bodu, který lze ztotožnit se středem snímku jako průsečíkem rámových značek (u správně seřízené komory M'=H', jinak je potřeba určit dx'a dy') > distorze objektivu - je udána výrobcem pro danou komoru či objektivu Prvky vnější orientace • Vnější orientace určuje polohu snímacího systému (středu promítání-O vstupní pupily) vzhledem k reálným souřadnicím. Určují i orientaci osy záběru v prostoru • Její prvky jsou většinou neznámé. • K prvkům vnější orientace patří šest následujících hodnot: Ground Point P X > tři souřadnice středu optického systému (O)v dané souřadné soustavě - Xo, Yo, Zo > tři úhly definující polohu osy záběru vůči souřadnicovým osám (směr osy záběru, sklon osy záběru a pootočení snímku - co,cp, k) rotace: 1. pncna 2. podélná 3. ve směru letu k - pootočení snímku ve vlastní rovině (rotace okolo osy z ... 3) cp - podélný sklon =^> sklon osy záběru od svislice ve směru letu (rotace okolo osy y... 2) co - příčný sklon =^> sklon osy záběru od svislice napříč směru letu (rotace okolo osy x ... 1) Určení prvků vnitřní orientace Kalibrační protokol - obsahuje informace o konstantě komory (f), přesné souřadnice rámových značek a radiální zkreslení. Přesné změření rozměrů fotografie (rámových značek), f je na okraji fotografie - potom lze zrekonstruovat polohu středu promítání. Údaje d rámových zrtačkáůri a raťUilrtfffl zkresleni kamery. K-EdiLiíbte dne Ů&. 12.2ÍH13 ;1 o o ZákJ.iiLJill LiJjjHď: Typ ham«y Čfetokafneť> Kioneásnila fcameíy LMK ti 152. tas Čito AŽ1 1 2 3 4 5 d t 8 K +11Ůj033 -109,964 +Ů.03Ů + 1tŮ,035 -1Q9J9ÍÔ +i tů.ůao ¥ +0,011 +0,024 +1tŮOf? -109J94S + fůa,93ft +11Ů022 -tos^sft 5ůUfáJntť ííi;Jl. yŕfinm rtawtí boťf atrii* ofirnaGe +ŮJ015 +0.O36 ♦qjůgů +cotů Piůseůk apcpi rettoi^eft Rí fCC +G.QŮB lí. jlIíjIm zí, i -a sil-a ni: £*nfin 0 20 40 30 10Ů 120 140 4 0 1 1 -2 -? -4 1 4 e; 1 1 2 2 5 t r 0 2 t Q ^3 -2 2 3 a 0 t 1 íí ■2 -t 3 2 přímé* Ů 1 t Ů -3 -t 3 3 Určení prvků vnější orientace Prvky vnější orientace se určují dodatečně (početně - analytické metody - a nebo empiricky), nejčastěji pomocí tzv. vlícovacích bodů V závislosti na konkrétních postupech orientace snímků (jeden snímek, překrývající se dvojice, blok snímků) je zapotřebí mít jistý minimální počet bodů, u nichž je známa poloha X, Y, resp. X, Y, Z. Vlícovací body mohu být značené či neznačené. Jejich přesné souřadnice se dříve zjišťovaly geodeticky, V současné době převládají metody GPS - určují souřadnice středu promítání. Tři úhly rotace se určují z měření IMS (inertial measuring unit). V některých případech je možné s dostatečnou přesností určit prvky vnější orientace přímo měřením GPS (v reálném čase) s přesností cca 15 cm. Určení prvků vnější orientace Vnější orientace modelu za pomoci vlícovacích bodů je tedy založena nejprve na procesu prostorového protínání zpět (a) - z vlícovacích bodů do modelu Po vypočtení prvků vnější orientace a obnovení modelu je potom možno určovat polohu všech ostatních bodů prostorovým protínáním vpřed (b). Relativní a absolutní orientace Především u zpracování snímků analytickými postupy na stereoplotrech se určení prvků vnější orientace provádí ve dvou krocích. • Relativní orientace - orientace stereoskopického páru v přístroji tak, aby vytvořil stereomodel v relativních souřadnicích - libovolně prostorově orientovaný, bez vazby na geodetické souřadnice. • Absolutní orientace - pootočení (rotace) a posun stereomodelu do geodetických souřadnic pomocí vlícovacích bodů. H' ■ŕ* O' z snímkové souřadnice modelové souřadnice 0P * OJ 0[Xq,YoZo] k Prvky vnitřní a vnější orientace pro pozemní a leteckou fotogrammetrii Základní transformační vztahy v FGM Transformace souřadnic - vzájemné zobrazení mezi 2 kartézskými souřadnicovými systémy a) rovinné transformace • shodnostní (posun, otočení) • podobnostní (posun, otočení, změna měřítka - 2 identické body), použití ve FM - např. převod zaměřených vlícovacích bodů, mezi místní a geodet, soustavou • afinní (posun, otočení, změna měřítka, tvaru - zkosení - 3 ident. body), např. převod měřených plošných souř. do soustavy snímkových souřadnic (pomocí rámových značek) • kolineární - středové zobrazení dvou rovinných souř. systémů - zachovává dvojpoměr, 4 ident. body, ve FM jednosnímková fotogrammetrie b) prostorové transformace Vztahy mezi souřadnými soustavami Fotogrammetrie řeší převod snímkových souřadnic objektu (x\ y\ z') na souřadnice geodetické (X, Y, Z). Tento převod zahrnuje obecně tři dílčí úkoly: 1) postupnou změnu orientace soustavy snímkových souřadnic (tzv. rotaci), 2) posunutí (tzv. translaci) počátku soustavy snímkových souřadnic 3) změnu měřítka Celou transformaci lze řešit postupně po krocích, které zahrnují převod snímku do ideální polohy (kolmý snímek), poté převod do soustavy modelových souřadnic a konečně převod souřadnic modelových na geodetické. Pootočení snímkového souřadného systému tak, aby tento byl rovnoběžný se souřadným systémem geodetickým. Modelové souřadnice x, y z jsou souřadnice rovnoběžné s reálným systémem X, Y, Z Získáme je rotací z původního systému snímkových souřadnic x', y', z'(f). Rotace je ve vztazích mezi těmito trojrozměrnými souřadnými soustavami vyjádřena tzv. rotační maticí M o rozměru 3x3: M = m12 m22 m23 vm3i m32 m33/ Rotace v trojrozměrném systému spočívá v trojím postupném pootočení vždy kolem jedné osy systému Nejprve se systém otočí o úhel - co kolem osy x, poté o úhel 9 kolem osy y a konečně o úhel k kolem osy z Souřadnice bodu A \/nové soustavě yp z1 pootočené o úhel co z původní soustavy y', z1: X yx - ycosú) + zsinco zx = -/sin# + zlcos*0 nove souřadnice t cos u> \ původní souřadnice y , z v maticovém zápisu: f 1 o 0 v 0 cos ú) sinú) 0 - sin co cos co yi a nebo zkráceně: y = M 2. Pootočení kolem osy y1 o úhel cp Souřadnice bodu A v dvakrát rotovaném systému x2y2 z2 budou ?2 \ \ x2 = xx cos^ + z1 sin^ z2 = —xľ sin^ + zx cos0 v maticovém zápisu: 2 x, ) f cos0 0 sin^ ]( x, y2 0 1 o v-sin0 0 cos^y a nebo zkráceně: x2 X^ = My yi \Z2j \ziJ yi Souřadnice bodu A v již třikrát rotovaném systému x, y, z : x = x2 cos K + y2 sin K y — -x2 sin K + y2 cos K z — z2 v maticovém zápisu: \ y COS K - sin K cos K 0 0 0 1 y2 a nebo zkráceně: "*2 x2 y = Mz y2 KZ2j Celý proces postupných rotací z původního systému souřadnic (x', y', z) do nového systému, který bude rovnoběžný se systémem geodetických souřadnic lze vyjádřit následovně: x2 y = MzMy yi - MZMyMx y KZJ \Z2j \zi) (X') y = M y a nebo zkráceně X = MX' Jednotlivé prvky matice představují tzv. směrové cosiny rotace a jsou určeny z následujících vztahů: m^ = COs((p) cos(k) m12= - cos(cp) sííi(k) m13 = sin((p) m21 = cos(co) síd(k) + sin(co) sin((p) cos(k) m22 = cos(co) cos(k) - sin(co) sin((p) síd(k) m23= -sin(co) cos(cp) m31 = sin(co) síd(k) - cos(co) sin((p) cos(k) m32= sin(co) cos(k) + cos(co) sin((p) síd(k) m33= cos(co) cos(cp) Rotační matice je maticí ortogonální, která má následující vlastnost: M~l =MT tedy inverzní matice se rovná matici transponované. Tato vlastnost je ve fotogrammetrii důležitá pro sestavení vztahu určujícího snímkové souřadnice bodu: x' = m^x + m12y + m13z X'=MTX anebo: y'= m21x + m22y + m23z z' = m31x + m32y + m33z Podmínka kolinearíty Bod na zemském povrchu, obraz tohoto bodu na snímku a střed promítání leží na jedné přímce. snímkové souřadnice libovolného bodu xp> yV snímkové souřadnice středu promítání * & y o> ť geodetické souřadnice středu promítání a je vektor ze středu promítání O do bodu p (na snímku) A je vektor z bodu O do bodu P (na zemském povrchu) Podmínka kolinearity: (J — k * A (k- měřítkové číslo) Ground Point P X Rovnice ko linearity Velikost vektoru a vyjádřená snímkovými souřadnicemi: Velikost vektoru A vyjádřená skutečnými souřadnicemi: X p X o a = y p-y o -f j A = Y -Y p o -Z ° J Vztah mezi snímkovými souřadnicemi libovolného bodu a skutečnými souřadnicemi tohoto bodu vyjádřený pomocí rotační matice M\ a = k-M A tj.maticově: í ' ' ^ xp x0 f*' -x0] y p - y'o = k-M ■ yp {-f J 1Zp -Zo) Rovnice kol i near i ty Výše uvedený vztah lze vyjádřit jako soustavu tří rovnic: x p -Xo =km[1(Xp-X0)+n\2(Yp-Y0)+n\£Zp -Z0)] yP -y0 =k[m21(Xp -X^+m^p -Y0)+m^{Zp-Z0)] -f =k[m31(Xp-X0)+m32(Yp -YQ)+r^3(Zp -Z0)] Podělíme-li první a druhou rovnici rovnicí třetí, obdržíme tzv. rovnice kolinearity. Tyto definují vztah mezi snímkovými a skutečnými souřadnicemi X p — X o f mn(Xp -X0) + mu(Yp-Y0) + mn(Zp -Z0) m3í(Xp-X0) + m32(Yp-Y0) + m33(Zp -Z0) ■ _ ■ rm2i(xp ~ Xo) + m22ÍYp -Y0) + m23(Zp -Z0) m31(Xp - X0) + m32(Yp -Y0) + m33(Zp -ZQ) Rovnice ko linea rity Analogicky lze kolineární vztah využít inverzně pro určení skutečných souřadnic bodu následovně: x = x + (z - z )mn^Xp -x°) + mi2(y p -y'o) + ^i3(-/) p ° p 0 m^x p - x 0) + m32(y p - y o) + m33(-f) Y =Y +(Z -Z )m2l(X^ -X^) + m22(};^ ->7^) + m23(-/) P ° P ° m3l(x P - x o) + m32(y P - y o) + m33(-f) Kolineárních rovnic je ve fotogrammetrii možné využít k určování prvků vnější orientace a dále především k transformaci souřadnic každého bodu na snímku do nové polohy vyjádřené v geodetickém souřadném systému - tedy např. k tvorbě ortofoto. Rovnice ko linea rity Je tam ale jeden problém: X = X + (Z - Z }mn(x r ~x o) + mn(y r ~ ľ o ) + mi3 ("/) 77 ° P ° m31 (x p - x o) + m32 (y p - y o) + m33 (-f) Y =Y +(Z -Z )m2l(X^ ~X°) + m22();^ ->7^) + m23(-/) P ° P ° m31(x p - x o) + m32(j p - y o) + m33(-f) Pro každý snímkový bod existuje nekonečné množství řešení (ze dvou měřených hodnot potřebujeme vypočítat 3 neznámé souřadnice) Z jednoho snímku nelze zrekonstruovat polohu 3D objektu z 2D souřadnic (snímkových) Co s tím? • potřebujeme informaci o souřadnici Z • nebo potřebujeme ještě jeden snímek Způsoby určení prvků vnější orientace V závislosti na počtu zpracovávaných snímků (jeden snímek či stereopár) a s tím souvisejícím potřebném počtu lícovacích bodů lze k fotogrammetrickým pracem využít následujících řešení, která využívají výše odvozených kolineárních rovnic: > zpětné promítání (space resection) - určení prvků vnější orientace samostatně pro jeden snímek > prostorové protínání vpřed (space forward intersection) - určení prvků vnější orientace společně pro dvojici překrývajících se snímků > blokové vyrovnání (bundle block adjustment) - určení prvků vnější orientace bloku snímků metodami aerotriangulace Ve fotogrammetrii existuje několik postupů k určení šesti neznámých prvků vnější orientace (X0, Y0, Z0, co,cp, k). Tyto postupy lze rozdělit na: 1. Početní - skládá se ze dvou kroků. Nejprve se provede relativní orientace, jejímž základem je změření tzv. vertikálních paralax na min. pěti bodech ve vyhodnocovacím přístroji. Poté následuje výpočet šesti neznámých prvků a tzv. absolutní orientace 2. Analytické - využívá se přímého vztahu mezi snímkovými a geodetickými souřadnicemi (základem je změření snímkových souřadnic). 3. Empirické - relativní orientace založená na postupném „ručním" odstraňování vertikálních paralax na orientačních bodech a následná absolutní orientace (posun, otočení a určení měřítka).