function [smoment,entropy,contrast,homogeneity]=GraidFeature(pmatrix)
%计算基于空间域的纹理特征
%input：
%pmatrix：归一化共生矩阵,包括 0度，45度，90度，135度
%output：
%smoment:二阶矩（能量）
%entropy：熵
%contrast:对比度
%homogeneity:同质性
[w,h]=size(pmatrix{1});%归一化共生矩阵的行，列大小
smoment=zeros(1,4);
entropy=zeros(1,4);
contrast=zeros(1,4);
homogeneity=zeros(1,4);
for i=1:w
for j=1:h
for k=1:4
smoment(1,k)=smoment(1,k)+pmatrix{k}(i,j).^2;%二阶矩
if(pmatrix{k}(i,j)>0)
entropy(1,k)=entropy(1,k)+pmatrix{k}(i,j)*log(pmatrix{k}(i,j));%熵
end
homogeneity(1,k)=homogeneity(1,k)+pmatrix{k}(i,j)/(1+i*j);%同质性
contrast(1,k)=contrast(1,k)+i*j*pmatrix{k}(i,j);%对比度
end
end
end
entropy=-entropy;
评论0