internal override double GeAssessment(AssesmentsSegment[] segments, IColorSheme colorSheme)
{
int k = segments.Length;
int n = segments[0].center.Length;
Matrix sumMatrix = new Matrix(n);
for (int i = 0; i < k; i++)
{
AssesmentsSegment current = segments[i];
sumMatrix += current.Length * GetMatrixOfCovariations(current, colorSheme);
}
return(sumMatrix.Determinant());
}
private Matrix GetMatrixOfCovariations(AssesmentsSegment segment, IColorSheme colorSheme)
{
Func <double, double, int, double> d = null;
if (colorSheme is GrayScaleColorSheme || colorSheme is RGBColorSheme)
{
d = (p1, p2, dim) => (p1 - p2) / 255d;
}
else if (colorSheme is LabColorSheme)
{
d = (p1, p2, dim) => dim == 0 ? (p1 - p2) / 100d : (p2 - p2) / 128d;
}
int n = segment.center.Length;
double[] averages = segment.center;
double[][] cov = ArrayMatrix.GetJaggedArray(n, n);
for (int k = 0; k < n; k++)
{
for (int p = 0; p < n; p++)
{
double v = 0;
for (int l = 0; l < segment.Length; l++)
{
//v += ((segment.points[l][k] - averages[k])) * ((segment.points[l][p] - averages[p]));
v += d(segment.points[l][k], averages[k], k) * d(segment.points[l][p], averages[p], p);
}
cov[k][p] = v / segment.Length;
}
}
return(new Matrix(cov));
}