.TH r.fillnulls 1grass "" "GRASS 6.4.4" "Grass User's Manual"
.SH NAME
\fI\fBr.fillnulls\fR\fR  - Fills no-data areas in raster maps using spline interpolation.
.SH KEYWORDS
raster, elevation, interpolation
.SH SYNOPSIS
\fBr.fillnulls\fR
.br
\fBr.fillnulls help\fR
.br
\fBr.fillnulls input\fR=\fIstring\fR \fBoutput\fR=\fIstring\fR  [\fBtension\fR=\fIfloat\fR]   [\fBsmooth\fR=\fIfloat\fR]  \fBmethod\fR=\fIstring\fR  [\-\-\fBoverwrite\fR]  [\-\-\fBverbose\fR]  [\-\-\fBquiet\fR] 
.SS Flags:
.IP "\fB\-\-overwrite\fR" 4m
.br
Allow output files to overwrite existing files
.IP "\fB\-\-verbose\fR" 4m
.br
Verbose module output
.IP "\fB\-\-quiet\fR" 4m
.br
Quiet module output
.PP
.SS Parameters:
.IP "\fBinput\fR=\fIstring\fR" 4m
.br
Name of input raster map in which to fill nulls
.IP "\fBoutput\fR=\fIstring\fR" 4m
.br
Name for output raster map with nulls filled by interpolation
.IP "\fBtension\fR=\fIfloat\fR" 4m
.br
Spline tension parameter
.br
Default: \fI40.\fR
.IP "\fBsmooth\fR=\fIfloat\fR" 4m
.br
Spline smoothing parameter
.br
Default: \fI0.1\fR
.IP "\fBmethod\fR=\fIstring\fR" 4m
.br
Interpolation method
.br
Options: \fIbilinear,bicubic,rst\fR
.br
Default: \fIrst\fR
.PP
.SH DESCRIPTION
\fIr.fillnulls\fR fills NULL pixels (no data areas) in input map and
stores filled map to a new output map. The fill data are interpolated 
from the no data area boundaries buffer using \fIv.surf.rst\fR or
\fIv.surf.bspline\fR interpolation.
.SH NOTES
Each area boundary buffer is set to three times the map resolution to get nominally
three points around the edge. This way the algorithm interpolates into the hole with
a trained slope and curvature at the edges, in order to avoid that such a flat plane
is generated in a hole.
.PP
During the interpolation following warning may occur when using the RST method:
.PP
Warning: strip exists with insufficient data
.br
Warning: taking too long to find points for interpolation--please change
the region to area where your points are
.PP
This warning is generated if large data holes exist within the surface.
As the idea of \fIr.fillnulls\fR is to fill such holes, the user may
ignore the warning. The interpolation will be continued. However, the user
may pay attention to below notes.
.SH NOTES
When using the default RST method, the algorithm is based on \fIv.surf.rst\fR
regularized splines with tension interpolation module which interpolates the
raster cell values for NULL data areas from the boundary values of the NULL
data area. An eventual raster MASK is respected during the NULL data area(s)
filling. The interpolated values are patched into the NULL data area(s) of
the input map and saved into a new raster map.
Otherwise, either the bilinear or bicubic method can be selected (based on
\fIv.surf.bspline\fR).
.SH WARNING
Depending on the shape of the NULL data area(s) problems may occur due to an
insufficient number of input cell values for the interpolation process. Most
problems will occur if a NULL data area reaches a large amount of the map
boundary. The user will have to carefully check the result using
\fIr.mapcalc\fR (generating a difference map to the
input map and applying the "differences" color table with \fIr.colors\fR)
and/or \fId.what.rast\fR to query individual cell values.
.SH EXAMPLE
In this example, the SRTM elevation map in the
North Carolina sample dataset location is filtered for outlier
elevation values; missing pixels are then re-interpolated to obtain
a complete elevation map:
\fC
.DS
.br
g.region rast=elev_srtm_30m \-p
.br
d.mon x0
.br
d.histogram elev_srtm_30m
.br
# remove too low elevations (esp. lakes)
.br
r.mapcalc "elev_srtm_30m_filt = if(elev_srtm_30m < 50.0, null(), elev_srtm_30m)"
.br
d.histogram elev_srtm_30m_filt
.br
d.rast elev_srtm_30m_filt
.br
r.fillnulls input=elev_srtm_30m_filt output=elev_srtm_30m_complete tension=20
.br
d.histogram elev_srtm_30m_complete
.br
d.rast elev_srtm_30m_complete
.br
.DE
\fR
.SH SEE ALSO
\fI
r.fill.dir, 
r.mapcalc, 
v.surf.bspline,
v.surf.rst
\fR
.SH REFERENCES
.PP
 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,
pp.481-492
.PP
Mitasova H., Mitas L.,  Brown W.M.,  D.P. Gerdes, I.
Kosinovsky, Baker, T.1995, Modeling spatially and temporally distributed
phenomena: New methods and tools for GRASS GIS. \fIInternational Journal of
GIS\fR, 9 (4), special issue on Integrating GIS and Environmental modeling,
433-446.
.PP
Mitasova H.
and Mitas L. 1993: Interpolation by Regularized Spline with Tension: I.
Theory and Implementation, \fIMathematical Geology\fR 25, 641-655.
.PP
Mitasova H.
and Hofierka L. 1993: Interpolation by Regularized Spline with Tension:
II. Application to Terrain Modeling and Surface Geometry Analysis,
\fIMathematical Geology\fR 25, 657-667.
.SH AUTHORS
r.fillnulls: Markus Neteler, University of Hannover and Fondazione Edmund Mach
.PP
and authors of v.surf.rst
.br
Improvement by Hamish Bowman, NZ
.PP
\fILast changed: $Date: 2014-03-15 15:28:47 +0100 (Sat, 15 Mar 2014) $\fR
.PP
Full index
.PP
© 2003-2014 GRASS Development Team