Python for Processing Digital Elevation Models

English Version | 中文版

This is a Python implementation of transforming, projecting, converting, and visualizing Digital Elevation Models as well as reading their elevation of given locations.

Table of Contents

• Terms and Abbreviations
• Functions
• DEM Data
  • ASTERGDEM
  • EUDEM
• Dependencies
• Pipeline of Processing ASTERGDEM
• Usage
  • Usage of ASTERGDEM
  • Usage of EUDEM
• Results
  • DEM Images of London
  • Elevation of London
• References

Terms and Abbreviations

Functions

• Transforming DEM from one GCS (e.g., WGS-84) to another GCS (e.g., OSGB-36).

• Projecting DEM from one GCS (e.g., OSGB-36) to the related PCS (e.g., BNG).

• Converting DEM from one PCS (e.g., LAEA) to the related GCS (e.g., ETRS-89).

• Visualizing DEM in PCS as a 2D image.

• Reading the elevation of given locations from DEM in GCS.

DEM Data

ASTERGDEM

• Update To access and download ASTERGDEM data set, see link.

• ASTERGDEMv2.0 download link: https://earthexplorer.usgs.gov/

• See Introduction of ASTGDEMv20.pdf for detailed introduction and user guide (recommended).

EUDEM

• EUDEMv1.1 download link: https://land.copernicus.eu/imagery-in-situ/eu-dem/eu-dem-v1.1

• See here for detailed introduction and user guide.

Dependencies

• gdal (osgeo) 1.11.3
• matplotlib 3.0.2
• numpy 1.15.4
• pandas 0.23.4

The following procedures (1-2) for installing GDAL are recapitulated from mothergeo.

1. Install GDAL Development Libraries and export environment variables for the compiler in Terminal.

$ sudo apt-get install libgdal-dev

$ export CPLUS_INCLUDE_PATH=/usr/include/gdal
$ export C_INCLUDE_PATH=/usr/include/gdal

2. Install GDAL Python Libraries in Terminal.

$ pip install GDAL

If it comes with the error: cpl_vsi_error.h: No such file or directory, please try the following procedures.

3. Check the required version of GDAL Python Libraries in Terminal.

$ gdal-config --version

4. Download the source file (e.g., gdal-1.11.3.tar.gz) of the related GDAL version (e.g., 1.11.3) from here.

5. Manually install GDAL Python Libraries in Terminal.

$ cd path/of/downloaded/gdal/package

~$ tar -xvzf gdal-{VERSION}.tar.gz
# E.g., ~$ tar -xvzf gdal-1.11.3.tar.gz

~$ cd extracted/gdal/folder
# E.g., ~$ cd gdal-1.11.3

~$ cd swig
~$ cd python
# The setup.py exists in this directory.

~$ python setup.py build_ext --include-dirs=/usr/include/gdal/
~$ python setup.py install

6. Run >>> from osgeo import gdal in Python. If no error occurs, the installation is accomplished.

Pipeline of Processing ASTERGDEM

Usage

The directory of source DEMs is shown as below.

└── Data/
    ├── DATA_ASTGDEMv20/
    │   └── EPSG4326_s/
    │       ├── ASTGTM2_N51W001_dem.tif
    │       └── ASTGTM2_N51E000_dem.tif
    └── DATA_EUDEMv11/
        └── EUDEMv11_EPSG3035.tif

Usage of ASTERGDEM

1. Download DEMs ASTGTM2_N51W001_dem.tif and ASTGTM2_N51E000_dem.tif from here.

2. Copy DEMs to the folder DATA/DATA_ASTGDEMv20/EPSG4326_s/.

3. To process ASTREGDEMv2.0 DEMs, run run_DEM_ASTGDEMv20.py or run_DEM_ASTGDEMv20_ipy.ipynb.

4. The log of running run_DEM_ASTGDEMv20.py can be seen in Log_run_DEM_ASTGDEMv20.txt.

Usage of EUDEM

1. Download DEM eu_dem_v11_E30N30.tif from here.

2. Copy DEM to the folder DATA/DATA_EUDEMv11/ and rename it as EUDEMv11_EPSG3035.tif.

3. To process EUDEMv1.1 DEM, run run_DEM_EUDEMv11.py or run_DEM_EUDEMv11_ipy.ipynb.

4. The log of running run_DEM_EUDEMv11.py can be seen in Log_run_DEM_EUDEMv11.txt.

Results

DEM Images of London

• The DEM Image of ASTERGDEMv2.0 in Pseudo Mercator PCS

• The DEM Image of ASTERGDEMv2.0 in BNG PCS

• The DEM Image of ASTERGDEMv2.0 in LAEA PCS

Elevation of London

• The elevation of 24 locations in London (meter)

References

[1] EPSG 4326 vs EPSG 3857 (projections, datums, coordinate systems, and more!)

[2] Coordinate systems and projections for beginners

[3] CSDN Blog

Please report an issue if you have any question about this repository, I will respond ASAP.

Please star this repository if you found its content useful. Thank you very much. ^_^