[1pt] PR: Clip WBD to 3DEP DEM domain

[mluck]

Creates a polygon of 3DEP DEM domain (to eliminate errors caused by stream networks with no DEM data in areas of HUCs that are outside of the U.S. border) and uses the polygon layer to clip the WBD and stream network (to a buffer inside the WBD). Addresses #748.

Additions

• data/usgs/acquire_and_preprocess_3dep_dems.py: Adds creation of 3DEP domain polygon by polygonizing all HUC6 3DEP DEMs and then dissolving them.
• src/gms/run_by_unit.sh: Adds 3DEP domain polygon .gpkg as input to src/clip_vectors_to_wbd.py

Changes

• src/clip_vectors_to_wbd.py: Clips WBD to 3DEP domain polygon and clips streams to a buffer inside the clipped WBD polygon.

Testing

1. Ran gms_run_unit.sh and gms_run_branch.sh on test HUCs 10060006, 04300109, 04100001 (specified in [8pt] US Border HUC branches are throwing error 61 #748) and 17010206.

Screenshots

Streams are clipped where the DEM data ends at the U.S. - Canada border

Streams are clipped to the buffer within the WBD

Notes

• This does not resolve the Exit status: 1 for HUC 04100001 (branch 1316000004, catchment in purple) which is caused by an unrelated issue

Todos

• Parallelize the HUC6 polygons creation: the current serialized procedure takes approximately 12 hours to create the 3DEP DEM domain polygon (but this only needs to be created if the DEM data change).

[RobHanna-NOAA]

New BED test, looked pretty good. One HUC failed outright (gms_run_unit) which is 010170104 and is a known issue.

95 branches failed of which 2 were exit status 1 and the other 93 were exit status 61. The two that were exit status of 1 are known issues. I spot checked some of the others they were for streams that are outside our stream buffer (which is 30 meters (3 * cellsize)) and that scenario is just fine. No further fix required.