~ubuntu-branches/ubuntu/vivid/grass/vivid-proposed

« back to all changes in this revision

Viewing changes to raster/r.surf.area/r.surf.area.html

  • Committer: Package Import Robot
  • Author(s): Bas Couwenberg
  • Date: 2015-02-20 23:12:08 UTC
  • mfrom: (8.2.6 experimental)
  • Revision ID: package-import@ubuntu.com-20150220231208-1u6qvqm84v430b10
Tags: 7.0.0-1~exp1
* New upstream release.
* Update python-ctypes-ternary.patch to use if/else instead of and/or.
* Drop check4dev patch, rely on upstream check.
* Add build dependency on libpq-dev to grass-dev for libpq-fe.h.
* Drop patches applied upstream, refresh remaining patches.
* Update symlinks for images switched from jpg to png.

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
<h2>DESCRIPTION</h2>
 
2
 
 
3
<em>r.surf.area</em> calculates area of regular 3D triangulated points
 
4
(centers of cells) in current region by adding areas of triangles.
 
5
Therefore, area of a flat surface will be reported as (<tt>rows + cols
 
6
-1) * (area of cell)</tt> less than area of flat region due to a half
 
7
row and half column missing around the perimeter.
 
8
 
 
9
<h2>NOTES</h2>
 
10
 
 
11
This calculation is heavily dependent on data resolution (think of it
 
12
as a fractal shoreline problem, the more resolution the more detail,
 
13
the more area, etc). This module uses the <b>current region
 
14
settings</b>, not the resolution of the raster map. This is especially
 
15
important for surfaces with <tt>NULL</tt> values and highly irregular
 
16
edges. The module does not [currently] attempt to correct for the
 
17
error introduced by this <em>edge effect</em>.
 
18
 
 
19
<p>This version actually calculates area twice for each triangle pair,
 
20
keeping a running minimum and maximum area depending on the direction
 
21
of the diagonal used.
 
22
 
 
23
<p>Reported totals are:
 
24
<ol>
 
25
<li>"Plan" area of <tt>NULL</tt> values within the current GRASS
 
26
region</li>
 
27
<li>"Plan" area within calculation region (<tt>rows-1 * cols-1 *
 
28
cellarea</tt>)</li>
 
29
<li>Average of the minimum and maximum calculated 3d triangle area
 
30
within this region</li>
 
31
<li>"Plan" area within current computational region (<tt>rows * cols *
 
32
cellarea</tt>)</li>
 
33
<li>Scaling of calculated area to current region</li>
 
34
</ol> 
 
35
 
 
36
<p><em>r.surf.area</em> works best when the surface being evaluated
 
37
extends to the edges of the current region and the cell resolution is
 
38
small. Surfaces which are especially long and thin and have highly
 
39
irregular boudaries will tend to have underestimated surface areas.
 
40
Setting a high cell resolution (small area) will greatly reduce this
 
41
impact, but will cause longer processing times.
 
42
 
 
43
<h2>EXAMPLES</h2>
 
44
 
 
45
<div class="code"><pre>
 
46
g.region -p raster=elevation
 
47
 
 
48
r.surf.area map=elevation units=hectares
 
49
Null value area ignored in calculation: 0.000000
 
50
Plan area used in calculation: 20221.510000
 
51
Surface area calculation(low, high, avg):
 
52
        20294.310421 20320.936368 20307.623395
 
53
Current region plan area: 20250.000000
 
54
Estimated region Surface Area: 20336.234719
 
55
</pre></div>
 
56
 
 
57
<h2>SEE ALSO</h2>
 
58
 
 
59
<em>
 
60
<a href="g.region.html">g.region</a>,
 
61
<a href="r.surf.idw.html">r.surf.idw</a>,
 
62
<a href="r.surf.fractal.html">r.surf.fractal</a>,
 
63
<a href="r.surf.gauss.html">r.surf.gauss</a>,
 
64
<a href="r.volume.html">r.volume</a>,
 
65
<a href="r.slope.aspect.html">r.slope.aspect</a>,
 
66
<a href="v.to.rast.html">v.to.rast</a>
 
67
</em>
 
68
 
 
69
<h2>AUTHOR</h2>
 
70
Bill Brown, USACERL December 21, 1994
 
71
<br>
 
72
Modified for floating point rasters and <tt>NULL</tt> values by Eric
 
73
G. Miller (October 17, 2000)
 
74
<br>
 
75
Updated for GRASS 7, and units option by Martin Landa, Czech Technical
 
76
University in Prague, Czech Republic (October 2011)
 
77
 
 
78
<p><i>Last changed: $Date: 2014-12-27 23:51:00 +0100 (Sat, 27 Dec 2014) $</i>