THESE/codes/wave/allcode/main_dwt.m

   1 % Main function
   2 % Usage :
   3 %        main
   4 % INPUT :
   5 %        Raw Lena image
   6 % OUTPUT :
   7 %        PSNR value of the denoised image
   8 %
   9 % Load clean image
  10 fid = fopen('barbara','r');
  11 s  = fread(fid,[512 512],'unsigned char');
  12 fclose(fid)
  13 %s = imread('zoe.bmp','bmp');
  14 %s = double(s);
  15 N = 512;
  16
  17 % Noise variance
  18 sigma_n = 25;
  19 n = sigma_n*randn(N);
  20
  21 % Add noise
  22 x = s + n;
  23
  24 % Run local adaptive image denoising algorithm
  25 y = denoising_dwt(x);
  26
  27 % Calculate the error
  28 err = s - y;
  29
  30 % Calculate the PSNR value
  31 PSNR = 20*log10(256/std(err(:)))