Mercurial > hg > smallbox
annotate Problems/Pierre_reconstruct.m @ 245:96d17e5dc5d3
Replaced the call to SMALL_ImgDenoise_DL_test_KSVDvsSPAMS with SMALL_ImgDenoise_DL_test_KSVDvsRLSDLAvsTwoStepMOD in Example 2 of SMALLBoxSetup.m.
author | luisf <luis.figueira@eecs.qmul.ac.uk> |
---|---|
date | Wed, 31 Oct 2012 11:53:07 +0000 |
parents | 8e660fd14774 |
children |
rev | line source |
---|---|
idamnjanovic@46 | 1 function reconstructed=Pierre_reconstruct(y, Problem) |
ivan@128 | 2 %% Pierre Villars Example - reconstruction function |
ivan@128 | 3 % |
ivan@128 | 4 % This example is based on the experiment suggested by Professor Pierre |
ivan@128 | 5 % Vandergheynst on the SMALL meeting in Villars. |
ivan@128 | 6 % The function is using sparse representation y in dictionary Problem.A |
ivan@128 | 7 % to reconstruct the patches of the target image. |
ivan@128 | 8 |
idamnjanovic@46 | 9 % |
idamnjanovic@46 | 10 % Centre for Digital Music, Queen Mary, University of London. |
idamnjanovic@46 | 11 % This file copyright 2009 Ivan Damnjanovic. |
idamnjanovic@46 | 12 % |
idamnjanovic@46 | 13 % This program is free software; you can redistribute it and/or |
idamnjanovic@46 | 14 % modify it under the terms of the GNU General Public License as |
idamnjanovic@46 | 15 % published by the Free Software Foundation; either version 2 of the |
idamnjanovic@46 | 16 % License, or (at your option) any later version. See the file |
idamnjanovic@46 | 17 % COPYING included with this distribution for more information. |
ivan@128 | 18 %% |
idamnjanovic@46 | 19 |
ivan@128 | 20 |
ivan@128 | 21 |
idamnjanovic@46 | 22 |
idamnjanovic@46 | 23 imout=Problem.A*y; |
idamnjanovic@46 | 24 |
idamnjanovic@46 | 25 % combine the patches into reconstructed image |
idamnjanovic@46 | 26 |
idamnjanovic@46 | 27 im=col2im(imout,Problem.blocksize,size(Problem.imageTrg),'disctint'); |
idamnjanovic@46 | 28 |
idamnjanovic@46 | 29 % bound the pixel values to [0,255] range |
idamnjanovic@46 | 30 im(im<0)=0; |
idamnjanovic@46 | 31 im(im>255)=255; |
idamnjanovic@46 | 32 |
idamnjanovic@46 | 33 %% output structure image+psnr %% |
idamnjanovic@46 | 34 reconstructed.image=im; |
idamnjanovic@46 | 35 reconstructed.psnr = 20*log10(Problem.maxval * sqrt(numel(Problem.imageTrg(:))) / norm(Problem.imageTrg(:)-im(:))); |
idamnjanovic@46 | 36 end |