10.03.2015 Views

sparse image representation via combined transforms - Convex ...

sparse image representation via combined transforms - Convex ...

sparse image representation via combined transforms - Convex ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

7.2. MODIFYING EDGELET DICTIONARY 131<br />

7.2 Modifying Edgelet Dictionary<br />

By inspection of Figure 6.2, one sees that the existing edgelet dictionary can be further<br />

improved. Indeed, in those figures, one sees that the wavelet component of the reconstruction<br />

is carrying a significant amount of the burden of edge <strong>representation</strong>. It seems that<br />

the terms in our edgelet dictionary do not have a width matched to the edge width, and<br />

that in consequence, many fine-scale wavelets are needed in the <strong>representation</strong>. In future<br />

experiments, we would try using edgelet features that have a finer width at fine scales and<br />

coarser width at coarse scales. This intuitive discussion matches some of the concerns in<br />

the paper [54].<br />

The idea in [54] is to construct a new tight frame intended for efficient <strong>representation</strong><br />

of 2-D objects with singularities along curves. The frame elements exhibit a range of<br />

dyadic positions, scales and angular orientations. The useful frame elements are highly<br />

directionally selective at fine scales. Moreover, the width of the frame elements scales with<br />

length according to the square of length.<br />

The frame construction combines ideas from ridgelet and wavelet analysis. One tool is<br />

the monoscale ridgelet transform. The other tool is the use of multiresolution filter banks.<br />

The frame coefficients are obtained by first separating the object into special multiresolution<br />

subbands f s , s ≥ 0 by applying filtering operation ∆ s f =Ψ 2s ∗ f for s ≥ 0, where<br />

Ψ 2s is built from frequencies in an annulus extending from radius 2 2s to 2 2s+2 . To the s-th<br />

subband one applies the monoscale ridgelet transform at scale s. Note that the monoscale<br />

ridgelet transform at scale index s is composed with multiresolution filtering near index 2s.<br />

The (s, 2s) pairing makes the useful frame elements highly anisotropic at fine scales.<br />

The frame gives rise to a near-optimal atomic decomposition of objects which have<br />

discontinuities along a closed C 2 curve. Simple thresholding of frame coefficients gives rise<br />

to new methods of approximation and smoothing that are highly anisotropic and provably<br />

optimal.<br />

We refer readers to the original paper cited at the beginning of this section for more<br />

details.<br />

7.3 Accelerating the Iterative Algorithm<br />

The idea of block coordinate relaxation can be found in [124, 125]. A proof of the convergence<br />

can be found in [135]. We happen to know an independent work in [130]. The idea

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!