comparison util/SMALL_ImgDeNoiseResult.m @ 77:62f20b91d870

add routines from sparco problems privite folder to {root}\util
some changes to ksvd vs rlsdla image denoising example

 m=size(SMALL.solver,2)+1;
 maxval=SMALL.Problem.maxval;
 im=SMALL.Problem.Original;
 imnoise=SMALL.Problem.Noisy;
 
-subplot(2, m, 1); imshow(im/maxval);
+subplot(2, m, 1); imagesc(im/maxval);colormap(gray);axis off; axis image;        % Set aspect ratio to obtain square pixels
 title('Original image');
 
-subplot(2,m,m+1); imshow(imnoise/maxval);
+subplot(2,m,m+1); imagesc(imnoise/maxval);axis off; axis image; 
 title(sprintf('Noisy image, PSNR = %.2fdB', SMALL.Problem.noisy_psnr ));
 
 for i=2:m
     
-    subplot(2, m, i); imagesc(SMALL.solver(i-1).reconstructed.Image/maxval);colormap(gray);
+    subplot(2, m, i); imagesc(SMALL.solver(i-1).reconstructed.Image/maxval);axis off; axis image; 
     title(sprintf('%s Denoised image, PSNR: %.2f dB in %.2f s',...
-        SMALL.DL(i-1).name, SMALL.solver(i-1).reconstructed.psnr, SMALL.solver(i-1).time ));
+        SMALL.DL(i-1).name, SMALL.solver(i-1).reconstructed.psnr, SMALL.solver(i-1).time ),'Interpreter','none');
     if strcmpi(SMALL.DL(i-1).name,'ksvds')
         D = kron(SMALL.Problem.basedict{2},SMALL.Problem.basedict{1})*SMALL.DL(i-1).D;
     else
         D = SMALL.DL(i-1).D;
     end
     dictimg = showdict(D,SMALL.Problem.blocksize,...
         round(sqrt(size(D,2))),round(sqrt(size(D,2))),'lines','highcontrast');
     
-    subplot(2,m,m+i);imagesc(dictimg);colormap(gray);
+    subplot(2,m,m+i);imagesc(dictimg);axis off; axis image; 
     title(sprintf('%s dictionary in %.2f s',...
-        SMALL.DL(i-1).name, SMALL.DL(i-1).time));
+        SMALL.DL(i-1).name, SMALL.DL(i-1).time),'Interpreter','none');
     
 end