#Pouzito z http://gama.fsv.cvut.cz/wiki/index.php/Zpracov%C3%A1n%C3%AD_obrazov%C3%BDch_dat #import dat grass /tmp/harta/cviko r.in.gdal input=/data/Landsat_7_ETM/vstup/189_25_020504_ms_vyrez.img output=tm location=tm exit grass /tmp/tm/PERMANENT #barevne kompozice d.mon x0 d.rast map=tm.1 d.rgb blue=tm.1 green=tm.2 red=tm.3 d.histogram tm.1 r.colors tm.1 col=grey.eq r.colors tm.2 col=grey.eq r.colors tm.3 col=grey.eq d.redraw d.rgb blue=tm.2 green=tm.3 red=tm.4 d.redraw d.rgb blue=tm.3 green=tm.5 red=tm.4 #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 (tm.6) ... 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 (tm.6) ... 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 (tm.6) ... Ukazuje půdní druh (jílovitá, hlinitá, písčitá). #ndvi r.mapcalc 'ndvi = 1.0 * (tm.4 - tm.3) / (tm.4 + tm.3)' r.colors ndvi col=rules < tb-ndvi.txt d.redraw d.rast ndvi d.rast.leg map=ndvi r.mapcalc 'temp1 = 100 * ndvi' r.reclass input=temp1 output=r_ndvi < reklas-ndvi.txt r.colors r_ndvi col=rules < tb-r_ndvi.txt d.rast r_ndvi d.rast.leg r_ndvi #nerizena klasifikace i.group group=tmx subgroup=tmx input=tm.1,tm.2,tm.3,tm.4,tm.5,tm.6 i.cluster group=tmx subgroup=tmx sigfile=tm_nerpri20 reportfile=tm_rep20.txt classes=20 i.maxlik group=tmx subgroup=tmx sigfile=tm_nerpri20 class=tm_nerkl20 reject=tm_nerza20 d.rast.leg tm_nerkl20