Replies: 1 comment 1 reply
-
Did you try out the example data? import os
import pkg_resources
from lidar import *
# identify the sample data directory of the package
package_name = 'lidar'
data_dir = pkg_resources.resource_filename(package_name, 'data/')
# use the sample dem. Change it to your own dem if needed
in_dem = os.path.join(data_dir, 'dem.tif')
# set the output directory
out_dir = os.getcwd()
# parameters for identifying sinks and delineating nested depressions
min_size = 1000 # minimum number of pixels as a depression
min_depth = 0.5 # minimum depth as a depression
interval = 0.3 # slicing interval for the level-set method
bool_shp = True # output shapefiles for each individual level
# extracting sinks based on user-defined minimum depression size
out_dem = os.path.join(out_dir, "median.tif")
in_dem = MedianFilter(in_dem, kernel_size=3, out_file=out_dem)
sink_path = ExtractSinks(in_dem, min_size, out_dir)
dep_id_path, dep_level_path = DelineateDepressions(sink_path,
min_size,
min_depth,
interval,
out_dir,
bool_shp)
print('Results are saved in: {}'.format(out_dir)) |
Beta Was this translation helpful? Give feedback.
1 reply
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
Hi, I was able to install the ArcGIS Pro Lidar package to extract sinks and depression hierarchies. However, I do not have an output depressions_info.csv file. Is this no longer an output for the latest lidar package or am I overlooking it somehow? I am utilizing bathymetry with negative elevations, so maybe that's a problem for the depressions_info.csv output file? I'm not sure if the previous reported issue may have caused a halt for this .csv output, but the information within the depressions_info.csv file (Table 1 in the supporting lidar documentation) is a central theme to a chapter of my PhD dissertation and I hope to publish and cite your work within the next year. I'm very thankful for your contributions so far!
Beta Was this translation helpful? Give feedback.
All reactions