Framework project: Digital Elevation Models (DEMs)

DEM process setup (Arctic DEM)

Introduction

A Digital Elevation Model (DEM) is a general representation of heights above or below a reference geoid (where the geoid is a mathematical shape of an orb, for instance Earth). A Digital Terrain Model (DTM) refers to the natural surface of the Earth, i.e. disregarding vegetation, buildings and other 3 dimensional (3D) human structures. A Digital Surface Model (DSM) instead includes all human 3D objects. Theoretically the conceptual definitions of DTM and DSM are more strict, with DEM being more widely defined. In reality the three concepts are mixed and note seldom used interchangeably. The Arctic DEM is defined as a DSM.

Overview

In this post you will add the ArcticDEM to your GeoImagine Framework. To be effective and cover the correct region, you have to start with defining a default region that fits the ArcticDEM. Then you have to link that default region to a project owned by your user (i.e. the system superuser). When you have created your own project and region you can import the ArcticDEM as an ancillary layer, and then tile it to one the predefined systems available in the Framework (i.e. the MODIS SIN grid, or an EASE-grid).

Arctic DEM

The Arctic DEM is an effort led by the United States National Geospatial-Intelligence Agency (NGA) and the National Science Foundation (NSF) to create a high resolution DEM of the Arctic land area. The DEM is created from optical stereo imagery, high-performance computing, and open source photogrammetry software. The DEM covers all land areas north of 60 degrees, all territory of Greenland, the State of Alaska in entirety, and the Kamchatka Peninsula of the Russian Federation.

Arctic Digital Elevation Model (DEM) in polar stereographic projection (orignal download format).

Create a default regions

To accommodate and work with the Arctic region in Karttur´s GeoImagine Framework, you need to create corresponding default regions. In my version I created three new default regions using the Framework command DefaultRegionFromCoords as defined in the json file below:

Create user tracts linked to arctic regions

The Framework is built in such a way that no user, including the super users, can use default regions for actual data processing. All processing can only be done vis-a-vis regions that are owned by the logged in user. Thus you have to create user regions (called “tracts”) that correspond to the default regions that you created above. The Framework process to use is ManageDefRegProj:

Data access and download

The Artic DEM is freely available from GRACE TELLUS. The data are available through ftp (ftp.data.pgc.umn.edu/ and then under the folder tree elev/dem/setsm/ArcticDEM). The mosaics at 100 m and coarser resolution are available as single file mosaics.

Karttur’s GeoImagine Framework does not include any process for downloading the Arctic DEM. To download the Arctic DEM, please use any FTP client (for example Filezilla).

At time of writing, the Arctic DEM version is 3.0, and the two files that I downloaded where the 100 and 500 m mosaics:

  • arcticdem_mosaic_100m_v3.0.tif
  • arcticdem_mosaic_500m_v3.0.tif

Import as ancillary data

DEM data is imported to the Framework as ancillary data, the process name is OrganizeAncillary. In the example import below both the 100 and 500 m versions are imported. Note that the project and tract for importing are set to the extended arctic region, as defined above.

One of the basic concepts behind Karttur´s GeoImagine Framework is to use tiling systems for processing the Earth’s surface piecewise. The main tiling systems in use are :

  • MODIS sinusoidal grid
  • EASE-grids (north polar, south polar and tropical)
  • MGRS (UTM based)

The json command file below applies the process LinkDefaultRegionTiles to identify all MODIS SIN tiles that overlap with the Arctic regions (defined in the json object above). To change to any other default system you just have to change the system parameter under userproject (e.g. to ease2n for EASE-grid north [EPSG:6931]).

If your Framework lack the MODIS tiles, you need to setup the regions, including the MODIS tiling system as described in the post on setting up regions [LINK NEEDED]. Setting the parameter defregmask to land limits the identified tiles to those that actually cover at least a single pixel of land surface.

To run the actual tiling below, you must have identified the tiles covering the region.

Tiling

At this stage your Framework should contain:

  • ArcticDEM (imported as ancillary layer)
  • Arctic default regions coinciding with the spatial extent of the arcticDEM
  • User project and tract built from the ArcticDEM fitted default region(s)
  • MODIS tiles coinciding with the ArcticDEM fitted default region(s)

If you have all of the above, you can tile the ArcticDEM data to the MODIS SIN grid system using the process TileAncillaryRegion.

The json above generates a shell script that you must run manually. The code snippet below shows the code for generating a single tile (h18v02 central Scandinavia) shown in the figure below.

mkdir -p /Volumes/karttur3tb/modis/NGA-NSF/tiles/dem/h18/v02/0;
/Library/Frameworks/GDAL.framework/Versions/2.1/Programs/gdalwarp -t_srs "+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs " -tr 463.313000000 463.313000000  -te 0.0 6671703.118 1111950.51966667 7783653.63766667  -r near  /Volumes/karttur3tb/ancillary/NGA-NSF/region/dem/arcticdemext/0/dem500_arcticdem_arcticdemext_0_v30.tif /Volumes/karttur3tb/modis/NGA-NSF/tiles/dem/h18/v02/0/dem500_arcticdem_h18v02_0_v30.tif;
Arctic Digital Elevation Model (DEM) projected to MODIS SIN grid h18v02 (central Scandinavia).

As you can see in the figure above there are quite a lot of no data in the Arctic DEM. Alternative DEMs and how to download them is described in the post Digital Elevation Models (DEMs).