matlab纹理特征提取代码

 function [t]=statxture(f,scale)
 if nargin==1
     scale(1:6)=1;
 else
     scale=scale(1:6)';
 end
 p=imhist(f);                  %p是256*1的列向量
 p=p./numel(f);
 L=length(p);
 [v,mu]=statmoments(p,3);
 %计算六个纹理特征
 t(1)=mu(1);                   %平均值
 t(2)=mu(2).^0.5;              %标准差
 varn=mu(2)/(L-1)^2;
 t(3)=1-1/(1+varn);            %平滑度首先为（0~1）区间通过除以（L-1）^2将变量标准化
 t(4)=mu(3)/(L-1)^2;           %三阶矩（通过除以（L-1）^2将变量标准化）
 t(5)=sum(p.^2);               %一致性
 t(6)=-sum(p.*(log2(p+eps)));  %熵
 T=[t(1) t(2) t(3) t(4) t(5) t(6)]
 %缩放值，默认为1
 t=t.*scale;
 end

 function [v,unv]=statmoments(p,n)
 Lp=length(p);
 if (Lp~=256)&(Lp~=65536)
     error('p must be a 256- or 65536-element vector.');
 end
 G=Lp-1;
 p=p/sum(p);p=p(:);
 z=0:G;
 z=z./G;
 m=z*p;
 z=z-m;
 v=zeros(1,n);
 v(1)=m;
 for j=2:n
     v(j)=(z.^j)*p;
 end
 if nargout>1
     unv=zeros(1,n);
     unv(1)=m.*G;
     for j=2:n
         unv(j)=((z*G).^j)*p
     end
 end
 end

f为输入的图片，必须是灰度图像。

因为纹理特征提取是灰度的纹理特征。

matlab中建立一个function的M文件

拷贝上面的代码运行就可以了。