%分析不同图像的纹理特征
close all;% 关闭所有图形窗口
clear all;%清除工作空间所有变量
clc;%清空命令行
% 步骤1,定义自相关函数zxcor(),建立zxcor.n文件
function[epsilon,eta,C]=zxcor(f,D,m,n);
% 自相关函数zxcor(),f为读入的图像数据,D为偏移距离【m,n】为图像尺寸,返回图像相关函数c的值
% epsilon和eta为图像自相关函数c的偏移量《偏移距离为人为设置,自己可根据需要设定》
for epsilon=1:D
for eta=1:D
temp=0;
fp=0;
for x=1:m
for y=1:n
if(x+epsilon-1)>m|(y+eta-1)>n
f1=0;
else
f1=f(x,y)*f(x+epsilon-1,y+eta-1);
end
temp=f1+temp;
fp=f(x,y)*f(x,y)+fp;
end
end
f2(epsilon,eta)=temp;
f3(epsilon,eta)=fp;
C(epsilon,eta)=f2(epsilon,eta)/f3(epsil
close all;% 关闭所有图形窗口
clear all;%清除工作空间所有变量
clc;%清空命令行
% 步骤1,定义自相关函数zxcor(),建立zxcor.n文件
function[epsilon,eta,C]=zxcor(f,D,m,n);
% 自相关函数zxcor(),f为读入的图像数据,D为偏移距离【m,n】为图像尺寸,返回图像相关函数c的值
% epsilon和eta为图像自相关函数c的偏移量《偏移距离为人为设置,自己可根据需要设定》
for epsilon=1:D
for eta=1:D
temp=0;
fp=0;
for x=1:m
for y=1:n
if(x+epsilon-1)>m|(y+eta-1)>n
f1=0;
else
f1=f(x,y)*f(x+epsilon-1,y+eta-1);
end
temp=f1+temp;
fp=f(x,y)*f(x,y)+fp;
end
end
f2(epsilon,eta)=temp;
f3(epsilon,eta)=fp;
C(epsilon,eta)=f2(epsilon,eta)/f3(epsil