## DESCRIPTION

r.topidx creates topographic index (wetness 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).

Input maps may have NULL values. For example, if you have a MASK for a watershed (basin map from r.water.outlet), the following command will create a masked elevation map (belev):

    r.mapcalc "belev = if(isnull(basin), basin, elev)"

r.stats -Anc prints out averaged statistics for topographic index.

## SEE ALSO

r.topmodel, r.water.outlet, r.mapcalc

## REFERENCE

Moore, I.D., R.B. Grayson, and A.R. Ladson, 1991. Digital terrain modeling: A review of hydrological, geomorphological, and biological applications. Hydrol. Processes 5:3-30.

## AUTHORS

Main algorithm sources are rewritten from GRIDATB.FOR.
Thanks to Keith Beven.

GRASS port by Huidae Cho
Hydro Laboratory, Kyungpook National University, South Korea

Last changed: $Date: 2011-11-08 03:29:50 -0800 (Tue, 08 Nov 2011) $