get energy from the dualtree transform
square the wavelet coefficients and apply the bias correction
get_en(dt, correct = "none", N = ncol(dt))
dt
: array of J x nx x ny x 6
complex wavelet coefficients, output of dtcwt(..., dec=FALSE)
correct
: type of correction, either "b"
or "b_bp"
, any other value results in no correction at all.N
: the smallest whole power of two larger than or equal to the dimensions of the input image, usually just ncol(dt)
an array of the same dimensions as dt
The bias correction matrix should correspond to the filter bank used in the transform, for details on the matrices see A
.
Nelson, J. D. B., A. J. Gibberd, C. Nafornita, and N. Kingsbury (2018) doi:10.1007/s11222-017-9784-0
A
Useful links