GRASS logo

NAME

r.statistics - Category or object oriented statistics.

KEYWORDS

raster

SYNOPSIS

r.statistics
r.statistics help
r.statistics [-c] base=string cover=string method=string [output=string] [--overwrite]

Flags:

-c
Cover values extracted from the category labels of the cover map
--overwrite
Force overwrite of output files

Parameters:

base=string
Base raster map
cover=string
Cover raster map
method=string
Method of object-based statistic
Options: distribution,average,mode,median,avedev,stddev,variance,skewness,kurtosis,min,max,sum
output=string
Resultant raster map (not used with 'distribution')

DESCRIPTION

r.statistics is a tool, where different class, resp. object oriented statistical analysis methods are possible.
The methods are the following (english - german): The calculations will be performed for areas with data of the covering-layers which belong in the base-layer o the same category.

Setting the -c flag the category lables of the covering raster layer will be used. This is nice to avoid the GRASS limitation to interger in raster maps because using category values floating point numbers can be stored.

The output-layer is a reclassified version of the base-layer with identical category values. The results of the calculations are stored in the category labels of the output-layer.

The described output is valid for all calculations with exception of distribution. Here the output is given to stdout. A file name eventually specified will be ignored. The result will be a table with three columns. In the 1. column are the category values of the base-layer (a), in the 2. column the associated value of the cover-Layers (b) and in the 3. column the percentage of b from a. Example:

1   124  23
1   201  47
1   273  30
2   101  5
2   152  16
2   167  60
2   187  19
     .
     .
     .
To transfer the values stored as category labels into cell values, r.mapcalc can be used ('@' operator).

EXAMPLE

Calculation of average elevation of each field in the Spearfish region:
r.statistics base=fields cover=elevation.dem out=elevstats method=average
r.cats elevstats
r.mapcalc "fieldelev=@elevstats"
r.univar fieldelev

REMARK

The source code of average, mode and median was taken from some GRASS raster modules.

BUGS

unkown

SEE ALSO

r.mode, r.median, r.average, r.mapcalc, r.neighbors

AUTHOR

Martin Schroeder, Geographisches Institut Heidelberg, Germany

Last changed: $Date: 2006/08/08 13:09:15 $


Main index - raster index - Full index