## Crop watershed-delineated DEMs from these retiles, using features labelled by the `index` attribute in a separate vector image
mkdir HUCs.d
NUMBER_OF_HUCS=739
for huc in $(seq 0 $(( ${NUMBER_OF_HUCS} - 1)) ); do gdalwarp -multi -cutline HUCs.shp -cl HUCs -cwhere "index=${huc}" -crop_to_cutline albers-warped.d.vrt HUCs.d/HUC${huc}.tif ; done
Consider re-including in the
README.md