图像中通常采用自相关函数作为纹理测度
自相关函数的定义为：
调用自定义函数 zxcor（）对砖墙面和大理石面纹理进行分析：
自定义函数 zxcor（）：
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(epsilon,eta); end end epsilon = 0:(D-1); eta = 0:(D-1); end
调用函数的测试代码如下：
close all;clear all;clc; f11 = imread('zhuanqiang.jpg'); f1 = rgb2gray(f11); f1 = double(f1); [m,n] = size(f1); D = 20; [epsilon1,eta1,C1]=zxcor(f1,D,m,n); f22 = imread('dalishi.jpg'); f2 = rgb2gray(f22); f2 = double(f2); [m,n] = size(f2); [epsilon2,eta2,C2]=zxcor(f2,20,m,n); figure; subplot(121);imshow(f11);title('砖墙'); subplot(122);imshow(f22);title('大理石'); figure; subplot(121);mesh(epsilon1,eta1,C1); xlabel('epsilon');ylabel('eta');title('砖墙'); subplot(122);mesh(epsilon2,eta2,C2); xlabel('epsilon');ylabel('eta');title('大理石'); 实验输出如下：
可以看到，自相关函数随着 
和 
的增加，砖墙的下降趋势比大理石的要快很多。可见砖墙表面的粗糙度更大。大理石面纹理的自相关函数下降幅度较小，说明其纹理表面相对光滑。