function?T=coocurrence（Image）

Gray=rgb2gray（Image）;
[MN]?=?size（Gray）;
%為了減少計算量，對原始圖像灰度級壓縮，將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
%計算四個共生矩陣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（