func_peaktransform.asv

一些个 人收集的小波变换的程序

function [out]=func_PeakTransform(inputimage)
[Row,Clom]=size(inputimage);
SUB_Row=Row/2;
SUB_Clom=Clom/2;
decompose_pic=inputimage;                     %  图像分为四块   
LL=decompose_pic(1:SUB_Row,1:SUB_Clom);      %  在矩阵左上方为低频分量--fi(x)*fi(y)
LH=decompose_pic(1:SUB_Row,SUB_Clom+1:Clom);    %  矩阵右上为--fi(x)*psi(y)
HL=decompose_pic(SUB_Row+1:Row,1:SUB_Clom);    %  矩阵左下为--psi(x)*fi(y)
HH=decompose_pic(SUB_Row+1:Row,SUB_Clom+1:Clom);  %  右下方为高频分量--psi(x)*psi(y)
%%%%%%%%%%%%%%%%%%%%%%【前向峰值变换】%%%%%%%%%%%%%%%%%%%%%%%%%%%
% PeakMapLL=LL*0;
% forwardPeakTransformLL=LL*0;
% for mm=1:SUB_Row
%     aLL=LL(mm,:);
%     for  row1=1:SUB_Clom-2
%         if(aLL(row1+1)>=aLL(row1)&aLL(row1+2)<=aLL(row1+1))
%            PeakMapLL(mm,row1+1)=aLL(row1+1);
%            forwardPeakTransformLL(mm,:)=func_forwardPeakTransform(LL(mm,:), PeakMapLL(mm,:));
%         end
%      end
% end
%%%%%%%%%【水平高频系数图像变换分析】%%%%%%%%%%%%
PeakMapLH=LH*0;
forwardPeakTransformLH=LH*0;
for mm=1:SUB_Row
    aLH=LH(mm,:);
    for  row1=1:SUB_Clom-2
        if(aLH(row1+1)>=aLH(row1)&aLH(row1+2)<=aLH(row1+1))
           PeakMapLH(mm,row1+1)=aLH(row1+1);
           forwardPeakTransformLH(mm,:)=func_forwardPeakTransform(LH(mm,:), PeakMapLH(mm,:));
        end
     end
end
%%%%%%%%%%【垂直高频系数图像变换分析】%%%%%%%%%%%%
PeakMapHL=HL*0;
forwardPeakTransformHL=HL*0;
for mm=1:SUB_Row
    aHL=HL(mm,:);
    for  row1=1:SUB_Clom-2
        if(aHL(row1+1)>=aHL(row1)&aHL(row1+2)<=aHL(row1+1))
           PeakMapHL(mm,row1+1)=aHL(row1+1);
           forwardPeakTransformHL(mm,:)=func_forwardPeakTransform(HL(mm,:), PeakMapHL(mm,:));
        end
    end
end
%%%%%%%%%%【斜线高频系数图像变换分析】%%%%%%%%%%%%
PeakMapHH=HH*0;
forwardPeakTransformHH=HH*0;
for mm=1:SUB_Row
    aHH=HH(mm,:);
    for  row1=1:SUB_Clom-2
        if(aHH(row1+1)>=aHH(row1)&aHH(row1+2)<=aHH(row1+1))
           PeakMapHH(mm,row1+1)=aHH(row1+1);
           forwardPeakTransformHH(mm,:)=func_forwardPeakTransform(HH(mm,:), PeakMapHH(mm,:));
        end
     end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%【后向峰值变换】%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% backwardPeakTransformLL=LL*0;
% for nn=1:SUB_Clom  
%    backwardPeakTransformLL(nn,:)=func_backwardPeakTransform(forwardPeakTransformLL(nn,:), PeakMapLL(nn,:));
% end   
%%%%%%%%%【水平高频系数图像变换分析】%%%%%%%%%%%%%%
backwardPeakTransformLH=LH*0;
for nn=1:SUB_Clom  
   backwardPeakTransformLH(nn,:)=func_backwardPeakTransform(forwardPeakTransformLH(nn,:), PeakMapLH(nn,:));
end   
%%%%%%%%%【垂直高频系数图像变换分析】%%%%%%%%%%%%%%
backwardPeakTransformHL=LH*0;
for nn=1:SUB_Clom  
   backwardPeakTransformHL(nn,:)=func_backwardPeakTransform(forwardPeakTransformHL(nn,:), PeakMapHL(nn,:));
end  
%%%%%%%%%【斜线高频系数图像变换分析】%%%%%%%%%%%%%%
backwardPeakTransformHH=LH*0;
for nn=1:SUB_Clom  
   backwardPeakTransformHH(nn,:)=func_backwardPeakTransform(forwardPeakTransformHH(nn,:), PeakMapHH(nn,:));
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
backwardPeakTransformout(1:SUB_Row,1:SUB_Clom)=LL;
backwardPeakTransformout(1:SUB_Row,SUB_Clom+1:Clom)=backwardPeakTransformLH;
backwardPeakTransformout(SUB_Row+1:Row,1:SUB_Clom)=backwardPeakTransformHL;
backwardPeakTransformout(SUB_Row+1:Row,SUB_Clom+1:Clom)=backwardPeakTransformHH;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
out=backwardPeakTransformout;