Коэффициент корреляции Пирсона

Я написал следующий код C#, чтобы найти коэффициент корреляции Пирсона между двумя изображениями. Полный исходный код находится здесь, в DotNetFiddle.

Исходный код корреляции:

public sealed class PearsonCorrelation 
{
    public static double GetSimilarityScore(double[,] p, double[,] q)
    {
        int Width = p.GetLength(0);
        int Height = p.GetLength(1);

        if (Width != q.GetLength(0) || Height != q.GetLength(1))
        {
            throw new ArgumentException("Input vectors must be of the same dimension.");
        }

        double pSum = 0, qSum = 0, pSumSq = 0, qSumSq = 0, productSum = 0;
        double pValue, qValue;

        for (int y = 0; y < Height; y++)
        {
            for (int x = 0; x < Width; x++)
            {
                pValue = p[y, x];
                qValue = q[y, x];

                pSum += pValue;
                qSum += qValue;
                pSumSq += pValue * pValue;
                qSumSq += qValue * qValue;
                productSum += pValue * qValue;
            }
        }

        double numerator = productSum - ((pSum * qSum) / (double)Height);
        double denominator = Math.Sqrt((pSumSq - (pSum * pSum) / (double)Height) * (qSumSq - (qSum * qSum) / (double)Height));

        return (denominator == 0) ? 0 : numerator / denominator;
    }
}

Результат:

введите описание изображения здесь

Одна и та же картинка загружается в две графические коробки.

Значение их коэффициента корреляции стало, -1,

Это правильный результат?

Если нет, что я должен сделать, чтобы исправить это?

1 ответ

Решение

Соотношение идентичных последовательностей должно быть 1

Кажется, у вас есть проблема в самом конце рутины; вместо

 double numerator = productSum - ((pSum * qSum) / (double)Height);
 double denominator = Math.Sqrt((pSumSq - (pSum * pSum) / (double)Height) * (qSumSq - (qSum * qSum) / (double)Height));

ты должен положить

 double n = ((double) Width) * Height;

 double numerator = productSum - ((pSum * qSum) / n);
 double denominator = 
   Math.Sqrt((pSumSq - (pSum * pSum) / n) * (qSumSq - (qSum * qSum) / n));

увидеть

https://en.wikipedia.org/wiki/Pearson_product-moment_correlation_coefficient

Другие вопросы по тегам