Difference between revisions of "HOWTO import SRTM elevation data"

From GRASS-Wiki
Jump to: navigation, search
m (See also)
m (Import of multiple SRTM tiles in one step)
Line 22: Line 22:
  
 
Original SRTM tiles:
 
Original SRTM tiles:
 +
<source lang=bash>
 
  # unzip all
 
  # unzip all
 
  for i in *.hgt.zip ; do unzip  $i ; done
 
  for i in *.hgt.zip ; do unzip  $i ; done
Line 28: Line 29:
 
  # import
 
  # import
 
  r.in.gdal input=srtm_mosaik output=srtm_mosaik
 
  r.in.gdal input=srtm_mosaik output=srtm_mosaik
 
+
</source>
 
CGIAR SRTM tiles:
 
CGIAR SRTM tiles:
 +
<source lang=bash>
 
  # create mosaik (optionally reproject on the fly with -t_srs)
 
  # create mosaik (optionally reproject on the fly with -t_srs)
 
  gdalwarp *.tif srtm_mosaik.tif
 
  gdalwarp *.tif srtm_mosaik.tif
 
  # import
 
  # import
 
  r.in.gdal input=srtm_mosaik output=srtm_mosaik
 
  r.in.gdal input=srtm_mosaik output=srtm_mosaik
 +
</source>
  
 
== Void filling of original SRTM tiles ==
 
== Void filling of original SRTM tiles ==

Revision as of 07:51, 23 June 2012

Background: Space Shuttle Radar Topography Mission

Download

Import of original SRTM tiles in HGT format

Use the r.in.srtm module for automated import into a Latitude-Longitude location.

See:

  • M.H. Bowman. Mapping freely available high resolution global elevation and vector data in GRASS. GRASS Newsletter, Vol.3, pp 7-10, June 2005. [ PDF ]

Import of CGIAR SRTM void-filled tiles

  • SRTM 90 data in GeoTIFF format are available from http://srtm.csi.cgiar.org/. These come in 5 degree by 5 degree tiles. Horizontal resolution: 3 arc second (90 m at equator). Coverage: world. Import with r.in.gdal

Import of multiple SRTM tiles in one step

To simplify the import of multiple SRTM tiles, the tiles can be mosaicked first with gdalwarp:

Original SRTM tiles:

 # unzip all
 for i in *.hgt.zip ; do unzip  $i ; done
 # create mosaik (optionally reproject on the fly with -t_srs)
 gdalwarp *.hgt srtm_mosaik.tif
 # import
 r.in.gdal input=srtm_mosaik output=srtm_mosaik

CGIAR SRTM tiles:

 # create mosaik (optionally reproject on the fly with -t_srs)
 gdalwarp *.tif srtm_mosaik.tif
 # import
 r.in.gdal input=srtm_mosaik output=srtm_mosaik

Void filling of original SRTM tiles

The module r.fillnulls can be used. See screenshot

See:

  • M. Neteler, 2005. SRTM and VMAP0 data in OGR and GRASS. GRASS Newsletter, Vol.3, pp 2-6, June 2005. [ PDF | Erratum 1/2006 (txt) ]

See also