SRTM

Z FreeGIS portál
Přejít na: navigace, hledání

SRTM (Shuttle Radar Topography Mission) je mezinárodní výzkumný program, zastřešený NGA a NASA, jehož cílem je získat údaje nadmořských výšek k vytvoření nejkompletnější digitální topografické databáze ve vysokém rozlišení. Mise SRTM sestávala ze speciálně upravených radarových systémů jež byly na palubě raketoplánu Endeavour během 11denní mise STS-99 v únoru 2000. Více informací zde.

Stažení dat

Data jsou ke stažení ve formátu GeoTIFF a HGT z adresy http://dds.cr.usgs.gov/srtm/.

Zpracování dat v GRASS GIS

Další informace na wiki projektu GRASS.

Import

Před importem je nutné vytvořit pro data lokaci - viz postup pro vytvoření lokace na základě geodat. V tomto případě jsou data, na základě kterých je lokace vytvořena, zároveň i do této lokace naimportována. Alternativně lze lokaci vytvořit na základě EPSG 4326.

Vytvoření lokace pro SRTM z příkazové řádky
grass70 -c EPSG:4326 /opt/grassdata/srtm

Import se provede standardním postupem - viz návod pro import rastrových dat.

Import dat SRTM z wxGUI

Data ve formátu HGT lze naimportovat pomocí specializovaného nástroje r.in.srtm.

Nastavení tabulky barev

Tabulku barev nastavíme pomocí modulu r.colors.

Nastavení tabulky barev z kontextového menu Layer Manageru
Nastavení tabulky barev pro data SRTM
Nastavení tabulky barev SRTM z příkazové řádky
r.colors map=SRTM_u01_n060w150 color=srtm

Předzpracování (no-data)

Vzhledem k použité technologii (interferometrická měření) obsahují data oblasti bez pokrytí (tzv. díry). Hodnoty bez informací jsou kódovány jako hodnota -32768 (viz obr níže). Tyto hodnoty (tzv. no-data) převedeme na hodnotu NULL, která se v systému GRASS používá pro vyjádření digitálních hodnot bez informace. Tuto operaci provedeme pomocí modulu r.null. Daná rastrová mapa musí být uložena v aktuální mapsetu.

Dotaz na hodnotu bez informací (tzv. no-data)
r.null map=SRTM_u01_n060w150 setnull=-32768
Dotaz na hodnotu bez informací po nahrazení hodnotou NULL

Vyplnění děr

Naším úkolem bude chybějící data dopočítat pomocí vhodné interpolace. Pro tento účel můžeme použít modul r.fillnulls, který používá metodu interpolace RST.

r.fillnulls input=SRTM_u01_n060w150T output=SRTM_u01_n060w150_fill method=rst

V případě neúspěchu (viz obr. níže)

Vyplnění děr (RST interpolace - neúspěch)

lze zvolit alternativní postup založený na kubické interpolaci (method=cubic).

r.fillnulls input=SRTM_u01_n060w150T output=SRTM_u01_n060w150_fill method=cubic
Vyplnění děr (kubická interpolace - úspěch)