%調用舉例：
%image=rgb2gray（imread（‘007.jpg‘））;
%f=coarseness（image5）
function?Fcrs?=?coarseness（?graypickmax?）%graphic為待處理的灰度圖像，2^kmax為最大窗口
[hw]=size（graypic）;?%獲取圖片大小
A=zeros（hw2^kmax）;?%平均灰度值矩陣A
%計算有效可計算范圍內每個點的2^k鄰域內的平均灰度值
for?i=2^（kmax-1）+1:h-2^（kmax-1）
????for?j=2^（kmax-1）+1:w-2^（kmax-1）
????????for?k=1:kmax
????????????A（ijk）=mean2（graypic（i-2^（k-1）:i+2^（k-1）-1j-2^（k-1）:j+2^（k-1）-1））;
????????end
????end
end
%對每個像素點計算在水平和垂直方向上不重疊窗口之間的Ak差
for?i=1+2^（kmax-1）:h-2^（kmax-1）
????for?j=1+2^（kmax-1）:w-2^（kmax-1）
????????for?k=1:kmax
????????????Eh（ijk）=abs（A（i+2^（k-1）jk）-A（i-2^（k-1）j））;
????????????Ev（ijk）=abs（A（ij+2^（k-1）k）-A（ij-2^（k-1）））;
????????end
????end
end
%對每個像素點計算使E達到最大值的k
for?i=2^（kmax-1）+1:h-2^（kmax-1）
????for?j=2^（kmax-1）+1:w-2^（kmax-1）
???

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

?????文件???????1200??2018-07-29?15:19??Tamura\coarseness.m

?????文件????????460??2018-07-29?15:39??Tamura\contrast.m

?????文件???????1171??2018-07-29?15:40??Tamura\directionality.m

?????文件???????3515??2018-07-29?15:40??Tamura\linelikeness.m

?????文件????????830??2018-07-29?15:41??Tamura\regularity.m

?????文件????????375??2018-07-29?15:41??Tamura\Tamura.m

?????文件???????7667??2018-07-27?18:40??Tamura\Tamura紋理特征.txt

?????目錄??????????0??2018-08-10?17:21??Tamura

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

????????????????15218????????????????????8