www.pudn.com > texturedxl.rar > texture.m

function [hfreq,absffreq,pm]=texture(I1) 
I=double(I1); 
[m,n]=size(I); 
m1=m-2; 
n1=n-2; 
A1=zeros(m1,n1); 
A2=zeros(m1,n1); 
A3=zeros(m1,n1); 
A4=zeros(m1,n1); 
for x=1:m1; 
     for y=1:n1; 
         A1(x,y)=abs(I(x+1,y)-I(x+1,y+2)); 
         A2(x,y)=abs(I(x,y+1)-I(x+2,y+1)); 
         A3(x,y)=abs(I(x,y+2)-I(x+2,y)); 
         A4(x,y)=abs(I(x,y)-I(x+2,y+2)); 
     end; 
end; 
A5=(A1+A2+A3+A4)./4; 
A=uint8(A5); 
[counts,X]=imhist(A); 
hfreq=counts./(m1*n1); 
ffreq=fft(hfreq,256); 
absffreq=abs(ffreq); 
pm=max(hfreq);