DIRECTIONAL COMPACTLY SUPPORTED BOX SPLINE TIGHT FRAMELETS
WITH SIMPLE STRUCTURE

arXiv:1708.08421v1 [cs.IT] 28 Aug 2017

BIN HAN, TAO LI, AND XIAOSHENG ZHUANG
Abstract. To effectively capture singularities in high-dimensional data and functions, multivariate
compactly supported tight framelets, having directionality and derived from refinable box splines, are
of particular interest in both theory and applications. The d-dimensional Haar refinable function χ[0,1]d
is a simple example of refinable box splines. For every dimension d ∈ N, in this paper we construct a
directional compactly supported d-dimensional Haar tight framelet such that all its high-pass filters in
its underlying tight framelet filter bank have only two nonzero coefficients with opposite signs and they
exhibit totally (3d − 1)/2 directions in dimension d. Furthermore, applying the projection method to
such directional Haar tight framelets, from every refinable box spline in every dimension, we construct
a directional compactly supported box spline tight framelet with simple structure such that all the
high-pass filters in its underlying tight framelet filter bank have only two nonzero coefficients with
opposite signs. Moreover, such compactly supported box spline tight framelets can achieve arbitrarily
high numbers of directions by using refinable box splines with increasing supports.

To capture singularities in many high-dimensional data such as images/videos, directional representations are of great importance in both theory and applications, for example, see curvelets and
shearlets in [2, 7] and tensor product complex tight framelets in [9, 15]. On the other hand, (refinable) box splines are widely used in both approximation theory and wavelet analysis. Motivated
by the interesting example of a two-dimensional directional Haar tight framelet constructed in [17]
which has impressive performance in parallel magnetic resonance imaging (pMRI), in this paper we
construct compactly supported tight framelets with directionality and very simple structures from
the Haar refinable functions and all refinable box splines in all dimensions. All the high-pass filters
in such directional tight framelets have only two nonzero coefficients with oppositive signs. Consequently, all of them naturally exhibit directionality and their associated fast framelet transforms can
be efficiently implemented through simple difference operations.
Let us first recall some definitions and notation. By l0 (Zd ) we denote the set of all finitely supported
d
d
d
sequences/filters
k∈Z : Z → C on Z . For a filter a ∈ l0 (Z ), its Fourier series is defined
P a = {a(k)}
to be b
a(ξ) := k∈Zd a(k)e−ik·ξ for ξ ∈ Rd , which is a 2πZd -periodic trigonometric polynomial. In
particular, by δ we denote the Dirac sequence such that δ(0) = 1 and δ(k) = 0 for all Zd \{0}. For
γ ∈ Zd , we also use the notation δ γ to stand for the sequence δ(· − γ), i.e., δ γ (γ) = 1 and δ γ (k) = 0
for all k ∈ Zd \{γ}. Note that c
δ γ (ξ) = e−iγ·ξ . For filters a, b1 , . . . , bs ∈ l0 (Zd ), we say that a filter
bank {a; b1 , . . . , bs } is a (d-dimensional dyadic) tight framelet filter bank if
s
X
b
a(ξ)b
bbℓ (ξ)bbℓ (ξ + πω) = δ(ω),
∀ ξ ∈ Rd , ω ∈ {0, 1}d .
(1)
a(ξ + πω) +
ℓ=1

In the spatial domain, the equations in (1) can be equivalently rewritten as
s X
X
X
bℓ (γ + 2k)bℓ (n + γ + 2k) = 2−d δ(n),
a(γ + 2k)a(n + γ + 2k) +
k∈Zd

ℓ=1

k∈Zd

∀ γ ∈ {0, 1}d

(2)

2010 Mathematics Subject Classification. 42C40, 42C15, 41A15, 65D07.
Key words and phrases. Directional tight framelets, tight framelet filter banks, Haar refinable functions, refinable
box splines, box spline filters, directionality.
Research of B. Han is supported in part by NSERC Canada under Grant RGP 228051. Research of T. Li is
supported by a Mitacs global link project. Research of X. Zhuang is supported in part by the Research Grants Council
of Hong Kong (Project No. CityU 11304414).
1

2

BIN HAN, TAO LI, AND XIAOSHENG ZHUANG

for all n ∈ Zd . Let φ, ψ1 , . . . , ψs ∈ L2 (Rd ). We say that {φ; ψ1 , . . . , ψs } is a (nonhomogeneous) tight
framelet in L2 (Rd ) if
kf k2L2 (Rd )

=

X

k∈Zd

2

|hf, φ(· − k)i| +

∞ X
s X
X

j=0 ℓ=1 k∈Zd

|hf, 2jd/2ψℓ (2j · −k)i|2 ,

∀ f ∈ L2 (Rd ).

(3)

Let {a; b1 , . . . , bs } be a (d-dimensional
dyadic) tight framelet filter bank with all filters a, b1 , . . . , bs ∈
P
d
l0 (Z ) such that b
a(0) = k∈Zd a(k) = 1. Then we can define compactly supported tempered distributions φ and ψ1 , . . . , ψs on Rd through
b :=
φ(ξ)

∞
Y
j=1

b
a(2−j ξ),

ξ ∈ Rd

b
and ψbℓ (ξ) = bbℓ (ξ/2)φ(ξ/2),

ξ ∈ Rd , ℓ = 1, . . . , s.

(4)

Then we must have φ, ψ1 , . . . , ψs ∈ L2 (Rd ) and {φ; ψ1 , . . . , ψs } is a (d-dimensional dyadic) tight
framelet in L2 (Rd ) (see [12, Lemma 2.1], [10, Corollary 12 and Theorem 17] and [5, 18]). See
[3, 4, 5, 10, 12, 14, 16, 18] and many references therein for extensive investigation on tight framelets
derived from refinable functions. The tempered distribution φ in (4) is called a refinable function
b =b
b
satisfying the refinement equation φ(ξ)
a(ξ/2)φ(ξ/2)
for ξ ∈ Rd with the refinement filter a.
This paper is motivated by the interesting paper [17], where a two-dimensional directional Haar
tight framelet has been constructed and applied with impressive performance to pMRI. Applying
finite linear combinations to the standard tensor product two-dimensional Haar wavelet, in [17, (3.5)],
the authors constructed a two-dimensional directional tight framelet filter bank {aH ; b1 , . . . , b6 }, where
aH is the two-dimensional Haar low-pass filter and



1
1
1



if
k
=
(0,
0),
if
k
=
(1,
0),
if k = (0, 0),
4
4
4
1
1
1
b1 (k) = − 4 if k = (1, 1),
b2 (k) = − 4 if k = (0, 1),
b3 (k) = − 4 if k = (0, 1),



0

0
otherwise,
0
otherwise,
otherwise,



1
1
1



if k = (0, 0),
if k = (1, 0),
if k = (0, 1),
4
4
4
1
1
1
b5 (k) = − 4 if k = (1, 1),
b6 (k) = − 4 if k = (1, 1),
b4 (k) = − 4 if k = (1, 0),




0
0
otherwise,
0
otherwise,
otherwise.

Here the d-dimensional Haar low-pass filter is given by
(
2−d if k ∈ {0, 1}d ,
H
a (k) =
0
otherwise.

(5)

Since each high-pass filter has only two nonzero coefficients with opposite signs, it naturally has directionality and very simple structures. To deal with problems such as video processing in dimensions
higher than two, it is very natural to ask
Q1. Is it possible to construct a directional Haar tight framelet for every dimension such that each
high-pass filter has only two nonzero coefficients with opposite signs?
A similar construction method/argument as in [17] will quickly run into difficulty, since there are
so many possible linear combinations even at the dimension three. Fortunately, adopting a geometric
viewpoint, we can positively answer the question Q1 completely as follows:
Theorem 1. Let aH be the d-dimensional Haar low-pass filter in (5). Define the high-pass filters
d
b1 , . . . , bs with s := 22 = 2d−1 (2d − 1) in the following way: 2−d (δ γ1 − δ γ2 ) for all undirected edges
with endpoints γ1 , γ2 ∈ {0, 1}d and γ1 6= γ2 . Then {aH ; b1 , . . . , bs } is a tight framelet filter bank such
that all the high-pass filters b1 , . . . , bs have directionality and exhibit totally 21 (3d − 1) directions in
dimension d. Define functions φ and ψ1 , . . . , ψs as in (4). Then {φ; ψ1 , . . . , ψs } is a d-dimensional
directional compactly supported Haar tight framelet in L2 (Rd ).

DIRECTIONAL BOX SPLINE TIGHT FRAMELETS

3

Proof. To prove that {aH ; b1 , . . . , bs } is a tight framelet filter bank, we have to check the conditions
in (2) for all γ ∈ {0, 1}d and n ∈ Zd . Since all the filters are supported inside {0, 1}d, it is trivial to
observe that all the filters aH , b1 , . . . , bs vanish at the position γ + 2k for all k ∈ Zd \{0}. Therefore,
the equations (2) become
aH (γ)aH (n + γ) +

s
X
ℓ=1

bℓ (γ)bℓ (n + γ) = 2−d δ(n),

γ ∈ {0, 1}d, n ∈ Zd .

(6)

Case 1: n = 0. Then |aH (γ)|2 = 2−2d . By the definition of thePhigh-pass filters, there are totally
2 − 1 filters whose P
supports contain the point γ. Consequently, sℓ=1 |bℓ (γ)|2 = 2−2d (2d − 1). Thus,
we have |aH (γ)|2 + sℓ=1 |bℓ (γ)|2 = 2−2d + 2−2d (2d − 1) = 2−d which proves (6) with n = 0.
Case 2: n 6= 0 and n + γ 6∈ {0, 1}d. For this case, since all the filters are supported inside {0, 1}d,
we trivially have aH (n + γ) = 0 and bℓ (n + γ) = 0 for all ℓ = 1, . . . , s. Hence, (6) is trivially true for
n 6= 0 and n + γ 6∈ {0, 1}d.
Case 3: n 6= 0 and n + γ ∈ {0, 1}d. Then n and n + γ are two distinct points in {0, 1}d . By
the definition of the high-pass filters, there exists exactly one integer j with 1 6 j 6 s such that
bj (γ)bj (n + γ) = −2−2d and bℓ (γ)bℓ (n + γ) = 0 for all ℓ ∈ {1, . . . , s}\{j}. Noting that aH (γ) =
aH (n + γ) = 2−d , we conclude
d

H

a

(γ)aH (n

+ γ) +

s
X
ℓ=1

bℓ (γ)bℓ (n + γ) = aH (γ)aH (n + γ) + bj (γ)bj (n + γ) = 2−2d − 2−2d = 0,

which proves (6) for n 6= 0 and n + γ ∈ {0, 1}.
Therefore, {aH ; b1 , . . . , bs } is a tight framelet filter bank. Since each high-pass filter has only two
nonzero coefficients with opposite signs, all the high-pass filters b1 , . . . , bs trivially have directionality.
We now count the total number of directions of all the high-pass filters. Note that the direction of
a high-pass filter 2−d (δ γ1 − δ γ2 ) can be represented by the vector v = γ1 − γ2 or v = γ2 − γ1 . Such
a direction vector v is unique if we additionally require that the first nonzero entry of v should be
positive. Note that v ∈ {−1, 0, 1}d \{0} with each entry of v belonging to {−1, 0, 1}. Let S be the
set of all the nonzero vectors v ∈ {−1, 0, 1}d such that the first nonzero entry of v is positive (i.e.,
1). Hence, any direction vector v of a high-pass filter belongs to S. Conversely, for every vector
v ∈ S, we can uniquely write v = γ1 − γ2 with γ1 , γ2 ∈ {0, 1}d and γ1 + γ2 ∈ {0, 1}d by separating the
positive and negative entries of v. Therefore, the vector v represents the direction of the high-pass
filter 2−d (δ γ1 − δ γ2 ). Hence, the total number of directions of all the high-pass filters is equal to the
cardinality of the set S. Consider the subset Sj whose elements are in S with the first nonzero entry
at the position j for j = 1, . . . , d. Clearly, the cardinality of Sj is 3d−j . Since S is the disjoint union
of S1 , . . . , Sd , we conclude that the cardinality of S is 3d−1 + 3d−2 + · · · + 3d−d = (3d − 1)/2.

The tight framelet filter bank in Theorem 1 with d = 1 is just the standard Haar orthogonal
wavelet filter bank and the case d = 2 recovers the directional Haar tight framelet filter bank
in [17]. We now provide an alternative algebraic proof to Theorem 1 from the viewpoint of discrete framelet transforms. Because all the filters in Theorem 1 are supported inside {0, 1}d , the
d-dimensional discrete framelet transform using the filter bank {a; b1 , . . . , bs } in Theorem 1 can be
simply implemented by applying the discrete framelet transform acting on data supported on each
disjoint {0, 1}d + 2k, k ∈ Zd , where {0, 1}d is the set of all vertices of the unit cube [0, 1]d . For
simplicity, we list the vertices as {v1 , . . . , v2d } = {0, 1}d and assume that the data value at the
point vj is xj ∈ R. Then all the high-pass filters in Theorem 1 are given by ±2−d (δ vj − δ vk ) with
1 6 j < k 6 2d . The framelet coefficient produced by this high-pass filter is simply ±2d/2 2−d (xj − xk )
(see [9] for discrete framelet transforms). The coefficient produced by the Haar low-pass filter aH
in Theorem 1 is 2d/2 2−d (x1 + · · · + x2d ). Hence, the total squared energy of all the framelet coeffi2 P
2
cients is 2d/2 2−d (x1 + · · · + x2d ) + 16j<k62d ±2d/2 2−d (xj − xk ) . Noting that (x1 + · · ·+ x2d )2 =

4

BIN HAN, TAO LI, AND XIAOSHENG ZHUANG

P
x21 + · · · + x22d + 16j<k62d 2xj xk and (xj − xk )2 = (x2j + x2k ) − 2xj xk , we conclude that the total
squared energy of all the framelet coefficients is
2

2
X 
2d/2 2−d (x1 + · · · + x2d ) +
± 2d/2 2−d (xj − xk )
16j<k62d

= 2−d (x21 + · · · + x22d ) + 2−d

X

16j<k62d

(x2j + x2k ) = x21 + · · · + x22d ,

which proves the energy preservation property of the discrete framelet transform. By [9, Theorem 2.4],
the filter bank {a; b1 , . . . , bs } in Theorem 1 must be a d-dimensional tight framelet filter bank.
One possible shortcoming of Theorem 1 is that all the Haar refinable functions χ[0,1]d are discontinuous and often introduce unpleasant block effects in data/image processing. Therefore, smooth
refinable functions and directional smooth tight framelets are often preferred in applications. Refinable box splines can be made arbitrarily smooth and are widely used in approximation theory and
wavelet analysis. This naturally leads us to ask
Q2. Can we construct directional compactly supported tight framelets in L2 (Rd ) with simple
structure from every refinable box spline in every dimension?
Using Theorem 1 and the projection method, we positively answer the question Q2 painlessly. To
do so, let us recall the definition of box splines which are closely linked to the projection method.
Let P be a d × n real-valued matrix of rank d with d 6 n. A box spline MP with the d × n direction
matrix P is defined to be
Y 1 − e−ik·ξ
d
MP (ξ) :=
,
ξ ∈ Rd ,
(7)
ik · ξ
k∈P

where k ∈ P means that k is a column vector of P and k goes through all the columns of P once
and only once. A box spline can be also defined through the projection method ([11, 8]). For an
integrable function f ∈ L1 (Rn ), we can define the projected function P f on Rd by
Pcf (ξ) := fb(P T ξ),

ξ ∈ Rd .

(8)

Since fb is continuous on Rn , the function Pcf is a well-defined continuous function on Rd . In the
spatial domain, the definition of the d-dimensional projected function P f in (8) can be equivalently
expressed as
Z
1
[P f ](x) = p
f dS,
x ∈ Rd ,
(9)
T
−1
det(P P ) P x
where S is the surface element on the superplane P −1 x := {y ∈ Rn : P y = x}. In fact, for
Q
1−e−ik·ξ
f ∈ L1 (Rn ), the projected function P f ∈ L1 (Rd ) (see [8, 11]). Note that χ
\
[0,1]n (ξ) =
k∈{0,1}n
ik·ξ
and
Y 1 − e−i(P k)·ξ
Y 1 − e−ik·(P Tξ)
T
d
\
=
=M
\
P χ[0,1]n (ξ) = χ
P (ξ).
[0,1]n (P ξ) =
T ξ)
ik
·
(P
i(P
k)
·
ξ
n
n
k∈{0,1}

k∈{0,1}

Hence, the box spline MP is nothing else but the projected function P χ[0,1]n of the n-dimensional
Haar function along the direction matrix P . See the book [1] for extensive study on box splines.
The projection method can be also applied to filters on Zn provided that P is a d × n integer
matrix. For an n-dimensional filter a ∈ l0 (Zn ), the projected filter P a ∈ l0 (Zd ) is defined by
X
Pca(ξ) := b
a(P Tξ),
ξ ∈ Rd , or equivalently, [P a](j) =
a(k),
j ∈ Zd ,
(10)
k∈P −1 j

where P −1j := {k ∈ Zn : P k = j}. Because b
a is a 2πZn -periodic trigonometric polynomial and
P is an integer matrix, Pca is a well-defined 2πZd -periodic trigonometric polynomial. If aH is the
n-dimensional Haar low-pass filter in (5), then we define aP := P aH to be the box spline refinement
filter/mask for the box spline MP with a d × n direction matrix. For an integer projection matrix

DIRECTIONAL BOX SPLINE TIGHT FRAMELETS

5

d
d
P , the box spline function MP in (7) is refinable: M
c
P (2ξ) = a
P (ξ)MP (ξ), since P χ[0,1]n = MP and
c
H
χ[0,1]n (ξ). Various compactly
\
the n-dimensional Haar function χ[0,1]n is refinable: χ
[0,1]n (2ξ) = a (ξ)\
supported tight framelets have been constructed from refinable box splines with integer direction
matrices in the literature, e.g., see Ron and Shen [19], Lai and Stöckler [16], Han [8, 11], and Fan, Ji
and Shen [6].
Suppose that P is a d × n integer matrix of rank d with d 6 n satisfying
P T(Zd \[2Zd ]) ⊆ Zn \[2Zn ].

(11)

For every tight framelet {φ; ψ1 , . . . , ψs } in L2 (Rn ) regardless whether {φ; ψ1 , . . . , ψs } has an associated
underlying filter bank or not, it is known in [11, Theorem 4] and [8, Theorem 2.3 and Corollary 5.3]
that {P φ; P ψ1, . . . , P ψs } must be a tight framelet in L2 (Rd ). Similarly, for every n-dimensional tight
framelet filter bank {a; b1 , . . . , bs }, then {P a; P b1, . . . , P bs } must be a d-dimensional tight framelet
filter bank. For completeness, let us recall the argument from [8, 11] here. By definition, {a; b1 , . . . , bs }
is an n-dimensional tight framelet filter bank if and only if
s
X
b
a(ζ)b
∀ ζ ∈ Rn , β ∈ {0, 1}n .
(12)
bbℓ (ζ)bbℓ (ζ + πβ) = δ(β),
a(ζ + πβ) +
ℓ=1

The condition in (11) is equivalent to saying that P Tω 6∈ 2Zn for all ω ∈ {0, 1}d\{0}, Consequently,
for every ω ∈ {0, 1}d\{0}, we must have P T ω ∈ Zn \[2Zn ]. Therefore, it is trivial to deduce from (12)
with ζ = P Tξ and β = P T ω that
s
X
T
T
T
bbℓ (P Tξ)bbℓ (P T ξ + πP T ω) = δ(ω),
∀ ξ ∈ Rd , ω ∈ {0, 1}d.
a(P ξ + πP ω) +
b
a(P ξ)b
ℓ=1

By the definition of the projected filters in (10), the above equations simply become
s
X
d
d
c
c
P
bℓ (ξ)P
P a(ξ)P a(ξ + πω) +
bℓ (ξ + πω) = δ(ω),
∀ ξ ∈ Rd , ω ∈ {0, 1}d .
ℓ=1

That is, {P a; P b1 , . . . , P bs } must be a d-dimensional tight framelet filter bank.
The condition in (11) is P
equivalent to saying ([8, Theorem 2.5]) that the filter aP has the sum rules
of order at least one, i.e., k∈Zd aP (γ + 2k) = 2−d for all γ ∈ {0, 1}d. If the condition in (11) fails,
then the box spline filter aP does not have any sum rules and therefore, no tight framelets can be
ever derived from the box spline MP , see [8, Theorem 2.5] for more detail. The projection method
is originally introduced in Han [13] to study optimal smoothness of interpolating refinable functions
and has been further developed in [12, 11] and other papers to study refinable functions, biorthogonal
wavelets, and framelets. See Han [11, 8] for some applications of the projection method in wavelet
analysis.
Note that if a high-pass filter b has only two nonzero coefficients with opposite signs, then either
P b = 0 or P b has only two nonzero coefficients with opposite signs. Applying the projection method
to the Haar tight framelets in Theorem 1, we have the following result positively answering Q2.
Theorem 2. Let P be a d ×
matrix of rank d with d 6 n such that (11) holds. Let
 n integer
2n
n−1 n
H
{a ; b1 , . . . , bs } with s := 2 = 2 (2 − 1) be the n-dimensional Haar tight framelet filter bank
constructed in Theorem 1. Then {P aH ; P b1 , . . . , P bs } is a d-dimensional tight framelet filter bank
with P aH being the box spline refinement filter aP such that all the high-pass filters have only two
nonzero coefficients with opposite signs. Define φ and ψ1 , . . . , ψs as in (4) with a and b1 , . . . , bs being
replaced by P aH and P b1 , . . . , P bs , respectively. Then {φ; ψ1 , . . . , ψs } is a d-dimensional directional
compactly supported tight framelet in L2 (Rd ) with φ being the box spline MP in (7) having the direction
matrix P .
Let {a; b1 , . . . , bs } be a d-dimensional tight framelet filter bank.
p If b1 = c1 b and b2 = c2 b for some
d
constants c1 , c2 ∈ C and b ∈ l0 (Z ), then it is trivial that {a; |c1 |2 + |c2 |2 b, b3 , . . . , bs } is a tight

6

BIN HAN, TAO LI, AND XIAOSHENG ZHUANG

framelet filter bank. That is, we can combine high-pass filters which are almost the same up to a
multiplicative constant. Hence, the number of high-pass filters in Theorem 2 can be reduced. Now
we provide a geometric construction for the box spline tight framelet filter bank in Theorem 2 but
with similar filters combined to reduce the number of filters. First we calculate the support supp(aP )
of the box spline filter aP , which must be the set P {0, 1}n ⊆ Zd . Now the set {0, 1}n of the vertices
of the unit cube [0, 1]n can be written as a disjoint union of the subsets P −1 k, k ∈ supp(aP ). Then
all the high-pass filters are constructed in the following
p way: For every pair of two distinct points
γ1 , γ2 ∈ supp(aP ), construct the high-pass filter 2−n (#P −1 γ1 )(#P −1 γ2 )(δ γ1 − δ γ2 ), where #S is
the cardinality of a set S. Clearly, there are totally m2 with m := #supp(aP ) number of high-pass
filters. We now provide two examples to illustrate this construction.
Example 1. Let P be the 2 × 3 integer matrix


1 0 −1
P =
.
0 1 −1

Then MP is the three-direction interpolating linear box spline with the refinement filter given by
−3
−iξ1
ac
)(1 + e−iξ2 )(1 + ei(ξ1 +ξ2 ) ).
P (ξ1 , ξ2 ) = 2 (1 + e

Note that supp(aP ) = {−1, 0, 1}2\{(−1, 1)T , (1, −1)T} with #supp(aP ) = 7 and P −1 (0, 0)T =
3
T
{(0, 0, 0)T, (1, 1, 1)T}, while P −1 γ contains
 only one point in Z for every γ ∈ supp(aP )\{(0, 0) }.
7
Consequently, there are totally 21 (by 2 = 21) high-pass filters given by
√
√
√
√
√
√
2
2
2
2
2
2
(δ (1,0) −δ),
(δ (1,1) −δ),
(δ (0,1) −δ),
(δ (−1,0) −δ),
(δ (−1,−1) −δ),
(δ (0,−1) −δ)
8
8
8
8
8
8
and all the rest 15 filters are given by choosing a pair of two distinct points from supp(aP )\{(0, 0)T}
(i.e., {(1, 0)T, (1, 1)T , (0, 1)T, (−1, 0)T, (−1, −1)T , (0, −1)T}) with value 1/8 at one point and −1/8
at the other point. This tight framelet filter bank {aP ; b1 , . . . , b21 } exhibits totally 6 directions in
dimension two. See Figure 1 for details.
(0, 1)

y
•

•(1, 1)

•

•(0, 0)

•(1, 0)x

•

•(0, −1)

(−1, 0)

(−1, −1)

Figure 1. Each edge connecting every two vertices indicates a high-pass filter with
coefficients of the same weight but opposite signs at its two endpoints. Weight for each
√
dashed blue edge (total 15) is 81 , while weight for each solid red edge (total 6) is 82 .
The total number of all the high-pass filters is 21 with 6 solid red edges and 15 dashed
blue edges. The total number of directions/slopes for the 21 high-pass filters is 6 with
angles 0◦ (5 edges), 26.6◦ (= arctan( 21 ), 2 edges), 45◦ (5 edges), 63.4◦ (= arctan(2), 2
edges), 90◦ (5 edges), and −45◦ (2 edges).
Example 2. Let P be the 2 × 4 integer matrix


1 0 −1 0
P =
.
0 1 0 −1

DIRECTIONAL BOX SPLINE TIGHT FRAMELETS

7

Then MP is the tensor product of the piecewise linear B-spline with the refinement filter aP given
−2
−iξ1 2
by ac
| |1 + e−iξ2 |2 . Note that supp(aP ) = {−1, 0, 1}2 with #supp(aP ) = 9 and
P (ξ1 , ξ2 ) = 2 |1 + e
P −1(0, 0)T = {(0, 0, 0, 0)T, (1, 0, 1, 0)T, (0, 1, 0, 1)T, (1, 1, 1, 1)T},

P −1(1, 0)T = {(1, 0, 0, 0)T, (1, 1, 0, 1)T},

P −1 (0, 1)T = {(1, 0, 0, 0)T, (1, 1, 0, 1)T},

P −1 (−1, 0)T = {(0, 0, 1, 0)T, (0, 1, 1, 1)T},

P −1 (0, −1)T = {(0, 0, 0, 1)T, (1, 0, 1, 1)T},

and P −1γ contains only one point in Z4 for every γ ∈ {(1, −1)T , (−1, 1)T, (1, 1)T , (−1, −1)T }. Let
S1 := {(1, −1)T , (−1, 1)T , (1, 1)T,(−1, −1)T } and S2 := {(1, 0)T , (0, 1)T, (−1, 0)T, (0, −1)T }. Consequently, there are total 36 (by 92 = 36) high-pass filters given by
√
√
1
2
2
(δ γ − δ) ∀ γ ∈ S1 ,
(δ γ − δ) ∀ γ ∈ S2 ,
(δ γ1 − δ γ2 ) ∀ γ1 ∈ S1 , γ2 ∈ S2 ,
8
8
16
1
1
(δ γ1 − δ γ2 ) ∀ γ1 6= γ2 , γ1 , γ2 ∈ S1 ,
(δ γ − δ γ2 ) ∀ γ1 6= γ2 , γ1, γ2 ∈ S2
16
8 1
with 4, 4, 16, 6, and 6 filters for each group, respectively. This tight framelet filter bank {aP ; b1 , . . . , b36 }
exhibits totally 8 directions in dimension two. See Figure 2 for details.
y
(−1, 1)

•

(−1, 0)

(0, 1)

•

(0, 0)

•

•

•

•

(−1, −1)

(0, −1)

(1, 1)

•

•

(1, 0)

x

•

(1, −1)

Figure 2. Each edge connecting every two vertices indicates a high-pass filter with
coefficients of the same weight but opposite signs at its two endpoints. Weight for each
1
small-dotted
brown arc-edge (total 6) is 16
, weight for each dotted blue edge (total 16)
√
2
is 16 , weight for √each dashed green edge (total 10) is 18 , and weight for each solid red
edge (total 4) is 82 . The total number of directions/slopes for the 36 high-pass filters
is 8 with angles 0◦ (9 edges), ±26.6◦ (= arctan( 12 ), 2 edges each), ±45◦ (5 edges each),
±63.4◦ (= arctan( 21 ), 2 edges each), and 90◦ (9 edges).
Note that if {a; b1 , . . . , bs } is a tight framelet filter bank, then {a; b1 (· − 2k1 ), . . . , bs (· − 2ks )} is also
a tight framelet filter bank for all k1 , . . . , ks ∈ Zd . Using this simple observation, the number of the
high-pass filters in Example 2 can be reduced from 36 to 30. Indeed, the two filters corresponding
to the edges {(1, −1)T , (1, 1)T} and {(−1, −1)T , (−1, 1)T} can be combined into one single filter; the
two filters for the edges {(0, 1)T, (1, 1)T} and {(0, −1)T , (1, −1)T} can be combined into one filter;
the two filters for the edges {(−1, 1)T, (0, 1)T } and {(−1, −1)T , (0, −1)T} can be combined into one
filter. We can also perform the same operation for filters with the vertical direction. Using this and
more complicated argument/technique, the number of high-pass filters in a box spline tight framelet
filter bank in Theorem 2 often can be reduced. We shall not address this issue in this paper. We
shall explore some applications of the directional tight framelets in Theorems 1 and 2 elsewhere.

8

BIN HAN, TAO LI, AND XIAOSHENG ZHUANG

References
[1] C. de Boor, K. Höllig, S. Riemenschneider. Box splines. Series in Appl. Math. Sci. 98. Springer-Verlag, New York,
1993.
[2] E. J. Candès and D. L. Donoho, New tight frames of curvelets and optimal representations of objects with
piecewise C 2 singularities. Comm. Pure Appl. Math. 57 (2004), 219–266.
[3] C. K. Chui, W. He, and J. Stöckler, Compactly supported tight and sibling frames with maximum vanishing
moments. Appl. Comput. Harmon. Anal. 13 (2002), 224–262.
[4] I. Daubechies, A. Grossmann, and Y. Meyer, Painless nonorthogonal expansions. J. Math. Phys. 27 (1986),
1271–1283.
[5] I. Daubechies, B. Han, A. Ron, Z. Shen, Framelets: MRA-based constructions of wavelet frames. Appl. Comput.
Harmon. Anal. 14 (2003), 1–46.
[6] Z. Fan, H. Ji, and Z. Shen, Dual Gramian analysis: duality principle and unitary extension principle. Math. Comp.
85 (2016), 239–270.
[7] K. Guo, G. Kutyniok, and D. Labate, Sparse multidimensional representations using anisotropic dilation and
shear operators, Wavelets and Splines (Athens, GA, 2005), Nashboro Press, Nashville, TN (2006) 189-201.
[8] B. Han, The projection method for multidimensional framelet and wavelet analysis. Math. Model. Nat. Phenom.
9 (2014), 83–110.
[9] B. Han, Properties of discrete framelet transforms. Math. Model. Nat. Phenom. 8 (2013), 18–47.
[10] B. Han, Nonhomogeneous wavelet systems in high dimensions. Appl. Comput. Harmon. Anal. 32 (2012), 169–196.
[11] B. Han, The porjection method in wavelet analysis, in Splines and Wavelets: Athens 2005, G. Chen and M.J. Lai
eds., (2006), 202–225.
[12] B. Han, Compactly supported tight wavelet frames and orthonormal wavelets of exponential decay with a general
dilation matrix. J. Comput. Appl. Math. 155 (2003), 43–67.
[13] B. Han, Analysis and construction of optimal multivariate biorthogonal wavelets with compact support. SIAM J.
Math. Anal. 31 (1999), 274–304.
[14] B. Han, On dual wavelet tight frames. Appl. Comput. Harmon. Anal. 4 (1997), 380–413.
[15] B. Han and Z. Zhao, Tensor product complex tight framelets with increasing directionality. SIAM J. Imaging Sci.
7 (2014), 997–1034.
[16] M. J. Lai and J. Stöckler, Construction of multivariate compactly supported tight wavelet frames. Appl. Comput.
Harmon. Anal. 21 (2006), 324–348.
[17] Y.-R. Li, R. H. Chan, L. Shen, Y.-C. Hsu, and W.-Y. Tseng, An adaptive directional Haar framelet-based
reconstruction algorithm for parallel magnetic resonance imaging. SIAM J. Imaging Sci. 9 (2016), 794–821.
[18] A. Ron and Z. Shen, Affine systems in L2 (Rd ): the analysis of the analysis operator. J. Funct. Anal. 148 (1997),
408–447.
[19] A. Ron and Z. Shen, Compactly supported tight affine spline frames in L2 (Rd ). Math. Comp. 67 (1998), 191–207.
Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton,
berta, Canada T6G 2G1.
bhan@ualberta.ca, tao8@ualberta.ca

Al-

Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon Tong,
Hong Kong.
xzhuang7@cityu.edu.hk.

