Коэффициент корреляции Пирсона
Я написал следующий код 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