comparison Problems/generateImageDenoiseProblem.m @ 44:2c59257d734c

(none)

     cd([pathstr1,FS,'data',FS,'images']);
     [filename,pathname] = uigetfile({'*.png;'},'Select an image');
     [pathstr, name, ext, versn] = fileparts(filename);
     data.name=name;
     im = imread(filename);
-    im = double(im);
+    %im = double(im);
 end;
+im = double(im);
 cd(TMPpath);
 
 %% check input parameters %%
 
 if ~ exist( 'blocksize', 'var' ) || isempty(blocksize),blocksize = 8;end

 
 %% set parameters %%
 
 x = imnoise;
 p = ndims(x);
-
+psnr=20*log10(maxval * sqrt(numel(im)) / norm(im(:)-imnoise(:)));
 if (p==2 && any(size(x)==1) && length(blocksize)==1)
     p = 1;
 end
 
 % blocksize %

 for i = 1:bsize:size(X,2)
     blockids = i : min(i+bsize-1,size(X,2));
     X(:,blockids) = remove_dc(X(:,blockids),'columns');
 end
 
+%  Noisy image blocks 
+xcol=im2col(x,blocksize,'sliding');
+[b1, dc] = remove_dc(xcol,'columns');
+
 %% output structure %%
 
 data.Original = im;
 data.Noisy = imnoise;
+data.noisy_psnr=psnr;
 data.b = X;
+data.b1=b1;
+data.b1dc=dc;
 data.m = size(X,1);
 data.n = size(X,2);
 data.p = dictsize;
 data.blocksize=blocksize;
 data.sigma = sigma;