c++ - GaussianBlur 不应该是对称的吗？

Question

我希望高斯模糊操作是对称的，但使用 OpenCV 2.4.11 GaussianBlur我遇到了不同。

这是一个例子。我申请 GaussianBlur图像，以及图像的翻转版本。我已经单独验证了 flip操作不会改变图像像素值(未显示)。当我将模糊图像翻转回来时，我希望它与原始图像的模糊相同，但是 diff显示出很多小差异(在 0.0 和 6.103515625e-005 之间)。我知道这很小，但它对我的后续处理有链式 react 。

高斯核是对称的，所以结果应该是一样的。这仅仅是实现中的舍入错误吗？

int main(int, char **)
{
    // e.g. 2008_005541.jpg from VOC2012 dataset
    char const * const filename = "...";
    float const sig_diff = 1.24899971f;

    cv::Mat image = cv::imread(filename, cv::IMREAD_GRAYSCALE);
    cv::Mat gray_fpt;
    image.convertTo(gray_fpt, cv::DataType<float>::type, 1, 0);
    GaussianBlur(gray_fpt, gray_fpt, cv::Size(), sig_diff, sig_diff);

    cv::Mat mirror;
    flip(image, mirror, 1);
    cv::Mat mirror_gray_fpt;
    mirror.convertTo(mirror_gray_fpt, cv::DataType<float>::type, 1, 0);
    GaussianBlur(mirror_gray_fpt, mirror_gray_fpt, cv::Size(), sig_diff, sig_diff);

    flip(mirror_gray_fpt, mirror_gray_fpt, 1);
    cv::Mat diff = abs(gray_fpt - mirror_gray_fpt);
    double minval, maxval;
    minMaxLoc(diff, &minval, &maxval);

    // minval = 0.0;
    // maxval = 6.103515625e-005;

    // easier to visualise the differences with this:
    normalize(diff, diff, 0.0, 1.0, cv::NORM_MINMAX, CV_32FC1);
    return 0;
}

编辑:我从 cv::DataType<float>::type 更改了类型至 cv::DataType<double>::type现在最大误差是 1.1368683772161603e-013，所以四舍五入似乎是问题所在。

Answer 1

将上面的代码更改为调用 gaussian_blur(下面)而不是 GaussianBlur 不会在我目前测试过的示例图像中产生任何差异。

由此可见，如果高斯运算的工作类型是double精度，那么float点精度的输出不会产生错误。这似乎是解决我的问题的好办法。

// Perform gaussian blur in double precision and convert back
void gaussian_blur(
    cv::Mat const &src, cv::Mat &dst,
    cv::Size ksize, double sigmaX, double sigmaY=0,
    int borderType=cv::BORDER_DEFAULT)
{
    cv::Mat src_dp;
    src.convertTo(src_dp, cv::DataType<double>::type, SIFT_FIXPT_SCALE, 0);

    cv::Mat dst_dp;
    GaussianBlur(src_dp, dst_dp, ksize, sigmaX, sigmaY, borderType);

    dst_dp.convertTo(dst, src.type(), SIFT_FIXPT_SCALE, 0);
}