SURFACE DEPRESSION MAPPING - LindaTheresB
SURFACE DEPRESSION MAPPING - LindaTheresB
SURFACE DEPRESSION MAPPING - LindaTheresB
net/publication/362413454
CITATIONS READS
0 100
1 author:
Linda Theres
SASTRA University
6 PUBLICATIONS 2 CITATIONS
SEE PROFILE
All content following this page was uploaded by Linda Theres on 02 August 2022.
▪ Next, we compute depressions, by subtracting focal statistics output from the fill output
in raster calculator
▪ Then we extract depression by using conditional statement in raster calculator. Say if
above step output name is ‘sink’, the conditional statement is Con (sink>0, 1). The
output is basically raster with depth values.
▪ Next two steps will group cells to form polygons. First using region group tool, give
abpve output as input and remaining values as default
▪ Then in the zonal statistics tool, input is the output of the above step and the original
DEM file as the input value raster.
▪ Then filter depressions by giving minimum depth conditions. We use a raster calculator.
If above output is ‘zonal_op’, the conditional statement is Con (zonal_op>2, 1). Note
here we have given 2 as minimum depth. This condition value can be chosen according
to the requirement.
▪ Then we convert the above output raster to a polygon using a conversion tool
▪ Then project the vector file to compute the area of the polygons derived. After
projection, create a field ‘area_sqkm’ in the attribute table, and using the ‘calculate
geometry’ option, compute the area of all the polygon features.
▪ Finally, we select the polygons by providing minimum area conditions and export data
as the final output and display as a surface depression map.
REFERENCES
https://github.com/giswqs/lidar
https://onlinelibrary.wiley.com/doi/10.1111/1752-1688.12689
https://joss.theoj.org/papers/10.21105/joss.02965