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

Z GeoWikiCZ
Přejít na: navigace, hledání
m (Metody zvýrazňování obrazu, histogram, barevné modely, skládání pásem: zmena nadpisu)
m (upozorneni)
 
(Není zobrazeno 15 mezilehlých verzí od 2 dalších uživatelů.)
Řádek 1: Řádek 1:
 +
<div class="boilerplate metadata" id="attention" style="border: 3px double #c32c2c; margin-bottom: 2em; padding: 0 .25em; background-color: #f96666;">
 +
<p style="font-size: 150%; text-align: center;"><br/>'''Tento text se vztahuje k akademickému roku 2006/2007. Aktuální stránky cvičení jsou dostupné [[153YZOD Zpracování obrazových dat - cvičení 3|zde]].'''<br/>&nbsp;</p>
 +
<div style="clear: both;"><br style="clear: both; display: none;" /></div>
 +
</div>
 +
{{ZOD}}
 +
[ [[Zpracování obrazových dat]] ]
 +
 
= Metody zvýrazňování obrazu, histogram, barevné modely, syntézy =
 
= Metody zvýrazňování obrazu, histogram, barevné modely, syntézy =
 
+
__TOC__
 
== osnova ==
 
== osnova ==
  
V&nbsp;první části tohoto kurzu budeme pracovat s&nbsp;tzv.&nbsp;histogramem
+
V&nbsp;první části tohoto kurzu budeme pracovat s&nbsp;tzv.&nbsp;''histogramem''
 
obrazových dat, poté si představíme
 
obrazových dat, poté si představíme
jednoduchou metodu zvýrazňování obrazu - hustotní řezy. Teoreticky budou popsány
+
jednoduchou metodu zvýrazňování obrazu - ''hustotní řezy''. Teoreticky budou popsány
i&nbsp;dva modely barev - RGB a&nbsp;IHS, v&nbsp;praktické rovině budeme pracovat
+
i&nbsp;dva ''modely barev'' - RGB a&nbsp;IHS, v&nbsp;praktické rovině budeme pracovat
s&nbsp;tzv.&nbsp;barevnými syntézami.
+
s&nbsp;tzv.&nbsp;''barevnými syntézami''.
  
 
== seznam použitých příkazů ==
 
== seznam použitých příkazů ==
  
* d.mon
+
* [http://grass.itc.it/grass60/manuals/html60_user/d.mon.html d.mon]
* d.histogram
+
* [http://grass.itc.it/grass60/manuals/html60_user/d.histogram.html d.histogram]
* r.stats
+
* [http://grass.itc.it/grass60/manuals/html60_user/r.stats.html r.stats]
* d.erase
+
* [http://grass.itc.it/grass60/manuals/html60_user/d.erase.html d.erase]
* r.info
+
* [http://grass.itc.it/grass60/manuals/html60_user/r.info.html r.info]
* r.colors
+
* [http://grass.itc.it/grass60/manuals/html60_user/r.colors.html r.colors]
* d.rast
+
* [http://grass.itc.it/grass60/manuals/html60_user/d.rast.html d.rast]
* d.redraw
+
* [http://grass.itc.it/grass60/manuals/html60_user/d.redraw.html d.redraw]
* g.region
+
* [http://grass.itc.it/grass60/manuals/html60_user/g.region.html g.region]
* d.rgb
+
* [http://grass.itc.it/grass60/manuals/html60_user/d.rgb.html d.rgb]
  
== odrazivost v jednotlivých kanálech LandSat-TM5 ==
+
== spekrální odrazivost v jednotlivých kanálech LandSat-TM5 ==
  
 
V&nbsp;průběhu minulého kurzu jsme si osvojili proces manipulace
 
V&nbsp;průběhu minulého kurzu jsme si osvojili proces manipulace
Řádek 30: Řádek 37:
  
 
Rastrové buňky
 
Rastrové buňky
jednotlivých kanálů družicové scény obsahují tzv.&nbsp;digitální hodnoty (DH). Tyto hodnoty
+
jednotlivých kanálů družicové scény obsahují tzv.&nbsp;''digitální hodnoty'' (DH). Tyto hodnoty
 
(již digitálně kódované) reprezentují v&nbsp;rámci daného spektrálního intervalu
 
(již digitálně kódované) reprezentují v&nbsp;rámci daného spektrálního intervalu
 
velikost skenerem zachyceného odraženého záření od zemského povrchu. Družicové snímky obsažené
 
velikost skenerem zachyceného odraženého záření od zemského povrchu. Družicové snímky obsažené
Řádek 46: Řádek 53:
 
pro vodu, zelenou vegetaci a&nbsp;písčitou půdu. V&nbsp;pravé části obrázku
 
pro vodu, zelenou vegetaci a&nbsp;písčitou půdu. V&nbsp;pravé části obrázku
 
je vidět dvoudimenzionální systém vytvořený
 
je vidět dvoudimenzionální systém vytvořený
třetím a&nbsp;čtvrtým kanálem družicové scény LandSat-TM5 (tzv.&nbsp;feature space)
+
třetím a&nbsp;čtvrtým kanálem družicové scény LandSat-TM5 (tzv.&nbsp;''feature space'')
 
s&nbsp;viditelnými shluky pixelů - náznaky jednotlivých tříd.
 
s&nbsp;viditelnými shluky pixelů - náznaky jednotlivých tříd.
  
Řádek 70: Řádek 77:
 
výsledek:''
 
výsledek:''
  
[[Soubor:ZOD-cv3-histogram-tm1_cast.png|frame|center|Obr č.3: Histogram prvního kanálu družicové scény LandSat-TM5 - vybraná část snímku]]"
+
[[Soubor:ZOD-cv3-histogram-tm1_cast.png|frame|center|Obr č.3: Histogram prvního kanálu družicové scény LandSat-TM5 - vybraná část snímku]]
  
 
V&nbsp;minulém kurzu jsme tutéž charakteristiku obrazu (tj.&nbsp;četnost DH) vyjádřili číselně, pouze
 
V&nbsp;minulém kurzu jsme tutéž charakteristiku obrazu (tj.&nbsp;četnost DH) vyjádřili číselně, pouze
 
pro osvěžení: vytiskneme celkový počet pixelů pro jednotlivé DH (pomocí systémové aplikace
 
pro osvěžení: vytiskneme celkový počet pixelů pro jednotlivé DH (pomocí systémové aplikace
<tt>awk</tt> omezíme výstup pouze na&nbsp;interval DH [55; 85]). Jako <u>ostřílení</u>
+
[http://www.gnu.org/software/gawk/gawk.html awk] omezíme výstup pouze na&nbsp;interval DH [55; 85]). Jako <u>ostřílení</u>
 
uživatelé dozajista víme jak na věc, pouze pro připomenutí:
 
uživatelé dozajista víme jak na věc, pouze pro připomenutí:
  
Řádek 112: Řádek 119:
 
  83 62548
 
  83 62548
 
  84 31241
 
  84 31241
 
----
 
  
 
Na&nbsp;následujícím obrázku jsou zobrazeny histogramy všech kanálů zkoumané družicové scény.
 
Na&nbsp;následujícím obrázku jsou zobrazeny histogramy všech kanálů zkoumané družicové scény.
Řádek 123: Řádek 128:
 
''Technická poznámka: Jak je vidno z&nbsp;výše uvedeného příkazu, na jeden řádek lze uvést hned několik příkazů oddělených středníkem.''
 
''Technická poznámka: Jak je vidno z&nbsp;výše uvedeného příkazu, na jeden řádek lze uvést hned několik příkazů oddělených středníkem.''
  
[[Soubor:ZOD-cv3-histogramy_tm_1.png|frame|center|Obr č.4: Histogramy jednotlivých kanálů družicové scény LandSat-TM5]]
+
[[Soubor:ZOD-cv3-histogramy_tm.png|frame|center|Obr č.4: Histogramy jednotlivých kanálů družicové scény LandSat-TM5]]
 
 
[[Soubor:ZOD-cv3-histogramy_tm_2.png|frame|center|Obr č.4: Histogramy jednotlivých kanálů družicové scény LandSat-TM5]]
 
  
 
== metoda hustotních řezů ==
 
== metoda hustotních řezů ==
Řádek 141: Řádek 144:
 
  #
 
  #
 
  GRASS:~&nbsp;>&nbsp;r.info tm1 | grep Range
 
  GRASS:~&nbsp;>&nbsp;r.info tm1 | grep Range
 
+
 
  |  Range of data:    min =  0 max = 255                                    |
 
  |  Range of data:    min =  0 max = 255                                    |
  
Řádek 162: Řádek 165:
 
Tuto tabulku můžete zadat ručně,
 
Tuto tabulku můžete zadat ručně,
 
či ji uložit do samostatného textové souboru (línější jedinci daný soubor pouze stáhnou na disk -
 
či ji uložit do samostatného textové souboru (línější jedinci daný soubor pouze stáhnou na disk -
[http://gama.fsv.cvut.cz/~martin/vyuka/zzoz/cv3/tb1.txt tb1.txt])
+
[http://gama.fsv.cvut.cz/~landa/geowikicz_data/zod/cv3/tb1.txt tb1.txt])
 
a&nbsp;nastavit tímto jednoduchým způsobem:
 
a&nbsp;nastavit tímto jednoduchým způsobem:
  
 
  #nastavení tabulky barev uložené v textovém souboru
 
  #nastavení tabulky barev uložené v textovém souboru
 +
#
 
  GRASS:~&nbsp;>&nbsp;r.colors map=tm1 color=rules &lt; tb1.txt
 
  GRASS:~&nbsp;>&nbsp;r.colors map=tm1 color=rules &lt; tb1.txt
  
Řádek 179: Řádek 183:
 
  GRASS:~&nbsp;>&nbsp;d.rast map=tm1
 
  GRASS:~&nbsp;>&nbsp;d.rast map=tm1
  
''Poznámka: Z&nbsp;"výchovných" důvodů jsem tentokrát uvedl úplný název parametru
+
''Poznámka: Z&nbsp;"výchovných" důvodů jsem tentokrát uvedl úplný název parametru - <tt>map</tt>. Většinou si z&nbsp;čiré lenosti příkazy zkracuji...''
- <tt>map</tt>. Většinou si z&nbsp;čiré lenosti příkazy zkracuji...''
 
  
 
[[Soubor:ZOD-cv3-tm1-tb1.png|center|frame|Obr č.6: První kanál družicové scény LandSat-TM5 s nastavenou tabulkou barev (tb1.txt)]]
 
[[Soubor:ZOD-cv3-tm1-tb1.png|center|frame|Obr č.6: První kanál družicové scény LandSat-TM5 s nastavenou tabulkou barev (tb1.txt)]]
 
   
 
   
 
Prostudujeme si blíže histogram snímku - četnosti jednotlivých DH a&nbsp;podle toho pozměníme
 
Prostudujeme si blíže histogram snímku - četnosti jednotlivých DH a&nbsp;podle toho pozměníme
naší tabulku barev, například takto ([http://gama.fsv.cvut.cz/~martin/vyuka/zzoz/cv3/tb2.txt tb2.txt]):
+
naší tabulku barev, například takto ([http://gama.fsv.cvut.cz/~landa/geowikicz_data/zod/cv3/tb2.txt tb2.txt]):
  
 
  0 0 0 0
 
  0 0 0 0
Řádek 205: Řádek 208:
 
A&nbsp;teď trochu teorie - popis barevných modelů '''RGB''' a&nbsp;'''IHS'''.
 
A&nbsp;teď trochu teorie - popis barevných modelů '''RGB''' a&nbsp;'''IHS'''.
 
Podobně jako geometrická
 
Podobně jako geometrická
data má i&nbsp;prostor barev svůj souřadnicový systém. Model RGB má aditivní charakter, kdy ostatní
+
data má i&nbsp;prostor barev svůj souřadnicový systém. Model RGB má ''aditivní'' charakter, kdy ostatní
 
barvy vznikají kombinací tří základních barev (červená, zelená a&nbsp;modrá). Druhý model IHS (HIS)
 
barvy vznikají kombinací tří základních barev (červená, zelená a&nbsp;modrá). Druhý model IHS (HIS)
má odlišnou podstatu, intenzita (intensity) odráží jas barvy, odstín (hue) koresponduje
+
má odlišnou podstatu, intenzita (''intensity'') odráží jas barvy, odstín (''hue'') koresponduje
s&nbsp;dominantní vlnovou délkou, třetí složkou je potom sytost (saturation) barvy.
+
s&nbsp;dominantní vlnovou délkou, třetí složkou je potom sytost (''saturation'') barvy.
 
Na&nbsp;obr.&nbsp;č.8 jsou zobrazeny
 
Na&nbsp;obr.&nbsp;č.8 jsou zobrazeny
 
oba výše zmíněné barevné modely.
 
oba výše zmíněné barevné modely.
Řádek 224: Řádek 227:
  
 
Nyní si ukážeme praktický postup při zobrazování
 
Nyní si ukážeme praktický postup při zobrazování
tzv.&nbsp;barevných syntéz - kompozitních snímků.
+
tzv.&nbsp;barevných syntéz - ''kompozitních snímků''.
 
Takový snímek vznikne přiřazením třech vybraných kanálů scény třem základním barvám modelu RGB
 
Takový snímek vznikne přiřazením třech vybraných kanálů scény třem základním barvám modelu RGB
 
- červené, zelené a&nbsp;modré (pracujeme-li s modelem RGB).
 
- červené, zelené a&nbsp;modré (pracujeme-li s modelem RGB).
  
Zobrazíme tedy barevnou kompozici v&nbsp;tzv.&nbsp;skutečných barvách
+
Zobrazíme tedy barevnou kompozici v&nbsp;tzv.&nbsp;''skutečných barvách''
 
(první kanál přiřadíme modré ose,
 
(první kanál přiřadíme modré ose,
 
druhý zelené a&nbsp;konečně třetí kanál červené ose RGB krychle). Takto složený snímek zobrazíme dvakrát,
 
druhý zelené a&nbsp;konečně třetí kanál červené ose RGB krychle). Takto složený snímek zobrazíme dvakrát,
Řádek 249: Řádek 252:
 
  #jak jednoduché a průhledné;-)
 
  #jak jednoduché a průhledné;-)
 
  #
 
  #
  GRASS:~&nbsp;>&nbsp;for map in `g.mlist rast pat=tm*`;do echo r.colors map=$map color=grey;done | head -n3 | sh
+
  GRASS:~&nbsp;>&nbsp;for map in `g.mlist rast pat=tm*`;do echo r.colors map=$map color=grey;done | \
 +
head -n3 | sh
  
 
A&nbsp;konečně zobrazíme danou barevnou syntézu v&nbsp;GRASS monitoru:
 
A&nbsp;konečně zobrazíme danou barevnou syntézu v&nbsp;GRASS monitoru:
Řádek 255: Řádek 259:
 
  GRASS:~&nbsp;>&nbsp;d.rgb blue=tm1 green=tm2 red=tm3
 
  GRASS:~&nbsp;>&nbsp;d.rgb blue=tm1 green=tm2 red=tm3
  
[[Soubo:ZOD-cv3-rgb_grey.png|frame|center|Obr č.9: Barevná syntéza ve skutečných barvách - tabulka barev grey]]
+
[[Soubor:ZOD-cv3-rgb_grey.png|frame|center|Obr č.9: Barevná syntéza ve skutečných barvách - tabulka barev grey]]
  
 
Výsledek se nám příliš nelíbí - změníme tedy tabulku barev a&nbsp;kompozitní snímek
 
Výsledek se nám příliš nelíbí - změníme tedy tabulku barev a&nbsp;kompozitní snímek
Řádek 272: Řádek 276:
 
[[Soubor:ZOD-cv3-rgb_grey.eq.png|center|frame|Obr č.10: Barevná syntéza ve skutečných barvách - tabulka barev grey.eq]]
 
[[Soubor:ZOD-cv3-rgb_grey.eq.png|center|frame|Obr č.10: Barevná syntéza ve skutečných barvách - tabulka barev grey.eq]]
  
Na závěr si vyzkoušíme další kombinace - o nepravých (falešných) barvách
+
Na závěr si vyzkoušíme další kombinace - o ''nepravých (falešných) barvách''
 
hovoříme v případě, kdy složkami výsledné kompozice nejsou kanály z viditelného oboru světla.
 
hovoříme v případě, kdy složkami výsledné kompozice nejsou kanály z viditelného oboru světla.
  
''Poznámka: Následující tabulka je převzata z dodatku GRASS příručky.''
+
''Poznámka: Následující tabulka je převzata z dodatku [http://grass.fsv.cvut.cz/wiki/index.php/GIS_GRASS_-_Praktick%C3%A1_rukov%C4%9B%C5%A5 GRASS příručky].''
  
<center>
+
{|align="center"
{|
 
|-
 
!kombinace !! !! popis
 
 
|-
 
|-
 
|123
 
|123
Řádek 323: Řádek 324:
  
 
[[Soubor:ZOD-cv3-rgb_354.png|frame|center|Obr č.12: Barevná syntéza v nepravých barvách - 354]]
 
[[Soubor:ZOD-cv3-rgb_354.png|frame|center|Obr č.12: Barevná syntéza v nepravých barvách - 354]]
 +
 +
 +
[ [[Zpracování obrazových dat]] ]

Aktuální verze z 19. 11. 2008, 19:54

[ Zpracování obrazových dat ]

Metody zvýrazňování obrazu, histogram, barevné modely, syntézy

osnova

V první části tohoto kurzu budeme pracovat s tzv. histogramem obrazových dat, poté si představíme jednoduchou metodu zvýrazňování obrazu - hustotní řezy. Teoreticky budou popsány i dva modely barev - RGB a IHS, v praktické rovině budeme pracovat s tzv. barevnými syntézami.

seznam použitých příkazů

spekrální odrazivost v jednotlivých kanálech LandSat-TM5

V průběhu minulého kurzu jsme si osvojili proces manipulace s tabulkami barev (v odborné literatuře se často používá termín look-up tables - LUT). Pro družicové snímky používáme většinou vyrovnané (ekvalizované) odstíny šedi (grey.eq).

Rastrové buňky jednotlivých kanálů družicové scény obsahují tzv. digitální hodnoty (DH). Tyto hodnoty (již digitálně kódované) reprezentují v rámci daného spektrálního intervalu velikost skenerem zachyceného odraženého záření od zemského povrchu. Družicové snímky obsažené v lokaci sevcech (LandSat-TM5) jsou kódovány jako 8-bitová data (tzn. že mohou obsahovat až 256 různých odstínů - hodnoty v rozsahu [0; 255]). Níže uvedený obrázek ukazuje spektrální rozlišení jednotlivých kanálů družicové scény LandSat-TM5.

Obr č.1: Spekrální odrazivost v jednotlivých kanálech družicové scény LandSat-TM5 (převzato z GRASS příručky)

První tři kanály tak pokrývají viditelné spektrum, zbylé kanály potom infračervený obor. Šestý kanál je v tepelném spektru - je označován jako panchromatický a od ostatních kanálů se navíc liší i geometrickým rozlišením - 120 m. Na obrázku č.1 je navíc zobrazena křivka odrazivosti pro vodu, zelenou vegetaci a písčitou půdu. V pravé části obrázku je vidět dvoudimenzionální systém vytvořený třetím a čtvrtým kanálem družicové scény LandSat-TM5 (tzv. feature space) s viditelnými shluky pixelů - náznaky jednotlivých tříd.

histogram

Grafické znázornění (v podobě sloupcového grafu) četností jednotlivých DH poskytuje tzv. histogram. Na ose X jsou naneseny jednotlivé DH ([0; 255]), na ose Y pak četnost jejich výskytu. Zobrazíme si tedy zkušebně histogram prvního kanálu družicové scény LandSat-TM5. Otevřeme GRASS monitor a aplikujeme příkaz d.histogram (pokud již máme v GRASS monitoru zobrazena nějaká data, je nanejvýš vhodné, před zobrazením samotného histogramu jeho obsah vymazat - d.erase):

GRASS:~ > d.mon x0
#
#zobrazení histogramu
#
GRASS:~ > d.histogram tm1
Obr č.2: Histogram prvního kanálu družicové scény LandSat-TM5

Poznámka: Nikoho jistě nepřekvapí, že i tento modul respektuje nastavení aktivního regionu. Pokud nemáme nastaven aktivní region na celý družicový snímek, můžeme dostat například takovýto výsledek:

Obr č.3: Histogram prvního kanálu družicové scény LandSat-TM5 - vybraná část snímku

V minulém kurzu jsme tutéž charakteristiku obrazu (tj. četnost DH) vyjádřili číselně, pouze pro osvěžení: vytiskneme celkový počet pixelů pro jednotlivé DH (pomocí systémové aplikace awk omezíme výstup pouze na interval DH [55; 85]). Jako ostřílení uživatelé dozajista víme jak na věc, pouze pro připomenutí:

#vytisknout četnost DH pouze pro interval hodnot [55; 85]
#
GRASS:~ > r.stats tm1 -c | awk '{if ($1>55 && $1<85) print $1,$2}'


r.stats:  100%
56 6
57 84
58 1582
59 9348
60 17800
61 63277
62 113072
63 110607
64 174193
65 220580
66 281778
67 189538
68 156562
69 205687
70 201391
71 140880
72 104717
73 165492
74 196217
75 203551
76 134829
77 129954
78 143600
79 85170
80 59787
81 56945
82 67504
83 62548
84 31241

Na následujícím obrázku jsou zobrazeny histogramy všech kanálů zkoumané družicové scény. V rámci procvičení a zkoumání rozdělení četnosti DH v jednotlivých kanálech si zobrazte postupně histogramy vybraných snímků, např.:

GRASS:~ > d.erase; d.histogram tm2

Technická poznámka: Jak je vidno z výše uvedeného příkazu, na jeden řádek lze uvést hned několik příkazů oddělených středníkem.

Obr č.4: Histogramy jednotlivých kanálů družicové scény LandSat-TM5

metoda hustotních řezů

Jednou ze základních (velmi jednoduchých) metod zvýraznění obrazu je metoda hustotních řezů (ang. density slicing). V podstatě jde pouze o nastavení vhodné tabulky barev - jednotlivým DH jsou přiřazeny dané barevné hodnoty. Interval DH kanálu je rozdělen na tři podintervaly a mezní hodnoty takto vzniklých intervalů přiřazeny třem základním barvám - modré, zelené a červené.

Tuto metodu si demonstrujeme například na prvním kanálu družicové scény, nejprve zjistíme rozsah DH (např. pomocí modulu r.info):

#zjištění intervalu hodnot (pomocí aplikace grep vybereme pouze řádek obsahující slovo "Range")
#
GRASS:~ > r.info tm1 | grep Range

|   Range of data:    min =  0 max = 255                                     |

Tento interval tedy rozdělíme na tři podintervaly s následujícími mezními hodnotami: 85 a 170. Na základě tohoto rozdělení nastavíme příslušnou tabulku barev:

0 255 0 0
85 255 255 0
170 0 255 0
255 0 0 255

Poznámka: Místo RGB trojice můžeme použít názvy standardních barev GRASSu:

0 red
85 yellow
170 green
255 blue

Tuto tabulku můžete zadat ručně, či ji uložit do samostatného textové souboru (línější jedinci daný soubor pouze stáhnou na disk - tb1.txt) a nastavit tímto jednoduchým způsobem:

#nastavení tabulky barev uložené v textovém souboru
#
GRASS:~ > r.colors map=tm1 color=rules < tb1.txt

Zobrazíme-li znovu histogram tohoto snímku, dostaneme stejný sloupcový graf - pouze v odlišném barevném znázornění korespondujícím s nově nastavenou tabulkou barev - viz obr. č.5.

Obr č.5: Histogram prvního kanálu družicové scény LandSat-TM5 (vlastní tabulka barev)

Zobrazíme-li konečně daný snímek v GRASS monitoru, zjistíme, že jeho barevná interpretace není zrovna příliš "šťastná".

GRASS:~ > d.rast map=tm1

Poznámka: Z "výchovných" důvodů jsem tentokrát uvedl úplný název parametru - map. Většinou si z čiré lenosti příkazy zkracuji...

Obr č.6: První kanál družicové scény LandSat-TM5 s nastavenou tabulkou barev (tb1.txt)

Prostudujeme si blíže histogram snímku - četnosti jednotlivých DH a podle toho pozměníme naší tabulku barev, například takto (tb2.txt):

0 0 0 0
50 255 0 0
80 255 255 0
120 0 255 0
255 0 0 255

Tabulku nastavíme již známým způsobem a snímek (aktuální obsah GRASS monitoru) překreslíme:

GRASS:~ > r.colors map=tm1 color=rules < tb2.txt
GRASS:~ > d.redraw
Obr č.7: První kanál družicové scény LandSat-TM5 s nastavenou tabulkou barev (tb2.txt)

barevné modely RGB a IHS

A teď trochu teorie - popis barevných modelů RGBIHS. Podobně jako geometrická data má i prostor barev svůj souřadnicový systém. Model RGB má aditivní charakter, kdy ostatní barvy vznikají kombinací tří základních barev (červená, zelená a modrá). Druhý model IHS (HIS) má odlišnou podstatu, intenzita (intensity) odráží jas barvy, odstín (hue) koresponduje s dominantní vlnovou délkou, třetí složkou je potom sytost (saturation) barvy. Na obr. č.8 jsou zobrazeny oba výše zmíněné barevné modely.

Obr č.8: Barevné modely RGB a IHS

Pixel v prostoru barev RGB je tedy dán specifickou polohou uvnitř barevné krychle, podobně v případě modelu IHS uvnitř hexagonu.

GRASS umožňuje snadnou konverzi mezi těmito dvěma modely barev - jde o moduly i.his.rgb. To pouze na okraj, prakticky si je představíme v jednom z následujících kurzů.

barevné syntézy

Nyní si ukážeme praktický postup při zobrazování tzv. barevných syntéz - kompozitních snímků. Takový snímek vznikne přiřazením třech vybraných kanálů scény třem základním barvám modelu RGB - červené, zelené a modré (pracujeme-li s modelem RGB).

Zobrazíme tedy barevnou kompozici v tzv. skutečných barvách (první kanál přiřadíme modré ose, druhý zelené a konečně třetí kanál červené ose RGB krychle). Takto složený snímek zobrazíme dvakrát, v prvním případě nastavíme podkladovým snímkům tabulku barev grey, v druhém případě vhodnější grey.eq. Rozdíl je na první pohled zřetelný - viz obr. č.9 a obr. č.10.

#čert nikdy nespí...
#
GRASS:~ > g.region tm
#
#tabulka barev grey
#
GRASS:~ > r.colors tm1 col=grey
GRASS:~ > r.colors tm2 col=grey
GRASS:~ > r.colors tm3 col=grey

Jelikož se považujeme za pokročilé uživatele OS GNU/Linux, jsme neskutečně líní a chceme "ohromit" naše blízké okolí, nastavíme tabulku barev pro první tři kanály hromadně(!):

#jak jednoduché a průhledné;-)
#
GRASS:~ > for map in `g.mlist rast pat=tm*`;do echo r.colors map=$map color=grey;done | \
head -n3 | sh

A konečně zobrazíme danou barevnou syntézu v GRASS monitoru:

GRASS:~ > d.rgb blue=tm1 green=tm2 red=tm3
Obr č.9: Barevná syntéza ve skutečných barvách - tabulka barev grey

Výsledek se nám příliš nelíbí - změníme tedy tabulku barev a kompozitní snímek překreslíme:

#tabulka barev grey.eq
#
GRASS:~ > r.colors tm1 col=grey.eq
GRASS:~ > r.colors tm2 col=grey.eq
GRASS:~ > r.colors tm3 col=grey.eq
#
#stačí překreslit obsah okna
#
GRASS:~ > d.redraw
Obr č.10: Barevná syntéza ve skutečných barvách - tabulka barev grey.eq

Na závěr si vyzkoušíme další kombinace - o nepravých (falešných) barvách hovoříme v případě, kdy složkami výsledné kompozice nejsou kanály z viditelného oboru světla.

Poznámka: Následující tabulka je převzata z dodatku GRASS příručky.

123 ... barevný snímek ve skutečných barvách - jsou použity první tři kanály pokrývající interval viditelného světla - viz obr. č.10.
234 ... Citlivý na zelenou vegetaci (znázorněno červenou barvou), jehličnaté lesy jsou v tmavější odstínu než listnaté. Silnice a vodní plochy jsou dobře rozeznatelné - viz obr. č.11.
243 ... Zelená vegetace je zelená, ale jehličnany nejsou tak dobře patrné jako na syntéze 234.
247 ... Asi nejlepší kombinace pro určení typu lesa. Dobře použitelné pro mapování těžebních ploch a silnic.
345 ... Obsahuje hlavní obrazové kanály (VIS, NIR, SWIR). Zelená vegetace je znázorněna zeleně, SWIR ukazuje stres na vegetaci, silnice už jsou rozeznatelné hůře.
347 ... Podobně jako kombinace 345, ale hodí se lépe pro zobrazení požářišť.
354 ... Působí jako barevný infračervený obrázek - viz obr. č.12.
374 ... Podobně jako 354.
457 ... Ukazuje půdní druh (jílovitá, hlinitá, písčitá).
Obr č.11: Barevná syntéza v nepravých barvách - 234
Obr č.12: Barevná syntéza v nepravých barvách - 354


[ Zpracování obrazových dat ]