Abstract:It has been universally recognized by the previous researches that depressions in DEM data are the major obstacle of determining hydrologic flow directions, so numerous methods were proposed to address this problem. However, when massive DEM data, such as SRTM and GTOPO30, become available through internet, the conventional approaches are too time- consuming and inefficient in dealing with such massive data for hydrologic analyses. A new innovative approach, unlike all existing methods that process DEM data straightforward without utilizing the topographic features implied in them, was proposed in this study to improve the DEM-processing efficiency. The new approach has two steps: first, classify the initial DEM data into eight groups according to the elevation values using the quantile classification method, and set each category cell values to its quantile and store them in a transient matrix. Second, scan the transient matrix from the minimum category to the maximum one, restore its initial value if its initial value is larger than or equal to its neighbors', else set its value to the minimum of its neighbors' if its quantile is larger than its neighbor's, and repeat this process until all the depressions are filled. Furthermore, unlike the traditional method, which have to scan all the DEM cells in each loop, the new one only scan the data needed to be further processed by storing their location information into two stacks. As a result, the total cells and scanning times are dramatically decreased due to stack's nature of "first in last out", and the location information of depressions is all stored in one of the two stacks, which facilitate the subsequent hydrologic analysis. DEMs of various sizes were used for testing and evaluating the new method. The results show the new method is thousands of times faster than the famous Jenson and Domingue's (1988) method and is over seven times faster than Planchon and Darboux's on average The new algorithm is fully detailed and pseudo-codes are provided in this paper.

A Novel Parallel Depression Removing Algorithm for Hydrology Analysis in Digital Elevation Models

Parallel Algorithms for Extracting Drainage Network from Digital Elevation Models

Parallel Priority-Flood Depression Filling For Trillion Cell Digital Elevation Models On Desktops Or Clusters

Parallel identification and filling of depressions in raster digital elevation models.

Extraction of Drainage Network from Grid Terrain Datasets Using Parallel Computing

A Parallel Algorithm for Erosion Slope Length Calculation Based on Grid Terrain Datasets

A Novel Method for Filling the Depressions in Massive DEM Data

Parallel Contributing Area Calculation with Granularity Control on Massive Grid Terrain Datasets

An Efficient Algorithm for Depression Filling and Flat-Surface Processing in Raster DEMs

Efficient Delineation of Nested Depression Hierarchy in Digital Elevation Models for Hydrological Analysis Using Level‐Set Method

Priority-Flood: An Optimal Depression-Filling and Watershed-Labeling Algorithm for Digital Elevation Models

Parallel Calculation and Efficiency Analysis for Neighborhood Statistic Algorithm in Digital Terrain Analysis

Improved Priority‐Flood Method for Depression Filling by Redundant Calculation Optimization in Local Micro‐relief Areas

Retaining Relative Height Information: an Enhanced Technique for Depression Treatment in Digital Elevation Models

Parallel assignment of flow directions over flat surfaces in massive digital elevation models

An Efficient Variant of the Priority-Flood Algorithm for Filling Depressions in Raster Digital Elevation Models

Methods for Processing Depression and Flat Areas in Extracting Drainage Networks Based on the DEM

A New Method for Treating the Depressions and Flat Areas in DEM for Large-Scale Hydrology and Climate Models

An Integrated Algorithm for Depression Filling and Assignment of Drainage Directions over Flat Surfaces in Digital Elevation Models

Another Fast and Simple DEM Depression-Filling Algorithm Based on Priority Queue Structure

An efficient method for applying a differential equation to deriving the spatial distribution of specific catchment area from gridded digital elevation models