153ZODH / 12. cvičení: Porovnání verzí
(→Zobrazení dat: avg) |
|||
| Řádek 35: | Řádek 35: | ||
d.vect map=us_natlas_urban type=boundary | 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 | 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í {{GrassPrikaz|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. | |||
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 | |||
---- | ---- | ||
Verze z 5. 1. 2009, 17:29
< 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.
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