A Multiresolution Image Completion Algorithm for Compressing Digital Color Images

<table><tr><td><table class="algorithm" id="alg2"><tr><td colspan="2">(1) Decompose <svg height="21.25" id="M47" style="vertical-align:-4.15506pt" version="1.1" viewbox="0 0 29.387501 21.25" width="29.387501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,16.013)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
<g transform="matrix(.012,-0,0,-.012,10.975,21.05)"><path d="M483 97q-42 -50 -88.5 -79.5t-68.5 -29.5q-37 0 -17 93l22 102h-2q-54 -79 -144 -149q-59 -46 -100 -46q-24 0 -43 29t-19 86q0 78 34.5 153t94.5 120q41 31 94 51.5t98 20.5q29 0 72 -9q26 -6 39 -6l2 -4q-30 -117 -67 -323q-8 -41 2 -41q16 0 79 58zM374 387
q-32 15 -73 15q-52 0 -83 -23q-48 -36 -78 -108.5t-30 -152.5q0 -33 8.5 -50.5t20.5 -17.5q31 0 107 79t99 132q15 40 29 126z" id="x1D44E"></path></g><g transform="matrix(.012,-0,0,-.012,10.975,7.85)"><path d="M394 607q0 -25 -15.5 -43t-39.5 -18q-19 0 -32 13t-13 35t16 41.5t40 19.5q20 0 32 -13.5t12 -34.5zM360 349l-57 -283q-38 -188 -122 -263q-72 -64 -114 -64q-16 0 -30 10.5t-14 22.5q0 24 21 45q15 13 26 6q30 -20 54 -20q41 0 65 95q37 147 86 415q6 34 6.5 42.5
t-5.5 11.5q-18 0 -79 -54l-14 25q48 48 96.5 79t71.5 31q29 0 9 -99z" id="x1D457"></path></g><g transform="matrix(.012,-0,0,-.012,16.008,7.85)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.012,-0,0,-.012,22.993,7.85)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g>
</svg> into a low pass image <svg height="21.25" id="M48" style="vertical-align:-4.15506pt" version="1.1" viewbox="0 0 17.6875 21.25" width="17.6875" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,16.013)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
<g transform="matrix(.012,-0,0,-.012,10.975,7.85)"><path d="M394 607q0 -25 -15.5 -43t-39.5 -18q-19 0 -32 13t-13 35t16 41.5t40 19.5q20 0 32 -13.5t12 -34.5zM360 349l-57 -283q-38 -188 -122 -263q-72 -64 -114 -64q-16 0 -30 10.5t-14 22.5q0 24 21 45q15 13 26 6q30 -20 54 -20q41 0 65 95q37 147 86 415q6 34 6.5 42.5
t-5.5 11.5q-18 0 -79 -54l-14 25q48 48 96.5 79t71.5 31q29 0 9 -99z" id="x1D457"></path></g><g transform="matrix(.012,-0,0,-.012,10.975,21.05)"><path d="M483 97q-42 -50 -88.5 -79.5t-68.5 -29.5q-37 0 -17 93l22 102h-2q-54 -79 -144 -149q-59 -46 -100 -46q-24 0 -43 29t-19 86q0 78 34.5 153t94.5 120q41 31 94 51.5t98 20.5q29 0 72 -9q26 -6 39 -6l2 -4q-30 -117 -67 -323q-8 -41 2 -41q16 0 79 58zM374 387
q-32 15 -73 15q-52 0 -83 -23q-48 -36 -78 -108.5t-30 -152.5q0 -33 8.5 -50.5t20.5 -17.5q31 0 107 79t99 132q15 40 29 126z" id="x1D44E"></path></g>
</svg> and a high pass image <svg height="23.0375" id="M49" style="vertical-align:-4.74467pt" version="1.1" viewbox="0 0 18.262501 23.0375" width="18.262501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.062)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
<g transform="matrix(.012,-0,0,-.012,10.975,7.85)"><path d="M394 607q0 -25 -15.5 -43t-39.5 -18q-19 0 -32 13t-13 35t16 41.5t40 19.5q20 0 32 -13.5t12 -34.5zM360 349l-57 -283q-38 -188 -122 -263q-72 -64 -114 -64q-16 0 -30 10.5t-14 22.5q0 24 21 45q15 13 26 6q30 -20 54 -20q41 0 65 95q37 147 86 415q6 34 6.5 42.5
t-5.5 11.5q-18 0 -79 -54l-14 25q48 48 96.5 79t71.5 31q29 0 9 -99z" id="x1D457"></path></g><g transform="matrix(.012,-0,0,-.012,10.975,22.838)"><path d="M530 686q-79 -330 -126 -595q-5 -24 4 -24q22 0 77 53l16 -24q-42 -49 -88 -78.5t-72 -29.5q-37 0 -19 83l21 99h-2q-62 -83 -153 -141q-65 -41 -97 -41q-28 0 -48 31.5t-20 91.5q0 71 32.5 144t88.5 118q36 30 94 52.5t94 22.5q38 0 68 -13l29 158q7 40 2 47.5t-38 7.5
h-35l1 26q39 4 76 13.5t58.5 17t27.5 7.5q16 0 9 -26zM387 375q-11 11 -38.5 20t-51.5 9q-45 0 -79 -23q-46 -32 -77.5 -107.5t-31.5 -143.5q0 -40 9.5 -58.5t23.5 -18.5q40 0 113.5 70t99.5 118z" id="x1D451"></path></g>
</svg> by </td></tr><tr><td colspan="2">  applying Laplacian Pyramid Scheme.</td></tr><tr><td colspan="2">(2) Use pseudo polar grid to compute <svg height="23.0375" id="M50" style="vertical-align:-4.74467pt" version="1.1" viewbox="0 0 40.512501 23.0375" width="40.512501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.062)"><path d="M619 482q0 -69 -40.5 -119t-96 -72.5t-118.5 -26.5h-44l-70 20l-31 -151q-14 -67 0.5 -83t89.5 -22l-5 -28h-287l8 28q65 7 81.5 22t29.5 83l79 398q12 56 0.5 70.5t-78.5 20.5l7 28h255q108 0 164 -43t56 -125zM524 478q0 141 -146 141q-25 0 -47 -8q-16 -6 -20.5 -13.5
t-10.5 -39.5l-43 -241q37 -13 83 -13q67 0 125.5 45t58.5 129z" id="x1D443"></path></g><g transform="matrix(.017,-0,0,-.017,14.647,17.062)"><path d="M170 255q0 -25 -16 -41t-40 -16t-40 16t-16 41t16 41.5t40 16.5t40 -16.5t16 -41.5z" id="x22C5"></path></g><g transform="matrix(.017,-0,0,-.017,22.297,17.062)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
<g transform="matrix(.012,-0,0,-.012,33.212,7.85)"><path d="M394 607q0 -25 -15.5 -43t-39.5 -18q-19 0 -32 13t-13 35t16 41.5t40 19.5q20 0 32 -13.5t12 -34.5zM360 349l-57 -283q-38 -188 -122 -263q-72 -64 -114 -64q-16 0 -30 10.5t-14 22.5q0 24 21 45q15 13 26 6q30 -20 54 -20q41 0 65 95q37 147 86 415q6 34 6.5 42.5
t-5.5 11.5q-18 0 -79 -54l-14 25q48 48 96.5 79t71.5 31q29 0 9 -99z" id="x1D457"></path></g><g transform="matrix(.012,-0,0,-.012,33.212,22.838)"><path d="M530 686q-79 -330 -126 -595q-5 -24 4 -24q22 0 77 53l16 -24q-42 -49 -88 -78.5t-72 -29.5q-37 0 -19 83l21 99h-2q-62 -83 -153 -141q-65 -41 -97 -41q-28 0 -48 31.5t-20 91.5q0 71 32.5 144t88.5 118q36 30 94 52.5t94 22.5q38 0 68 -13l29 158q7 40 2 47.5t-38 7.5
h-35l1 26q39 4 76 13.5t58.5 17t27.5 7.5q16 0 9 -26zM387 375q-11 11 -38.5 20t-51.5 9q-45 0 -79 -23q-46 -32 -77.5 -107.5t-31.5 -143.5q0 -40 9.5 -58.5t23.5 -18.5q40 0 113.5 70t99.5 118z" id="x1D451"></path></g>
</svg><i>.  </i></td></tr><tr><td colspan="2">(3) Band pass filter the output matrix <svg height="23.0375" id="M51" style="vertical-align:-4.74467pt" version="1.1" viewbox="0 0 40.512501 23.0375" width="40.512501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.062)"><path d="M619 482q0 -69 -40.5 -119t-96 -72.5t-118.5 -26.5h-44l-70 20l-31 -151q-14 -67 0.5 -83t89.5 -22l-5 -28h-287l8 28q65 7 81.5 22t29.5 83l79 398q12 56 0.5 70.5t-78.5 20.5l7 28h255q108 0 164 -43t56 -125zM524 478q0 141 -146 141q-25 0 -47 -8q-16 -6 -20.5 -13.5
t-10.5 -39.5l-43 -241q37 -13 83 -13q67 0 125.5 45t58.5 129z" id="x1D443"></path></g><g transform="matrix(.017,-0,0,-.017,14.647,17.062)"><path d="M170 255q0 -25 -16 -41t-40 -16t-40 16t-16 41t16 41.5t40 16.5t40 -16.5t16 -41.5z" id="x22C5"></path></g><g transform="matrix(.017,-0,0,-.017,22.297,17.062)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
<g transform="matrix(.012,-0,0,-.012,33.212,7.85)"><path d="M394 607q0 -25 -15.5 -43t-39.5 -18q-19 0 -32 13t-13 35t16 41.5t40 19.5q20 0 32 -13.5t12 -34.5zM360 349l-57 -283q-38 -188 -122 -263q-72 -64 -114 -64q-16 0 -30 10.5t-14 22.5q0 24 21 45q15 13 26 6q30 -20 54 -20q41 0 65 95q37 147 86 415q6 34 6.5 42.5
t-5.5 11.5q-18 0 -79 -54l-14 25q48 48 96.5 79t71.5 31q29 0 9 -99z" id="x1D457"></path></g><g transform="matrix(.012,-0,0,-.012,33.212,22.838)"><path d="M530 686q-79 -330 -126 -595q-5 -24 4 -24q22 0 77 53l16 -24q-42 -49 -88 -78.5t-72 -29.5q-37 0 -19 83l21 99h-2q-62 -83 -153 -141q-65 -41 -97 -41q-28 0 -48 31.5t-20 91.5q0 71 32.5 144t88.5 118q36 30 94 52.5t94 22.5q38 0 68 -13l29 158q7 40 2 47.5t-38 7.5
h-35l1 26q39 4 76 13.5t58.5 17t27.5 7.5q16 0 9 -26zM387 375q-11 11 -38.5 20t-51.5 9q-45 0 -79 -23q-46 -32 -77.5 -107.5t-31.5 -143.5q0 -40 9.5 -58.5t23.5 -18.5q40 0 113.5 70t99.5 118z" id="x1D451"></path></g>
</svg></td></tr><tr><td colspan="2">(4) Re-assemble the Cartesian sampled values </td></tr><tr><td colspan="2">(5) Take the 2D IFFT.                                                     </td></tr></table></td></tr></table>

Construction of discrete shearlet transform.

Journal of Applied Mathematics

alg2

Algorithm 2

Algorithm 2: A Multiresolution Image Completion Algorithm for Compressing Digital Color Images 