纹理增强matlab代码,纹理特征描述之自相关函数法 纹理粗糙性与自相关函数的扩展成正比 matlab代码实现...

图像中通常采用自相关函数作为纹理测度

自相关函数的定义为：

​

调用自定义函数 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('大理石');

实验输出如下：

​

​

可以看到，自相关函数随着

  和 

 的增加，砖墙的下降趋势比大理石的要快很多。可见砖墙表面的粗糙度更大。大理石面纹理的自相关函数下降幅度较小，说明其纹理表面相对光滑。