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).
[dMask, dDepths] = findDepressions(heights)
Returns both the logical mask of where depressions exist and the depths of those depressions.
See chapter 4 of my dissertation.
I wrote the function for this paper specifically.
I try to understand this algorithm. Can you give me some papers from which you take inspiration to write this algorithm ?
Inspired by: extrema.m, extrema2.m