For example, if one raster map layer were specified, the output would look like:
1:1350000.00 2:4940000.00 3:8870000.00If 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.00Within 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.
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.00The 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.
Last changed: $Date$