将 2 个直方图与卡方进行比较

Question

我想比较 2 个直方图，它们有 2 个维度。 为此，我想使用 卡方-Metric。 我的比较器看起来像这个函数：

double Histogram::compareHistogram(Histogram *hist){
    double result=0;
    double a=0;
    double b=0;
    for (int y=0 ; y < bins_1 ; y++) {
         for (int x=0 ; x < bins_2 ; x++) {
             a=getHistogramValue(x,y)-hist->getHistogramValue(x,y);
             b=getHistogramValue(x,y)+hist->getHistogramValue(x,y);
             if(fabs(b)>0.0){
               result+=a*a/b;
             }
         }
    }
    return result;
}

我将结果与 OpenCv 的 cv::compareHist() 函数的结果进行了比较，结果有所不同。我不知道为什么。

在比较直方图之前，我使用 MINMAX-Norm 对直方图进行规范。 我将我的标准化直方图与 openCV 的标准化直方图进行了比较，它们是相等的。 所以我认为，问题出在我的 compareHist 函数中。 但是在哪里？

最好的问候，

Answer 1

来自OpenCV的相关部分源码如下：

if( method == CV_COMP_CHISQR )
{
    for( j = 0; j < len; j++ )
    {
        double a = h1[j] - h2[j];
        double b = h1[j];
        if( fabs(b) > DBL_EPSILON )
            result += a*a/b;
    }
}

所以你可以看到你的代码的不同是这一行

b=getHistogramValue(x,y)+hist->getHistogramValue(x,y);

应该是

b=getHistogramValue(x,y);