NAME¶

r.stats - Generates area statistics for raster map layers.

KEYWORDS¶

raster, statistics

SYNOPSIS¶

r.stats r.stats help r.stats [-1AacplgxrnNCi] input=name[,name,...] [ output=name] [fs= character|space|tab] [nv=string] [nsteps= integer] [--verbose] [--quiet]

Flags:¶

-1

 

One cell (range) per line

-A

 

Print averaged values instead of intervals

-a

 

Print area totals

-c

 

Print cell counts

-p

 

Print APPROXIMATE percents (total percent may not be 100%)

-l

 

Print category labels

-g

 

Print grid coordinates (east and north)

-x

 

Print x and y (column and row)

-r

 

Print raw indexes of fp ranges (fp maps only)

-n

 

Suppress reporting of any NULLs

-N

 

Suppress reporting of NULLs when all values are NULL

-C

 

Report for cats fp ranges (fp maps only)

-i

 

Read fp map as integer (use map's quant rules)

--verbose

 

Verbose module output

--quiet

 

Quiet module output

Parameters:¶

input=name[,name,...]

 

Name of input raster map(s)

output=name

 

Name for output file (if omitted or "-" output to stdout)

fs=character|space|tab

 

Field separator

 

Output field separator

 

Default: space

nv=string

 

String representing no data cell value

 

Default: *

nsteps=integer

 

Number of fp subranges to collect stats from

 

Default: 255

DESCRIPTION¶

r.stats calculates the area present in each of the categories of user-selected raster map layer(s). Area statistics are given in units of square meters and/or cell counts. This analysis uses the current geographic region and mask settings. Output can be sent to a file in the user's current working directory. If a single map layer is specified on the command line, a list of areas in square meters (assuming the map's coordinate system is in meters) for each category in the raster map layer will be printed. (If the -c option is chosen, areas will be stated in number of cells.) If multiple raster map layers are specified on the command line, a cross-tabulation table of areas for each combination of categories in the map layers will be printed. For example, if one raster map layer were specified, the output would look like:
1:1350000.00
2:4940000.00
3:8870000.00 If three raster map layers a, b, and c, were specified, the output would look like:
0:0:0:8027500.00
0:1:0:1152500.00
1:0:0:164227500.00
1:0:1:2177500.00
1:1:0:140092500.00
1:1:1:3355000.00
2:0:0:31277500.00
2:0:1:2490000.00
2:1:0:24207500.00
2:1:1:1752500.00
3:0:0:17140000.00
3:1:0:11270000.00
3:1:1:2500.00 Within each grouping, the first field represents the category value of map layer a, the second represents the category values associated with map layer b, the third represents category values for map layer c, and the last field gives the area in square meters for the particular combination of these three map layers' categories. For example, above, combination 3,1,1 covered 2500 square meters. Fields are separated by colons.

NOTES¶

r.stats works in the current geographic region with the current mask. If a nicely formatted output is desired, pipe the output into a command which can create columnar output. For example, the command: r.stats input=a,b,c | pr -3 | cat -s will create a three-column output 1:4:4:10000.00 2:1:5:290000.00 2:4:5:2090000.00 1:4:5:1340000.00 2:2:5:350000.00 3:1:2:450000.00 2:1:1:1090000.00 2:4:1:700000.00 3:1:3:5280000.00 2:1:3:410000.00 2:4:3:10000.00 3:1:5:3140000.00 The output from r.stats on more than one map layer is sorted. Note that the user has only the option of printing out cell statistics in terms of cell counts and/or area totals. Users wishing to use different units than are available here should use the GRASS program r.report.

SEE ALSO¶

g.region, r.coin, r.describe, r.report, r.statistics, r.univar

AUTHOR¶

Michael Shapiro, U.S. Army Construction Engineering Research Laboratory Last changed: $Date: 2008-05-16 21:09:06 +0200 (Fri, 16 May 2008) $ Full index © 2003-2014 GRASS Development Team