%計算Tamura紋理特征
close?all
clear?all
clc
t0=cputime;
I?=?imread（‘D1.gif‘）;
[NxNy]?=?size（I）;
Ng=256;
G=double（I）;
%計算粗糙度（coarseness）
Sbest=zeros（NxNy）;
E0h=zeros（NxNy）;
E0v=zeros（NxNy）;
E1h=zeros（NxNy）;
E1v=zeros（NxNy）;
E2h=zeros（NxNy）;
E2v=zeros（NxNy）;
E3h=zeros（NxNy）;
E3v=zeros（NxNy）;
E4h=zeros（NxNy）;
E4v=zeros（NxNy）;
E5h=zeros（NxNy）;
E5v=zeros（NxNy）;
flag=0;

for?i=1:Nx
????for?j=2:Ny
????????E0h（ij）=G（ij）-G（ij-1）;
????end
end
E0h=E0h/2;

for?i=1:Nx-1
????for?j=1:Ny
????????E0v（ij）=G（ij）-G（i+1j）;
????end
end
E0v=E0v/2;

%圖片大小必須大于4*4才能計算E1h、E1v
if?（Nx<4||Ny<4）
????flag=1;
end
if（flag==0）
????for?i=1:Nx-1
????????for?j=3:Ny-1
????????????E1h（ij）=sum（sum（G（i:i+1j:j+1）））-sum（sum（G（i:i+1j-2:j-1）））;
????????end
????end
????for?i=2:Nx-2
????????for?j=2:Ny
????????????E1v（ij）=sum（sum（G（i-1:ij-1:j）））-sum（sum（G（i+1:i+2j-1:j）））;
????????end
????end
????E1h=E1h/4;
????E1v=E1v/4;
end


%圖片大小必須大于8*8才能計算E2h、E2v
if?（Nx<8||Ny<8）
????flag=1;
end
if（flag==0）
????for?i=2:Nx-2
????????for?j=5:Ny-3
????????????E2h（ij）=sum（sum（G（i-1:i+2j:j+3）））-sum（sum（G（i-1:i+2j-4:j-1）））;
????????end
????end
????for?i=4:Nx-4
????????for?j=3:Ny-1
????????????E2v（ij）=sum（sum（G（i-3:ij-2:j+1）））-sum（sum（G（i+1:i+4j-2:j+1）））;
????????end
????end
????E2h=E2h/16;
????E2v=E2v/16;
end


%圖片大小必須大于16*16才能計算E3h、E3v
if?（Nx<16||Ny<16）
????flag=1
end
if（flag==0）
????for?i=4:Nx-4
????????for?j=9:Ny-7
????????????E3h（ij）=sum（sum（G（i-3:i+4j:j+7）））-sum（sum（G（i-3:i+4j-8:j-1）））;
????????end
????end
????for?i=8:Nx-8
????????for?j=5:Ny-3
????????????E3v（ij）=sum（sum（G（i-7:ij-4:j+3）））-sum（sum（G（i+1:i+8j-4:j+3）））;
????????end
????end
????E3h=E3h/64;
????E3v=E3v/64;
end


%圖片大小必須大于32*32才能計算E4h、E4v?
if?（Nx<32||Ny<32）
????flag=1;
end
if（flag==0）
????for?i=8:Nx-8
????????for?j=17:Ny-15
????????????E4h（ij）=sum（sum（G（i-7:i+8j:j+15）））-sum（sum（G（i-7:i+8j-16:j-1）））;
????????end
????end
????for?i=16:Nx-16
????????for?j=9:Ny-7
????????????E4v（ij）=sum（sum（G（i-15:ij-8:j+7）））-sum（sum（G（i+1:i+16j-8:j+7）））;
????????end
????end
????E4h=E4h/256;
????E4v=E4v/256;
end


%圖片大小必須大于64*64才能計算E5h、E5v
if?（Nx<64||Ny<64）
????flag=1;
end
if（flag==0）
????for?i=16:Nx-16
????????for?j=33:Ny-31
????????????E5h（ij）=sum（sum（G（i-15:i+16j:j+31）））-sum（sum（G（i-15:i+16j-32:j-31）））;
????????end
????end
????for?i=32:Nx-32
????????for?j=17:Ny-15
????????????E5v（ij）=sum（sum（G（i-31:ij-16:j+15）））-sum（sum（G（i+1:i+32j-16:j+15）））;
????????end
????end
????E5h=E5h/1024;
????E5v=E5v/1024;
end

for?i=1:Nx
????for?j=1:Ny
????????[maxvindex]=max（[E0h（ij）E0v（ij）E1h（ij）E1v（ij）E2h（ij）E2v（ij）E3h（ij）E3v（ij）E4h（ij）E4v（ij）E5h（ij）E5v（ij）]）;
????????k=floor（（index+1）/2）;
????????Sbest（ij）=2.^k;
????end
end
Fcoarseness=sum（sum（Sbest））/（Nx*Ny）;
%計算對比度
[countsgraylevels]=imhist（I）;
PI=counts/（Nx*Ny）;
averagevalue=sum（graylevels.*PI）;
u4=sum（（graylevels-repmat（averagevalue[2561]））.^4.*PI）;
standarddeviation=sum（（graylevels-repmat（averagevalue[2561]））.^2.*PI）;
alpha4=u4/standarddeviation^2;
Fcontrast=

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

?????文件???????1626??2008-10-26?20:50??Roughness\ReadMe.txt

?????文件????????450??2008-10-26?20:50??Roughness\Resource.h

?????文件????????246??2008-10-26?20:50??Roughness\Roughness.clw

?????文件??????12922??2008-10-27?14:56??Roughness\Roughness.cpp

?????文件???????4394??2008-10-26?20:50??Roughness\Roughness.dsp

?????文件????????526??2008-10-26?20:50??Roughness\Roughness.dsw

?????文件????????532??2008-10-27?14:49??Roughness\Roughness.h

?????文件??????50176??2008-10-27?15:10??Roughness\Roughness.ncb

?????文件??????49664??2008-10-27?15:10??Roughness\Roughness.opt

?????文件???????1983??2008-10-27?15:10??Roughness\Roughness.plg

?????文件???????1457??2008-10-26?20:50??Roughness\Roughness.rc

?????文件????????892??2008-10-26?19:49??Roughness\Roughness.sln

????..A..H.?????21504??2008-10-27?15:10??Roughness\Roughness.suo

?????文件????????296??2008-10-26?20:50??Roughness\StdAfx.cpp

?????文件???????1106??2008-10-26?20:50??Roughness\StdAfx.h

?????目錄??????????0??2008-11-07?18:36??Roughness

?????文件?????282020??2006-11-18?12:14??brodatz紋理庫\D1.gif

?????文件?????275134??2006-11-18?12:25??brodatz紋理庫\D10.gif

?????文件?????372412??2006-11-18?12:17??brodatz紋理庫\D2.gif

?????文件?????335990??2006-11-18?12:17??brodatz紋理庫\D3.gif

?????文件?????415253??2006-11-18?12:21??brodatz紋理庫\D4.gif

?????文件?????351599??2006-11-18?12:20??brodatz紋理庫\D5.gif

?????文件?????325555??2006-11-18?12:22??brodatz紋理庫\D6.gif

?????文件?????305841??2006-11-18?12:23??brodatz紋理庫\D7.gif

?????文件?????242409??2006-11-18?12:23??brodatz紋理庫\D8.gif

?????文件?????429522??2006-11-18?12:25??brodatz紋理庫\D9.gif

?????文件???????4704??2008-12-15?11:19??brodatz紋理庫\Tamura.m

????..A.SH.?????59392??2008-11-16?18:18??brodatz紋理庫\Thumbs.db

?????目錄??????????0??2008-12-16?15:52??brodatz紋理庫

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

............此處省略2個文件信息