aboutsummaryrefslogtreecommitdiffstats
path: root/libcolour-template.c
diff options
context:
space:
mode:
Diffstat (limited to 'libcolour-template.c')
-rw-r--r--libcolour-template.c731
1 files changed, 731 insertions, 0 deletions
diff --git a/libcolour-template.c b/libcolour-template.c
new file mode 100644
index 0000000..3862fce
--- /dev/null
+++ b/libcolour-template.c
@@ -0,0 +1,731 @@
+/* See LICENSE file for copyright and license details. */
+
+#if defined(__GNUC__) && !defined(__clang__)
+# pragma GCC diagnostic ignored "-Wfloat-equal"
+#elif defined(__clang__)
+# pragma clang diagnostic ignored "-Wcomma"
+# pragma clang diagnostic ignored "-Wfloat-equal"
+# pragma clang diagnostic ignored "-Wvla"
+# pragma clang diagnostic ignored "-Wtautological-compare"
+# pragma clang diagnostic ignored "-Wcovered-switch-default"
+# pragma clang diagnostic ignored "-Wfloat-conversion"
+# pragma clang diagnostic ignored "-Wconversion"
+# pragma clang diagnostic ignored "-Wdouble-promotion"
+# pragma clang diagnostic ignored "-Wswitch-enum"
+# pragma clang diagnostic ignored "-Wcast-align"
+#endif
+
+
+#define MARSHAL_VERSION 0
+
+
+TYPE
+libcolour_srgb_encode(TYPE t)
+{
+ TYPE sign = 1;
+ if (t < 0) {
+ t = -t;
+ sign = -1;
+ }
+ t = t <= D(0.0031306684425217108) ? D(12.92) * t :
+ D(1.055) * xpow(t, 1 / D(2.4)) - D(0.055);
+ return t * sign;
+}
+
+TYPE
+libcolour_srgb_decode(TYPE t)
+{
+ TYPE sign = 1;
+ if (t < 0) {
+ t = -t;
+ sign = -1;
+ }
+ t = t <= D(0.0031306684425217108) * D(12.92) ? t / D(12.92) :
+ xpow((t + D(0.055)) / D(1.055), D(2.4));
+ return t * sign;
+}
+
+
+int
+libcolour_delta_e(const libcolour_colour_t *a, const libcolour_colour_t *b, TYPE *e)
+{
+ libcolour_cielab_t u, v;
+ u.model = v.model = LIBCOLOUR_CIELAB;
+ if (libcolour_convert(a, (libcolour_colour_t *)&u))
+ return -1;
+ if (libcolour_convert(b, (libcolour_colour_t *)&v))
+ return -1;
+ u.L -= v.L, u.L *= u.L;
+ u.a -= v.a, u.a *= u.a;
+ u.b -= v.b, u.b *= u.b;
+ *e = xsqrt(u.L + u.a + u.b);
+ return 0;
+}
+
+
+#define eliminate(M, n, m) eliminate_(n, m, &(M))
+static int
+eliminate_(size_t n, size_t m, TYPE (*Mp)[n][m])
+{
+#define M (*Mp)
+ size_t r1, r2, c;
+ TYPE d;
+ TYPE *R1;
+ TYPE *R2;
+ for (r1 = 0; r1 < n; r1++) {
+ R1 = M[r1];
+ if (R1[r1] == 0) {
+ for (r2 = r1 + 1; r2 < n; r2++)
+ if (M[r2][r1] != 0)
+ break;
+ if (r2 == n) {
+ errno = EDOM;
+ return -1;
+ }
+ R2 = M[r2];
+ for (c = 0; c < m; c++)
+ d = R1[c], R1[c] = R2[c], R2[c] = d;
+ }
+ d = R1[r1];
+ for (c = 0; c < m; c++)
+ R1[c] /= d;
+ for (r2 = r1 + 1; r2 < n; r2++) {
+ R2 = M[r2];
+ d = R2[r1];
+ for (c = 0; c < m; c++)
+ R2[c] -= R1[c] * d;
+ }
+ }
+ for (r1 = n; --r1;) {
+ R1 = M[r1];
+ for (r2 = r1; r2--;) {
+ R2 = M[r2];
+ d = R2[r1];
+ for (c = 0; c < m; c++)
+ R2[c] -= R1[c] * d;
+ }
+ }
+ return 0;
+#undef M
+}
+
+
+int
+libcolour_proper(libcolour_colour_t *colour)
+{
+ libcolour_colour_t r, g, b, w;
+ TYPE m[3][4];
+ switch (colour->model) {
+ case LIBCOLOUR_CIELUV:
+ colour->cieluv.white.model = LIBCOLOUR_CIEXYZ;
+ break;
+ case LIBCOLOUR_CIELCHUV:
+ colour->cielchuv.white.model = LIBCOLOUR_CIEXYZ;
+ if (!colour->cielchuv.one_revolution)
+ colour->cielchuv.one_revolution = D(360.);
+ break;
+ case LIBCOLOUR_RGB:
+ colour->rgb.red.model = LIBCOLOUR_CIEXYY;
+ colour->rgb.green.model = LIBCOLOUR_CIEXYY;
+ colour->rgb.blue.model = LIBCOLOUR_CIEXYY;
+ colour->rgb.white.model = LIBCOLOUR_CIEXYY;
+ colour->rgb.red.Y = 1;
+ colour->rgb.green.Y = 1;
+ colour->rgb.blue.Y = 1;
+ r.model = g.model = b.model = w.model = LIBCOLOUR_CIEXYZ;
+ libcolour_convert((const libcolour_colour_t *)&colour->rgb.red, &r);
+ libcolour_convert((const libcolour_colour_t *)&colour->rgb.green, &g);
+ libcolour_convert((const libcolour_colour_t *)&colour->rgb.blue, &b);
+ libcolour_convert((const libcolour_colour_t *)&colour->rgb.white, &w);
+ m[0][0] = r.ciexyz.X, m[0][1] = g.ciexyz.X, m[0][2] = b.ciexyz.X, m[0][3] = w.ciexyz.X;
+ m[1][0] = r.ciexyz.Y, m[1][1] = g.ciexyz.Y, m[1][2] = b.ciexyz.Y, m[1][3] = w.ciexyz.Y;
+ m[2][0] = r.ciexyz.Z, m[2][1] = g.ciexyz.Z, m[2][2] = b.ciexyz.Z, m[2][3] = w.ciexyz.Z;
+ if (eliminate(m, 3, 4))
+ return -1;
+ colour->rgb.red.Y = m[0][3];
+ colour->rgb.green.Y = m[1][3];
+ colour->rgb.blue.Y = m[2][3];
+ break;
+ default:
+ break;
+ }
+ return 0;
+}
+
+
+static int
+get_primaries(libcolour_rgb_t *cs)
+{
+ libcolour_colour_t r, g, b, w;
+ TYPE Sr, Sg, Sb;
+ TYPE M[3][6];
+
+ r.model = g.model = b.model = w.model = LIBCOLOUR_CIEXYZ;
+ r.ciexyz.Y = 1;
+ g.ciexyz.Y = 1;
+ b.ciexyz.Y = 1;
+
+ Sr = cs->M[1][0] * cs->white_r;
+ Sg = cs->M[1][1] * cs->white_g;
+ Sb = cs->M[1][2] * cs->white_b;
+
+ r.ciexyz.X = cs->M[0][0] / Sr, r.ciexyz.Z = cs->M[2][0] / Sr;
+ g.ciexyz.X = cs->M[0][1] / Sg, g.ciexyz.Z = cs->M[2][1] / Sg;
+ b.ciexyz.X = cs->M[0][2] / Sb, b.ciexyz.Z = cs->M[2][2] / Sb;
+
+ M[0][0] = r.ciexyz.X, M[0][1] = g.ciexyz.X, M[0][2] = b.ciexyz.X, M[0][3] = 1, M[0][4] = 0, M[0][5] = 0;
+ M[1][0] = r.ciexyz.Y, M[1][1] = g.ciexyz.Y, M[1][2] = b.ciexyz.Y, M[1][3] = 0, M[1][4] = 1, M[1][5] = 0;
+ M[2][0] = r.ciexyz.Z, M[2][1] = g.ciexyz.Z, M[2][2] = b.ciexyz.Z, M[2][3] = 0, M[2][4] = 0, M[2][5] = 1;
+
+ if (eliminate(M, 3, 6))
+ return -1;
+
+ memcpy(M[0], M[0] + 3, 3 * sizeof(TYPE)), M[0][3] = Sr;
+ memcpy(M[1], M[1] + 3, 3 * sizeof(TYPE)), M[1][3] = Sg;
+ memcpy(M[2], M[2] + 3, 3 * sizeof(TYPE)), M[2][3] = Sb;
+
+ if (eliminate(M, 3, 4))
+ return -1;
+
+ w.ciexyz.X = M[0][3];
+ w.ciexyz.Y = M[1][3];
+ w.ciexyz.Z = M[2][3];
+
+ cs->red.model = LIBCOLOUR_CIEXYY;
+ cs->green.model = LIBCOLOUR_CIEXYY;
+ cs->blue.model = LIBCOLOUR_CIEXYY;
+ cs->white.model = LIBCOLOUR_CIEXYY;
+
+ libcolour_convert(&r, (libcolour_colour_t *)&cs->red);
+ libcolour_convert(&g, (libcolour_colour_t *)&cs->green);
+ libcolour_convert(&b, (libcolour_colour_t *)&cs->blue);
+ libcolour_convert(&w, (libcolour_colour_t *)&cs->white);
+
+ return 0;
+}
+
+
+static int
+get_matrices(libcolour_rgb_t *cs)
+{
+ libcolour_colour_t r, g, b, w;
+ TYPE M[3][6];
+ TYPE Sr, Sg, Sb;
+
+ r.model = g.model = b.model = LIBCOLOUR_CIEXYZ;
+ w.model = LIBCOLOUR_CIEXYY;
+ libcolour_convert((const libcolour_colour_t *)&cs->red, &w), w.ciexyy.Y = 1, libcolour_convert(&w, &r);
+ libcolour_convert((const libcolour_colour_t *)&cs->green, &w), w.ciexyy.Y = 1, libcolour_convert(&w, &g);
+ libcolour_convert((const libcolour_colour_t *)&cs->blue, &w), w.ciexyy.Y = 1, libcolour_convert(&w, &b);
+ w.model = LIBCOLOUR_CIEXYZ;
+ libcolour_convert((const libcolour_colour_t *)&cs->white, &w);
+
+ M[0][0] = r.ciexyz.X, M[0][1] = g.ciexyz.X, M[0][2] = b.ciexyz.X, M[0][3] = 1, M[0][4] = 0, M[0][5] = 0;
+ M[1][0] = r.ciexyz.Y, M[1][1] = g.ciexyz.Y, M[1][2] = b.ciexyz.Y, M[1][3] = 0, M[1][4] = 1, M[1][5] = 0;
+ M[2][0] = r.ciexyz.Z, M[2][1] = g.ciexyz.Z, M[2][2] = b.ciexyz.Z, M[2][3] = 0, M[2][4] = 0, M[2][5] = 1;
+
+ if (eliminate(M, 3, 6))
+ return -1;
+
+ Sr = M[0][3] * w.ciexyz.X + M[0][4] * w.ciexyz.Y + M[0][5] * w.ciexyz.Z;
+ Sg = M[1][3] * w.ciexyz.X + M[1][4] * w.ciexyz.Y + M[1][5] * w.ciexyz.Z;
+ Sb = M[2][3] * w.ciexyz.X + M[2][4] * w.ciexyz.Y + M[2][5] * w.ciexyz.Z;
+
+ Sr /= cs->white_r;
+ Sg /= cs->white_g;
+ Sb /= cs->white_b;
+
+ cs->M[0][0] = Sr * r.ciexyz.X, cs->M[0][1] = Sg * g.ciexyz.X, cs->M[0][2] = Sb * b.ciexyz.X;
+ cs->M[1][0] = Sr * r.ciexyz.Y, cs->M[1][1] = Sg * g.ciexyz.Y, cs->M[1][2] = Sb * b.ciexyz.Y;
+ cs->M[2][0] = Sr * r.ciexyz.Z, cs->M[2][1] = Sg * g.ciexyz.Z, cs->M[2][2] = Sb * b.ciexyz.Z;
+
+ memcpy(M[0], cs->M[0], 3 * sizeof(TYPE)), M[0][3] = 1, M[0][4] = 0, M[0][5] = 0;
+ memcpy(M[1], cs->M[1], 3 * sizeof(TYPE)), M[1][3] = 0, M[1][4] = 1, M[1][5] = 0;
+ memcpy(M[2], cs->M[2], 3 * sizeof(TYPE)), M[2][3] = 0, M[2][4] = 0, M[2][5] = 1;
+
+ if (eliminate(M, 3, 6))
+ return -1;
+
+ memcpy(cs->Minv[0], M[0] + 3, 3 * sizeof(TYPE));
+ memcpy(cs->Minv[1], M[1] + 3, 3 * sizeof(TYPE));
+ memcpy(cs->Minv[2], M[2] + 3, 3 * sizeof(TYPE));
+
+ return 0;
+}
+
+
+#define invert(Minv, M, n) invert_(n, &(Minv), &(M))
+static int
+invert_(size_t n, TYPE (*Minvp)[n][n], TYPE (*Mp)[n][n])
+{
+#define Minv (*Minvp)
+#define M (*Mp)
+#define J (*Jp)
+ TYPE J[n][2 * n] = alloca(sizeof(TYPE[n][2 * n]));
+ size_t i;
+ for (i = 0; i < n; i++) {
+ memcpy(J[i], M[i], n * sizeof(TYPE));
+ memset(J[i] + n, 0, n * sizeof(TYPE));
+ J[i][n + i] = 1;
+ }
+ if (eliminate(J, n, 2 * n))
+ return -1;
+ for (i = 0; i < n; i++)
+ memcpy(Minv[i], J[i] + n, n * sizeof(TYPE));
+ return 0;
+#undef J
+#undef M
+#undef Minv
+}
+
+
+static TYPE
+transfer_function_l_star(TYPE t)
+{
+ TYPE sign = 1;
+ if (t < 0) {
+ t = -t;
+ sign = -1;
+ }
+ t = t > D(216.) / D(24389.) ? D(1.16) * xcbrt(t) - D(0.16) : t * D(24389.) / D(2700.);
+ return t * sign;
+}
+
+static TYPE
+invtransfer_function_l_star(TYPE t)
+{
+ TYPE sign = 1;
+ if (t < 0) {
+ t = -t;
+ sign = -1;
+ }
+ t = t > D(0.08) ? (((D(1000000.) * t + D(480000.)) * t +
+ D(76800.)) * t + D(4096.)) / D(1560896.) :
+ t * D(2700.) / D(24389.);
+ return t * sign;
+}
+
+static TYPE
+transfer_function_oetf_hlg(TYPE t)
+{
+ TYPE sign = 1;
+ if (t < 0) {
+ t = -t;
+ sign = -1;
+ }
+ t = 12 * t <= 1 ? xsqrt(3 * t) :
+ D(0.17883277) * xlog(t - D(0.02372241)) + D(1.004293468902569985701234145381);
+ return t * sign;
+}
+
+static TYPE
+invtransfer_function_oetf_hlg(TYPE t)
+{
+ TYPE sign = 1;
+ if (t < 0) {
+ t = -t;
+ sign = -1;
+ }
+ t = t <= D(0.5) ? t * t / 3 :
+ xexp(t - D(1.004293468902569985701234145381)) / D(0.17883277) + D(0.02372241);
+ return t * sign;
+}
+
+
+static void
+get_transfer_function(libcolour_colour_t *cs)
+{
+ if (cs->model == LIBCOLOUR_RGB) {
+ switch (cs->rgb.colour_space) {
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2:
+ cs->rgb.to_encoded_red = cs->rgb.to_encoded_green = cs->rgb.to_encoded_blue = transfer_function_l_star;
+ cs->rgb.to_decoded_red = cs->rgb.to_decoded_green = cs->rgb.to_decoded_blue = invtransfer_function_l_star;
+ break;
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_PQ:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_PQ:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_PQ:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_HLG:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_HLG:
+ /* TODO http://www.itu.int/dms_pubrec/itu-r/rec/bt/R-REC-BT.2100-0-201607-I!!PDF-E.pdf */
+ break;
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_HLG:
+ cs->rgb.to_encoded_red = cs->rgb.to_encoded_green = cs->rgb.to_encoded_blue = transfer_function_oetf_hlg;
+ cs->rgb.to_decoded_red = cs->rgb.to_decoded_green = cs->rgb.to_decoded_blue = invtransfer_function_oetf_hlg;
+ break;
+ default:
+ break;
+ }
+ }
+}
+
+
+int
+libcolour_get_rgb_colour_space(libcolour_rgb_t *cs, libcolour_rgb_colour_space_t space)
+{
+#define XYY(XVALUE, YVALUE) (libcolour_ciexyy_t){.model = LIBCOLOUR_CIEXYY, .x = D(XVALUE), .y = D(YVALUE), .Y = 1}
+
+ switch (space) {
+ case LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MEASUREMENTS:
+ if (get_matrices(cs))
+ return -1;
+ if (cs->encoding_type == LIBCOLOUR_ENCODING_TYPE_REGULAR)
+ cs->transitioninv = cs->transition * cs->slope;
+ return 0;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MATRIX:
+ if (invert(cs->Minv, cs->M, 3) || get_primaries(cs))
+ return -1;
+ if (cs->encoding_type == LIBCOLOUR_ENCODING_TYPE_REGULAR)
+ cs->transitioninv = cs->transition * cs->slope;
+ return 0;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_INV_MATRIX:
+ if (invert(cs->M, cs->Minv, 3) || get_primaries(cs))
+ return -1;
+ if (cs->encoding_type == LIBCOLOUR_ENCODING_TYPE_REGULAR)
+ cs->transitioninv = cs->transition * cs->slope;
+ return 0;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_SRGB:
+ cs->red = XYY(0.6400, 0.3300);
+ cs->green = XYY(0.3000, 0.6000);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->offset = D(0.055);
+ cs->gamma = D(2.4);
+ cs->slope = D(12.92);
+ cs->transition = D(0.0031306684425217108);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ADOBE_RGB:
+ cs->red = XYY(0.6400, 0.3300);
+ cs->green = XYY(0.2100, 0.7100);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_APPLE_RGB:
+ cs->red = XYY(0.6250, 0.3400);
+ cs->green = XYY(0.2800, 0.5950);
+ cs->blue = XYY(0.1550, 0.0700);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(1.8);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_BEST_RGB:
+ cs->red = XYY(0.7347, 0.2653);
+ cs->green = XYY(0.2150, 0.7750);
+ cs->blue = XYY(0.1300, 0.0350);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_BETA_RGB:
+ cs->red = XYY(0.6888, 0.3112);
+ cs->green = XYY(0.1986, 0.7551);
+ cs->blue = XYY(0.1265, 0.0352);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_BRUCE_RGB:
+ cs->red = XYY(0.6400, 0.3300);
+ cs->green = XYY(0.2800, 0.6500);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_CIE_RGB:
+ cs->red = XYY(0.7350, 0.2650);
+ cs->green = XYY(0.2740, 0.7170);
+ cs->blue = XYY(0.1670, 0.0090);
+ cs->white = LIBCOLOUR_ILLUMINANT_E;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_COLORMATCH_RGB:
+ cs->red = XYY(0.6300, 0.3400);
+ cs->green = XYY(0.2950, 0.6050);
+ cs->blue = XYY(0.1500, 0.0750);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(1.8);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_D65:
+ cs->red = XYY(0.6380, 0.3200);
+ cs->green = XYY(0.2650, 0.6900);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.6);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_THEATER:
+ cs->red = XYY(0.6380, 0.3200);
+ cs->green = XYY(0.2650, 0.6900);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = XYY(0.314, 0.351);
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.6);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_DON_RGB_4:
+ cs->red = XYY(0.6960, 0.3000);
+ cs->green = XYY(0.2150, 0.7650);
+ cs->blue = XYY(0.1300, 0.0350);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB:
+ cs->red = XYY(0.6700, 0.3300);
+ cs->green = XYY(0.2100, 0.7100);
+ cs->blue = XYY(0.1400, 0.0800);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(1.8);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2:
+ cs->red = XYY(0.6700, 0.3300);
+ cs->green = XYY(0.2100, 0.7100);
+ cs->blue = XYY(0.1400, 0.0800);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_CUSTOM;
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_EKTA_SPACE_PS5:
+ cs->red = XYY(0.6950, 0.3050);
+ cs->green = XYY(0.2600, 0.7000);
+ cs->blue = XYY(0.1100, 0.0050);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.2);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_625_LINE:
+ cs->red = XYY(0.6400, 0.3300);
+ cs->green = XYY(0.2900, 0.6000);
+ cs->blue = XYY(0.1500, 0.0060);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_525_LINE:
+ cs->red = XYY(0.6300, 0.3400);
+ cs->green = XYY(0.3100, 0.5950);
+ cs->blue = XYY(0.1550, 0.0700);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_709:
+ cs->red = XYY(0.6300, 0.3300);
+ cs->green = XYY(0.3000, 0.6000);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2020:
+ cs->red = XYY(0.7080, 0.2920);
+ cs->green = XYY(0.1700, 0.7970);
+ cs->blue = XYY(0.1310, 0.0460);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_PQ:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_PQ:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_PQ:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_HLG:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_HLG:
+ case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_HLG:
+ cs->red = XYY(0.7080, 0.2920);
+ cs->green = XYY(0.1700, 0.7970);
+ cs->blue = XYY(0.1310, 0.0460);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_CUSTOM;
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_LIGHTROOM_RGB:
+ cs->red = XYY(0.7347, 0.2653);
+ cs->green = XYY(0.1596, 0.8404);
+ cs->blue = XYY(0.0366, 0.0001);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_LINEAR;
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_NTSC_RGB:
+ cs->red = XYY(0.6700, 0.3300);
+ cs->green = XYY(0.2100, 0.7100);
+ cs->blue = XYY(0.1400, 0.0800);
+ cs->white = LIBCOLOUR_ILLUMINANT_C;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_PAL_SECAM_RGB:
+ cs->red = XYY(0.6400, 0.3300);
+ cs->green = XYY(0.2900, 0.6000);
+ cs->blue = XYY(0.1500, 0.0600);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_PROPHOTO_RGB:
+ cs->red = XYY(0.7347, 0.2653);
+ cs->green = XYY(0.1596, 0.8404);
+ cs->blue = XYY(0.0366, 0.0001);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(1.8);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_SGI_RGB:
+ cs->red = XYY(0.6250, 0.3400);
+ cs->green = XYY(0.2800, 0.5950);
+ cs->blue = XYY(0.1550, 0.0700);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(1.47);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_240M_RGB:
+ cs->red = XYY(0.6300, 0.3400);
+ cs->green = XYY(0.3100, 0.5950);
+ cs->blue = XYY(0.1550, 0.0700);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.1115721957735072);
+ cs->slope = D(4.0);
+ cs->transition = D(0.022821585552393633);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_C_RGB:
+ cs->red = XYY(0.6300, 0.3400);
+ cs->green = XYY(0.3100, 0.5950);
+ cs->blue = XYY(0.1550, 0.0700);
+ cs->white = LIBCOLOUR_ILLUMINANT_D65;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_REGULAR;
+ cs->gamma = 1 / D(0.45);
+ cs->offset = D(0.09929682680944);
+ cs->slope = D(4.5);
+ cs->transition = D(0.018053968510807);
+ break;
+
+ case LIBCOLOUR_RGB_COLOUR_SPACE_WIDE_GAMUT_RGB:
+ cs->red = XYY(0.7350, 0.2650);
+ cs->green = XYY(0.1150, 0.8260);
+ cs->blue = XYY(0.1570, 0.0180);
+ cs->white = LIBCOLOUR_ILLUMINANT_D50;
+ cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
+ cs->gamma = D(2.19921875);
+ break;
+
+ default:
+ errno = EINVAL;
+ return -1;
+ }
+ if (cs->encoding_type == LIBCOLOUR_ENCODING_TYPE_REGULAR)
+ cs->transitioninv = cs->transition * cs->slope;
+ cs->colour_space = space;
+ cs->white_r = cs->white_g = cs->white_b = 1;
+ if (get_matrices(cs) || libcolour_proper((libcolour_colour_t *)cs))
+ return -1;
+ get_transfer_function((libcolour_colour_t *)cs);
+ return 0;
+
+#undef XYY
+}
+
+
+size_t
+libcolour_marshal(const libcolour_colour_t *colour, void *buf)
+{
+ if (buf)
+ *(int *)buf = MARSHAL_VERSION;
+ switch (colour->model) {
+#define X(C, T, N)\
+ case C:\
+ if (buf)\
+ memcpy((char *)buf + sizeof(int), colour, sizeof(T));\
+ return sizeof(int) + sizeof(T);
+ LIST_MODELS(X)
+#undef X
+ default:
+ errno = EINVAL;
+ return 0;
+ }
+}
+
+
+size_t
+libcolour_unmarshal(libcolour_colour_t *colour, const void *buf)
+{
+ enum libcolour_model model;
+ size_t r;
+ int ver;
+ ver = *(const int *)buf;
+ if (ver != MARSHAL_VERSION) {
+ errno = EINVAL;
+ return 0;
+ }
+ model = *(const enum libcolour_model *)((const char *)buf + sizeof(int));
+ switch (model) {
+#define X(C, T, N)\
+ case C:\
+ if (colour)\
+ memcpy(colour, (const char *)buf + sizeof(int), sizeof(T));\
+ r = sizeof(int) + sizeof(T);\
+ break;
+ LIST_MODELS(X)
+#undef X
+ default:
+ errno = EINVAL;
+ return 0;
+ }
+ if (colour) {
+ if (colour->model == LIBCOLOUR_RGB) {
+ colour->rgb.to_encoded_red = colour->rgb.to_encoded_green = colour->rgb.to_encoded_blue = NULL;
+ colour->rgb.to_decoded_red = colour->rgb.to_decoded_green = colour->rgb.to_decoded_blue = NULL;
+ }
+ get_transfer_function(colour);
+ }
+ return r;
+}
generated by cgit v1.2.3-70-g09d2 (git 2.46.0) at 2024-11-10 00:29:37 +0000