Dual-tree complex 2D discrete wavelet transform (DWT).
cplxdual2D(x, J, Faf, af)icplxdual2D(w, J, Fsf, sf)
For the analysis of x
, the output is
wavelet
coefficients indexed by [[j]][[i]][[d1]][[d2]]
, where
For the synthesis of
w
, the output is
output signal.
2D array.
number of stages.
first stage analysis filters for tree i.
analysis filters for the remaining stages on tree i.
wavelet coefficients.
last stage synthesis filters for tree i.
synthesis filters for the preceeding stages.
Matlab: S. Cai, K. Li and I. Selesnick; R port: B. Whitcher
FSfarras
, farras
, afb2D
,
sfb2D
.
if (FALSE) {
## EXAMPLE: cplxdual2D
x = matrix(rnorm(32*32), 32, 32)
J = 5
Faf = FSfarras()$af
Fsf = FSfarras()$sf
af = dualfilt1()$af
sf = dualfilt1()$sf
w = cplxdual2D(x, J, Faf, af)
y = icplxdual2D(w, J, Fsf, sf)
err = x - y
max(abs(err))
}
Run the code above in your browser using DataLab