Interpolation/
extrapsmooth.pro
similar to extrapolate but could to the job in a better way
because the extrapolated values are smoothed...
takes more time than extrapolate.
extrapolate data where mskin is equal 0 by filling
step by step the coastline points with the mean value of the 8 neighbors.
extrapsmooth
Interpolation
result = extrapsmooth(in, mskin, X_PERIODIC=X_PERIODIC, MINVAL=scalar, MAXVAL=scalar, GE0=scalar 0 or 1, _EXTRA=_EXTRA)
Return value
the extrapolated array with no more masked values
Parameters
in
in
required
type: 2d array
data to be extrapolate
mskin
in
required
type: 2d array or -1
a 2D array, the land-sea mask of the output data (1 on ocean, 0 on land)
put -1 if input data are not masked
Keywords
X_PERIODIC
MINVAL
type: scalar
default: not used
to specify a minimum value to the extrapolated values
MAXVAL
type: scalar
default: not used
to specify a maximum value to the extrapolated values
GE0
type: scalar 0 or 1
default: 0
put 1 to force the extrapolated values to be larger than 0, same as using minval=0.
_EXTRA
to be able to call extrapsmooth with _extra keyword
Examples
IDL> a=extrapsmooth(dist(jpi,jpj),tmask[*,*,0],/x_periodic)
IDL> tvplus, a
compare to extrapolate result:
IDL> b=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic)
IDL> tvplus, b, window = 1
Version history
Version
$Id: extrapsmooth.pro 371 2008-08-07 09:32:02Z pinsard $
History
January 2007: Sebastien Masson (smasson@lodyc.jussieu.fr)
Known issues
Restrictions
You cannot specify the number of iterations done in the extrapolation process
Produced by IDLdoc 2.0.