Question

我期望高斯模糊操作是对称的，但是使用OpenCV 2.4.11 var orderedSubList = subList.OrderBy ( s => masterList.IndexOf(masterList.FirstOrDefault(m => m.RollNumber == s.RollNumber)) );我会得到差异。

这是一个例子。我将GaussianBlur应用于图像，并应用于图像的翻转版本。我已单独验证GaussianBlur操作不会更改图像像素值（未显示）。当我将模糊图像翻转回来时，我预计它会与原始模糊图像相同，但flip显示出很多小差异（介于0.0和6.103515625e-005之间）。我知道这很小，但它在我后续的处理过程中会产生连锁反应。

高斯核是对称的，因此结果应该是相同的。这只是实现中的舍入错误吗？

diff

编辑：我将类型从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，现在最大错误为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);
}

GaussianBlur不应该是对称的吗？

1 个答案: