%**************************************************************************
%???????????????????圖像檢索——紋理特征
%基于共生矩陣紋理特征提取，d=1θ=0°45°90°135°共四個矩陣
%所用圖像灰度級均為256
%參考《基于顏色空間和紋理特征的圖像檢索》
%function?:?T=Texture（Image）?
%Image????:?輸入圖像數(shù)據(jù)
%T????????:?返回八維紋理特征行向量
%**************************************************************************
function?T?=?Texture（Image）
%?Image?=?imread（‘E:\1\3.tiff‘）;
%?[MNO]?=?size（Image）;
M?=?256;?
N?=?256;

%--------------------------------------------------------------------------
%1.將各顏色分量轉(zhuǎn)化為灰度
%--------------------------------------------------------------------------
Gray?=?double（0.3*Image（::1）+0.59*Image（::2）+0.11*Image（::3））;

%--------------------------------------------------------------------------
%2.為了減少計算量，對原始圖像灰度級壓縮，將Gray量化成16級
%--------------------------------------------------------------------------
for?i?=?1:M
????for?j?=?1:N
????????for?n?=?1:256/16
????????????if?（n-1）*16<=Gray（ij）&Gray（ij）<=（n-1）*16+15
????????????????Gray（ij）?=?n-1;
????????????end
????????end
????end
end

%--------------------------------------------------------------------------
%3.計算四個共生矩陣P取距離為1，角度分別為04590135
%--------------------------------------------------------------------------
P?=?zeros（16164）;
for?m?=?1:16
????for?n?=?1:16
????????for?i?=?1:M
????????????for?j?=?1:N
????????????????if?j????????????????????P（mn1）?=?P（mn1）+1;
????????????????????P（nm1）?=?P（mn1）;
????????????????end
????????????????if?i>1&j????????????????????P（mn2）?=?P（mn2）+1;
????????????????????P（nm2）?=?P（mn2）;
????????????????end
????????????????if?i????????????????????P（mn3）?=?P（mn3）+1;
????????????????????P（nm3）?=?P（mn3）;
????????????????end
????????????????if?i????????????????????P（mn4）?=?P（mn4）+1;
????????????????????P（nm4）?=?P（mn4）;
????????????????end
????????????end
????????end
????????if?m==n
????????????P（mn:）?=?P（mn:）*2;
????????end
????end
end

%%---------------------------------------------------------
%?對共生矩陣歸一化
%%---------------------------------------------------------
for?n?=?1:4
????P（::n）?=?P（::n）/sum（sum（P（::n）））;
end

%--------------------------------------------------------------------------
%4.對共生矩陣計算能量、熵、慣性矩、相關(guān)4個紋理參數(shù)
%--------------------------------------------------------------------------
H?=?zeros（14）;
I?=?H;
Ux?=?H;??????Uy?=?H;
deltaX=?H;??deltaY?=?H;
C?=H;
for?n?=?1:4
????E（n）?=?sum（sum（P（::n）.^2））;?%%能量
????for?i?=?1:16
????????for?j?=?1:16
????????????if?P（ijn）~=0
????????????????H（n）?=?-P（ijn）*log（P（ijn））+H（n）;?%%熵
????????????end
????????????I（n）?=?（i-j）^2*P（ijn）+I（n）;??%%慣性矩
???????????
????????????Ux（n）?=?i*P（ijn）+Ux（n）;?%相關(guān)性中μx
????????????Uy（n）?=?j*P（ijn）+Uy（n）;?%相關(guān)性中μy
????????end
????end
end
for?n?=?1:4

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3879??2008-01-17?11:39??紋理特征的提取\Texture.m

?????目錄??????????0??2010-04-22?09:47??紋理特征的提取

-----------?---------??----------?-----??----

?????????????????4097????????????????????3