swr: Spherical Wavelet Reconstruction of 'swd' Object

Description Usage Arguments Details Value References See Also Examples

Description

This function performs spherical wavelet reconstruction.

Usage

1

Arguments

swd

an object of class ‘swd’

Details

This function performs spherical wavelet reconstruction.

Value

recon

the spherical wavelet reconstruction

References

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.

See Also

sbf, swd, swthresh

Examples

 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="")

SeoncheolPark/SphereWaveu documentation built on May 29, 2019, 12:04 a.m.