gusucode.com > 红外图像增强及目标检测演示界面matlab源码程序 > code/TrousTransform.m

    function C=TrousTransform(F,N)
%trous wavelet transformation use B-样条小波
%f为一幅灰度图像，W为f的trous小波变换高频部分，C为相应的低频
%边界处理采用镜像对称的方法
[row col]=size(F);
h=[1/256 1/64 3/128 1/64 1/256;
   1/64  1/16 3/32  1/16 1/64 ;
   3/128 3/32 9/64 3/32 3/128;
   1/64  1/16 3/32  1/16 1/64 ;
   1/256 1/64 3/128 1/64 1/256];

%W=f;
%C=f;%初始化
C=zeros(row,col,N);
for i=1:N
    C(1:row,1:col,i)=conv2(F,h,'same');
    F=C(1:row,1:col,i);
end
%for x=1:row 
 %   for y=1:col
  %      C(x,y)=0;
   %     for i=-2:2
    %        for j=-2:2
     %           C(x,y)=C(x,y)+h(i+3,j+3)*tft(f,x+i,y+j);
     %      end
     %    end
     % end
     %end