public static MyMatrix GenerateBradfordMatrix(ColorProfile from, ColorProfile to)
{
var ConeResponse = GetConeResponseMatrix();
var InverseConeResponse = GetInverseConeResponseMatrix();
double fromX;
double fromY;
double fromZ;
double toX;
double toY;
double toZ;
double m;
if (from.White_Y < 1e-5)
{
fromX = 0;
fromY = 0;
fromZ = 0;
}
else
{
fromY = 100;
m = fromY / from.White_Y;
fromX = from.White_X * m;
fromZ = from.White_Z * m;
}
if (to.White_Y < 1e-5)
{
toX = 0;
toY = 0;
toZ = 0;
}
else
{
toY = 100;
m = toY / to.White_Y;
toX = to.White_X * m;
toZ = to.White_Z * m;
}
var v1 = ConeResponse.MultiplyByVectorRight(fromX, fromY, fromZ);
var v2 = ConeResponse.MultiplyByVectorRight(toX, toY, toZ);
var middleMatrix = new MyMatrix();
middleMatrix[0, 0] = v2.X / v1.X;
middleMatrix[1, 1] = v2.Y / v1.Y;
middleMatrix[2, 2] = v2.Z / v1.Z;
var result = InverseConeResponse * middleMatrix * ConeResponse;
return(result);
}
public static void ConvertImage(SimpleColor[,] imageSource, SimpleColor[,] imageDest, ColorProfile from, ColorProfile to)
{
var bradfordMatrix = MyMatrix.GenerateBradfordMatrix(from, to);
Parallel.For(0, imageSource.GetLength(0), i =>
{
for (int j = 0; j < imageSource.GetLength(1); j++)
{
DoubleColor color = new DoubleColor(imageSource[i, j]);
var c1 = ConvertColorToXYZ(color, from);
var c2 = ConvertColorFromXYZ(c1, to, bradfordMatrix);
imageDest[i, j] = c2.ToSimpleColor();
}
});
}
private static MyMatrix GetInverseConeResponseMatrix()
{
var bm = new MyMatrix();
bm[0, 0] = 0.9870;
bm[0, 1] = -0.1471;
bm[0, 2] = 0.1600;
bm[1, 0] = 0.4323;
bm[1, 1] = 0.5184;
bm[1, 2] = 0.0493;
bm[2, 0] = -0.0085;
bm[2, 1] = 0.0400;
bm[2, 2] = 0.9685;
return(bm);
}
private static MyMatrix GetConeResponseMatrix()
{
var bm = new MyMatrix();
bm[0, 0] = 0.8951;
bm[0, 1] = 0.2664;
bm[0, 2] = -0.1614;
bm[1, 0] = -0.7502;
bm[1, 1] = 1.7135;
bm[1, 2] = 0.0367;
bm[2, 0] = 0.0389;
bm[2, 1] = -0.0685;
bm[2, 2] = 1.0296;
return(bm);
}
public static MyMatrix operator *(MyMatrix m1, MyMatrix m2)
{
MyMatrix result = new MyMatrix();
for (int i = 0; i < size; i++)
{
for (int j = 0; j < size; j++)
{
result[i, j] = 0;
for (int k = 0; k < size; k++)
{
result[i, j] += m1[i, k] * m2[k, j];
}
}
}
return(result);
}
public MyMatrix Invert()
{
MyMatrix inv = new MyMatrix();
double A = array[1, 1] * array[2, 2] - array[1, 2] * array[2, 1];
double B = array[1, 2] * array[2, 0] - array[1, 0] * array[2, 2];
double C = array[1, 0] * array[2, 1] - array[1, 1] * array[2, 0];
double det = array[0, 0] * A + array[0, 1] * B + array[0, 2] * C;
if (Math.Abs(det) < 1e-5)
{
return(null);
}
inv[0, 0] = A / det;
inv[1, 0] = B / det;
inv[2, 0] = C / det;
inv[0, 1] = (array[0, 2] * array[2, 1] - array[0, 1] * array[2, 2]) / det;
inv[1, 1] = (array[0, 0] * array[2, 2] - array[0, 2] * array[2, 0]) / det;
inv[2, 1] = (array[0, 1] * array[2, 0] - array[0, 0] * array[2, 1]) / det;
inv[0, 2] = (array[0, 1] * array[1, 2] - array[0, 2] * array[1, 1]) / det;
inv[1, 2] = (array[0, 2] * array[1, 0] - array[0, 0] * array[1, 2]) / det;
inv[2, 2] = (array[0, 0] * array[1, 1] - array[0, 1] * array[1, 0]) / det;
return(inv);
}
public static DoubleColor ConvertColorFromXYZ(DoubleColor color, ColorProfileEnum from, ColorProfileEnum to)
{
ColorProfile cp = null;
switch (to)
{
case ColorProfileEnum.sRGB:
{
cp = ColorProfileFactory.GetFactory().sRBGcolorProfile;
}
break;
case ColorProfileEnum.adobeRGB:
{
cp = ColorProfileFactory.GetFactory().adobeRBGcolorProfile;
}
break;
case ColorProfileEnum.appleRGB:
{
cp = ColorProfileFactory.GetFactory().appleRBGcolorProfile;
}
break;
case ColorProfileEnum.CIE_RGB:
{
cp = ColorProfileFactory.GetFactory().CIE_RBGcolorProfile;
}
break;
case ColorProfileEnum.WideGamut:
{
cp = ColorProfileFactory.GetFactory().WideGamutcolorProfile;
}
break;
case ColorProfileEnum.PAL:
{
cp = ColorProfileFactory.GetFactory().PALcolorProfile;
}
break;
default:
break;
}
//color.InverseGammaCorrection(cp.Gamma);
MyMatrix bradford = MyMatrix.GenerateBradfordMatrix(from, to);
if (bradford != null)
{
(color.R, color.G, color.B) = bradford.MultiplyByVectorRight(color.R, color.G, color.B);
}
var v = cp.XYZtoRGB.MultiplyByVectorRight(color.R, color.G, color.B);
DoubleColor XYZ = new DoubleColor(v.X, v.Y, v.Z);
XYZ.GammaCorrection(cp.Gamma);
return(XYZ);
//return new DoubleColor(-1, -1, -1);
}
public static DoubleColor ConvertColorFromXYZ(DoubleColor color, ColorProfile dest, MyMatrix bradfordMatrix = null)
{
if (bradfordMatrix != null)
{
(color.R, color.G, color.B) = bradfordMatrix.MultiplyByVectorRight(color.R, color.G, color.B);
}
var v = dest.XYZtoRGB.MultiplyByVectorRight(color.R, color.G, color.B);
DoubleColor XYZ = new DoubleColor(v.X, v.Y, v.Z);
XYZ.GammaCorrection(dest.Gamma);
return(XYZ);
}
public static MyMatrix GenerateBradfordMatrix(ColorProfileEnum from, ColorProfileEnum to)
{
//from is D65
if (from == ColorProfileEnum.sRGB || from == ColorProfileEnum.adobeRGB || from == ColorProfileEnum.appleRGB || from == ColorProfileEnum.PAL)
{
//to is D65
if (to == ColorProfileEnum.sRGB || to == ColorProfileEnum.adobeRGB || to == ColorProfileEnum.appleRGB || to == ColorProfileEnum.PAL)
{
return(null);
}
//to is D50
else if (to == ColorProfileEnum.WideGamut)
{
var bm = new MyMatrix();
bm[0, 0] = 1.0478;
bm[0, 1] = 0.0229;
bm[0, 2] = -0.0501;
bm[1, 0] = 0.0295;
bm[1, 1] = 0.9905;
bm[1, 2] = -0.0171;
bm[2, 0] = -0.0092;
bm[2, 1] = 0.0150;
bm[2, 2] = 0.7521;
return(bm);
}
//to is E
else if (to == ColorProfileEnum.CIE_RGB)
{
var bm = new MyMatrix();
bm[0, 0] = 1.0503;
bm[0, 1] = 0.0271;
bm[0, 2] = -0.0233;
bm[1, 0] = 0.0391;
bm[1, 1] = 0.9730;
bm[1, 2] = -0.0093;
bm[2, 0] = -0.0024;
bm[2, 1] = 0.0026;
bm[2, 2] = 0.9181;
return(bm);
}
}
//from is D50
else if (from == ColorProfileEnum.WideGamut)
{
//to is D65
if (to == ColorProfileEnum.sRGB || to == ColorProfileEnum.adobeRGB || to == ColorProfileEnum.appleRGB || to == ColorProfileEnum.PAL)
{
var bm = new MyMatrix();
bm[0, 0] = 0.9556;
bm[0, 1] = -0.0230;
bm[0, 2] = 0.0632;
bm[1, 0] = -0.0283;
bm[1, 1] = 1.0099;
bm[1, 2] = 0.0210;
bm[2, 0] = 0.0123;
bm[2, 1] = -0.0205;
bm[2, 2] = 1.3299;
return(bm);
}
//to is D50
else if (to == ColorProfileEnum.WideGamut)
{
return(null);
}
//to is E
else if (to == ColorProfileEnum.CIE_RGB)
{
var bm = new MyMatrix();
bm[0, 0] = 1.0026;
bm[0, 1] = 0.0036;
bm[0, 2] = 0.0360;
bm[1, 0] = 0.0097;
bm[1, 1] = 0.9819;
bm[1, 2] = 0.0106;
bm[2, 0] = 0.0089;
bm[2, 1] = -0.0161;
bm[2, 2] = 1.2209;
return(bm);
}
}
//from is E
else if (from == ColorProfileEnum.CIE_RGB)
{
//to is D65
if (to == ColorProfileEnum.sRGB || to == ColorProfileEnum.adobeRGB || to == ColorProfileEnum.appleRGB || to == ColorProfileEnum.PAL)
{
var bm = new MyMatrix();
bm[0, 0] = 0.9532;
bm[0, 1] = -0.0266;
bm[0, 2] = 0.0239;
bm[1, 0] = -0.0382;
bm[1, 1] = 1.0288;
bm[1, 2] = 0.0094;
bm[2, 0] = 0.0026;
bm[2, 1] = -0.0030;
bm[2, 2] = 1.0893;
return(bm);
}
//to is D50
else if (to == ColorProfileEnum.WideGamut)
{
var bm = new MyMatrix();
bm[0, 0] = 0.9978;
bm[0, 1] = -0.0042;
bm[0, 2] = -0.0294;
bm[1, 0] = -0.0098;
bm[1, 1] = 1.0183;
bm[1, 2] = -0.0085;
bm[2, 0] = -0.0074;
bm[2, 1] = 0.0134;
bm[2, 2] = 0.8192;
return(bm);
}
//to is E
else if (to == ColorProfileEnum.CIE_RGB)
{
return(null);
}
}
return(null);
}