Terrain and hydrological analysis based on LiDAR-derived digital elevation models (DEM)
Project description
lidar
Author: Qiusheng Wu (https://wetlands.io | wqs@binghamton.edu)
lidar is a toolset for terrain and hydrological analysis using digital elevation models (DEMs). It is particularly useful for analyzing high-resolution topographic data, such as DEMs derived from Light Detection and Ranging (LiDAR) data.
GitHub repo: https://github.com/giswqs/lidar
Documentation: https://lidar.readthedocs.io.
Free software: MIT license
Features
Extracting depressions from DEMs (see filling.py).
Filtering out small artifact depressions based on user-specified minimum depression size (see filling.py).
Generating refined DEMs with small depressions filled but larger depressions kept intact (see filling.py).
Delineating depression nested hierarchy using the level-set method (see slicing.py).
Computing topological and geometric properties of depressions, including size, volume, mean depth, maximum depth, lowest elevation, and spill elevation (see slicing.py).
Exporting depression properties as a csv file (see slicing.py).
Using It
Install the Python package using the following command:
pip install lidarAnd use:
import os
import pkg_resources
import lidar
# identify the sample data directory of the package
package_name = 'lidar'
data_dir = pkg_resources.resource_filename(package_name, 'data/')
# use the sample dem. Change it to your own dem if needed
in_dem = os.path.join(data_dir, 'dem.tif')
# set output directory. By default, use the temp directory under user's home directory
out_dir = os.path.join(os.path.expanduser("~"), "temp")
# parameters for identifying sinks and delineating nested depressions
min_size = 1000 # minimum number of pixels as a depression
min_depth = 0.5 # minimum depth as a depression
interval = 0.3 # slicing interval for the level-set method
bool_shp = False # output shapefiles for each individual level
# extracting sinks based on user-defined minimum depression size
sink_path = lidar.ExtractSinks(in_dem, min_size, out_dir)
dep_id_path, dep_level_path = lidar.DelineateDepressions(sink_path, min_size, min_depth, interval, out_dir, bool_shp)Check the example.py for more details.
Examples
The images below show working examples of the level set method for delineating nested depressions in the Cottonwood Lake Study Area (CLSA), North Dakota. More test datasets (e.g., the Pipestem watershed in the Prairie Pothole Region of North Dakota) can be downloaded from http://gishub.org/2018-JAWRA-Data
The following example was conducted on a 64-bit Linux machine with a quad-core Intel i7-7700 CPU and 16 GB RAM. The average running time of the algorithm for this DEM was 0.75 seconds.
Publications
The level-set algorithm used in the lidar package has been published in the following articles.
Wu, Q., Lane, C.R., Wang, L., Vanderhoof, M.K., Christensen, J.R., & Liu, H. (2018). Efficient Delineation of Nested Depression Hierarchy in Digital Elevation Models for Hydrological Analysis Using Level-Set Methods. Journal of the American Water Resources Association. (forthcoming)
Credits
This algorithms are built on richdem, numpy, scipy, scikit-image, and pygdal.
This package was created with Cookiecutter and the audreyr/cookiecutter-pypackage project template.
History
0.1.5 (2018-05-15)
0.1.4 (2018-05-15)
0.1.3 (2018-05-15)
0.1.2 (2018-05-15)
0.1.1 (2018-05-14)
0.1.0 (2018-05-14)
First release on PyPI.
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
File details
Details for the file lidar-0.1.5.tar.gz.
File metadata
- Download URL: lidar-0.1.5.tar.gz
- Upload date:
- Size: 675.1 kB
- Tags: Source
- Uploaded using Trusted Publishing? No
File hashes
| Algorithm | Hash digest | |
|---|---|---|
| SHA256 | 866b2c1d1a066b07e4b92f4d167225d0bef4be9df262773ff02f84e7b56d094f |
|
| MD5 | 378c38f9079bfca3b316c6c75f157a86 |
|
| BLAKE2b-256 | 2c757351879e75aefe25ee948cd15e4fd7c0ff3036441665f61c7c8717eac232 |
