- Abstract A new method for extracting drainage systems from Digital Elevation Models (DEMs) is presented. The main algorithm of the proposed method performs a skeletonization process of the set of elevations in the DEM and produces a skeleton of flow paths. An enumeration algorithm performs the removal of loops from the initial flow path. A preprocess for filling depressions is described as is the necessary postprocessing for determining the drainage network through depressions. The new method does not suffer from any of the maladies of former methods described in the literature, such as flow cutoffs, loops of flow, and basin flooding. The new method is tested on several real-world DEMs and produced connected, complete, and loopless networks.