[GRASS-SVN] r53216 - grass/branches/releasebranch_6_4/raster/r.surf.contour
svn_grass at osgeo.org
svn_grass at osgeo.org
Tue Sep 18 13:39:02 PDT 2012
Author: neteler
Date: 2012-09-18 13:39:02 -0700 (Tue, 18 Sep 2012)
New Revision: 53216
Modified:
grass/branches/releasebranch_6_4/raster/r.surf.contour/description.html
Log:
new example
Modified: grass/branches/releasebranch_6_4/raster/r.surf.contour/description.html
===================================================================
--- grass/branches/releasebranch_6_4/raster/r.surf.contour/description.html 2012-09-18 20:32:10 UTC (rev 53215)
+++ grass/branches/releasebranch_6_4/raster/r.surf.contour/description.html 2012-09-18 20:39:02 UTC (rev 53216)
@@ -60,7 +60,6 @@
the output map. Run <em> <a href="r.slope.aspect.html">r.slope.aspect</a>
</em> on <em>r.surf.contour</em> results to locate potential anomalies.
-
<p>
The running of <em>r.surf.contour</em> is very sensitive to the resolution of
rasterized vector map. If multiple contour lines go through the same raster,
@@ -69,6 +68,37 @@
raster cells). Since a flood fill algorithm is used, the program's running
time will grow exponentially with the distance between contour lines.
+<h2>EXAMPLE</h2>
+
+Example to create contour lines from elevation model, then recreating DEM
+from these contour lines along with differences analysis (North Carolina
+sample data set):
+
+<div class="code"><pre>
+g.region rast=elevation -p
+
+# get minimum elevation value
+r.univar elevation
+
+# generate vector contour lines
+r.contour input=elevation output=contours_5m step=5 minlevel=50
+
+# rasterize contour lines
+v.info -c contours_5m
+v.to.rast input=contours_5m output=contours_5m use=attr column=level
+
+# generate DEM from rasterized contour lines
+r.surf.contour input=contours_5m output=elevation_from_cont5m
+
+# calculate difference map
+r.mapcalc "diff = elevation - elevation_from_cont5m"
+r.colors diff color=differences
+
+# analyze differences statistically
+r.univar diff
+</pre></div>
+
+
<h2>SEE ALSO</h2>
<em><a href="r.mapcalc.html">r.mapcalc</a></em>,
More information about the grass-commit
mailing list