Geografic-Information-System/raster/r.param.scale/r.param.scale.html

<h2>DESCRIPTION</h2>

<b>r.param.scale</b> extracts terrain parameters from a DEM. Uses a
multi-scale approach by fitting a bivariate quadratic polynomial to a given
window size using least squares.

<p>The module calculates the following parameters (terminology is from Wood,
1996 with related terminology used in other GRASS modules listed in
brackets):
<ul>
<li>
<i>elev</i>: Generalised elevation value (for resampling purposes at different
scale)</li>

<li>
<i>slope</i>: Magnitude of maximum gradient (steepest slope angle)</li>

<li>
<i>aspect</i>: Direction of maximum gradient (steepest slope direction=flow direction)</li>

<li>
<i>profc</i> or profile curvature (curvature intersecting with the plane
defined by <i>Z</i> axis and maximum gradient direction). Positive values
describe convex profile curvature, negative values concave profile
curvature.</li>

<li>
<i>planc</i> or plan curvature (horizontal curvature, intersecting with
the <i>XY</i> plane)</li>

<li>
<i>longc</i> or longitudinal curvature (profile curvature intersecting
with the plane defined by the surface normal and maximum gradient direction)</li>

<li>
<i>crosc</i> or cross-sectional curvature (tangential curvature intersecting
with the plane defined by the surface normal and a tangent to the contour
- perpendicular to maximum gradient direction)</li>

<li>
<i>maxic</i> or maximum curvature (can be in any direction)</li>

<li>
<i>minic</i> or minimum curvature (in direction perpendicular to the direction
of of maximum curvature)</li>

<!--
<li>
<i>meanc</i> or mean curvature (average of maximum and minimum curvatures).</li>
-->

<li>
<i>feature</i>: Morphometric features: peaks, ridges, passes, channels, pits and planes</li>
</ul>

<h2>NOTES</h2>
In <i>r.param.scale</i> the direction of maximum gradient (considered
downslope) is stored as (West is 0 degree, East is +/- 180 degree):

<ul>
<li>
0..+180 degree from West to North to East</li>

<li>
0..-180 degree from West to South to East</li>
</ul>

Note that the aspect map is calculated differently from
<a href="r.slope.aspect.html">r.slope.aspect</a>.

<h2>Still to do</h2>

Fix bug when `constrain through central cell' option selected. Create color
tables for all output files (presently only on features).

<h2>EXAMPLE</h2>

The next commands will create a geomorphological map of the Spearfish region:

<div class="code"><pre>
g.region rast=elevation.10m -p
r.param.scale in=elevation.10m output=morphology param=feature size=9
</pre></div>

<p><center>
<img src="r_param_scale_morph.jpg" alt="r.param.scale generated geomorphological map"><br>
<i>Geomorphological map of a subregion in the Spearfish (SD) area</i>
</center>


<h2>SEE ALSO</h2>
<!-- not ported to GRASS 6 due to non-GPLness of numerical recipes.
<i><a href="d.param.scale.html">d.param.scale</a></i>
-->
<p>Java Code in
<a href="http://www.geog.le.ac.uk/jwo/research/LandSerf">LandSerf</a>
that implements the same procedure

<h2>REFERENCE</h2>

Wood, J. (1996): The Geomorphological characterisation of Digital Elevation
Models. Diss., Department of Geography, University of Leicester, U.K.

<br>online at:
<br><a href="http://www.soi.city.ac.uk/~jwo/phd/">http://www.soi.city.ac.uk/~jwo/phd/</a>

<h2>AUTHOR</h2>

<address>
<a href="MAILTO:jwo@le.ac.uk">jwo@le.ac.uk</a>
- <a href="http://www.geog.le.ac.uk/assist/index.html">ASSIST's home</a></address>

<p>Update to FP 3/2002: L. Potrich, M. Neteler, S. Menegon (ITC-irst)

<p><i>Last changed: $Date$</i>