153ZODH / 12. cvičení
< Stránky předmětu • Předchozí cvičení
Osnova
Pokročilé analýzy jsou založeny na monitorování zemského povrchu a atmosféry v různých prostorových rozlišení a časových intervalech. Zaměříme se na analýzu časové řady dat MODIS (Moderate Resolution Imaging Spectroradiometer, satelit Terra a Aqua NASA). K dispozici jsou následující data:
- Aqua (1:30h)
- Terra (10:30h)
- Aqua (13:30h)
- Terra (22:30h)
Seznam příkazů
Zobrazení dat
Nejprve přidáme mapset obsahující data MODIS do cesty.
g.mapsets add=modis20021st g.list rast mapset=modis20021st
Data zobrazíme.
# noc, 22:30h, teplota ve stupních Celsia d.rast.leg map=terra_1st_night_20020921 d.vect map=us_natlas_hydrogp type=boundary color=blue d.vect map=us_natlas_urban type=boundary d.vect map=us_natlas_urban diplay=attrb attrcol=NANE type=centroid where="AREA >= 0.002" lsize=12
# den, 10:30h d.rast.leg map=terra_1st_night_20020922 d.vect map=us_natlas_hydrogp type=boundary color=blue d.vect map=us_natlas_urban type=boundary d.vect map=us_natlas_urban diplay=attrb attrcol=NANE type=centroid where="AREA >= 0.002" lsize=12
Diferenci LST noc-den vypočteme pomocí r.mapcalc.
r.mapcalc "diff_den_noc = terra_1st_day20020922 - terra_1st_day_20020921" d.rast.leg map=diff_den_noc d.vect map=us_natlas_hydrogp type=boundary color=black d.vect map=us_natlas_urban type=boundary color=brown d.vect map=us_natlas_urban diplay=attrb attrcol=NANE type=centroid where="AREA >= 0.002" lcolor=black lsize=12
Závěr: Na rozdíl od zastavěných ploch se teplotně stabilní se jeví vodní plochy.
Indikátory časových řad
Data využijeme k výpočtu řady indikátorů časových řad. Data jsou již předzpracována, tj. byly odstraněny digitální hodnoty odpovídající teplotě mraků. Během výpočtu musíme s NULL hodnotami počítat, např. při určování průměrné teploty.
Vypočteme denní LST průměry v měsíci září 2002, pomocí modulu r.series.
for d in `seq 1 30` ; do
DAY=`echo $d | awk '{printf "%02d\n", $1}'`
LIST=`g.mlist type=rast pattern="*lst*200209${DAY}" separator=","`
echo "$LIST"
r.series -n input=$LIST output=lst_200209${DAY}_avg method=average
d.rast.leg map=1st_200209${DAY}_avg
d.vect map=us_natlas_urban type=boundary
done