cpt.2d(x, crystal.names=NULL, basis="block.cpt.2d", cost.fun="threshold", n.levels=NULL, taper="poly2", dct.type=2, boundary="periodic", n.taper=NULL, scale=NULL, thresh=NULL, p=1, prob=.75) icpt.2d(x)
icpt.2d returns an image if x is an object of class cpt.2d.
The algorithms for the taper functions are given in the S+WAVELETS User's Manual, in the section "Cosine Packet Algorithms". They are discussed in greater depth in Wickerhauser (1994).
xx <- phone-mean(phone) par(mfrow=c(1,2)) image(xx) bb2 <- crystal.names("block.dct.2d", 2) cc2 <- cpt.2d(xx, crystal.names=bb2, taper="poly2") thresh <- rep(3.7, 16) ss2 <- shrink(cc2, thresh) zz <- reconstruct(ss2) image(zz)