findDepressions

Finds the area coverage and depth of localized depressions on a surface.

This function finds depressions on a 2D surface. It was intended to model the formation of liquid pools on a simulated geological surface (lakes pooling in depressions).

Usage:
[dMask, dDepths] = findDepressions(heights)

Returns both the logical mask of where depressions exist and the depths of those depressions.