Luca Delucchi
021fa9d879
r.topidx: added examples and figures, thanks to Jerry Huang and Google Code-in
|
7 gadi atpakaļ | |
|---|---|---|
| .. | ||
| Makefile | 16 gadi atpakaļ | |
| README | 17 gadi atpakaļ | |
| arc.to.gridatb.pl | 13 gadi atpakaļ | |
| file_io.c | 11 gadi atpakaļ | |
| global.h | 11 gadi atpakaļ | |
| gridatb.to.arc.pl | 13 gadi atpakaļ | |
| main.c | 8 gadi atpakaļ | |
| r.topidx.html | 7 gadi atpakaļ | |
| r_topidx.png | 7 gadi atpakaļ | |
| topidx.c | 10 gadi atpakaļ | |
README
Main algorithm sources are rewritten from GRIDATB.FOR.
GRIDATB.FOR is available from http://www.es.lancs.ac.uk/hfdg/topmodel.html.
Thanks to Keith Beven.
[r.topidx]
creates topographic index, ln(a/tan(beta)), map from elevation map.
where
a: the area of the hillslope per unit contour length that drains
through any point
tan(beta): the local surface topographic slope
(delta vertical) / (delta horizontal)
Elevation map need not to be rectangular if study basin was already determined.
Non-interesting region may have NULL values. Study basin map can be taken
from r.water.outlet then the following r.mapcalc statements will create
study basin extracted elevation map.
r.mapcalc "belev = if(isnull(basin), basin, elev)"
GRIDATB.FOR is available from http://www.es.lancs.ac.uk/hfdg/topmodel.html.
Thanks to Keith Beven.
[r.topidx]
creates topographic index, ln(a/tan(beta)), map from elevation map.
where
a: the area of the hillslope per unit contour length that drains
through any point
tan(beta): the local surface topographic slope
(delta vertical) / (delta horizontal)
Elevation map need not to be rectangular if study basin was already determined.
Non-interesting region may have NULL values. Study basin map can be taken
from r.water.outlet then the following r.mapcalc statements will create
study basin extracted elevation map.
r.mapcalc "belev = if(isnull(basin), basin, elev)"