|
@@ -299,30 +299,31 @@ Using the <b>where</b> parameter, the interpolation can be limited to
|
|
|
use only a subset of the input vectors.
|
|
|
|
|
|
<p>
|
|
|
-Spearfish example (we simulate randomly distributed elevation
|
|
|
-measures):
|
|
|
+North Carolina example (we simulate randomly distributed elevation
|
|
|
+measures which we interpolate to a gap-free elevation surface):
|
|
|
|
|
|
<div class="code"><pre>
|
|
|
-g.region raster=elevation.10m -p
|
|
|
-# random elevation extraction
|
|
|
-r.random elevation.10m vector_output=elevrand n=200
|
|
|
+g.region raster=elevation -p
|
|
|
+# random elevation extraction of 500 samplings
|
|
|
+r.random elevation vector_output=elevrand n=500
|
|
|
v.info -c elevrand
|
|
|
v.db.select elevrand
|
|
|
|
|
|
# interpolation based on all points
|
|
|
v.surf.rst elevrand zcol=value elevation=elev_full
|
|
|
-r.colors elev_full rast=elevation.10m
|
|
|
+# apply the color table of the original raster map
|
|
|
+r.colors elev_full raster=elevation
|
|
|
d.rast elev_full
|
|
|
d.vect elevrand
|
|
|
|
|
|
# interpolation based on subset of points (only those over 1300m/asl)
|
|
|
v.surf.rst elevrand zcol=value elevation=elev_partial where="value > 1300"
|
|
|
-r.colors elev_partial rast=elevation.10m
|
|
|
+r.colors elev_partial raster=elevation
|
|
|
d.rast elev_partial
|
|
|
d.vect elevrand where="value > 1300"
|
|
|
</pre></div>
|
|
|
|
|
|
-<h2> REFERENCES</h2>
|
|
|
+<h2>REFERENCES</h2>
|
|
|
|
|
|
<ul>
|
|
|
<li><a href="http://www4.ncsu.edu/~hmitaso/gmslab/papers/IEEEGRSL2005.pdf">
|