contourlet代码

图像去噪方面用到contourlet变换，源代码

function y = backsamp（y）
% BACKSAMP    Backsampling the subband images of the directional filter bank
%
%   y = backsamp（y）
%
% Input and output are cell vector of dyadic length
%
% This function is called at the end of the DFBDEC to obtain subband images 
% with overall sampling as diagonal matrices
%
% See also: DFBDEC

% Number of decomposition tree levels
n = log2（length（y））;

if （n ~= round（n）） | （n < 1）
    error（‘Input must be a cell vector of dyadic length‘）;
end

if n == 1
    % One level the decomposition filterbank shoud be Q1r
    % Undo the last resampling （Q1r = R2 * D1 * R3）
    for k = 1:2
	y{k} = resamp（y{k} 4）;
	y{k}（: 1:2:end） = resamp（y{k}（: 1:2:end） 1）;
	y{k}（: 2:2:end） = resamp（y{k}（: 2:2:end） 1）;	
    end    
    
elseif n > 2    
    N = 2^（n-1）;
    
    for k = 1:2^（n-2）
	shift = 2*k - （2^（n-2） + 1）;
	
	% The first half channels
	y{2*k-1} = resamp（y{2*k-1} 3 shift）;
	y{2*k} = resamp（y{2*k} 3 shift）;
	
	% The second half channels
	y{2*k-1+N} = resamp（y{2*k-1+N} 1 shift）;
	y{2*k+N} = resamp（y{2*k+N} 1 shift）;	
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2390  2003-11-12 23:50  ChangeLog
     文件        5127  2003-11-05 00:10  computescale.m
     文件        3537  2003-11-11 10:45  Contents.m
     文件         706  2003-11-11 11:41  COPYING
     文件        3437  2003-11-10 12:22  decdemo.m
     文件        2472  2003-11-13 00:17  denoisedemo.m
     文件        2223  2003-11-05 22:44  dfbdec.m
     文件        1768  2003-11-05 23:54  dfbdec_l.m
     文件        1748  2003-11-05 08:11  dfbimage.m
     文件        2012  2003-11-10 00:52  dfbrec.m
     文件        1736  2003-11-06 00:18  dfbrec_l.m
     文件        2729  2003-11-06 10:39  dfilters.m
     文件         755  2003-04-10 10:30  dup.m
     文件        1103  2003-04-10 10:30  efilter2.m
     文件        1861  2003-11-05 22:25  extend2.m
     文件        2098  2003-04-10 10:30  fbdec.m
     文件        1507  2003-11-06 10:41  fbdec_l.m
     文件        2360  2003-04-10 10:30  fbrec.m
     文件        1426  2003-11-06 00:23  fbrec_l.m
     文件         485  2003-04-10 10:30  ffilters.m
     文件         856  2003-04-10 10:30  ld2quin.m
     文件         662  2003-11-05 08:17  ldfilter.m
     文件         722  2003-09-13 02:30  lpdec.m
     文件         935  2003-09-13 02:48  lprec.m
     文件         882  2003-11-05 08:19  mctrans.m
     文件         775  2003-04-10 10:30  modulate2.m
     文件        5507  2003-11-10 18:07  nlademo.m
     文件        1238  2003-11-11 10:26  nlademo2.m
     文件        1254  2003-11-12 23:43  pdfb2vec.m
     文件         931  2003-11-13 00:25  pdfb_nest.m
     文件        1254  2003-11-13 00:09  pdfb_tr.m
............此处省略34个文件信息