在C ++中使用imfiltermatlab

如何解决在C ++中使用imfiltermatlab

我有一个带过滤器的matlab代码。现在我想将其转换为c ++。我读了这个here 但是我不能用openCV3.3实现filterEngine（不再有filterEngine了。）

这是matlab中的过滤器：

imfilter(im2double(image),k,'same','replicate');

那是我的内核k：

k = fspecial('disk',3)

此行的结果是一个7x7矩阵，值介于0-1之间。

因此滤镜将获得输入图像。该值介于0和1之间（im2double）。因此，我必须将我的8位图像转换为此范围，我认为它的工作方式如下：

cv::normalize(inImage,outImage,1,cv::NORM_MINMAX);

然后使用内核过滤图像。 最后，我想再次输出一个8位图像。

所以我想将两个给定的Matlab行转换为C ++。

我在链接中读到，您不能简单地使用filter2D。 但是我该怎么办呢？有人有主意吗？

++++编辑++++ （测试zteffi的代码段）

我已经测试了代码。到目前为止，结果是不同的。

内核“磁盘”的结构略有不同（更多零）。其中的值也不同于Matlab。

Matlabkernel：

0   0.000280919186665314    0.0110250278796606  0.0171905962763014  0.0110250278796606  0.000280919186665314    0
0.000280919186665314    0.0245167426505023  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0245167426505023  0.000280919186665314
0.0110250278796606  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0110250278796606
0.0171905962763014  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0171905962763014
0.0110250278796606  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0110250278796606
0.000280919186665314    0.0245167426505023  0.0353677651315323  0.0353677651315323  0.0353677651315323  0.0245167426505023  0.000280919186665314
0   0.000280919186665314    0.0110250278796606  0.0171905962763014  0.0110250278796606  0.000280919186665314    0

c ++内核：

[0,0.0303030303030303,0;
 0,0;
 0.0303030303030303,0.0303030303030303;
 0.0303030303030303,0.0303030303030303;
 0,0]

我不知道如何在matlab中计算值。

以下行应在tha矩阵中创建值的“圆”。 但是在行维度中，边缘上有3 * 2个零，而在列维度中有2 * 2个零。但是每个角应该只有一个零。

Mat kernel = cv::getStructuringElement(cv::MORPH_ELLIPSE,cv::Size(radius * 2 + 1,radius * 2 + 1));

解决方法

fspecial('disk',3)的等价项为cv::getStructuringElement(cv::MORPH_ELLIPSE,cv::Size(radius * 2 + 1,radius * 2 + 1))，除以内核之和，因此总和等于1。所以类似：

Mat m = cv::imread("lena.jpg");
int radius = 3;
Mat kernel = cv::getStructuringElement(cv::MORPH_ELLIPSE,radius * 2 + 1));
double sum = cv::sum(cv::sum(kernel))[0]; // use cv::sum twice to get sum of all elements
kernel.convertTo(kernel,CV_64F,1 / sum);
cv::filter2D(m,m,kernel);
m.convertTo(m,CV_8U);

除非您需要显示图像，否则无需将图像转换为范围。在这种情况下，您实际上也不需要将其转换为double并将以下两行交换为：

cv::filter2D(m,CV_8U,kernel);

,

通过在编辑器中单击函数名称上的鼠标右键，然后在菜单中选择“打开'函数名称'”项，可以看到matlab对函数的实现。它将显示该功能的matlab代码，您可以将其准确地移植到C ++。