function?OUT=GrayExitMonentR（IN）
%?灰度共生矩陣為p2;
%?二階矩?f1;
%?相關(guān)度?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
?%復(fù)制一個共生矩陣
p2=p1;
%?figureimshow（p2）;


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

%計算二階矩?f1
f1=p2.^2;
f1=sum（f1（:））;

%?計算相關(guān)度?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