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.

4.3. MINIMUM l 1 NORM SOLUTION 85<br />

l 0 norm is a nonconvex function while l 1 norm is convex.<br />

Some inspiring applications of the minimum l 1 norm phenomenon are given in [55]:<br />

• Recovering missing segments of a bandlimited signal. Suppose s is a bandlimited<br />

signal with frequency components limited in B. Buts is missing in a time interval T .<br />

Given B, if|T | is small enough, then the minimizer of ‖s − s ′ ‖ 1 among all s ′ whose<br />

frequency components are limited in B is exactly s.<br />

• Recovery of a “<strong>sparse</strong>” wide-band signal from narrow-band measurements. Suppose s<br />

is a signal that is “<strong>sparse</strong>” in the time domain. As we know, a time-<strong>sparse</strong> signal must<br />

occupy a wide band in the frequency domain. Suppose N t ≥|s|. Suppose we can only<br />

observe a bandlimited fraction of s, r = P B s, here P B functioning like a bandpass<br />

filter. We can perfectly recover s by finding the minimizer of ‖r−P B s ′ ‖ 1 among all the<br />

s ′ satisfying |s ′ |≤N t . This phenomenon has application in several branches of applied<br />

science, where instrumental limitations make the available observation bandlimited.<br />

On the other hand, if we consider the symmetry between the time domain and the<br />

frequency domain, this is a dual of the previous case. We simply switch the positions<br />

of the time domain and the frequency domain.<br />

Another recent advance in theory [43] has given an interesting result. Suppose that the<br />

overcomplete dictionary we consider is a combination of two complementary orthonormal<br />

bases. By “complementary” we mean that the maximum absolute value of the inner product<br />

of any two elements (one from each basis) is upper bounded by a small value. Suppose the<br />

observed signal y is made by a small number of atoms in the dictionary. Solving the<br />

minimum l 1 norm problem will give us the same solution as solving the minimum l 0 norm<br />

problem. More specifically, if a dictionary is made by two bases—Dirac basis and Fourier<br />

basis—and if the observation y is made by fewer than √ N/2 atoms from the dictionary,<br />

where N is the size of the signal, then the minimum l 1 norm decomposition is the same as<br />

the minimum l 0 norm decomposition.<br />

‖x‖ 0 denotes the number of nonzero elements in the vector x. Thel 0 norm is generally<br />

regarded as the measure of sparsity. So the minimum l 0 norm decomposition is generally<br />

regarded as the <strong>sparse</strong>st decomposition. Note the minimum l 0 norm problem is a combinatorial<br />

problem and in general is NP hard. But the minimum l 1 norm problem is a convex<br />

optimization problem and can be solved by some polynomial-time optimization methods,

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

Saved successfully!

Ooh no, something went wrong!