Description Usage Arguments Details Value References See Also Examples
This function performs spherical wavelet reconstruction.
1 |
swd |
an object of class ‘swd’ |
This function performs spherical wavelet reconstruction.
recon |
the spherical wavelet reconstruction |
Oh, H-S. and Li, T-H. (2004) Estimation of global temperature fields from scattered observations by a spherical-wavelet-based spatially adaptive method. Journal of the Royal Statistical Society Ser. B, 66, 221–238.
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 | ### Observations of year 1967
data(temperature)
names(temperature)
# Temperatures on 939 weather stations of year 1967
temp67 <- temperature$obs[temperature$year == 1967]
# Locations of 939 weather stations
latlon <- temperature$latlon[temperature$year == 1967, ]
### Network design by BUD
data(netlab)
### Bandwidth for Poisson kernel
eta <- c(0.961, 0.923, 0.852, 0.723, 0.506)
### SBF representation of the observations by pls
out.pls <- sbf(obs=temp67, latlon=latlon, netlab=netlab, eta=eta,
method="pls", grid.size=c(50, 100), lambda=0.89)
### Decomposition
out.dpls <- swd(out.pls)
### Thresholding
out.univ <- swthresh(out.dpls, policy="universal", by.level=TRUE,
type="hard", nthresh=4)
### Reconstruction
out.rec <- swr(out.univ)
sw.plot(z=out.rec, type="recon", xlab="", ylab="")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.