sw.plot {SpherWave} | R Documentation |
This function performs plotting of observation, network design, field, SW coefficients, decomposition or reconstruction result.
sw.plot(sw = NULL, z = NULL, latlon = NULL, latlim = NULL, lonlim = NULL, type = "field", nlevel = 256, pch = NULL, cex = NULL, ...)
sw |
sbf or swd object |
z |
observations, network design labels or reconstruction |
latlon |
grid points of observation sites in degree |
latlim |
range of latitudes in degree |
lonlim |
range of longitudes in degree |
type |
specifies the type "obs", "network", "field", "swcoeff", "decom" or "recon" |
nlevel |
number of color levels used in legend strip |
pch |
either an integer specifying a symbol or a single character to be used as the default in plotting points |
cex |
a numerical value giving the amount by which plotting text and symbols should be scaled relative to the default |
... |
the usual arguments to the image function or plot function |
This function plots spherical wavelet results. Possible types are
`"obs"' for observations
`"network"' for network design
`"field"' for field
`"swcoeff"' for spherical wavelet coefficients
`"decom"' for decomposition results
`"recon"' for reconstruction result.
For `sw', sbf
or swd
object must be provided. For sbf
object, type `"obs"', `"network"', `"field"'
are possible whereas all types are possible for swd
object. Or specify `z' and `latlon'
without `sw'.
### 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, ] ### Draw the temperature data sw.plot(z=temp67, latlon=latlon, type="obs") ### Network design by BUD data(netlab) sw.plot(z=netlab, latlon=latlon, type="network") ### SBF representation of the observations #eta <- c(0.961,0.923,0.852,0.723,0.506) #out.pls <- sbf(obs=temp67, latlon=latlon, netlab=netlab, eta=eta, # method="pls", grid.size=c(50, 100), lambda=0.89) # observation #sw.plot(out.pls, type="obs") # network design #sw.plot(out.pls, type="network") # field #sw.plot(out.pls, type="field") ### Decomposition #out.dpls <- swd(out.pls) # observation #sw.plot(out.dpls, type="obs") # network design #sw.plot(out.dpls, type="network") # SBF representation of the observations #sw.plot(out.dpls, type="field") # sw coefficient #sw.plot(out.dpls, type="swcoeff") # decomposition result #sw.plot(out.dpls, type="decom") # Thresholding #out.univ <- swthresh(out.dpls, policy="universal", by.level=TRUE, # type="hard", nthresh=4) #par(oma=c(0,0,3.5,0)) #sw.plot(out.univ, type="decom") #mtext("Decomposition & Threshold", side = 3, outer = TRUE, # cex = 1.2, line = 1) # Reconstruction #out.rec <- swr(out.univ) #sw.plot(z=out.rec, type="recon", xlab="", ylab="")