HOWTO import SRTM elevation data: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
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 14: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