An orthonormal wavelet system in R, d ∈ N, is a countable collection of functions {ψ j,k}, j ∈ Z, k ∈ Z, ` = 1, . . . , L, of the form ψ j,k(x) = | deta|−j/2ψ`(a−jx− k) ≡ (Daj Tk ψ)(x) that is an orthonormal basis for L2(Rd), where a ∈ GLd(R) is an expanding matrix. The first such system to be discovered (almost one hundred years ago) is the Haar system for which L = d = 1, ψ1(x) = ψ(x) = χ[0,1...