.TH r.surf.area 1grass "" "GRASS 7.6.0" "Grass User's Manual" .SH NAME \fI\fBr.surf.area\fR\fR \- Prints estimation of surface area for raster map. .SH KEYWORDS raster, surface, statistics, area estimation .SH SYNOPSIS \fBr.surf.area\fR .br \fBr.surf.area \-\-help\fR .br \fBr.surf.area\fR \fBmap\fR=\fIname\fR [\fBvscale\fR=\fIfloat\fR] [\fBunits\fR=\fIstring\fR] [\-\-\fBhelp\fR] [\-\-\fBverbose\fR] [\-\-\fBquiet\fR] [\-\-\fBui\fR] .SS Flags: .IP "\fB\-\-help\fR" 4m .br Print usage summary .IP "\fB\-\-verbose\fR" 4m .br Verbose module output .IP "\fB\-\-quiet\fR" 4m .br Quiet module output .IP "\fB\-\-ui\fR" 4m .br Force launching GUI dialog .SS Parameters: .IP "\fBmap\fR=\fIname\fR \fB[required]\fR" 4m .br Name of raster map .IP "\fBvscale\fR=\fIfloat\fR" 4m .br Vertical scale .br Default: \fI1.0\fR .IP "\fBunits\fR=\fIstring\fR" 4m .br Output units .br Default: square map units .br Options: \fImiles, feet, meters, kilometers, acres, hectares\fR .SH DESCRIPTION \fIr.surf.area\fR calculates area of regular 3D triangulated points (centers of cells) in current region by adding areas of triangles. Therefore, area of a flat surface will be reported as (rows + cols \-1) * (area of cell) less than area of flat region due to a half row and half column missing around the perimeter. .SH NOTES This calculation is heavily dependent on data resolution (think of it as a fractal shoreline problem, the more resolution the more detail, the more area, etc). This module uses the \fBcurrent region settings\fR, not the resolution of the raster map. This is especially important for surfaces with NULL values and highly irregular edges. The module does not [currently] attempt to correct for the error introduced by this \fIedge effect\fR. .PP This version actually calculates area twice for each triangle pair, keeping a running minimum and maximum area depending on the direction of the diagonal used. .PP Reported totals are: .IP .IP \fB1\fR \(dqPlan\(dq area of NULL values within the current GRASS region .IP \fB2\fR \(dqPlan\(dq area within calculation region (rows\-1 * cols\-1 * cellarea) .IP \fB3\fR Average of the minimum and maximum calculated 3d triangle area within this region .IP \fB4\fR \(dqPlan\(dq area within current computational region (rows * cols * cellarea) .IP \fB5\fR Scaling of calculated area to current region .PP .PP \fIr.surf.area\fR works best when the surface being evaluated extends to the edges of the current region and the cell resolution is small. Surfaces which are especially long and thin and have highly irregular boundaries will tend to have underestimated surface areas. Setting a high cell resolution (small area) will greatly reduce this impact, but will cause longer processing times. .SH EXAMPLES .br .nf \fC g.region \-p raster=elevation r.surf.area map=elevation units=hectares Null value area ignored in calculation: 0.000000 Plan area used in calculation: 20221.510000 Surface area calculation(low, high, avg): 20294.310421 20320.936368 20307.623395 Current region plan area: 20250.000000 Estimated region Surface Area: 20336.234719 \fR .fi .SH SEE ALSO \fI g.region, r.surf.idw, r.surf.fractal, r.surf.gauss, r.volume, r.slope.aspect, v.to.rast \fR .SH AUTHOR Bill Brown, USACERL December 21, 1994 .br Modified for floating point rasters and NULL values by Eric G. Miller (October 17, 2000) .br Updated for GRASS 7, and units option by Martin Landa, Czech Technical University in Prague, Czech Republic (October 2011) .PP \fILast changed: $Date: 2016\-09\-19 11:37:30 +0200 (Mon, 19 Sep 2016) $\fR .SH SOURCE CODE .PP Available at: r.surf.area source code (history) .PP Main index | Raster index | Topics index | Keywords index | Graphical index | Full index .PP © 2003\-2019 GRASS Development Team, GRASS GIS 7.6.0 Reference Manual