DoG (Difference of Gaussian)角点检测

转载原文地址:http://blog.csdn.net/abcjennifer/article/details/7639488


DoG (Difference of Gaussian)是灰度图像增强和角点检测的方法,其做法较简单,证明较复杂,具体讲解如下:

Difference of Gaussian(DOG)是高斯函数的差分。我们已经知道可以通过将图像与高斯函数进行卷积得到一幅图像的低通滤波结果,即去噪过程,这里的Gaussian和高斯低通滤波器的高斯一样,是一个函数,即为正态分布函数。


那么difference of Gaussian 即高斯函数差分是两幅高斯图像的差,

一维表示:

二维表示:


具体到图像处理来讲,就是将两幅图像在不同参数下的高斯滤波结果相减,得到DoG图。


1. 处理一幅图像在不同参数下的DoG

[cpp]  view plain copy
  1. A = Process(Im, 0.3, 0.4, x);  
  2. B = Process(Im, 0.6, 0.7, x);  
  3. a = getExtrema(A, B, C, thresh);  
其中,

[cpp]  view plain copy
  1. function [ out_img ] = Process( img, sig1, sig2, size )  

是求图像DoG的结果,两个高斯平滑参数分别为sig1和sig2,结果如下:

[cpp]  view plain copy
  1. A = Process(Im, 0.3, 0.4, x);  

DoG (Difference of Gaussian)角点检测_第1张图片

[cpp]  view plain copy
  1. B = Process(Im, 0.6, 0.7, x);  
DoG (Difference of Gaussian)角点检测_第2张图片

[cpp]  view plain copy
  1. C = Process(Im, 0.7, 0.8, x);  
DoG (Difference of Gaussian)角点检测_第3张图片

2. 根据DOG求角点

Theory:DOG三维图中的最大值和最小值点是角点

DoG (Difference of Gaussian)角点检测_第4张图片

X标记当前像素点,绿色的圈标记邻接像素点,用这个方式,最多检测26个像素点。X被标记为特征点,如果它是所有邻接像素点的最大值或最小值点。

因此在上一步计算出的A,B,C三个DOG图中求图B中是极值的点,并标记(max:1;min:-1)

[cpp]  view plain copy
  1. a = getExtrema(A, B, C, thresh);  
  2. figure;  
  3. imshow(a, [-1 1]);  

[cpp]  view plain copy
  1.   
结果:

黑色为极小值,白色为极大值

DoG (Difference of Gaussian)角点检测_第5张图片

最后在原图上予以显示:

DoG (Difference of Gaussian)角点检测_第6张图片


就得到了一幅图的DOG角点检测结果。


代码下载链接:http://download.csdn.net/detail/abcjennifer/4363144

你可能感兴趣的:(matlab,图像处理,Dog,角点检测)