Digital Elevation Model (DEM) Processing

From AWF-Wiki
(Difference between revisions)
Jump to: navigation, search
(→‎Loading and reprojection of DEM data)
(→‎Clipping and Resampling of DEM data)
Line 27: Line 27:
 
## Open the raster layer properties by right-clicking the clipped DEM map in the [[TOC]] and selecting {{button|text=Properties}}, or by simply double clicking.  
 
## Open the raster layer properties by right-clicking the clipped DEM map in the [[TOC]] and selecting {{button|text=Properties}}, or by simply double clicking.  
 
## Select the {{button|text=Style}} tab. Under {{button|text=Render type}} select {{button|text=Singleband pseudocolor}}.
 
## Select the {{button|text=Style}} tab. Under {{button|text=Render type}} select {{button|text=Singleband pseudocolor}}.
−
## In the {{button|text=Generate new color map}} section select {{button|text=Spectral}}. In the {{button|text=Mode}} pulldown menu select {{button|text=Continuous}}. Set the number of {{button|text=Classes}} to 10 and click {{button|text=Classify}}.  
+
## In the {{button|text=Generate new color map}} section select {{button|text=Spectral}}. In the {{button|text=Mode}} pulldown menu select {{button|text=Continuous}}. Set the number of {{button|text=Classes}} to 20 and click {{button|text=Classify}}.  
 
## To finish, click {{button|text=Apply}}. If your're content with you settings, click {{button|text=OK}} to close the window (figure '''D''').  
 
## To finish, click {{button|text=Apply}}. If your're content with you settings, click {{button|text=OK}} to close the window (figure '''D''').  
−
# SAGA tools for resampling
 
−
## If not open, activate the processing toolbox by selecting {{mitem|text=Processing --> Toolbox}}.
 
−
## In the toolbox, browse for {{mitem|text= SAGA --> Grid-Spline --> Multilevel B-Spline Interpolation (from grid)}} or simply type {{typed|text=Multilevel}} into the search field to find the tool. Open it with a double-click.
 
−
## In the appearing menu, select the following options:
 
−
##* The clipped layer from the {{button|text=Grid}} pulldown menu.
 
−
##* For the {{button|text=Output extent}}, select the clipped layer by clicking {{button|text=...}} and selecting {{button|text=Use layer/canvas extent}}.
 
−
##* Set the {{button|text=Cell size}} to 10.
 
−
##* In the {{button|text=Grid}} input field at the bottom of the menu, select the directory and name for the output file.
 
−
##* Choose an output name and directory.
 
−
##* Leave the rest of the menu points as they are and click {{button|text=Run}}. The resampled layer will be added to the [[QGIS GUI|canvas]].
 
  
 
{| class="wikitable"
 
{| class="wikitable"
Line 44: Line 34:
 
| style="border: 0pt" | [[file:RemSens_Exercise05_04.png|thumb|left|450px|'''Figure D:''' Clipped DEM map displayed in pseudocolor]]
 
| style="border: 0pt" | [[file:RemSens_Exercise05_04.png|thumb|left|450px|'''Figure D:''' Clipped DEM map displayed in pseudocolor]]
 
|}
 
|}
−
 
  
 
==3D Layer styling==
 
==3D Layer styling==

Revision as of 11:34, 12 April 2018

In this article you will learn how to pre-process and enhance digital elevation model (DEM) data, as well as how to create a 3D-Model

Loading and reprojection of DEM data

  1. Add a ditial elevation model as raster layer n51_e009_e0010_merge.tif to a new QGIS project (from geodata\raster\).
  1. Reprojecting layers (warp)
  • In the search engine of the Processing Toolbox, type warp and select Warp (reproject) under Projections of GDAL.
  • Under the Parameter tab, select the dem file as input layer.
  • Select EPSG:4326 (WGS84 / Long Lat) under the Source SRS tab.
  • Select EPSG:32632 (WGS84 /UTM Zon 32 N) under the Target SRS tab.
  • Enter 30 as the output file resolution in meters.
  • Select the Resampling method bilinear.
  • Specify the Path and name of the output file.
  • Click on Run to execute the algorithm.

Qgis warp srtm.png

Clipping and Resampling of DEM data

  1. Create a boundary polygon
    1. Add the raster map geodata/raster/s2/Subset_S2A_MSIL2A_20170619T
    2. Select Vector --> Research tools --> Polygon from layer extent. Choose the Sentinel-2 map as input, select an output path and file name and confirm with OK. Confirm the following prompt to let the layer be added to the canvas (figure C).
  2. Buffer the boundary polygon
    • Choose Vector --> Geoprocessing tools --> Buffer and select the boundary polygon created before as input. Select a reasonable name and output directory under Output shapefile. Confirm with OK and let the output layer be added to the canvas. Place it below the boundary polygon in the TOC to see the effect of buffering.
  3. Clip the DEM map
    • Open the clipping tool under Raster --> Extractions --> Clipper. Select the reprojected DEM file as input and select the output directory and file name as usual. Select the Mask layer radio button to select the buffered boundary layer as blueprint for clipping. Don't forget to enter an output path and file name. Confirm with OK and let the result be added to the map.
  4. Map visualization
    1. Open the raster layer properties by right-clicking the clipped DEM map in the TOC and selecting Properties, or by simply double clicking.
    2. Select the Style tab. Under Render type select Singleband pseudocolor.
    3. In the Generate new color map section select Spectral. In the Mode pulldown menu select Continuous. Set the number of Classes to 20 and click Classify.
    4. To finish, click Apply. If your're content with you settings, click OK to close the window (figure D).
Figure C: Boundary polygon produced by polygon to layer
Figure D: Clipped DEM map displayed in pseudocolor

3D Layer styling

  1. Add the raster layer dem_10m.tif
  2. Right click on the layer name and then click Duplicate from the menu
  3. Open the layer settings menu. Go to Style and change the render type to Hillshade and press OK
  4. Right click on the duplicated file dem_10m copy and open the layer settings menu
  5. Change the render type to Singleband Pseudocolor and select a new color ramp for elevation
  6. Change the Blending mode under Color rendering to Burn
  7. Change the layer transparency to 30%
Figure E: Hillshade
Figure F: Pseudocolor topography style
Figure G: Colorized hillshade
Retrieved from "http://wiki.awf.forst.uni-goettingen.de/wiki/index.php?title=Digital_Elevation_Model_(DEM)_Processing&oldid=14291"
Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export