灰度梯度共生矩阵的提取方式，matlab实现的

灰度共生矩阵和灰度梯度共生矩阵的提取方式，是比较重要的纹理特征提取方法，用matlab实现的

function OUT=GrayExitMonentR（IN）
% 灰度共生矩阵为p2;
% 二阶矩 f1;
% 相关度 f2;
% 墒 f3;
% 对比度 f4;
% 逆差矩 f5;
% 和方差 f6;
% close all;
IN= imread（‘F:\picture\oil.jpg‘）;
IN=double（IN）;
IN=rgb2gray（IN）;
gray=256;




  
[RC]=size（IN）;
%把图象降为32像素
% IN=uint8（IN/8）;
% max（IN（:））


p1=zeros（gray）;
for M=1:R
    for N=1:（C-1）
        p1（IN（MN）+1IN（MN+1）+1）= p1（IN（MN）+1IN（MN+1）+1）+1;
        p1（IN（MN+1）+1IN（MN）+1）= p1（IN（MN+1）+1IN（MN）+1）+1;
    end
end
 %复制一个共生矩阵
p2=p1;
% figureimshow（p2）;


%归一化
p2=double（p2）/（M*N）;

%计算二阶矩 f1
f1=p2.^2;
f1=sum（f1（:））;

% 计算相关度 f2
colour=0:gray-1;
m1=sum（p2）*colour（:）;
m2=m1;
pp=sum（p2）;
v1=（（colour-m1）.^2）*pp（:）;
 v2=v1;
[mn]=meshgrid（colourcolour）;
pp=m.*n.*p2;
f2=（sum（pp（:））-m1*m2）/（v1+eps）*（v2+eps）;
%计算墒 f3
p

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1786  2009-09-16 14:03  texture_extraction\GrayExitMonentR.m

     文件       3105  2006-07-07 20:59  texture_extraction\GrayGradinet.m

     文件         86  2008-04-11 14:26  texture_extraction\read.txt

     文件       1729  2009-09-16 13:51  texture_extraction\GrayExitMonentR.asv

     目录          0  2009-10-09 20:27  texture_extraction

----------- ---------  ---------- -----  ----

                 6706                    5