Retrieving watershed boundaries from the National Hydrography Dataset (NHD)
Uncomment the following line to install leafmap if needed.
# %pip install -U leafmap
import leafmap
point_geometry = {"x": -114.452985, "y": 36.13323}
gdf = leafmap.get_wbd(point_geometry, digit=8, return_geometry=True)
gdf.explore()
bbox_geometry = {"xmin": -115.0954, "ymin": 36.1000, "xmax": -114.6658, "ymax": 36.1986}
gdf = leafmap.get_wbd(bbox_geometry, digit=8, return_geometry=True)
gdf.explore()
huc_id = "10160002"
gdf = leafmap.get_wbd(searchText=huc_id, digit=8)
gdf.explore()
gdf = leafmap.get_wbd(searchText=huc_id, digit=10)
gdf.explore()