THESE/codes/wave/allcode/main_dtdwt.m

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