/usr/share/doc/grass-doc/html/r.resamp.rst.html is in grass-doc 6.4.3-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<title>GRASS GIS manual: r.resamp.rst</title>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<link rel="stylesheet" href="grassdocs.css" type="text/css">
</head>
<body bgcolor="white">
<img src="grass_logo.png" alt="GRASS logo"><hr align=center size=6 noshade>
<h2>NAME</h2>
<em><b>r.resamp.rst</b></em> - Reinterpolates and optionally computes topographic analysis from input raster map to a new raster map (possibly with different resolution) using regularized spline with tension and smoothing.
<h2>KEYWORDS</h2>
raster, resample
<h2>SYNOPSIS</h2>
<b>r.resamp.rst</b><br>
<b>r.resamp.rst help</b><br>
<b>r.resamp.rst</b> [-<b>td</b>] <b>input</b>=<em>name</em> <b>ew_res</b>=<em>float</em> <b>ns_res</b>=<em>float</em> [<b>elev</b>=<em>string</em>] [<b>slope</b>=<em>string</em>] [<b>aspect</b>=<em>string</em>] [<b>pcurv</b>=<em>string</em>] [<b>tcurv</b>=<em>string</em>] [<b>mcurv</b>=<em>string</em>] [<b>smooth</b>=<em>string</em>] [<b>maskmap</b>=<em>string</em>] [<b>overlap</b>=<em>integer</em>] [<b>zmult</b>=<em>float</em>] [<b>tension</b>=<em>float</em>] [<b>theta</b>=<em>float</em>] [<b>scalex</b>=<em>float</em>] [--<b>overwrite</b>] [--<b>verbose</b>] [--<b>quiet</b>]
<h3>Flags:</h3>
<DL>
<DT><b>-t</b></DT>
<DD>Use dnorm independent tension</DD>
<DT><b>-d</b></DT>
<DD>Output partial derivatives instead of topographic parameters</DD>
<DT><b>--overwrite</b></DT>
<DD>Allow output files to overwrite existing files</DD>
<DT><b>--verbose</b></DT>
<DD>Verbose module output</DD>
<DT><b>--quiet</b></DT>
<DD>Quiet module output</DD>
</DL>
<h3>Parameters:</h3>
<DL>
<DT><b>input</b>=<em>name</em></DT>
<DD>Name of input raster map</DD>
<DT><b>ew_res</b>=<em>float</em></DT>
<DD>Desired east-west resolution</DD>
<DT><b>ns_res</b>=<em>float</em></DT>
<DD>Desired north-south resolution</DD>
<DT><b>elev</b>=<em>string</em></DT>
<DD>Output z-file (elevation) map</DD>
<DT><b>slope</b>=<em>string</em></DT>
<DD>Output slope map (or fx)</DD>
<DT><b>aspect</b>=<em>string</em></DT>
<DD>Output aspect map (or fy)</DD>
<DT><b>pcurv</b>=<em>string</em></DT>
<DD>Output profile curvature map (or fxx)</DD>
<DT><b>tcurv</b>=<em>string</em></DT>
<DD>Output tangential curvature map (or fyy)</DD>
<DT><b>mcurv</b>=<em>string</em></DT>
<DD>Output mean curvature map (or fxy)</DD>
<DT><b>smooth</b>=<em>string</em></DT>
<DD>Name of raster map containing smoothing</DD>
<DT><b>maskmap</b>=<em>string</em></DT>
<DD>Name of raster map to be used as mask</DD>
<DT><b>overlap</b>=<em>integer</em></DT>
<DD>Rows/columns overlap for segmentation</DD>
<DD>Default: <em>3</em></DD>
<DT><b>zmult</b>=<em>float</em></DT>
<DD>Multiplier for z-values</DD>
<DD>Default: <em>1.0</em></DD>
<DT><b>tension</b>=<em>float</em></DT>
<DD>Spline tension value</DD>
<DD>Default: <em>40.</em></DD>
<DT><b>theta</b>=<em>float</em></DT>
<DD>Anisotropy angle (in degrees)</DD>
<DT><b>scalex</b>=<em>float</em></DT>
<DD>Anisotropy scaling factor</DD>
</DL>
<h2>DESCRIPTION</h2>
<i>r.resamp.rst</i> reinterpolates the values a from given raster map (named
<i>input</i>) to a new raster map (named <i>elev</i>).
This module is intended for reinterpolation of continuous data
to a different resolution rather than for interpolation from scattered data
(use the <em>v.surf.*</em> modules for that purpose).
Reinterpolation (resampling) is done to higher, same or lower resolution
specified by the <i>ew_res</i> and <i>ns_res</i> parameters.
<p>
All resulting raster maps are created using the settings of the current
region (which may be different from that of the <i>input</i> raster map).
<p>
Optionally, and simultaneously with interpolation, topographic parameters
are computed from an input raster map containing z-values of elevation/depth: slope,
aspect, profile curvature (measured in the direction of steepest
slope), tangential curvature (measured in the direction of a tangent to
contour line) and/or mean curvature are computed from and saved as raster maps
as specified by the options <i>slope, aspect, pcurv, tcurv, mcurv</i> respectively.
<p>
If the <i>-d</i> flag is set the program outputs partial derivatives fx, fy,
fxx, fxy, and fyy instead of slope, aspect and curvatures.
<p>
For noisy data it is possible to define spatially variable smoothing
by providing a raster map named by the <i>smooth</i> option containing
smoothing parameters.
With the smoothing parameter set to zero (<i>smooth</i> is not given or
contains zero data), the resulting surface passes exactly through the data
points.
<p>
The user can also define a raster map (named with <i>maskmap</i>) which will
be used as a mask. The interpolation is skipped for cells which have zero
or NULL value in the mask.
<p>
Zero values will be assigned to these cells in all output raster maps.
<!-- is this still true???? -->
<p>
The <i>zmult</i> parameter allows the user to rescale the z-values which may
be useful, e.g., for transformation of elevations given in feet to meters,
so that the proper values of slopes and curvatures can be computed.
The default value is 1.
<p>
A regularized spline with tension method is used for the interpolation.
The <i>tension</i> parameter tunes the character of the resulting surface
from thin plate to membrane. Higher values of tension parameter reduce the
overshoots that can appear in surfaces with rapid change of gradient.
<p>
The <i>-t</i> flag can be set to use "dnorm independent tension".
<!-- but what is "dnorm independent tension" ??? -->
<p>
The interpolation is performed for overlapping rectangular segments.
The user can define the width of overlap (in number of cells) with
the <i>overlap</i> option. The default value is 3.
<br>
<h2>NOTES</h2>
<i>r.resamp.rst</i> uses regularized spline with tension for interpolation
(as described in Mitasova and Mitas, 1993).
<p>
The region is temporarily changed while writing output files with desired
resolution. Topographic parameters are computed in the same way as in
the <i>v.surf.rst</i> module. (See also Mitasova and Hofierka, 1993)
<p>
The raster map used with the <i>smooth</i> option should contain variable
smoothing parameters. These can be derived from errors, slope, etc. using
the <i>r.mapcalc</i> module.
<p>
The program gives warning when significant overshoots appear and higher
tension should be used. However, with tension set too high the resulting surface
changes its behavior to a membrane (rubber sheet stretched over the data
points resulting in a peak or pit in each given point and everywhere else
the surface goes rapidly to trend). Smoothing can be used to reduce
the overshoots. When overshoots occur the resulting <i>elev</i> file will
have white color in the locations of overshoots since the color table for
the output file is the same as colortable for raster input file.
<p>
The program checks the numerical stability of the algorithm by computation
of values at given points, and prints the maximum difference found into the
history file of raster map <i>elev</i> (view with <i>r.info</i>).
An increase in tension is suggested if the difference is unacceptable.
For computations with smoothing set to 0 this difference should be 0.
With a smoothing parameter greater than zero the surface will not pass through
the data points exactly, and the higher the parameter the closer the surface
will be to the trend.
<p>
The program writes the values of parameters used in computation into
the comment part of the <i>elev</i> map history file. Additionally the
following values are also written to assist in the evaluation of results
and choosing of suitable parameters:
<ul>
<li>minimum and maximum z values in the data file (zmin_data, zmax_data) and
in the interpolated raster map (zmin_int, zmax_int),</li>
<li>maximum difference between the given and interpolated z value at a given
point (errtotal),</li>
<li>rescaling parameter used for normalization (dnorm), which influences the
tension.</li>
</ul>
<p>
The program gives a warning when the user wants to interpolate outside
the region given by the <i>input</i> raster map's header data. Zooming into the
area where the points are is suggested in this case.
<p>
When a mask is used, the program uses all points in the given region for
interpolation, including those in the area which is masked out, to ensure
proper interpolation along the border of the mask. It therefore does not
mask out the data points; if this is desirable, it must be done outside
<i>r.resamp.rst</i> before processing.
<h2>SEE ALSO</h2>
<a href="g.region.html">g.region</a>,
<a href="r.info.html">r.info</a>,
<a href="r.resample.html">r.resample</a>,
<a href="r.mapcalc.html">r.mapcalc</a>,
<a href="r.surf.contour.html">r.surf.contour</a>,
<a href="v.surf.rst.html">v.surf.rst</a>
<h2>AUTHORS</h2>
<i>Original version of program (in FORTRAN):</i>
<br>Lubos Mitas, NCSA, University of Illinois at Urbana Champaign, Il
<br>Helena Mitasova, US Army CERL, Champaign, Illinois
<p><i>Modified program (translated to C, adapted for GRASS , segmentation
procedure):</i>
<br>Irina Kosinovsky, US Army CERL.
<br>Dave Gerdes, US Army CERL.
<h2>REFERENCES</h2>
Mitas, L., Mitasova, H., 1999, Spatial Interpolation. In: P.Longley, M.F.
Goodchild, D.J. Maguire, D.W.Rhind (Eds.), Geographical Information Systems:
Principles, Techniques, Management and Applications, Wiley, 481-492.
<p>Mitasova, H. and Mitas, L., 1993. Interpolation by regularized spline
with tension: I. Theory and implementation, Mathematical Geology No.25
p.641-656.
<p>Mitasova, H. and Hofierka, L., 1993. Interpolation by regularized spline
with tension: II. Application to terrain modeling and surface geometry
analysis, Mathematical Geology No.25 p.657-667.
<p>Talmi, A. and Gilat, G., 1977. Method for smooth approximation
of data, Journal of Computational Physics , 23, pp 93-123.
<p>Wahba, G., 1990. Spline models for observational data, CNMS-NSF Regional
Conference series in applied mathematics, 59, SIAM, Philadelphia, Pennsylvania.
<p><i>Last changed: $Date: 2008-05-16 12:09:06 -0700 (Fri, 16 May 2008) $</i>
<HR>
<P><a href="index.html">Main index</a> - <a href="raster.html">raster index</a> - <a href="full_index.html">Full index</a></P>
<P>© 2003-2013 <a href="http://grass.osgeo.org">GRASS Development Team</a></p>
</body>
</html>
|