OpenCV实现Gabor滤波(2)

最新推荐文章于 2024-08-08 18:04:14 发布

大唐游子

最新推荐文章于 2024-08-08 18:04:14 发布

阅读量1.2w

点赞数 4

分类专栏：图像处理

本文链接：https://blog.csdn.net/lichengyu/article/details/24534245

版权

本文针对上一版Gabor滤波代码进行优化，解决了filter2D函数ddepth参数导致的误差问题，将ddepth设为CV_64。调整了Gabor核参数，包括增大kernel_size至69x69，并根据Mian Zhou的论文建议，设置了v的范围。通过这些改动，滤波器能更好地适应不同图像大小，特别是捕捉细节信息。附带展示了修改后的滤波结果和Gabor核阵列。

摘要由CSDN通过智能技术生成

对上一版的Gabor滤波代码进行了修改：

1) 将filter2D函数的ddepth参数值改为CV_64, 因为输入图像的类型是CV_8U, 而gabor kernel的类型是CV_64, 如果ddepth为默认值-1，计算的结果会有误差(为什么?)。

2) 修改了garbor kernel的参数，包括kernel_size 和 v的范围，都是为了适应不同图像的大小。v 越小，Gabor函数的宽度越小，越能刻画细节信息，适应于较小的图像。按照Mian Zhou的学位论文里的推荐，v的范围选择为-1到3，而kernel_size的确定公式为：

，对应v的kernel_size大小分别为19x19, 25x25, 35x35, 49x49, 69x69.

我将所有kenel_size统一选为69x69。构建出的gabor bank如下图：

修改过的代码如下：

#include <opencv2/core/core.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <cmath>
#include <iostream>

using namespace cv;
using namespace std;

const double PI = 3.14159265;

// ref: http://blog.csdn.net/watkinsong/article/details/7876361
Mat getMyGabor(int width, int height, int U, int V, double Kmax, double f,
			   double sigma, int ktype, const string kernel_name)
{
	//CV_ASSERT(width % 2 == 0 && height % 2 == 0);
	//CV_ASSERT(ktype == CV_32F || ktype == CV_64F);

	int half_width = width / 2;
	int half_height = height / 2;
	double Qu = PI*U/8;
	double sqsigma = sigma*sigma;
	double Kv = Kmax/pow(f,V);
	double postmean = exp(-sqsigma/2);

	Mat kernel_re(width, height, ktype);
	Mat kern