Data is available in an archive on Hydroshare. ➡️Hydroshare
Data in this folder will support an example of producing an FLDPLN library.
- Processor:
Intel(R) Xeon(R) CPU E5-2630 (2 processors) - Installed RAM:
128 GB DDR-3 1333MHz (10 of 16) - System type:
64-bit operating system, x64-based processor - Edition:
Windows Server 2022 Standard - Version:
21H2 - OS Build:
20348.1787
- Processors:
12th Gen Intel(R) Core(TM) i7-12700H 2.3 GHz - Installed RAM:
32 GB - System type:
64-bit operating system, x64-based processor - Edition:
Windows 10 Education - Version:
21H2 - OS build:
19044.3086 - Experience:
Windows Feature Experience Pack
- Version:
R2022b Update 5 - Type:
64-bit (win64) - Parallel Tools:
Parallel Computation Tools from Matlab ADD-Ons
- Version:
ArcGIS Pro 3.1.2
- Version:
ArcGIS Pro 3.1 ➡️HydroTool
- 10-meter Resolution ➡️USGS National Map
- Request to Dr. Jude Kastens at the University of Kansas
- Request to Dr. Jude Kastens at the University of Kansas
- The Amite River Basin in Louisiana (LA.) is an example of preprocessing.
- To initiate the study, download the HUC-8 Watershed Boundary Data (WBD) and the DEM of the study area.
- The HUC-8 refers to a shapefile that outlines the watershed boundary of the United States (U.S.).
- The DEM refers to the elevation data of the United States, sourced from the United States Geological Survey (USGS).
- To obtain DEM data from the USGS, users can upload their own shapefile to the website.
- However, keep in mind that there may be unexpected errors, as mentioned on the website.
To distinguish the designated study area, the HUC-8, which is the Continental United States (CONUS) WBD as the type of polygon, is necessary.
[Step 1] Open the HUC-8 shapefile in ArcGIS Pro.
[Step 2] Right-click the shapefile and go to the attribute table.
[Step 3] Delete data in the attribute table that does not relate to the target study area.
[Step 4] Right-click the shapefile, click the Data, and choose Export Features to save.
Save Ex. Amite_basin.shp
[Result] You will only be provided with the shapefile for the Amite River watershed boundary.
If the study area covers more than one piece of DEM, it is necessary to combine them through processing.
When conducting a large-scale study of FIM, it is common to combine DEMs.
However, this process is not necessary if the study area is already located on a single DEM.
[Step 1] In the Data Management Tool in ArcGIS Pro, go to the Mosaic to New Raster.
[Step 2] Type the information into Input Rasters, Output Location, Extension, and Number of Bands.
The other options are not considered.
The Amite River Basin covers 4 pieces of raw DEM, thus the 4 DEMs should be placed into the Input Rasters.
Save Ex. Combined_DEM.tif
[Result] Once the 4 pieces of the downloaded DEM are combined, you will have a single, unified DEM.
By applying the created shapefile to mask the DEM, users can obtain a masked DEM for the study area only.
The square box DEM, which is a raw elevation model from the USGS, contains data on elevations that include unnecessary locations.
The following steps are showing how to mask the study area DEM by cutting the unnecessary part of the DEM off.
[Step 1] Open the extracted shapefile (Amite_basin.shp) from [1] and the DEM (Combined_DEM.tif) from [2].
[Step 2] At the Spatial Analysis Tools, go to the Extract by Mask.
[Step 3] The DEM and extracted shapefile should be placed into the input raster and feature mask, respectively.
[Step 4] Make sure to include the file extension as TIF when naming the output file.
Save Ex. Masked_Amite_DEM.tif
[Result] It appears that the masked DEM and extracted shapefile are matched.
[Step 1] Go to the Terrain Preprocessing in ArcHydro Tool.
[Step 2] In the DEM Manipulation, click the Fill Sinks.
[Step 3] In the Fill Sinks, Type the name and save it in TIF format.
Save Ex. Fil.tif
[Result] You will see the Filled DEM map in the main GIS window.
[Step 1] Go to the Terrain Preprocessing in ArcHydro Tool.
[Step 2] Click the Flow Direction.
[Step 3] Input is the output from Fill DEM (Fil.tif).
[Step 4] Type the name and save it in TIF format.
Save Ex. Fdr.tif
[Result] You will see the Flow Direction map in the main GIS window.
[Step 1] Go to the Terrain Preprocessing in ArcHydro Tool.
[Step 2] Click the Flow Accumulation.
[Step 3] Input is the output from Flow Direction (Fdr.tif).
[Step 4] Type the name and save it in TIF format.
Save Ex. Fac.tif
[Result] You will see the Flow Accumulation map in the main GIS window.
[Step 1] Go to the ArcGIS pro, Catalog pane.
[Step 2] In the Catalog pane, click the Computer tab and go to the file location you saved.
Ex. Masked_Amite_DEM.tif, Fil.tif, Fdr.tif, Fac.tif
[Step 3] Right-click on each TIF file and click the export to a different format.
[Step 4] Convert all TIF files to ESRI BIL format.
[Step 5] During the conversion of files, other options are not considered.
Only for Fac.tif converting, type the -9999 in NoData Value.
After completing this preprocessing, you should have 8 items listed below.
- Masked_Amite_DEM.tif and Masked_Amite_DEM.bil
- Fil.tif and Fil.bil
- Fdr.tif and Fdr.bil
- Fac.tif and Fac.bil
-
Home directory: 15 - 20 line
> dr1 = 'Tuscaloosa\'; > dr0 = ['C:\Users\jsong33\Desktop\Work\FIM\FLDPLN\StudyAreas\Tuscaloosa\',dr1]; > if(~exist(dr0,'dir')) > dr0 = ['C:\Users\jsong33\Desktop\Work\FIM\FLDPLN\StudyAreas\Tuscaloosa\',dr1]; > end -
Iteration step size (in vertical DEM units): 27 - 32 line
> %dh = 0.01; > %dh = 0.1; > %dh = 0.25; > %dh = 0.5; > %dh = 1; -
Identify synthetic stream network: 54 - 58 line
> mi2px = 5280*.3048/pxsz; % horizontal units assumed to be meters (otherwise modify as needed) > seglen = round(5*mi2px); % 5 miles maximum segment length > sqmi2px = (5280^2)*(.3048^2)/(pxsz^2); % horizontal units assumed to be meters (otherwise modify as needed) > facthr = round(25*sqmi2px); % 25 sq. mi minimum catchment size (min for segment break) > strthr = round(70*sqmi2px); % 70 sq. mi minimum catchment size (min for defining a stream) -
Seedpts: 69 - 70 line
> create_segdb_facthr_maxlen(strthr,facthr,seglen,fdrf,facf,segf,matf) % With no seedpts > create_segdb_facthr_maxlen_seedpts(shpf,facthr,seglen,fdrf,facf,segf,matf) % With seedpts -
Parallelization: 187 - 212 line
* Worker Pool > if(1>0) > g = gcp('nocreate'); > if(~isempty(g)) > p = gcp; > else > p = parpool; > end > for j = 1:num > pf(j) = parfeval(p,@fldpln_model_v5ram,1,seg_list(j),inp,fildat,fdrdat,fldmn,fldmx,dh,mxht,1,bg); > end > % wait(pf) > pfresults = cell(num,1); > for idx = 1:num > % fetchNext blocks until next results are available. > [completedIdx,value] = fetchNext(pf); > pfresults{completedIdx} = value; > fprintf('Segment %d completed\n', completedIdx); > end > % delete(gcp('nocreate')) > end > return; * Simple One-segment-at-a-time 1. Active lines below 2. Change 1>0 to 1<0 in the IF statement from Worker Pool > % parfor j = 1:num > % fldpln_model_v5ram(seg_list(j),inp,fildat,fdrdat,fldmn,fldmx,dh,mxht,1,bg); > % end > % return;