Nonsubsampled Contourlet Transform 非下采样contourlet相关去噪

非下采样contourlet相关去噪，论文程序，包括非下采样contourlet变换工具箱，已通过仿真

clear all
clc
pfilt =‘maxflat‘;
dfilt = ‘dmaxflat7‘;
%nlevs = [0 0 4 4 5];    % Number of levels for DFB at each pyramidal level
nlevs = [123];
th = 3;                     % lead to 3*sigma threshold denoising
rho = 3;                    % noise level

% Test image: the usual suspect...
im = imread（‘camera.tif‘）;
im = double（im）;
%im = double（im）/ 256;
[mn]=size（im）;
f= wgn（mn10）;
var=std（f（:））;
nim=im+n;
nim=nim/256;
% Generate noisy image. 
sig = std（im（:））;
sigma = sig / rho;
sig = std（im（:））;
sigma = sig / rho;
%nim = im + sigma * randn（size（im））;
nn=nim;
%nim=[12;34];

%%%%% Contourlet Denoising %%%%%
% Contourlet transform
y = nsctdec（nimnlevsdfilt pfilt）;
%y1=y{11};% 64x64  double array （低频系数）
%y2=y{12};%64x128  64x128 
%y3=y{13};%128x128x4 
%y4=y{14};%128x256x4   256x128x4 
          %y5=y{15};%[32x128 double]x8  [128x32 double]x8
          %y6=y{16};%[32x256 double]x16  [256x32 double]x16

%相关去噪

for x=2:3
EC=zeros（1length（nlevs）+1）;
EP=zeros（1length（nlevs）+1）;
for k=1:2.^（x-1）
    nimm=y{1x}{1k};
    [mn]=size（nimm）;
nnimm（::）=zeros（m+2n+2）;
nnimm（11）=0;nnimm（1end）=0;nnimm（end1）=0;nnimm（endend）=0;
        nnimm（12:end-1）=nimm（1:）;%复制边界象素点
        nnimm（2:end-11）=nimm（:1）;
        nnimm（m+22:end-1）=nimm（m:）;
        nnimm（2:end-1n+2）=nimm（:n）;
for i=1:m
    for j=1:n
      nnimm（i+1j+1）=nimm（ij）;
    end
end
for i=1:m
    for j=1:n
       Ec1=0;Ec2=0;
        for a=1:3 %窗口大小为3
            for b=1:3
                %c=nnim（i+a-1j+b-1）;
                c1=nnimm（i+a-1j+b-1）;
                Ec1=c1+Ec1;
               % c2=（nnim（i+a-1j+b-1）-Ec1/9）^2;
                %Ec2=c2+Ec2;
               % noise（ij）=sum（（nnim（i+a-1j+b-1）-sum（nnim（i+a-1j+b-1））/（a*b））^2）;
                
            end
        end
        c2=（nimm（ij）-Ec1/9）^2;
        Ec2=c2+Ec2;
       noise（ij）=Ec2; 
    end
end
    for i=1:m
        for j=1:n      
  %相关量 
 while 1
  P_coeffs{1x}{1k}= y{1x}{1k}.* y{1x+1}{1k};
  EC_coeffs{1x}{1k}=y{1x}{1k}.^2;
  %归一化 
   EC（1k）=EC（1k）+EC_coeffs{1x}{1k}（ij）;
   EP_coeffs{1x}{1k}= P_coeffs{1x}{1k}.^2;
   EP（1k）=EP（1k）+EP_coeffs{1x}{1k}（ij）;
     Pav_coeffs{1x}{1k}= sqrt（abs（P_coeffs{1x}{1k}）.*sqrt（EP（1k）/EC（1k）））;
     
      if（abs（ Pav_coeffs{1x}{1k}（ij））>abs（y{1x}{1k}（ij）））
          %ynew{1x}{1k}（ij）=Pav_coeffs{1x}{1k}（ij）;
          ynew{1x}{1k}（ij）=sgn（y{1x}{1k}（ij））.*max（abs（y{1x}{1k}（ij）-noise（ij）/var）0）;
          y{1x}{1k}（ij）=0;
          P_coeffs{1x}{1k}（ij）=0;
      
      end
   if（abs（EC_coeffs{1x}{1k}（ij））-min（noise（:））>1e-20）% do while
     break
    end
  end
        end
       disp（i）;
       disp（j）;
    end
end
       disp（x）;
       disp（k）;
end
ynew{11}=y{11};ynew{14}=y{14};
m = nsctrec（ynew dfilt pfilt）;

%XIAOBO
y = pdfbdec（nim ‘9-7‘‘pkva‘ zeros（length（nlevs） 1））;
[c s] = pdfb2vec（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        336  2008-11-09 23:23  Nonsubsampled Contourlet Transform\code\noscdenoise.m

     文件       2969  2004-10-17 16:42  Nonsubsampled Contourlet Transform\code\nsct\atrousc.c

     文件       7168  2004-10-17 16:45  Nonsubsampled Contourlet Transform\code\nsct\atrousc.dll

     文件       8664  2005-01-24 21:51  Nonsubsampled Contourlet Transform\code\nsct\atrousc.mexmac

     文件       1234  2005-11-03 09:10  Nonsubsampled Contourlet Transform\code\nsct\atrousdec.m

     文件      19406  2004-09-09 15:01  Nonsubsampled Contourlet Transform\code\nsct\atrousfilters.m

     文件       1012  2004-10-17 16:45  Nonsubsampled Contourlet Transform\code\nsct\atrousrec.m

     文件       3214  2005-01-21 15:59  Nonsubsampled Contourlet Transform\code\nsct\decdemo.m

     文件       1495  2004-12-22 22:48  Nonsubsampled Contourlet Transform\code\nsct\dfbdecdemo.m

     文件      13234  2004-10-18 12:28  Nonsubsampled Contourlet Transform\code\nsct\dfilters.m

     文件       2607  2004-12-22 21:26  Nonsubsampled Contourlet Transform\code\nsct\dmaxflat.m

     文件       1103  2003-04-10 11:30  Nonsubsampled Contourlet Transform\code\nsct\efilter2.m

     文件       1861  2003-11-06 00:25  Nonsubsampled Contourlet Transform\code\nsct\extend2.m

     文件        856  2003-04-10 11:30  Nonsubsampled Contourlet Transform\code\nsct\ld2quin.m

     文件        662  2003-11-05 10:17  Nonsubsampled Contourlet Transform\code\nsct\ldfilter.m

     文件        854  2003-11-05 10:19  Nonsubsampled Contourlet Transform\code\nsct\mctrans.m

     文件        775  2003-04-10 11:30  Nonsubsampled Contourlet Transform\code\nsct\modulate2.m

     文件       4662  2005-11-03 10:47  Nonsubsampled Contourlet Transform\code\nsct\nsctdec.m

     文件       3257  2005-11-03 10:33  Nonsubsampled Contourlet Transform\code\nsct\nsctrec.m

     文件       5046  2005-01-20 23:13  Nonsubsampled Contourlet Transform\code\nsct\nsdfbdec.m

     文件       4931  2005-01-20 23:14  Nonsubsampled Contourlet Transform\code\nsct\nsdfbrec.m

     文件        889  2005-12-19 11:12  Nonsubsampled Contourlet Transform\code\nsct\nsfbdec.m

     文件        821  2004-10-21 07:11  Nonsubsampled Contourlet Transform\code\nsct\nsfbrec.m

     文件       2424  2005-01-20 23:51  Nonsubsampled Contourlet Transform\code\nsct\nssfbdec.m

     文件       2785  2005-01-20 23:50  Nonsubsampled Contourlet Transform\code\nsct\nssfbrec.m

     文件       1552  2004-08-07 18:43  Nonsubsampled Contourlet Transform\code\nsct\parafilters.m

     文件       1161  2003-11-07 19:58  Nonsubsampled Contourlet Transform\code\nsct\qupz.m

     文件       1557  2005-01-20 23:45  Nonsubsampled Contourlet Transform\code\nsct\README.txt

     文件       1981  2003-04-10 11:30  Nonsubsampled Contourlet Transform\code\nsct\resampz.m

     文件       1001  2005-01-20 23:27  Nonsubsampled Contourlet Transform\code\nsct\shownsct.m

............此处省略64个文件信息