www.pudn.com > wenlituxiangjiansuo.rar > TextureDTCWT.m

clear; 
clc; 
level = 3; 
FV_e =[]; 
FV_std = []; 
 
image = imread('22.jpg'); 
image = RGB2GRAY(image); 
[Yl,Yh] = dtwavexfm2(image,level,'near_sym_b','qshift_b'); 
 
sz = size(Yl); 
FV_e(1) = sum(abs(Yl(:)))/sz(1)*sz(2); 
FV_std(1) = std(abs(Yl(:))); 
 
matric = []; 
for i=1:level 
    for j=1:6 
        matric = Yh{i,1}(:,:,j); 
        sz = size(matric); 
        FV_e((i-1)*6+1+j) = sum(abs(matric(:)))/(sz(1)*sz(2)); 
        FV_std((i-1)*6+1+j) = std(abs(matric(:))); 
    end 
end 
 
FV=[FV_e,FV_std];