← Back to branch summary

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

« back to all changes in this revision

Viewing changes to raster/r.slope.aspect/description.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:
raster/r.slope.aspect/description.html
1
 
<h2>DESCRIPTION</h2>
2
 
 
3
 
<b>r.slope.aspect </b>generates raster maps of slope, aspect, curvatures and
4
 
first and second order partial derivatives from a raster map of true
5
 
elevation values. The user must specify the input <i>elevation</i> file name
6
 
and at least one output file name. The user can also specify the
7
 
<i>format</i> for slope (degrees, percent; default=degrees), and the 
8
 
<i>zfactor</i>: multiplicative factor to convert elevation units to meters;
9
 
(default 1.0).
10
 
<p>
11
 
The <i>elevation</i> input raster map specified by the user must contain true
12
 
elevation values, <b>not</b> rescaled or categorized data. If the elevation
13
 
values are in feet or other units than meters (with a conversion factor
14
 
<i>meters:</i>, defined in PROJ_UNITS), they must be converted to meters using
15
 
the parameter <i>zfactor</i>.
16
 
<p>
17
 
The <i>aspect</i> output raster map indicates the direction that slopes are
18
 
facing. The aspect categories represent the number degrees of east. Category
19
 
and color table files are also generated for the aspect map layer. The aspect
20
 
categories represent the number degrees of east and they increase
21
 
counterclockwise: 90deg is North, 180 is West, 270 is South 360 is East. The
22
 
aspect value 0 is used to indicate undefined aspect in flat areas with slope=0.
23
 
<p>
24
 
The <i>slope</i> output raster map contains slope values, stated in degrees of
25
 
inclination from the horizontal if <i>format</i>=degrees option (the default)
26
 
is chosen, and in percent rise if <i>format</i>=percent option is chosen.
27
 
Category and color table files are generated.
28
 
<p>
29
 
Profile and tangential curvatures are the curvatures in the direction of
30
 
steepest slope and in the direction of the contour tangent respectively. The
31
 
curvatures are expressed as 1/metres, e.g. a curvature of 0.05 corresponds to a
32
 
radius of curvature of 20m. Convex form values are positive and concave form values
33
 
are negative.
34
 
 
35
 
<p>
36
 
<table width="100%" border="0">
37
 
 <tr valign="baseline">
38
 
  <td>
39
 
  <center>
40
 
    <img src="dem.png" border="1">
41
 
    <p>
42
 
    Example DEM
43
 
    <br><br>
44
 
  </center>
45
 
  </td>
46
 
  <td>
47
 
  </td>
48
 
 </tr>
49
 
 
50
 
 <tr valign="baseline">
51
 
  <td>
52
 
  <center>
53
 
    <img src="slope.png" border="1">
54
 
    <p>
55
 
    Slope (degree) from example DEM
56
 
    <br><br>
57
 
  </center>
58
 
  </td>
59
 
  <td>
60
 
  <center>
61
 
    <img src="aspect.png" border="1">
62
 
    <p>
63
 
    Aspect (degree) from example DEM
64
 
    <br><br>
65
 
  </center>
66
 
  </td>
67
 
 </tr>
68
 
 
69
 
 <tr valign="baseline">
70
 
  <td>
71
 
  <center>
72
 
    <img src="tcurv.png" border="1">
73
 
    <p>
74
 
    Tangential curvature (m<sup>-1</sup>) from example DEM
75
 
    <br><br>
76
 
  </center>
77
 
  </td>
78
 
  <td>
79
 
  <center>
80
 
    <img src="pcurv.png" border="1">
81
 
    <p>
82
 
    Profile curvature (m<sup>-1</sup>) from example DEM
83
 
    <br><br>
84
 
  </center>
85
 
  </td>
86
 
  <td>
87
 
  </td>
88
 
 </tr>
89
 
</table>
90
 
 
91
 
<p>
92
 
For some applications, the user will wish to use a reclassified raster map
93
 
of slope that groups slope values into ranges of slope. This can be done using
94
 
<i><a href="r.reclass.html">r.reclass</a></i>. An example of a useful
95
 
reclassification is given below:
96
 
<div class="code"><pre>          category      range   category labels
97
 
                     (in degrees)    (in percent)
98
 
 
99
 
             1         0-  1             0-  2%
100
 
             2         2-  3             3-  5%
101
 
             3         4-  5             6- 10%
102
 
             4         6-  8            11- 15%
103
 
             5         9- 11            16- 20%
104
 
             6        12- 14            21- 25%
105
 
             7        15- 90            26% and higher
106
 
 
107
 
     The following color table works well with the above
108
 
     reclassification.
109
 
 
110
 
          category   red   green   blue
111
 
 
112
 
             0       179    179     179
113
 
             1         0    102       0
114
 
             2         0    153       0
115
 
             3       128    153       0
116
 
             4       204    179       0
117
 
             5       128     51      51
118
 
             6       255      0       0
119
 
             7         0      0       0</pre></div>
120
 
 
121
 
<h2>NOTES</h2>
122
 
 
123
 
To ensure that the raster elevation map layer is not inappropriately resampled,
124
 
the settings for the current region are modified slightly (for the execution
125
 
of the program only): the resolution is set to match the resolution of
126
 
the elevation map and the edges of the region (i.e. the north, south, east
127
 
and west) are shifted, if necessary, to line up along edges of the nearest
128
 
cells in the elevation map. If the user really wants the elevation map
129
 
resampled to the current region resolution, the -a flag should be specified.
130
 
 
131
 
<p>
132
 
The current mask is ignored.
133
 
 
134
 
<p>
135
 
The algorithm used to determine slope and aspect uses a 3x3 neighborhood
136
 
around each cell in the elevation file. Thus, it is not possible to determine
137
 
slope and aspect for the cells adjacent to the edges in the elevation map
138
 
layer. These cells are assigned a "zero slope" value (category 0) in both
139
 
the slope and aspect raster map layers.
140
 
 
141
 
<p>
142
 
Horn's formula is used to find the first order derivatives in x and y directions.
143
 
 
144
 
<p>
145
 
Only when using integer elevation models, the aspect is biased in 0,
146
 
45, 90, 180, 225, 270, 315, and 360 directions; i.e., the distribution
147
 
of aspect categories is very uneven, with peaks at 0, 45,..., 360 categories.
148
 
When working with floating point elevation models, no such aspect bias occurs.
149
 
 
150
 
<p>
151
 
Because most cells with a very small slope end up having category 0,
152
 
45, ..., 360, it is sometimes possible to reduce the bias in these directions
153
 
by filtering out the aspect in areas where the terrain is almost flat. A new
154
 
option <i>min_slp_allowed</i> was added to specify the minimum slope for which
155
 
aspect is computed. The aspect for all cells with slope &lt;
156
 
<i>min_slp_allowed</i> is set to <b>null</b>.
157
 
 
158
 
 
159
 
<h2>REFERENCE</h2>
160
 
 
161
 
<ul>
162
 
<li> Horn, B. K. P. (1981). <i>Hill Shading and the Reflectance Map</i>, Proceedings
163
 
of the IEEE, 69(1):14-47.
164
 
<li> Mitasova, H. (1985). <i>Cartographic aspects of computer surface modeling. PhD thesis.</i>
165
 
Slovak Technical University , Bratislava
166
 
<li> Hofierka, J., Mitasova, H., Neteler, M., 2009. <i>Geomorphometry in GRASS GIS.</i>
167
 
In: Hengl, T. and Reuter, H.I. (Eds), <i>Geomorphometry: Concepts, Software, Applications. </i>
168
 
Developments in Soil Science, vol. 33, Elsevier, 387-410 pp,
169
 
<a href="http://www.geomorphometry.org">http://www.geomorphometry.org</a>
170
 
</ul>
171
 
 
172
 
<h2>SEE ALSO</h2>
173
 
<em><a href="r.mapcalc.html">r.mapcalc</a></em>,
174
 
<em><a href="r.neighbors.html">r.neighbors</a></em>,
175
 
<em><a href="r.reclass.html">r.reclass</a></em>,
176
 
<em><a href="r.rescale.html">r.rescale</a></em>
177
 
 
178
 
<h2>AUTHORS</h2>
179
 
Michael Shapiro, U.S.Army Construction Engineering Research Laboratory<br>
180
 
Olga Waupotitsch, U.S.Army Construction Engineering Research Laboratory
181
 
 
182
 
<p>
183
 
<i>Last changed: $Date: 2011-11-08 12:29:50 +0100 (Tue, 08 Nov 2011) $</i>