if not os.path.isfile(elevfile): print(('\nerror: file {} does not exist!'.format(elevfile))) print('double-check that the path is correct\n') raise # create an instance of the delineator and supply the path to the source files delineator = NHDPlusDelineator(VAAfile, flowfile, cfile, elevfile) # longitude, latitude of the point to delineate (the delineator looks for the # closest flowline to this point) longitude = -76.6056 latitude = 38.5839 # location to place the output (put it inside the existing HUC8 directory) gageoutput = '{}/01594670'.format(output) # file name plot of the output plot = '{}/hunting_watershed'.format(gageoutput) # extracts the catchments and flowlines for the gage's watershed and merge # the shapes together to make a boundary file delineator.delineate_watershed(longitude, latitude, output = gageoutput, plot = plot)
# paths to the different source files for the model data flowfile = '{}/flowlines'.format(output) # HUC8 flowline shapefile cfile = '{}/catchments'.format(output) # HUC8 catchment shapefile VAAfile = '{}/flowlineVAAs'.format(output) # NHDPlus value added attributes elevfile = '{}/elevations.tif'.format(output) # NED raster file watershed = '{}/delineated'.format(output) # directory for delineated files # create an instance of the delineator and supply the path to the source files delineator = NHDPlusDelineator(VAAfile, flowfile, cfile, elevfile) # longitude, latitude of the point to delineate (the delineator looks for the # closest flowline to this point) longitude = -76.6056 latitude = 38.5839 # extracts the catchments and flowlines for the gage's watershed and merge # the shapes together to make a boundary file gagewatershed = '{}/01594670'.format(output) delineator.delineate_watershed(longitude, latitude, output = gagewatershed) # make a plot of the watershed point = longitude, latitude plot = 'hunting_watershed' delineator.plot_delineated_watershed(point = point, output = plot)
if not os.path.isfile(elevfile): print('\nerror: file {} does not exist!'.format(elevfile)) print('double-check that the path is correct\n') raise # create an instance of the delineator and supply the path to the source files delineator = NHDPlusDelineator(VAAfile, flowfile, cfile, elevfile) # longitude, latitude of the point to delineate (the delineator looks for the # closest flowline to this point) longitude = -76.6056 latitude = 38.5839 # location to place the output (put it inside the existing HUC8 directory) gageoutput = '{}/01594670'.format(output) # file name plot of the output plot = '{}/hunting_watershed'.format(gageoutput) # extracts the catchments and flowlines for the gage's watershed and merge # the shapes together to make a boundary file delineator.delineate_watershed(longitude, latitude, output = gageoutput, plot = plot)