This GRASS GIS addon is the iterative version of r.accumulate -r. It is already integrated into r.accumulate as the default algorithm (-r flag for the recursive version) and deprecated r.lfp.
For some reason, Mendeley's data upload did not work successfully and my dataset page lists no data files at all. Sorry for the inconvenience. Please use the following files for your review:
create_test_location.shgrassdata.zip(5.3 GB)analysis.zip(397 MB)
System requirements: GRASS GIS 7.8.3 on Linux; this addon was not tested on Windows
- Download
create_test_location.sh,grassdata.zip, andanalysis.zipin the same directory - Unzip
grassdata.zipandanalysis.zip - Run
grass78grass78 -c test --exec ./create_test_location.sh
Five mapsets and 1,616 maps will be created in the test location:
PERMANENT: Not usedlfp-GA-input: Input rasters and vectors for Georgia (305 maps)lfp-GA-hdd: Results for Georgia from the Linux HDD system (202 maps); you will have access to the maps inlfp-GA-inputlfp-GA-ssd: Results for Georgia from the Linux SSD system (202 maps); you will have access to the maps inlfp-GA-inputlfp-GA-archydropro: Results for Georgia from the Windows system (101 maps); you will have access to the maps inlfp-GA-inputlfp-TX-input: Input rasters and vectors for Texas (305 maps)lfp-TX-hdd: Results for Texas from the Linux HDD system (200 maps); you will have access to the maps inlfp-TX-inputlfp-TX-ssd: Results for Texas from the Linux SSD system (200 maps); you will have access to the maps inlfp-TX-inputlfp-TX-archydropro: Results for Texas from the Windows system (101 maps); you will have access to the maps inlfp-TX-input
Both Georgia and Texas analyses used the same map names except for the following close-to-raw input maps:
lfp-GA-inputmapsetned_1_ga: USGS National Elevation Dataset (NED) raster for Georgia (not clipped to the boundary)statep010_ga: Mask raster for Georgia (clipped to the boundary)
lfp-TX-inputmapsetned_1_tx: USGS NED raster for Texas (not clipped to the boundary)statep010_tx: Mask raster for Texas (clipped to the boundary)
Inputs in the lfp-*-input mapset
drain: Flow direction raster created byr.watershedaccum: Flow accumulation raster created byr.accumulate, but it is only an intermediate output and not neededoutlets100: 100 random outlet points vector for the batch experimentoutlets100_1throughoutlets100_100: The same 100 outlets, but each point individually in one vector map for the non-batch experimentbasins100_subbasin_1throughbasins100_subbasin_100: Each subwatershed in an individual raster map for maskingbasins100_archydropro_1throughbasins100_archydropro_100: Each subwatershed in an individual vector map for setting the computational region; the above 100 subbasin rasters have the extent of the entire state
Outputs in both lfp-*-hdd and lfp-*-ssd mapsets
lfp100_recursive: 100 longest flow paths generated byr.accumulate -rfrom the batch mode experiment- Missing in the Texas mapsets because
r.accumulate -rfailed
- Missing in the Texas mapsets because
lfp100_recursive_1throughlfp100_recursive_100: 100 individual longest flow paths generated byr.accumulate -rfrom the non-batch mode experimentlfp100_recursive_51is missing in the Texas mapsets becauser.accumulate -rfailed
lfp100_iterative: 100 longest flow paths generated by thisr.accumulate-iterativemodulelfp100_iterative_1throughlfp100_iterative_100: 100 individual longest flow paths generated by thisr.accumulate-iterativemodule
Outputs in the lfp-*-archydropro mapset; imported from ArcGIS Pro geodatabases
lfp100_archydropro: 100 longest flow paths generated by the Longest Flow Path tool in the Arc Hydro toolbox for ArcGIS Prolfp100_archydropro_1throughlfp100_archydropro_100: 100 individual longest flow paths generated by the Longest Flow Path tool