DEM: Difference between revisions

From OCAD Wiki - English
Jump to navigation Jump to search
Line 51: Line 51:
[[File:DemImportWizard5.png]]
[[File:DemImportWizard5.png]]


In this example OCAD transfrom the DEM data from ''US State Plane (Feet)'' to the metric UTM Zone 18 North.
In this example OCAD transform the DEM data from ''US State Plane (Feet)'' to the metric UTM Zone 18 North.


Optional it is possible to convert also the height values from Feet to Meters.
Optional it is possible to convert also the height values from Feet to Meters.

Revision as of 14:15, 8 February 2016

A DEM (Digital Elevation Model) contains points with elevation data. DEM Data are based on LIDAR (Light Detection and Ranging) technology measurement, also known as Airborne Laser Scanning. There are DEM with point data arranged in a regular grid with a constant distance between the points. This distance is called cell size. Other DEM contain data points arranged irregularly (cloud-model).

Read more about this topic: http://en.wikipedia.org/wiki/Digital_elevation_model

DEM.jpg

DEM Import Data Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

The DEM Import Wizard is able to import files with regular and irregular DEM data. Supported DEM data formats are:

  • ESRI ASCII Grid (*.asc)
  • Raw data ASCII XYZ file (*.xyz)
  • ASCII Grid XYZ file (*.xyz)
  • LAS file (*.las)
  • compressed LAS file (*.laz or *.zlas)
  • SRTM file (*.hgt)

Start Import DEM Wizard

Select Import from the DEM menu. The DEM Import Wizard appears. Otherwise you can also drag and drop an DEM file from Windows File Explorer onto the OCAD drawing area.

DemImportWizard1.png

First use Add button to add at least 1 DEM to the DEM Import Wizard dialog. You have three different options:

  • Choose a las files which often contain DTM and DSM data
  • Choose only DTM data
  • Choose DTM and DSM data

Choose las Files

  • After adding a las file you can select the item and click the LAS Info button. The LAS File Header dialog opens wich shows the information from the file header like the Number of Points by Return or the Extent.

DemImportWizard2.png

Choose DTM Data

Add only the DTM file(s).

Please note that not all functions (e.g. vegetation height) are available if only DTM data are available.

DemImportWizard3.png

Choose DTM and DSM Data

Add the DTM file(s) on the left side and the DSM file(s) on the right.

DemImportWizard4.png

If the DSM have an other extent than the DTM then OCAD cuts to the DSM to the extent of the DTM or the choosen extent.

Coordinate System

Usually the DEM data are in the same coordinate system as the OCAD map. So no transformation is necessary.

OCAD supports the transformation during the DEM import. Click the Choose button to select the coordinate system of the DEM data.

DemImportWizard5.png

In this example OCAD transform the DEM data from US State Plane (Feet) to the metric UTM Zone 18 North.

Optional it is possible to convert also the height values from Feet to Meters.

Extent

Click the Next button. OCAD imports the file(s). Depends on the data this step can takes some minutes.

OCAD changes to the Settings dialog.

Settings

  1. Analyze the DEM to get some information about the DEM like extent, cell size etc. The DEM Import dialog extends after analyzing a LAS file and provides some additional settings.
  2. If the source file is a regular grid (data type of import files = grid), the Cell size box is set to read only. OCAD sets the same cell size for the imported DEM. For irregular DEM data source (data type of import files = raw) the cell size can be set in the 'cell size' box. For these DEM's OCAD interpolates a regular grid with the specified cell size during the import.
  3. At the end of the import procedure the imported DEM must be saved in the OCAD DEM file format (*.ocdDem) and it is loaded to the OCAD map.
  4. To see the extent of the loaded DEM activate the Show Frame menu item.


hint The Cell size range is between 0.01 and 650 m.
Hint.jpg The DEM import limit depends on the amount of memory available.

Coordinate System:
DEM data are sometimes in another coordinate system than the map. OCAD provides a new function to transform DEM data during the DEM import.

In this example bellow the map file with the loaded orthophoto has the UTM Zone 18 North coordinate system, but the Lidar data are only available in US State Plane New York Long Island (US Feet).

Click the Choose button in the DEM Import dialog to choose the coordinate system of the DEM file and check the option Convert height values from feet to meters. Click Analyze.

OCAD transforms all data points during the analysis procedure.


DemImportLas.PNG

LAS data import:

  • Click the Choose DTM button to create a terrain model from LAS data. This selects ground classification and the Last return option.
  • Click the Choose DSM button to create a surface model from LAS data. This selects all classifications and the First return option.

Hint.jpg During LAS data import a classification image and a intensity image are created and loaded as background maps. A difference image between the first and last return is created if the All returns option in the DEM Import dialog is choosen.
The classification tiff image has the following color index:

  • 0 Created, never classified (light gray)
  • 1 unclassified (red)
  • 2 Ground (yellow)
  • 3 Low Vegetation (light green)
  • 4 Medium Vegetation (green)
  • 5 High Vegetation (dark green)
  • 6 Building (Generic) (gray)
  • 7 Low point (noise) (gray)
  • 9 Water (blue)
  • 10 Bridge (brown)
  • 11 Road (brown)
  • 12 Overlap pointd (gray)


Hint.jpg A Radar beam splits as it hits objects. The result are multiple returns. LIDAR can collect up to four returns. The difference between first and last return can show object height. The last return doesn’t always reach the ground.

DemImportLas MultipleReturn.jpg
Source: Lohani, Bharat. Airborne Altimetric LiDAR: Principle, Data Collection, Processing and Applications.

Hint.jpg Compressed LAS file (*.laz): OCAD uses the laszip.exe tool form Martin Isenburg to decompress the laz file.

SRTM hgt data import:

  • This is a world wide available DEM data set from the Shuttle Radar Topography Mission (SRTM).
  • The data files and documentation: [1]

Hint.jpg SRTM hgt data import requires that a coordinate system set in the OCAD map file.

Hint.jpg If the ocd file is untiteled OCAD exports the DEM files in the temporary user directory (C:\Users\XXX\AppData\Roaming\OCAD\OCAD 12\Tmp).

Open Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.This function is available in OCAD Starter.This function is available in OCAD Course Setting.

Open an OCAD DEM file (*ocdDem).


Show Frame Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.This function is available in OCAD Starter.This function is available in OCAD Course Setting.

Shows blue rectangle with the extent of the loaded DEM.


Resize Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

Resize the loaded OCAD DEM file (make a subset) and save it as a new OCAD DEM file.


Info Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.This function is available in OCAD Starter.This function is available in OCAD Course Setting.

Shows information about OCAD DEM file.

Hint.jpg When moving the mouse cursor over the file name then the file name with path appears.


Close Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.This function is available in OCAD Starter.This function is available in OCAD Course Setting.

Close OCAD DEM file.


Merge DEM Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

miniatur

Choose Merge from DEM menu to merge two different DEMs.


  • DEM 1: Choose the first DEM.
  • DEM 2: Choose the second DEM.


Hint.jpg The two DEMs must have the same cell size.

Calculate DEM Difference Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.


CalculateDEMDifference.PNG
  1. Choose Calculate DEM Difference from DEM menu.
  2. The Calculate DEM Difference dialog box appears.
  3. Add Upper DEM = DSM data file
  4. Add Lower DEM = DEM data file
  5. Click OK.


To visualize the DEM difference choose Classify Vegetation Height

Hint.jpg The extent of the Upper DEM and Lower DEM can be different. OCAD takes the overlapping area for the new DEM.

Hint.jpg The cell size of Upper DEM and Lower DEM can be different. OCAD takes the cell size of Upper DEM for the new DEM.

Example

This is an example to show what can result from the Calculate DEM Difference function.

This is a DTM (Digital Terrain Model) with a cell size of 5m shown as hypsometric map:

DTM5mHypso.PNG

The next picture shows the DSM (Digital Surface Model) of the same area with a cell size of 2m as hypsometric map. The buildings (northern part) and forest (south western part) are slightly visible.

DSM2mHypso.PNG

The next picture shows a Difference DEM with a cell size of 2m shown as raster background map after calculating the DEM Difference. In addition, heights were colored using the Classify Vegetation Height function.

ClassifiedDifferenceDEM.PNG

  • The area with no difference of the DTM and the DSM are displayed white.
  • A height difference up to 15m appears red.
  • The greater the difference, the greener an area appears.

When moving the mouse cursor over the map the difference is shown in the Status Bar.

Data source: Test data Wabern from swisstopo.

Create Contour Lines Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

  • Choose Create Contour Lines from DEM menu.
  • The Create Contour Lines dialog box appears.

CreateContourLines.png

This function calculates contour lines based on the loaded DEM.

This function can take a lot of time. OCAD starts to calculate the contour lines from the lowest elevation (minimun contour). During the calculation process the current elevation is shown in the left status bar.

  • Define 1-3 contour intervals (for example 1m, 5m, 25m) and a pre-selected line symbol (according to the first three line symbols in the settings) for each type appears (can be changed).
Hint.jpg Contour interval values can be entered manually or chosen from the list.
  • Specify the minimum (lowest) and maximum (highest) contour for the calculation.
  • Activate the option Split DEM into tiles to speed up this calculate process. In this case OCAD splits the dem into small tiles and calculates tile by tile. Finally the contour lines are also splites at the tile borders. It is not visible on the screen but you see it when selecting a contour line. Use the Merge Contour Lines By Selected Symbols function to merge the contour lines if needed.

ContourLinesSplitted.PNG

The selected contour line is cutted at the tile borders.

ContourLines.PNG

Create Hypsometric Map Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

miniatur
This function calculates a grayscale or colored hypsometric map as GeoTIFF file.
Optionally it is directly loaded as background map.


Two different types of hypsometric maps can be edited:

  • Grayscale hypsometric map
  • Colored hypsometric map

HypsometricMap.PNG

Create Hill Shading Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

miniatur
This function calculates a shaded relief picture (hill shading).


There are two calculation methods available:
  • Slope shading is optimized to see outlines of features like paths in a slope.
  • Slope shading combined with oblique light shading is the recommended method if the hill shading should be used as a background relief of a map.
    Optionally the calculated hill shading is directly loaded as background map.
Aside from the chosen method, there is to define the Resolution and the Interpolation mode (if chosen).
The default interpolation mode is Bicubic, but there are 7 other algorithm.
Additionaly an Azimuth and a Declination of the light source has to be set. Standard settings are 315° (north-west) and 45°.
An Exaggeration factor of 4 is pre-selected and can be altered.

The Preview... button allows you get a first impression of the current setting for your hillshading and can be use to optimize it.

Click and move with the left mouse button allows to pan the current view.

The Preview... button is only enabeld when importing a dem grid file or an ocddem is already loaded. This function is disabled when importing las/laz files.


hintTo detect point and line objects like paths or watercourses from DEM we recommend using the same resolution as the DEM. To create a relief and if the DEM cell size isn’t high then we recommend to set a smaller resolution.


hintThe default export file format is JPEG and creates an 8 bit JPEG with in grayscale and a world file with the geo reference. If you decide to export the file in TIFF-format, only with resolution option ‘DEM cell size’ then OCAD creates an 8 bit grayscale tiff with color palette. Otherwise a 24 bit RGB tiff. HillShading.PNG

Error Message Bitmap is too big

The size of the hill shading image is limited. If this error message appears please enter a bigger value for the interpolation.

Calculate Slope Gradient Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

CalculateSlopeGradient.PNG

Choose Calculate Slope Gradient from DEM menu.
The Calculate Slope Gradient dialog box appears.


Select one of two different methods:

  • Continous (<x° = grayscale / >x° = black)
  • Black/White (<x° = white / >x° = black)

The resulting picture can be used to identify cliffs an rock faces. The result can sometimes be significantly improved with a slight adjustment of the gradient (between 42-45 degrees).

Slope gradient also shows paths or relief features independent from a azimut like the Hill Shading.

SlopeGradient.PNG

Classify Vegetation Height Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

miniatur

Choose Classify Vegetation Height from DEM menu.


There are two different options to show vegetation height classification:

  • Gray scale classification with options: Linear, Quadratic negative, Quadratic positive
  • Colored classification: Define classes with a height and color range
- Split a class into two classes by clicking the Split class button
- Remove a class by clicking the Remove class button
- Load the settings from a text file by clicking the Load button
- Save the settings to a text file by clicking the Save button
- Reset the classes and colors to the default settings by clicking the Reset classes to default button

VegetationHeightGrayAndColor.PNG


Merge Contour Lines By Selected Symbols Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

  • Select the contour line symbols
  • Choose Merge Contour Lines By Selected Symbols from DEM menu.
  • The Merge Contour Lines By Selected Symbols dialog box appears.

MergeContourLinesBySelectedSymbols.png

Create Profile Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.This function is available in OCAD Starter.This function is available in OCAD Course Setting.

ProfileExample.PNG

Find more information about this function on the DEM Profile page.

Export Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

  • Choose Export from DEM menu.
  • The DEM Export dialog box appears.

DEMExport.PNG

The function exports the loaded DEM file to the following formats:

  • ESRI ASCII Grid
  • ASCII Grid XYZ
  • OCAD 11 DEM

Select Create tiles for large DEMs to create tiles fron 1001x1001 points.

The OCAD 12 DEM's cannot be opened in OCAD 11. This export function creates a OCAD 11 compatible DEM.

ocdLas File Space.PNGThis function is available in OCAD Professional.This function is available in OCAD Orienteering.

Learn more about the ocdLas File functions on the ocdLas File page.



Back to Main Page

Previous Chapter: Printing Maps

Next Chapter: GPS