Note: This document is for an older version of GRASS GIS that will be discontinued soon. You should upgrade, and read the current manual page.
g.region raster=elev_lid792_1m -p r.slope.aspect elevation=elev_lid792_1m slope=slope aspect=aspect r.flow elevation=elev_lid792_1m flowaccumulation=flowacc # exponents m=1.3 and n=1.2 # multiply flowaccumulation by cell area/resolution to get contributing area per unit width r.mapcalc "sflowtopo = pow(flowacc * 1.,1.3) * pow(sin(slope),1.2)" # Compute sediment flow by combining the rainfall, soil and land cover factors # with the topographic sediment transport factor. # We use a constant value of 270. for rainfall intensity factor. r.mapcalc "sedflow = 270. * soils_Kfactor * cfactorbare_1m * sflowtopo" r.divergence magnitude=sedflow direction=aspect output=erosion_deposition # set suitable color table r.colors map=erosion_deposition rules=- << EOF 0% 100 0 100 #dark magenta -100 magenta -10 red -1 orange -0.1 yellow 0 200 255 200 #light green 0.1 cyan 1 aqua 10 blue 100 0 0 100 #dark blue 100% black EOF
Available at: r.divergence source code (history)
Latest change: Tuesday Nov 05 09:35:47 2024 in commit: f00067506253c89598cffadc3c7fa07cff7233fd
Main index | Raster index | Topics index | Keywords index | Graphical index | Full index
© 2003-2024 GRASS Development Team, GRASS GIS 8.3.3dev Reference Manual