sparse image representation via combined transforms - Convex ...

sparse image representation via combined transforms - Convex ... sparse image representation via combined transforms - Convex ...

convexoptimization.com
from convexoptimization.com More from this publisher
10.03.2015 Views

138 APPENDIX A. DIRECT EDGELET TRANSFORM (a) Huo image (b) Largest−200 coeff. (c) Largest−400 coeff. 10 20 30 40 50 60 20 40 60 10 20 30 40 50 60 20 40 60 10 20 30 40 50 60 20 40 60 (d) Edgelet trans. coeff.(largest−10000/26752) 1.4 1.2 1 0.8 0.6 0.4 0.2 0 0 5000 10000 10 20 30 40 50 60 (e) Largest−800 coeff. 20 40 60 10 20 30 40 50 60 (f) Largest−1600 coeff. 20 40 60 Figure A.1: Edgelet transform of the Chinese character “Huo”: (a) is the original; (d) is the sorted coefficients; (b), (c), (e) and (f) are reconstructions based on the largest 200, 400, 800, 1600 coefficients, respectively. The reason to do this is that both of them show some patchy patterns in the original images. The filtered images show more obvious edge features. See Figure A.2 (b) and Figure A.4 (b). A.3 Details . A direct way of computing edgelet coefficients is explained. It is direct in the sense that every single edgelet coefficient is computed by a direct method; the coefficient is a weighted sum of the pixel values (intensities) that the edgelet trespasses. For an N × N image, the order of the complexity of the direct method is O(N 3 log N). Section A.3.1 gives the definition of the direct edgelet transform. Section A.3.2 calculates the cardinality of the edgelet system. Section A.3.3 specifies the ordering. Section A.3.4

A.3. DETAILS 139 (a) Stick image (b) Filtered stick image (c) Edgelet trans. coeff.(largest−10000/362496) 5 20 20 4 40 60 40 60 3 80 80 2 100 120 20 40 60 80 100 120 100 120 20 40 60 80 100 120 1 0 0 5000 10000 (d) Edgelets assoc. with largest−100 coeff. (e) Largest−300 coeff. (f) Largest−500 coeff. 20 20 20 40 40 40 60 60 60 80 80 80 100 100 100 120 20 40 60 80 100 120 120 20 40 60 80 100 120 120 20 40 60 80 100 120 Figure A.2: Edgelet transform of the sticky image: (a) is the original; (b) is the filtered image; (c) is the sorted coefficients; (d), (e) and (f) are the reconstructions based on the largest 100, 300 and 500 coefficients, respectively. explains how to compute a single edgelet coefficient. Section A.3.5 shows how to do the adjoint transform. A.3.1 Definition In this subsection, we define what a direct edgelet transform is. For continuous functions, the edgelet transform is defined in [50]. But in modern computing practice, an image is digitalized: it is a discrete function, or a matrix. There are many possibilities to realize the edgelet transform for a matrix, due to many ways to interpolate a matrix as a continuous 2-D function. We explain a way of interpolation. Suppose the image is a squared image having N rows and N columns. The total number of pixels is N 2 . We assume N is a power of 2, N =2 n . Let the matrix I = I(i, j) 1≤i,j≤N ∈ R N×N correspond to the digital image,

A.3. DETAILS 139<br />

(a) Stick <strong>image</strong><br />

(b) Filtered stick <strong>image</strong><br />

(c) Edgelet trans. coeff.(largest−10000/362496)<br />

5<br />

20<br />

20<br />

4<br />

40<br />

60<br />

40<br />

60<br />

3<br />

80<br />

80<br />

2<br />

100<br />

120<br />

20 40 60 80 100 120<br />

100<br />

120<br />

20 40 60 80 100 120<br />

1<br />

0<br />

0 5000 10000<br />

(d) Edgelets assoc. with largest−100 coeff.<br />

(e) Largest−300 coeff.<br />

(f) Largest−500 coeff.<br />

20<br />

20<br />

20<br />

40<br />

40<br />

40<br />

60<br />

60<br />

60<br />

80<br />

80<br />

80<br />

100<br />

100<br />

100<br />

120<br />

20 40 60 80 100 120<br />

120<br />

20 40 60 80 100 120<br />

120<br />

20 40 60 80 100 120<br />

Figure A.2: Edgelet transform of the sticky <strong>image</strong>: (a) is the original; (b) is the filtered<br />

<strong>image</strong>; (c) is the sorted coefficients; (d), (e) and (f) are the reconstructions based on the<br />

largest 100, 300 and 500 coefficients, respectively.<br />

explains how to compute a single edgelet coefficient. Section A.3.5 shows how to do the<br />

adjoint transform.<br />

A.3.1<br />

Definition<br />

In this subsection, we define what a direct edgelet transform is.<br />

For continuous functions, the edgelet transform is defined in [50]. But in modern computing<br />

practice, an <strong>image</strong> is digitalized: it is a discrete function, or a matrix. There are<br />

many possibilities to realize the edgelet transform for a matrix, due to many ways to interpolate<br />

a matrix as a continuous 2-D function.<br />

We explain a way of interpolation. Suppose the <strong>image</strong> is a squared <strong>image</strong> having N<br />

rows and N columns. The total number of pixels is N 2 . We assume N is a power of<br />

2, N =2 n . Let the matrix I = I(i, j) 1≤i,j≤N ∈ R N×N correspond to the digital <strong>image</strong>,

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

Saved successfully!

Ooh no, something went wrong!