[GRASS-SVN] r42648 - grass/trunk/raster/r.resamp.bspline

svn_grass at osgeo.org svn_grass at osgeo.org
Wed Jun 23 13:05:55 EDT 2010


Author: mmetz
Date: 2010-06-23 17:05:55 +0000 (Wed, 23 Jun 2010)
New Revision: 42648

Modified:
   grass/trunk/raster/r.resamp.bspline/main.c
   grass/trunk/raster/r.resamp.bspline/r.resamp.bspline.html
Log:
update manual and option description

Modified: grass/trunk/raster/r.resamp.bspline/main.c
===================================================================
--- grass/trunk/raster/r.resamp.bspline/main.c	2010-06-23 14:02:35 UTC (rev 42647)
+++ grass/trunk/raster/r.resamp.bspline/main.c	2010-06-23 17:05:55 UTC (rev 42648)
@@ -89,7 +89,7 @@
     stepE_opt->type = TYPE_DOUBLE;
     stepE_opt->required = NO;
     stepE_opt->description =
-	_("Length of each spline step in the east-west direction. Default: ewres.");
+	_("Length of each spline step in the east-west direction. Default: 1.5 * ewres.");
     stepE_opt->guisection = _("Settings");
 
     stepN_opt = G_define_option();
@@ -97,7 +97,7 @@
     stepN_opt->type = TYPE_DOUBLE;
     stepN_opt->required = NO;
     stepN_opt->description =
-	_("Length of each spline step in the north-south direction. Default: nsres.");
+	_("Length of each spline step in the north-south direction. Default: 1.5 * nsres.");
     stepN_opt->guisection = _("Settings");
 
     method_opt = G_define_option();

Modified: grass/trunk/raster/r.resamp.bspline/r.resamp.bspline.html
===================================================================
--- grass/trunk/raster/r.resamp.bspline/r.resamp.bspline.html	2010-06-23 14:02:35 UTC (rev 42647)
+++ grass/trunk/raster/r.resamp.bspline/r.resamp.bspline.html	2010-06-23 17:05:55 UTC (rev 42648)
@@ -13,27 +13,29 @@
 <p>
 Spline step values <b><i>se</i></b> for the east-west direction and
 <b><i>sn</i></b> for the north-south direction should not be smaller than
-the east-west and north-south resolutions of the input map. For very large
-areas with missing values (NULL cells), larger spline step values may be
-required, but most of the time the defaults (1.5 x resolution) should be fine.
+the east-west and north-south resolutions of the input map. For a raster
+map without NULL cells, 1 * resolution can be used, but check for
+undershoots and overshoots. For very large areas with missing values
+(NULL cells), larger spline step values may be required, but most of the
+time the defaults (1.5 x resolution) should be fine.
 <p>
 The Tykhonov regularization parameter ("<b><i>lambda</i></b>") acts to
 smooth the interpolation. With a small <b><i>lambda</i></b>, the
 interpolated surface closely follows observation points; a larger value
 will produce a smoother interpolation. Reasonable values are 0.0001,
-0.001, 0.005, 0.01, 0.02, 0.05, 0.1. For seamless NULL cell interpolation,
-a small value is required and default is set to 0.001.
+0.001, 0.005, 0.01, 0.02, 0.05, 0.1 (needs more testing). For seamless
+NULL cell interpolation, a small value is required and default is set to 0.005.
 <p>
 From a theoretical perspective, the interpolating procedure takes place in two
-parts: the first is an estimate of the linear coefficients of a spline function
-is derived from the observation points using a least squares regression; the
+parts: the first is an estimate of the linear coefficients of a spline function;
+these are derived from the observation points using a least squares regression; the
 second is the computation of the interpolated surface (or interpolated vector
 points). As used here, the splines are 2D piece-wise non-zero polynomial
-functions calculated within a limited, 2D area. The length of each spline step
+functions calculated within a limited 2D area. The length of each spline step
 is defined by <b><i>se</i></b> for the east-west direction and
-<b><i>sn</i></b> for the north-south direction. For optimum performance, the
+<b><i>sn</i></b> for the north-south direction. For optimal performance, the
 spline step values should be no less than the east-west and north-south
-resolutions of the input map. Each non NULL cell observation is modeled as a
+resolutions of the input map. Each non-NULL cell observation is modeled as a
 linear function of the non-zero splines in the area around the observation.
 The least squares regression predicts the the coefficients of these linear functions.
 Regularization avoids the need to have one one observation and one coefficient



More information about the grass-commit mailing list