Technická podpora - zpracování rastrových obrázků. Přesuňte se do jihozápadního rohu

- Lineární algebra

Redukce kvadratického tvaru na hlavní osy

Dříve jsme zvažovali problém snížení reálného


q(x)= \sum_(n=1)^(n) \sum_(j=1)^(n) a_(ij)x_ix_j=x^TAx


n proměnných do


\widetilde(q)(y)=\lambda_1y_1^2+ \lambda_2y_2^2+ \ldots+ \lambda_ny_n^2


pomocí nedegenerované lineární změny proměnných x=Sy. K vyřešení tohoto problému jsme použili .


Zvažme jiný přístup k řešení. Lineární nedegenerovaná změna proměnných x=Sy s ortogonální maticí S~(S^(-1)=S^T) bude nazývána ortogonální změnou proměnných (neboli ortogonální transformací proměnných).


Pojďme formulovat problém redukce kvadratického tvaru na hlavní osy: je potřeba najít ortogonální změnu proměnných x=Sy (S^(-1)=S^T), čímž se kvadratická forma (9.23) dostane do kanonické formy (9.24).


K řešení použijeme následující geometrický význam úlohy. Pojďme počítat proměnné x_1, x_2,\ldots,x_n souřadnice vektoru \boldsymbol(x) n-rozměrného euklidovského prostoru \mathbb(E) na ortonormálním základě (\boldsymbol(e))= (\boldsymbol(e)_1,\ldots,\boldsymbol(e)_n), a matice A kvadratického tvaru (9.23) je maticí některých lineární transformace \mathcal(A)\dvojtečka \mathbb(E)\to \mathbb(E) na stejném základě. Navíc je tato transformace samoadjungovaná, protože její matice je symetrická: A^T=A. Kvadratická forma (9.23) může být reprezentována jako skalární součin


q(\boldsymbol(x))= \bigl\langle \mathcal(A)(\boldsymbol(x)), \boldsymbol(x)\bigr\rangle= \bigl\langle \boldsymbol(x), \mathcal(A )(\boldsymbol(x))\bigr\rangle.


Ortogonální změna proměnných x=Sy odpovídá přechodu z jedné ortonormální báze do druhé. Nechť S je skutečně přechodová matice z ortonormálního základu (\boldsymbol(e)) do ortonormálního základu (\boldsymbol(s))= (\boldsymbol(s)_1,\ldots,\boldsymbol(s)_n), tj. (\boldsymbol(s))= (\boldsymbol(e))S a S^(-1)=S^T. Potom souřadnice x vektoru \boldsymbol(x) v základu (\boldsymbol(e)) a souřadnice y stejného vektoru v základu (\boldsymbol(s)) souvisí vzorcem (8.11): x= Sy .


Problém redukování kvadratické formy na hlavní osy lze tedy formulovat následovně: je třeba najít v prostoru \mathbb(E) základ, ve kterém má samoadjungovaná transformační matice \mathcal(A) úhlopříčku formulář. Podle věty 9.10 je nutné zvolit ortonormální bázi z vlastních vektorů samoadjungované transformace. V tomto případě se přechodová matice S na kanonickou bázi ukáže jako ortogonální: S^T=S^(-1) .


Zformulujme tento výsledek pro kvadratickou formu.

Věta (9.12) o redukci kvadratické formy na hlavní osy

Reálnou kvadratickou formu (9.23) lze redukovat na kanonickou formu (9.24) pomocí ortogonální transformace proměnných x=Sy, kde - vlastní čísla matrice A.


Následek. Kvadratická forma (9.23) je kladně definitní (nezáporně definitní) právě tehdy, když jsou všechna vlastní čísla její matice kladná (nezáporná).


Poznámky 9.10


1. Při lineární nedegenerované změně proměnných se matice kvadratického tvaru mění podle vzorce (6.10): A"=S^TAS. Pro ortogonální matici S má tento vzorec tvar A"=S^(- 1)AS, který se shoduje se vzorcem (9.4) změny v lineární transformační matici při změně báze.


2. Najít kanonická forma(9.24) stačí určit všechny kořeny \lambda_1, \lambda_2,\ldots,\lambda_m(mezi kterými mohou být stejné) (rovnice) \det(A-\lambda E)=0, kde E je matice identity.


3. Důsledek věty 9.12 lze použít k analýze znaménka kvadratické formy:


– pokud jsou všechna vlastní čísla kladná (záporná), pak je kvadratická forma kladně (záporná) definitní;


– pokud jsou všechna vlastní čísla nezáporná (nekladná), pak je kvadratická forma nezáporná (nekladná) určitá;


– pokud existují vlastní hodnoty různých znaků, pak je kvadratická forma neurčitá (střídavá).


4. Výsledky formulované v odstavci 3 komentářů lze použít k ověření dostatečných a nutné podmínky druhého řádu v problému vyhledávání bezpodmínečný extrém funkcí. Chcete-li to provést, musíte najít vlastní čísla \lambda_1, \lambda_2,\ldots,\lambda_m \dfrac(d^2f(x))(dx^Tdx) v každém z pevných bodů x^(\ast) funkce f(x)=f(x_1,\ldots,x_n).


Pokud jsou všechna vlastní čísla kladná: \lambda_i>0,~ i=1,\ldots,n, pak v bodě x^(\ast) je lokální minimum;


– pokud jsou všechna vlastní čísla záporná: \lambda_i<0,~ i=1,\ldots,n , pak v bodě x^(\ast) je lokální maximum;


– pokud jsou všechna vlastní čísla nezáporná: \lambda_i\geqslant0,~ i=1,\ldots,n, pak v bodě x^(\ast) může být lokální minimum;


– pokud jsou všechna vlastní čísla nekladná: \lambda_i\leqslant0,~ i=1,\ldots,n, pak v bodě x^(\ast) může být lokální maximum;


– pokud vlastní hodnoty \lambda_i,~ i=1,\ldots,n, různá znaménka, pak v bodě x^(\ast) není žádný extrém;


– pokud jsou všechna vlastní čísla nula: \lambda_i=0,~ i=1,\ldots,n, pak je nutný další výzkum.


5. Problém redukce kvadratické formy na hlavní osy je řešen pomocí algoritmu pro redukci samoadjungované transformace na diagonální formu. V tomto případě je nalezen diagonální tvar matice kvadratické formy a ortogonální matice S změny proměnných x=Sy, čímž se kvadratická forma dostává do kanonického tvaru (k hlavním osám).

Příklad 9.7. Určete znaménko kvadratického tvaru tří proměnných


q(x)= x_1^2+2x_1x_2+2x_1x_3+x_2^2+2x_2x_3+x_3^2


a najít ortogonální změnu proměnných x=Sy, čímž se kvadratická forma dostane do kanonické formy (k hlavním osám).

Řešení. Sestavíme matici kvadratického tvaru: A=\začátek(pmatice) 1&1&1\\ 1&1&1\\ 1&1&1 \end(pmatrix). V příkladu 9.6 byly nalezeny vlastní hodnoty této matice: \lambda_(1,2)=0, \lambda_3=3. Všechna vlastní čísla jsou nezáporná, takže kvadratická forma je nezáporná a definitní (viz bod 4 v poznámkách 9.10).


Byla nalezena ortogonální matice


S=\begin(pmatrix) \dfrac(\sqrt(2))(2)& \dfrac(\sqrt(6))(6)& \dfrac(\sqrt(3))(3)\\ 0&-\ dfrac(\sqrt(6))(3)& \dfrac(\sqrt(3))(3)\\ -\dfrac(\sqrt(2))(2)& \dfrac(\sqrt(6))( 6)& \dfrac(\sqrt(3))(3) \end(pmatrix)\!,


redukující matici A na diagonální tvar \Lambda= \operatorname(diag) (0,0,3). Zapíšeme požadovanou ortogonální změnu proměnných x=Sy:


x_1= \frac(\sqrt(2))(2)\,y_1+ \frac(\sqrt(6))(6)\,y_2+ \frac(\sqrt(3))(3)\,y_3;\quad x_2= -\frac(\sqrt(6))(3)\,y_2+ \frac(\sqrt(3))(3)\,y_3;\quad x_3= -\frac(\sqrt(2))(2 )\,y_1+ \frac(\sqrt(6))(6)\,y_2+ \frac(\sqrt(3))(3)\,y_3.


a kvadratická forma v kanonické formě: \widetilde(q)(y)= 3y_3^2.

Příklad 9.8. Najděte lokální extrémy funkce dvou proměnných pomocí matic


f(x)=3x_1^5+x_1^4-5x_1^3-2x_1^2x_2+x_2^2.


Řešení. V kroku 1 byl nalezen gradient funkce a z nutné podmínky pro extrém prvního řádu tři stacionární body:


x^0= \begin(pmatrix)0&0 \end(pmatrix)^T,\qquad x^1=\begin(pmatrix) 1&1 \end(pmatrix)^T,\qquad x^2=\begin(pmatrix) - 1&1\konec(pmatice)^T.


Hessenská matice má tvar

\frac(df(x))(dx^Tdx)= \begin(pmatrix) 60x_1^3+12x_1^2-30x_1-4x_2&-4x_1\\-4x_1&2 \end(pmatrix)\!.


Pojďme najít vlastní čísla Hessovy matice v každém stacionárním bodě:


\frac(df(x^0))(dx^Tdx)= \begin(pmatrix)0&0\\ 0&2 \end(pmatrix)\!;\quad \frac(df(x^1))(dx^Tdx ) = \begin(pmatrix)38&-4\\ -4&2 \end(pmatrix)\!,\quad \frac(df(x^2))(dx^Tdx)= \begin(pmatrix) -22&4\\4&2\ konec (pmatrix)


a použijte odstavec 4 komentáře 9.10.

Na místě x^0=\begin(pmatrix)0\\0 \end(pmatrix) hessenská matice má tvar \begin(pmatrix) 0&0\\ 0&2\end(pmatrix). Z rov. \begin(vmatrix) -\lambda&0\\ 0&2-\lambda\end(vmatrix)=0 najdeme \lambda_1=0, \lambda_2=2 . Protože všechna vlastní čísla jsou nezáporná, může být v bodě x^0 lokální minimum a pro konečný závěr je zapotřebí další výzkum (viz příklad 6.13).


Na místě x^1=\začátek(pmatice)1\\1 \konec(pmatice) hessenská matice má tvar \begin(pmatrix) 38&-4\\ -4&2 \end(pmatrix). Z rov. \begin(vmatrix) 38-\lambda&-4\\ -4&2-\lambda\end(vmatrix)=0 nebo \lambda^2-40 \lambda+60=0 dostaneme \lambda_(1,2)= 20\pm2\sqrt(85). Protože všechna vlastní čísla jsou kladná, pak v bodě x^1 existuje lokální minimum funkce.


Na místě x^2=\začátek(pmatice)-1\\1 \end(pmatrix) hessenská matice má tvar \begin(pmatrix) -22&4\\ 4&2 \end(pmatrix). Z rov. \begin(vmatrix) -22-\lambda&4\\ 4&2-\lambda\end(vmatrix)=0 nebo \lambda^2+40 \lambda-60=0 dostaneme \lambda_(1,2)=-10\pm4\sqrt(10). Protože vlastní čísla mají různá znaménka, v bodě x^2 není žádný extrém.

Často máme papírovou mapu oblasti a chceme tuto mapu přidat do našeho projektu GIS. Podívejme se, jak vytvořit georeferenční snímek z naskenované nebo vyfotografované mapy na příkladu mapy rezervace Kvitucha Gora.

Ve výše uvedeném příkladu je vše provedeno v QGIS. Při práci budou použity následující moduly: Rastrová vazba, QuickMapServices, GeoSearch. Tyto pluginy je třeba nainstalovat a aktivovat, můžete si přečíst více o instalaci modulů. Moduly QuickMapServices a GeoSearch vyžadují ke své činnosti připojení k internetu.


Dalším krokem je nalezení oblasti zájmu na základní mapě. Po pečlivém prozkoumání naskenované mapy na ní najdeme název osady - „vesnice Milcha“.


Když známe název obce, můžeme ji najít pomocí jednoho z modulů „GeoSearch“, „osmSearch“ nebo „OSM place search“.


Po zmenšení mapy na místo zájmu přistoupíme přímo k propojení mapy. Pro georeferencování rastrových obrázků má QGIS vestavěný modul „Raster Referencing“ (Georeferencer). Modul se spouští ze sekce menu „Rastr“ - „Vazba rastru“.


Modul Georeferencer se otevře v novém okně.


Pomocí tlačítka „Otevřít rastr“ nebo kombinace kláves +přidat obrázek, na který budeme odkazovat.
V horní části okna se zobrazí obrázek a ve spodní části tabulka popisující kotevní body.


Dále je třeba vybrat body na základní mapě a snímek, ke kterému bude snímek georeferencován. Obvykle se jedná o křižovatky a odbočky silnic, mosty a další objekty, které jsou jasně viditelné na základní mapě a připojeném obrázku.

Zvětšujeme rozsah základní mapy k prvnímu kotevnímu bodu. Také zvětšíme obrázek, který je ukotven k vybranému kotevnímu bodu. Po najetí na kotevní bod v okně modulu klikněte na tlačítko „Přidat bod“ a klikněte myší na vybraný bod. Poté se otevře formulář pro zadání souřadnic. Souřadnice lze zadat buď pomocí vstupních polí, nebo je zachytit z mapy. Pokud máme souřadnice bodů získané např. pomocí GPS navigátoru, můžeme je zadat do příslušných polí. Chcete-li získat souřadnice ze základní mapy, klikněte na tlačítko „Z mapy“.


Po kliknutí na tlačítko „Z mapy“ se automaticky otevře hlavní okno QGIS. V něm kurzor myši vypadá jako bílý kříž. Vyberte kotevní bod na základní mapě a stiskněte levé tlačítko myši.


Po kliknutí se automaticky vrátíme do okna modulu rastrové vazby. Hodnoty souřadnic bodu se objeví ve vstupním formuláři. Vyplněné hodnoty mají souřadnicový systém projektu se základní mapou.


Po kliknutí se bod přidá do tabulky s popisem kotevních bodů. Tímto způsobem přidáme co nejvíce kotevních bodů. Je vhodné umístit body rovnoměrně přes propojený obrázek. Čím více je zdrojový obraz zkreslený, tím více kotevních bodů je zapotřebí. Minimální počet kotevních bodů jsou 3.


Dále nastavte parametry transformace. Chcete-li to provést, klepněte na ozubené kolo na panelu nástrojů. V okně, které se otevře, nastavte následující požadované hodnoty: typ transformace, metoda interpolace, cílový souřadnicový systém, cílový rastr. Zbývající parametry jsou volitelné a lze je ponechat s výchozími hodnotami.

Kvalita uchopení závisí na počtu bodů uchopení a volbě metody transformace. Můžete si přečíst více o transformačních metodách.


Jedním z klíčových bodů je správné určení cílového souřadnicového systému. Pokud jste zadali souřadnice získané pomocí GPS navigátoru, označte souřadnicový systém zadaný v nastavení GPS navigátoru, nejčastěji je to WGS 84 (EPSG:4326). Pokud jsme vzali souřadnice z mapy, pak označujeme souřadnicový systém pracovního projektu. V našem případě se jedná o WGS 84 / Pseudo Mercator (EPSG:3857), který je „nativní“ pro takové mapové služby jako OpenStreetMap, ArcGIS Online a mnoho dalších.


Po nastavení parametrů transformace zahájíme proces vazby kliknutím na zelený trojúhelník na panelu nástrojů nebo výběrem příslušné položky v nabídce „Soubor“. V důsledku rastrové vazby bude získán soubor ve formátu GeoTIFF.

Pokud jste v okně transformačních parametrů zaškrtli volbu „Otevřít výsledek QGIS“, pak po dokončení procesu vazby bude výsledná vrstva přidána do pracovního projektu nad základní mapu.

Důležitou nuancí je, že v důsledku činnosti modulu má výsledný rastr v parametrech transformace zadaný souřadnicový systém, ale neobsahuje informaci o tom, o jakou projekci rastru se přesně jedná. Z tohoto důvodu může být přítomen v seznamu vrstev, ale nezobrazen na mapě. V tomto případě musíte přejít do „vlastností hladiny“ a zadat požadovaný souřadnicový systém ručně.


Jakmile explicitně určíte správný souřadnicový systém, bude obrázek umístěn na správné místo.


Úpravou průhlednosti můžeme skrýt černá pole podél okrajů propojeného obrázku, která vznikla transformací.


Správnost vazby můžeme zkontrolovat i zadáním průhlednosti vrstvy na 50 %.


Vrcholem může být záznam informací o projekci přímo do georeferencovaného obrazového souboru ve formátu GeoTIFF. Chcete-li to provést, vyberte nabídku „Rastr“ - „Projekce“ - „Přiřadit projekci“. V okně, které se otevře, vyberte náš georeferenční soubor a souřadnicový systém. Tento nástroj můžeme použít i pro katalogy s rastry zaškrtnutím příslušného políčka a v případě potřeby ručně nastavit parametry programu gdalwarp. Můžete si přečíst více o parametrech gdalwarp.


Po zadání informací o souřadnicovém systému do souboru georeferencovaného obrázku jej budeme moci otevřít v jiných GIS projektech a GIS programech bez dalších manipulací ručním zadáním souřadnicového systému.

Můžeme také oříznout propojený obrázek na oblast zájmu, jak je popsáno, a stáhnout jej do tabletu nebo smartphonu, jak je popsáno.

Rastrová data se obvykle získávají skenováním papírových map nebo zpracováním leteckých snímků či satelitních snímků. Datové sady získané skenováním map obvykle nejsou georeferencovány (ať už v rámci sady, ani jako samostatný soubor). Někdy jsou informace o poloze poskytované v souborech leteckých nebo satelitních snímků nepřesné a taková data nelze správně kombinovat s jinými dostupnými prostorovými daty. Chcete-li tedy použít některé rastrové datové sady s ostatními prostorovými daty, musíte tyto rastrové datové sady georeferencovat na souřadnicový systém mapy. Souřadnicový systém mapy je specifikován pomocí mapové projekce (metoda promítání zakřiveného povrchu Země do roviny).

Georeferencování rastrových dat se provádí určením umístění pomocí mapových souřadnic a přiřazením souřadnicového systému k datovému rámci. Prostorově odkazovaná rastrová data vám umožňují prohlížet, dotazovat se a analyzovat data spolu s dalšími geografickými daty. Panel nástrojů Prostorová reference Umožňuje georeferencovat rastrové datové sady a vrstvy rastrových datových sad (které mohou mít rastrovou funkcionalitu), obrazové služby a rastrové produkty.

V zásadě se k nastavení prostorové reference rastrové datové sady používají následující kroky:

Chcete-li vidět ukázku georeferencování rastrové datové sady, podívejte se na video o georeferencování.

Uchopte rastr pomocí řídicích bodů

Georeferencování rastrových dat v zásadě probíhá pomocí existujících prostorových dat (cílových dat), jako je georeferenční rastr nebo vektorová třída prvků, která má požadovaný mapový souřadnicový systém. Tento proces zahrnuje definování řady řídicích bodů – známých souřadnic x,y – které spojují známá místa v rastrové datové sadě s odpovídajícími místy v georeferencovaných datech (cílových datech). Kontrolní body jsou místa, která lze snadno definovat v rastrových datových sadách a která mají přesné souřadnice. Identifikátory polohy mohou zahrnovat například křižovatky silnic nebo vodních toků, ústí řek, skalní výchozy, rohy ulic nebo náměstí a křižovatky lesních ochranných pásů.

Kontrolní body se používají ke konstrukci polynomiální transformace, která posune rastrový soubor dat do geograficky správného umístění. Spojení mezi referenčním bodem v rastrové datové sadě (point from) a odpovídajícím bodem v již zarovnaných datech (point to) se nazývá odkaz.

Níže uvedený příklad ukazuje referenční bod „do“ (žlutý kříž) umístěný na vektorové datové vrstvě na křižovatce ulic a související referenční bod „od“ (zelený kříž) umístěný na rastru. Spojení je znázorněno modrou čarou spojující oba body.

Počet požadovaných odkazů závisí na složitosti transformace, kterou plánujete použít k propojení rastrové datové sady s mapou. Přidání dalších odkazů nemusí nutně zlepšit přesnost transformace. Pokud je to možné, je nejlepší rozmístit vztahy rovnoměrně v celé rastrové datové sadě, spíše než je soustředit na jedno místo. Nejlepší výsledky obvykle poskytuje jeden odkaz v každém rohu rastru a několik ve středu.

Obecně platí, že čím větší je oblast překrytí mezi rastrovým datovým souborem a cílovými daty, tím lepší je výsledek georeferencování, protože kontrolní body, ze kterých je konstruován, mohou být rovnoměrněji rozloženy. Pokud například vaše cílová data pokrývají čtvrtinu plochy rastrové datové sady, budou kontrolní body soustředěny pouze v této oblasti. Cokoli, co nespadá do oblasti překrytí, bude tedy s největší pravděpodobností zaklapnuto nesprávně.

Mějte na paměti, že přesnost georeferencovaných dat nebude větší než přesnost dat, ze kterých georeferencujete. Pro snížení chyb by mělo být georeferencování prováděno s použitím dat s nejvyšším rozlišením v největším měřítku.

Rastrová konverze

Jakmile vytvoříte dostatek vztahů, můžete přichytit – nebo transformovat – sadu rastrových dat tak, aby odpovídala souřadnicím cílových dat. Chcete-li dosáhnout maximální shody souřadnic pro každou buňku, můžete použít polynomiální, spline, fiting nebo podobnostní transformace.

Polynomiální transformace využívá konstrukci polynomu založeného na referenčních bodech a algoritmu proložení metodou nejmenších čtverců (LSF). Tato metoda je optimalizována pro celkovou přesnost, ale nezaručuje místní přesnost. Polynomiální transformace používá dva vzorce: jeden pro výpočet souřadnic x ze vstupních (x,y) souřadnic umístění a druhý pro výpočet souřadnic y. Cílem nejmenších čtverců je vypočítat obecný vzorec, který platí pro všechny body, obvykle mírným posunutím pozic referenčních bodů. Počet nekorelovaných kontrolních bodů požadovaných pro tuto metodu je 1 pro nultý řád, 3 pro první řád, 6 pro druhý řád a 10 pro třetí řád. Polynomy nižšího řádu mají tendenci mít chybu náhodného typu, zatímco polynomy vyššího řádu mají tendenci mít chybu extrapolace.

Polynomiální transformace prvního řádu se běžně používá pro georeferencování obrazu. Níže je uvedena rovnice pro transformaci rastrové datové sady pomocí afinní (prvního řádu) polynomiální transformace. Šest parametrů určuje transformaci rastrových řádků a sloupců na mapové souřadnice.


K zkreslení vašich dat bude použit polynom nultého řádu. To se často používá v situaci, kdy jsou vaše data již georeferencována, ale malý posun vaše data lépe zarovná. K provedení polynomického posunu nulového řádu na vašich datech je potřeba pouze jeden odkaz. Doporučuje se provést několik připojení a poté vybrat to, které vypadá nejpřesněji.

Použijte transformaci prvního řádu (afinní) k posunutí, změně měřítka nebo otočení rastrové datové sady. To obvykle umožňuje zachovat rovné linie původního rastru v transformovaném rastrovém souboru dat. Čtverce a obdélníky rastrové datové sady jsou tedy typicky zkrouceny do rovnoběžníků spolu s náhodnými změnami v měřítku a úhlové orientaci.

Pomocí pouhých tří odkazů může matematický výraz použitý v transformaci prvního řádu přesně propojit každý rastrový bod s cílovými daty. Použití více než tří odkazů přináší deformace nebo rezidua, která jsou distribuována napříč všemi odkazy. Měli byste však zadat více než tři odkazy, protože pokud je jeden z nich pozičně nesprávný, bude to mít mnohem větší dopad na proces převodu. Proto i přes nárůst matematické chyby transformace s rostoucím počtem spojení se zvýší i celková přesnost transformace.

Čím vyšší je řád transformace, tím složitější zkreslení lze korigovat. Transformace vyšší než třetího řádu se však používají velmi zřídka. Transformace vyššího řádu vyžadují více připojení, a tedy delší dobu zpracování. Obecně platí, že pokud je třeba rastrovou datovou sadu roztáhnout, otočit nebo změnit její velikost, měli byste použít transformaci prvního řádu. Pro složitější transformace se používají transformace druhého nebo třetího řádu.


Spline transformace je skutečná metoda pryžového listu optimalizovaná tak, aby poskytovala místní, ale ne globální přesnost přichycení. Je založen na splajnové funkci, po částech polynomu, který poskytuje plynulé, spojité přechody mezi sousedními polynomy. Zajišťuje přesné zarovnání zdrojových a cílových referenčních bodů, ale u odstraněných pixelů není zaručeno, že budou přesně zarovnány. Tato transformace je užitečná, pokud je správné umístění kotevních bodů nanejvýš důležité. Přidání dalších řídicích bodů může zvýšit celkovou přesnost transformace spline. Vyžaduje alespoň 10 kontrolních bodů.

Transformace fitováním je vhodná pro zachování jak globální přesnosti (nejmenších čtverců, LSF), tak lokální přesnosti. Je postaven na algoritmu, který optimalizuje polynomiální transformaci a metody interpolace triangulační sítě (TIN). Polynomiální transformace se provádí pomocí dvou sad řídicích bodů, které je umístí co nejpřesněji, s přihlédnutím k technice interpolace TIN. Pro přizpůsobení hranic jsou zapotřebí alespoň tři referenční body.

Projektivní transformace může zachovat linie tak, aby zůstaly rovné. V tomto případě se paralelní čáry mohou ukázat jako neparalelní. Projektivní transformace je užitečná zejména pro zkreslené snímky, naskenované mapy a satelitní snímky Landsat nebo Digital Globe. K provedení projektivní transformace jsou zapotřebí alespoň čtyři připojení. Pokud existují čtyři propojení, bude chyba RMS nula. Pokud je připojení více, bude chyba RMS mírně nad nulou.

Podobnostní transformace je transformace prvního řádu, která se snaží zachovat tvar původního rastru. Střední kvadratická chyba má tendenci být vyšší než u jiných polynomických transformací, protože zachování tvaru má vyšší prioritu než přesné přizpůsobení.

Interpretace střední kvadratické chyby

Když se odvodí obecný vzorec a aplikuje se na kontrolní body, změří se chyba – reziduum. Tato chyba je rozdíl mezi skutečnou a vypočítanou polohou "bodu k". Celková chyba se vypočítá z druhé odmocniny (RMS) součtu chyb všech odkazů. Tato hodnota popisuje stupeň koordinace transformace mezi všemi referenčními body (vazbami). Pokud je chyba příliš velká, můžete ji odstranit nebo přidat kontrolní body, abyste ji snížili.

Ačkoli je chyba RMS dobrým odhadem přesnosti transformace, malé množství chyb RMS by nemělo být zaměňováno s přesností mapování. Výsledek transformace může mít například značné nepřesnosti v důsledku nesprávného výběru kontrolních bodů. Čím správnější kontrolní body použijete, tím přesnější bude výsledek převodu původních dat. Proložení hranic a splajnové transformace obvykle produkují chybu RMS přibližně nula nebo nula; to však neznamená, že je obraz přesně odkazován v prostoru.

Dopředný posun zobrazuje chybu ve stejných jednotkách, jaké jsou určeny prostorovým odkazem datového rámce. Reverzní offsety zobrazují chybu v pixelech. Zkreslení vpřed a vzad je mírou přesnosti a zobrazuje se v pixelech. Hodnoty offsetu bližší nule jsou považovány za přesnější.

Přepočet rastrové datové sady

Když transformujete, promítáte nebo převzorkováváte rastrovou datovou sadu, převádíte ji z jedné projekce na jinou nebo měníte velikost buňky, provádíte geometrickou transformaci. Geometrická transformace je proces změny geometrie rastrové datové sady z jednoho souřadnicového prostoru do druhého. Geometrické transformace zahrnují pryžovou fólii (běžně používanou při prostorových odkazech), projekci (používá informace o projekci k transformaci dat do jiné projekce), translaci (přesun všech souřadnic stejně), rotaci (otočení všech souřadnic o stejný úhel) a změnu velikosti rastru. buňky datové sady.

Po provedení geometrické transformace se středy buněk původního rastru zřídka shodují se středy buněk konečného rastru; hodnoty však musí být vycentrovány.

Můžete si myslet, že každá buňka rastrové datové sady je transformována do nového umístění na mapě, ale ve skutečnosti tomu tak není. Během georeferencování se pomocí souřadnic mapy vypočítá matice prázdných buněk. Každé prázdné buňce je pak přiřazena hodnota na základě výsledku přepočtu.

Tři nejčastěji používané metody převzorkování jsou: nejbližší soused, bilineární interpolace a kubická konvoluce. Všechny tyto metody přiřazují hodnotu každé prázdné buňce na základě výsledků zpracování buněk neověřené rastrové datové sady.

Výpočet nejbližšího sousedství je nejrychlejší metodou převzorkování a dobře se hodí pro zpracování kategorických nebo aktuálních dat, protože nemění hodnoty vstupních buněk. Poté, co jsou středy buněk výstupních dat porovnány se středy buněk vstupního rastru, metoda nejbližšího souseda určí středovou polohu nejbližší buňky vstupního rastru a přiřadí tuto hodnotu buňce výstupního rastru.

Při použití metody nejbližšího souseda se hodnoty buněk vstupního rastru nezmění. Hodnota 2 ve vstupním rastru zůstane vždy hodnotou 2 ve výstupním rastru; nezmění se na 2, 2 nebo 3. Protože výstupní hodnoty zůstávají stejné, měla by být použita metoda Nearest Neighbor s nominálními nebo ordinálními daty, kde každá hodnota představuje třídu, funkční jednotku nebo klasifikaci – to může být kategoriální údaje, jako jsou typy půdy, lesy nebo využití půdy.

Bilineární interpolace využívá hodnoty čtyř nejbližších středů vstupních buněk k výpočtu jedné hodnoty ve výstupním rastru. Hodnota výstupní buňky je vážený průměr těchto čtyř vstupních buněk, upravený podle jejich vzdálenosti od středu odpovídající buňky ve vstupním rastru. Ve srovnání s metodou nejbližšího souseda tato metoda převzorkování vytváří hladší obraz.

Protože hodnoty výstupních buněk se počítají v závislosti na jejich relativní poloze a hodnotách vstupních buněk, bilineární interpolace je vhodná pro data, ve kterých poloha vůči známému bodu nebo jevu určuje hodnotu buňky - tzn. pro souvislé plochy. Nadmořské výšky, sklony, intenzita letištního hluku nebo salinita podzemní vody v okolí ústí řeky jsou jevy, které jsou reprezentovány souvislými plochami a nejlépe je změnit měřítko pomocí metody bilineární interpolace.

Kubická konvoluce je podobná bilineární interpolaci s tím rozdílem, že vážený průměr se vypočítává pomocí hodnot ve středech 16 nejbližších buněk. Kubická konvoluce vytváří ostřejší obraz než bilineární interpolace, protože do procesu výpočtu výstupní hodnoty je zapojeno více buněk. Proto se tato metoda převzorkování často používá ke zpracování snímků, jako jsou satelitní a letecké snímky.

Bilineární interpolace a kubická konvoluce by se neměly používat s kategorickými daty, protože kategorie nebudou ve výstupu zachovány. Všechny tři metody však lze použít se spojitými daty, přičemž nejbližší okolí produkuje větší pixely, bilineární interpolace vytváří hladší povrch a kubická konvoluce vytváří nejostřejší obraz. Při zahájení velkého projektu převzorkování se doporučuje vytvořit prototyp a vyzkoušet několik metod převzorkování, abyste určili metodu, která nejlépe vyhovuje vašim datům.

Je nutné transformovat rastr?

Transformaci rastrové datové sady po georeferencování můžete uložit pomocí příkazu Transformace na panelu nástrojů Prostorová reference nebo pomocí nástroje Transform. Kromě toho můžete uložit informace o převodu do dalšího souboru pomocí příkazu panelu nástrojů Prostorová reference.

Transformace umožňují vytvořit novou rastrovou datovou sadu, která je georeferencována na mapové souřadnice. Můžete jej uložit ve formátu BIL, BIP, BMP, BSQ, DAT, GIF, GRID, IMG, JPEG, JPEG 2000, PNG nebo TIFF. ArcGIS nevyžaduje konverzi rastrové datové sady pro její zobrazení s jinými prostorovými daty; měli byste to však udělat, pokud jej zamýšlíte analyzovat nebo jej chcete použít v jiných programech, které nerozpoznají prostorovou referenci uloženou v samostatném georeferenčním souboru.

Při aktualizaci prostorové reference se transformace uloží do samostatného souboru – nevytvoří se nová rastrová datová sada, jako je tomu při ukládání transformace rastrové datové sady. U rastrové datové sady ve formě souboru, jako je TIFF, bude transformace obvykle zapsána do externího souboru XML s příponou .AUX.XML. Pokud je rastrová datová sada nezpracovaným obrázkem, jako je BMP, a transformace je afinní, bude zapsána do světového souboru. Pokud je rastrová datová sada uložena v geodatabázi, při spuštění příkazu Aktualizace prostorového odkazu, transformace bude také uložena v databázi jako další informace týkající se této rastrové datové sady. Aktualizace vrstvy rastru, obrazové služby nebo mozaiky aktualizuje vrstvu v mapovém dokumentu, ale neuloží změněnou prostorovou referenci do zdrojového souboru.

Níže uvedená tabulka ukazuje, jak jsou jednotlivé typy cílů uloženy.

Prostorová reference různých rastrů

Typ datVýsledek

Rastrová datová sada

Tým Aktualizace prostorového odkazu aktualizuje rastrovou datovou sadu.

Rastrová vrstva

Tým Aktualizace prostorového odkazu Aktualizuje rastrovou datovou vrstvu, zdrojové rastry nejsou ovlivněny.

Obslužná vrstva obrazu

Služba obrázků na serveru není aktualizována. Po provedení příkazu Aktualizace prostorového odkazu

Rastrový produkt

Rastrový produkt neaktualizuje podkladové soubory rastrových datových sad. Po provedení příkazu Aktualizace prostorového odkazu Můžete buď uložit mapový dokument (.mxd) nebo vytvořit soubor vrstvy (.lyr), abyste si uložili práci s georeferencováním.

Rastr s funkcí

Funkce rastru neaktualizuje podkladové rastrové soubory. Po provedení příkazu Aktualizace prostorového odkazu Můžete buď uložit mapový dokument (.mxd) nebo vytvořit soubor vrstvy (.lyr), abyste si uložili práci s georeferencováním.




Nahoru