153YZOD Zpracování obrazových dat 2006 - 8. cvičení: Porovnání verzí

Z GeoWikiCZ
(doplneni dalsich sekci)
Řádek 8: Řádek 8:


== metoda hlavních a obecných komponent ==
== metoda hlavních a obecných komponent ==
Multispekrální obrazové kanály velmi často obsahují korelace a to díky podobné spekrální odezvě sledovaných objektů nebo mírně se překrývající funkci filtru spekrálního skeneru. Metoda hlavních komponent (''Principal Component Transformation'') provede transformaci těchto dat do nového spekrálního prostoru bez korelací mezi jednotlivými kanály. To ve svém důsledku vede ke koncentaci obrazové informace do menšího počtu kanálů.
Jak ukazuje obr. č.1 a 2, metoda je založena na transformaci původního multispekrálního datasetu do nového spekrálního souřadnicového systému.
[[Soubor:ZOD-cv8-pca1.png|frame|center|Obr. č.1: Redukce dat příznakového prostoru dvou kanálů satelitního snímku, amplituda ''A'' se redukuje na ''A''' - převzato z GRASS příručky]]
[[Soubor:ZOD-cv8-pca2.png|frame|center|Obr. č.2: Zobrazení bodových dat jako standardizovaných datových vektorů ''x<sub>1</sub>'' až ''x<sub>5</sub>'' s&nbsp;prvními dvěma přidruženými ortogonálními komponentami - vektory ''PCA<sub>1</sub>'' a&nbsp;''PCA<sub>2</sub>'' v&nbsp;polárním zobrazení (vlevo) a&nbsp;bodová data v&nbsp;kartézském souřadnicovém systému (vpravo) - převzato z&nbsp;GRASS příručky]]
První kanál po&nbsp;transformaci (označme jej <font style="background-color: rgb(200,200,200)">PC1</font>) obsahuje maximální stupeň možné variance původních obrazových vrstev. Druhý kanál (<font style="background-color: rgb(200,200,200)">PC2</font>) pokrývá maximální varianci, která není obsažena v&nbsp;prvním kanálu. Další PC kanály obsahují zbytkové variance. Přičemž počet PC kanálů odpovídá počtu vstupních pásem. Poslední PC kanály nicméně obsahují pouze nekorelovaný šum. Tuto metodu tak lze použít pro zmenšení nároku na&nbsp;úložný prostor, obrazová informace je koncentrována v&nbsp;menším počtu kanálů.
Metoda hlavních komponent je v&nbsp;GRASSu dostupná pomocí modulu <tt>i.pca</tt>. Modul vyžaduje seznam vstupních obrazových vrstev (nejméně dvě), prefix pro výstupní kanály, volitelně je možné provést změnu rozsahu DH (parametr <tt>rescale</tt>).
GRASS&nbsp;>&nbsp;g.region rast=tm1
#transformace hlavních komponent v GRASSu
#
GRASS&nbsp;>&nbsp;i.pca in=tm1,tm2,tm3,tm4,tm5,tm7 out=pc
Před dalším zpracováním nastavíme tabulku barev (hromadně pro všechny PC kanály):
#hromadné nastavení tabulky barev pro kanály PC
GRASS&nbsp;>&nbsp;for pc in `g.mlist pat=pc*`;do \
r.colors map=$pc color=grey.eq;\
done
[[Soubor:ZOD-cv8-pc1-6.png|frame|center|Obr. č.3: Přehled PC kanálů]]


== transformace RGB - IHS ==
== transformace RGB - IHS ==

Verze z 12. 12. 2005, 20:07

Metoda hlavních komponent, transformace RGB - IHS

< Zpracování obrazových dat

osnona

seznam příkazů

metoda hlavních a obecných komponent

Multispekrální obrazové kanály velmi často obsahují korelace a to díky podobné spekrální odezvě sledovaných objektů nebo mírně se překrývající funkci filtru spekrálního skeneru. Metoda hlavních komponent (Principal Component Transformation) provede transformaci těchto dat do nového spekrálního prostoru bez korelací mezi jednotlivými kanály. To ve svém důsledku vede ke koncentaci obrazové informace do menšího počtu kanálů.

Jak ukazuje obr. č.1 a 2, metoda je založena na transformaci původního multispekrálního datasetu do nového spekrálního souřadnicového systému.

Obr. č.1: Redukce dat příznakového prostoru dvou kanálů satelitního snímku, amplituda A se redukuje na A' - převzato z GRASS příručky
Obr. č.2: Zobrazení bodových dat jako standardizovaných datových vektorů x1x5 s prvními dvěma přidruženými ortogonálními komponentami - vektory PCA1PCA2 v polárním zobrazení (vlevo) a bodová data v kartézském souřadnicovém systému (vpravo) - převzato z GRASS příručky

První kanál po transformaci (označme jej PC1) obsahuje maximální stupeň možné variance původních obrazových vrstev. Druhý kanál (PC2) pokrývá maximální varianci, která není obsažena v prvním kanálu. Další PC kanály obsahují zbytkové variance. Přičemž počet PC kanálů odpovídá počtu vstupních pásem. Poslední PC kanály nicméně obsahují pouze nekorelovaný šum. Tuto metodu tak lze použít pro zmenšení nároku na úložný prostor, obrazová informace je koncentrována v menším počtu kanálů.

Metoda hlavních komponent je v GRASSu dostupná pomocí modulu i.pca. Modul vyžaduje seznam vstupních obrazových vrstev (nejméně dvě), prefix pro výstupní kanály, volitelně je možné provést změnu rozsahu DH (parametr rescale).

GRASS > g.region rast=tm1
#transformace hlavních komponent v GRASSu
#
GRASS > i.pca in=tm1,tm2,tm3,tm4,tm5,tm7 out=pc

Před dalším zpracováním nastavíme tabulku barev (hromadně pro všechny PC kanály):

#hromadné nastavení tabulky barev pro kanály PC
GRASS > for pc in `g.mlist pat=pc*`;do \
r.colors map=$pc color=grey.eq;\
done
Obr. č.3: Přehled PC kanálů

transformace RGB - IHS