/** * Copyright © 2016 Mattias Andrée * * This program is free software: you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by * the Free Software Foundation, either version 3 of the License, or * (at your option) any later version. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with this program. If not, see . */ #include "libcolour.h" #include #include #include #include #include #include #define MARSHAL_VERSION 0 #define WASDIV0(x) (isinf(x) || isnan(x)) static void to_rgb(const libcolour_colour_t* restrict from, libcolour_rgb_t* restrict to); static void to_srgb(const libcolour_colour_t* restrict from, libcolour_srgb_t* restrict to); static void to_ciexyy(const libcolour_colour_t* restrict from, libcolour_ciexyy_t* restrict to); static void to_ciexyz(const libcolour_colour_t* restrict from, libcolour_ciexyz_t* restrict to); static void to_cielab(const libcolour_colour_t* restrict from, libcolour_cielab_t* restrict to); static void to_cieluv(const libcolour_colour_t* restrict from, libcolour_cieluv_t* restrict to); static void to_cielchuv(const libcolour_colour_t* restrict from, libcolour_cielchuv_t* restrict to); static void to_yiq(const libcolour_colour_t* restrict from, libcolour_yiq_t* restrict to); static void to_ydbdr(const libcolour_colour_t* restrict from, libcolour_ydbdr_t* restrict to); static void to_yuv(const libcolour_colour_t* restrict from, libcolour_yuv_t* restrict to); static void to_ypbpr(const libcolour_colour_t* restrict from, libcolour_ypbpr_t* restrict to); static void to_ycgco(const libcolour_colour_t* restrict from, libcolour_ycgco_t* restrict to); static void to_cie1960ucs(const libcolour_colour_t* restrict from, libcolour_cie1960ucs_t* restrict to); static void to_cieuvw(const libcolour_colour_t* restrict from, libcolour_cieuvw_t* restrict to); static void ciexyz_to_rgb(const libcolour_ciexyz_t* restrict from, libcolour_rgb_t* restrict to) { double X = from->X, Y = from->Y, Z = from->Z; to->R = to->Minv[0][0] * X + to->Minv[0][1] * Y + to->Minv[0][2] * Z; to->G = to->Minv[1][0] * X + to->Minv[1][1] * Y + to->Minv[1][2] * Z; to->B = to->Minv[2][0] * X + to->Minv[2][1] * Y + to->Minv[2][2] * Z; } static void to_rgb(const libcolour_colour_t* restrict from, libcolour_rgb_t* restrict to) { int have_gamma = 0, with_gamma = to->with_gamma; libcolour_ciexyz_t tmp; switch (from->model) { case LIBCOLOUR_CIEXYZ: ciexyz_to_rgb(&from->ciexyz, to); break; case LIBCOLOUR_RGB: if (!memcmp(from->rgb.M, to->M, sizeof(double[3][3]))) { have_gamma = from->rgb.with_gamma; *to = from->rgb; break; } /* fall through */ default: tmp.model = LIBCOLOUR_CIEXYZ; to_ciexyz(from, &tmp); ciexyz_to_rgb(&tmp, to); break; } if (have_gamma != with_gamma) { if (with_gamma) { switch (to->encoding_type) { case LIBCOLOUR_ENCODING_TYPE_LINEAR: break; case LIBCOLOUR_ENCODING_TYPE_SIMPLE: to->R = pow(to->R, 1 / to->gamma); to->G = pow(to->G, 1 / to->gamma); to->B = pow(to->B, 1 / to->gamma); break; case LIBCOLOUR_ENCODING_TYPE_REGULAR: to->R <= to->transition ? to->slope * to->R : (1 + to->offset) * pow(to->R, 1 / to->gamma) - to->offset; to->G <= to->transition ? to->slope * to->G : (1 + to->offset) * pow(to->G, 1 / to->gamma) - to->offset; to->B <= to->transition ? to->slope * to->B : (1 + to->offset) * pow(to->B, 1 / to->gamma) - to->offset; break; case LIBCOLOUR_ENCODING_TYPE_CUSTOM: to->R = (to->to_encoded_red)(to->R); to->G = (to->to_encoded_green)(to->G); to->B = (to->to_encoded_blue)(to->B); break; default: fprintf(stderr, "libcolour: invalid encoding type\n"); abort(); } } else { switch (to->encoding_type) { case LIBCOLOUR_ENCODING_TYPE_LINEAR: break; case LIBCOLOUR_ENCODING_TYPE_SIMPLE: to->R = pow(to->R, to->gamma); to->G = pow(to->G, to->gamma); to->B = pow(to->B, to->gamma); break; case LIBCOLOUR_ENCODING_TYPE_REGULAR: to->R <= to->transitioninv ? to->R * to->slope : pow((to->R + to->offset) / (1 + to->offset), to->gamma); to->G <= to->transitioninv ? to->G * to->slope : pow((to->G + to->offset) / (1 + to->offset), to->gamma); to->B <= to->transitioninv ? to->B * to->slope : pow((to->B + to->offset) / (1 + to->offset), to->gamma); break; case LIBCOLOUR_ENCODING_TYPE_CUSTOM: to->R = (to->to_decoded_red)(to->R); to->G = (to->to_decoded_green)(to->G); to->B = (to->to_decoded_blue)(to->B); break; default: fprintf(stderr, "libcolour: invalid encoding type\n"); abort(); } } } } static void ciexyz_to_srgb(const libcolour_ciexyz_t* restrict from, libcolour_srgb_t* restrict to) { double X = from->X, Y = from->Y, Z = from->Z; to->R = 3.240450 * X + -1.537140 * Y + -0.4985320 * Z; to->G = -0.969266 * X + 1.876010 * Y + 0.0415561 * Z; to->B = 0.0556434 * X + -0.204026 * Y + 1.0572300 * Z; } static void srgb_to_srgb(const libcolour_srgb_t* restrict from, libcolour_srgb_t* restrict to) { if (from->with_gamma == to->with_gamma) { *to = *from; } else if (to->with_gamma) { to->R = libcolour_srgb_encode(from->R); to->G = libcolour_srgb_encode(from->G); to->B = libcolour_srgb_encode(from->B); } else { to->R = libcolour_srgb_decode(from->R); to->G = libcolour_srgb_decode(from->G); to->B = libcolour_srgb_decode(from->B); } } static void yiq_to_srgb(const libcolour_yiq_t* restrict from, libcolour_srgb_t* restrict to) { double Y = from->Y, I = from->I, Q = from->Q; to->R = Y + I * 0.956 + Q * 0.621; to->G = Y - I * 0.272 - Q * 0.647; to->B = Y - I * 1.106 + Q * 1.703; } static void ydbdr_to_srgb(const libcolour_ydbdr_t* restrict from, libcolour_srgb_t* restrict to) { double Y = from->Y, Db = from->Db, Dr = from->Dr; to->R = Y + Db * 0.000092303716148 - Dr * 0.525912630661865; to->G = Y - Db * 0.129132898890509 + Dr * 0.267899328207599; to->B = Y + Db * 0.664679059978955 - Dr * 0.000079202543533; } static void ypbpr_to_srgb(const libcolour_ypbpr_t* restrict from, libcolour_srgb_t* restrict to) { double Y = from->Y, Pb = from->Pb, Pr = from->Pr; to->R = Pr + Y; to->B = Pb + Y; to->G = (Y - to->R * 0.2126 - to->B * 0.0722) / 0.7152; } static void ycgco_to_srgb(const libcolour_ycgco_t* restrict from, libcolour_srgb_t* restrict to) { double Y = from->Y, Cg = from->Cg, Co = from->Co; to->R = Y - Cg + Co; to->G = Y + Cg; to->B = Y - Cg - Co; } static void other_to_srgb(const libcolour_colour_t* restrict from, libcolour_srgb_t* restrict to) { libcolour_ciexyz_t tmp; tmp.model = LIBCOLOUR_CIEXYZ; to_ciexyz(from, &tmp); ciexyz_to_srgb(&tmp, to); } static void to_srgb(const libcolour_colour_t* restrict from, libcolour_srgb_t* restrict to) { libcolour_srgb_t tmp; switch (from->model) { case LIBCOLOUR_CIEXYZ: ciexyz_to_srgb(&from->ciexyz, to); break; case LIBCOLOUR_SRGB: srgb_to_srgb(&from->srgb, to); return; case LIBCOLOUR_YIQ: yiq_to_srgb(&from->yiq, to); break; case LIBCOLOUR_YDBDR: ydbdr_to_srgb(&from->ydbdr, to); break; case LIBCOLOUR_YPBPR: ypbpr_to_srgb(&from->ypbpr, to); break; case LIBCOLOUR_YCGCO: ycgco_to_srgb(&from->ycgco, to); break; default: other_to_srgb(from, to); break; } if (to->with_gamma) { tmp = *to; tmp.with_gamma = 0; srgb_to_srgb(&tmp, to); } } static void ciexyz_to_ciexyy(const libcolour_ciexyz_t* restrict from, libcolour_ciexyy_t* restrict to) { double X = from->X, Y = from->Y, Z = from->Z; double s = X + Y + Z; to->x = X / s; to->y = Y / s; to->Y = Y; if (WASDIV0(to->x) || WASDIV0(to->y)) to->x = to->y = 0; } static void other_to_ciexyy(const libcolour_colour_t* restrict from, libcolour_ciexyy_t* restrict to) { libcolour_ciexyz_t tmp; tmp.model = LIBCOLOUR_CIEXYZ; to_ciexyz(from, &tmp); ciexyz_to_ciexyy(&tmp, to); } static void srgb_to_ciexyy(const libcolour_srgb_t* restrict from, libcolour_ciexyy_t* restrict to) { libcolour_srgb_t tmp; tmp.model = LIBCOLOUR_SRGB; if (from->with_gamma) { tmp.with_gamma = 0; to_srgb((const libcolour_colour_t*)from, &tmp); } else { tmp = *from; } if (tmp.R == 0 && tmp.G == 0 && tmp.B == 0) { to->x = 0.31272660439158; to->y = 0.32902315240275; to->Y = 0; } else { other_to_ciexyy((const libcolour_colour_t*)&tmp, to); } } static void to_ciexyy(const libcolour_colour_t* restrict from, libcolour_ciexyy_t* restrict to) { switch (from->model) { case LIBCOLOUR_CIEXYY: *to = from->ciexyy; return; case LIBCOLOUR_CIEXYZ: ciexyz_to_ciexyy(&from->ciexyz, to); return; case LIBCOLOUR_SRGB: srgb_to_ciexyy(&from->srgb, to); return; default: other_to_ciexyy(from, to); return; } } static void rgb_to_ciexyz(const libcolour_rgb_t* restrict from, libcolour_ciexyz_t* restrict to) { libcolour_rgb_t tmp; double R, G, B; if (from->with_gamma) { tmp = *from; tmp.with_gamma = 0; to_rgb((const libcolour_colour_t*)from, &tmp); R = tmp.R, G = tmp.G, B = tmp.B; } else { R = from->R, G = from->G, B = from->B; } to->X = from->M[0][0] * R + from->M[0][1] * G + from->M[0][2] * B; to->Y = from->M[1][0] * R + from->M[1][1] * G + from->M[1][2] * B; to->Z = from->M[2][0] * R + from->M[2][1] * G + from->M[2][2] * B; } static void srgb_to_ciexyz(const libcolour_srgb_t* restrict from, libcolour_ciexyz_t* restrict to) { libcolour_srgb_t tmp; double R, G, B; if (from->with_gamma) { tmp.model = LIBCOLOUR_SRGB; tmp.with_gamma = 0; to_srgb((const libcolour_colour_t*)from, &tmp); R = tmp.R, G = tmp.G, B = tmp.B; } else { R = from->R, G = from->G, B = from->B; } to->X = 0.4124564 * R + 0.3575761 * G + 0.1804375 * B; to->Y = 0.2126729 * R + 0.7151522 * G + 0.0721750 * B; to->Z = 0.0193339 * R + 0.1191920 * G + 0.9503041 * B; } static void ciexyy_to_ciexyz(const libcolour_ciexyy_t* restrict from, libcolour_ciexyz_t* restrict to) { double x = from->x, y = from->y, Y = from->Y; double Yy = Y / y; if (WASDIV0(Yy)) { to->X = to->Y = to->Z = Y; } else { to->X = x * Yy; to->Y = Y; to->Z = (1 - x - y) * Yy; } } static inline double cielab_finv(double t) { return (t > 6. / 29.) ? t * t * t : (t - 4. / 29.) * 108 / 841; } static void cielab_to_ciexyz(const libcolour_cielab_t* restrict from, libcolour_ciexyz_t* restrict to) { double L = from->L, a = from->a, b = from->b; to->Y = (L + 16) / 116; to->X = to->Y + a / 500; to->Z = to->Y - b / 200; to->X = cielab_finv(to->X) * 0.95047; to->Y = cielab_finv(to->Y) * 1.08883; to->Z = cielab_finv(to->Z); } static void cieluv_to_ciexyz(const libcolour_cieluv_t* restrict from, libcolour_ciexyz_t* restrict to) { double X = from->white.X, Y = from->white.Y, L = from->L, L13 = from->L * 13; double t = X + 15 * Y + 3 * from->white.Z; double u = from->u / L13 + 4 * X / t; double v = from->v / L13 + 9 * Y / t; if (L <= 8) { Y *= L * 27 / 24389; } else { L = (L + 16) / 116; Y *= L * L * L; } to->X = 2.25 * Y * u / v; to->Y = Y; to->Z = Y * (3 / v - 0.75 * u / v - 5); } static void cielchuv_to_ciexyz(const libcolour_cielchuv_t* restrict from, libcolour_ciexyz_t* restrict to) { libcolour_cieluv_t tmp; tmp.model = LIBCOLOUR_CIELUV; tmp.white = from->white; to_cieluv((const libcolour_colour_t*)from, &tmp); cieluv_to_ciexyz(&tmp, to); } static void yuv_to_ciexyz(const libcolour_yuv_t* restrict from, libcolour_ciexyz_t* restrict to) { libcolour_ydbdr_t tmp; tmp.model = LIBCOLOUR_YDBDR; to_ydbdr((const libcolour_colour_t*)from, &tmp); to_ciexyz((const libcolour_colour_t*)&tmp, to); } static void cie1960ucs_to_ciexyz(const libcolour_cie1960ucs_t* restrict from, libcolour_ciexyz_t* restrict to) { double u = from->u, v = from->v, Y = from->Y; to->X = 1.5 * Y * u / v; to->Y = Y; to->Z = (4 * Y - Y * u - 10 * Y * v) / (2 * v); } static void cieuvw_to_ciexyz(const libcolour_cieuvw_t* restrict from, libcolour_ciexyz_t* restrict to) { libcolour_cie1960ucs_t tmp; tmp.model = LIBCOLOUR_CIE1960UCS; to_cie1960ucs((const libcolour_colour_t*)from, &tmp); cie1960ucs_to_ciexyz(&tmp, to); } static void other_to_ciexyz(const libcolour_colour_t* restrict from, libcolour_ciexyz_t* restrict to) { libcolour_srgb_t tmp; tmp.model = LIBCOLOUR_SRGB; tmp.with_gamma = 0; to_srgb(from, &tmp); srgb_to_ciexyz(&tmp, to); } static void to_ciexyz(const libcolour_colour_t* restrict from, libcolour_ciexyz_t* restrict to) { switch (from->model) { case LIBCOLOUR_RGB: rgb_to_ciexyz(&from->rgb, to); return; case LIBCOLOUR_SRGB: srgb_to_ciexyz(&from->srgb, to); return; case LIBCOLOUR_CIEXYY: ciexyy_to_ciexyz(&from->ciexyy, to); return; case LIBCOLOUR_CIEXYZ: *to = from->ciexyz; return; case LIBCOLOUR_CIELAB: cielab_to_ciexyz(&from->cielab, to); return; case LIBCOLOUR_CIELUV: cieluv_to_ciexyz(&from->cieluv, to); return; case LIBCOLOUR_CIELCHUV: cielchuv_to_ciexyz(&from->cielchuv, to); return; case LIBCOLOUR_YUV: yuv_to_ciexyz(&from->yuv, to); return; case LIBCOLOUR_CIE1960UCS: cie1960ucs_to_ciexyz(&from->cie1960ucs, to); return; case LIBCOLOUR_CIEUVW: cieuvw_to_ciexyz(&from->cieuvw, to); return; default: other_to_ciexyz(from, to); return; } } static inline double cielab_f(double t) { return (t > 216. / 24389.) ? cbrt(t) : t * 841. / 108. + 4. / 29.; } static void ciexyz_to_cielab(const libcolour_ciexyz_t* restrict from, libcolour_cielab_t* restrict to) { double X = from->X, Y = from->Y, Z = from->Z; X /= 0.95047; Y /= 1.08883; Y = cielab_f(Y); to->L = 116 * Y - 16; to->a = 500 * (cielab_f(X) - Y); to->b = 200 * (Y - cielab_f(Z)); } static void other_to_cielab(const libcolour_colour_t* restrict from, libcolour_cielab_t* restrict to) { libcolour_ciexyz_t tmp; tmp.model = LIBCOLOUR_CIEXYZ; to_ciexyz(from, &tmp); ciexyz_to_cielab(&tmp, to); } static void to_cielab(const libcolour_colour_t* restrict from, libcolour_cielab_t* restrict to) { switch (from->model) { case LIBCOLOUR_CIEXYZ: ciexyz_to_cielab(&from->ciexyz, to); return; case LIBCOLOUR_CIELAB: *to = from->cielab; return; default: other_to_cielab(from, to); return; } } static void ciexyz_to_cieluv(const libcolour_ciexyz_t* restrict from, libcolour_cieluv_t* restrict to) { double t, u, v, y, y2; t = to->white.X + 15 * to->white.Y + 3 * to->white.Z; u = 4 * to->white.X / t; v = 9 * to->white.Y / t; t = from->X + 15 * from->Y + 3 * from->Z; u = 4 * from->X / t - u; v = 9 * from->Y / t - v; y = from->Y / to->white.Y; y2 = y * 24389; y = y2 <= 216 ? y2 / 27 : cbrt(y) * 116 - 16; to->L = y; y *= 13; to->u = u * y; to->v = v * y; } static void cielchuv_to_cieluv(const libcolour_cielchuv_t* restrict from, libcolour_cieluv_t* restrict to) { libcolour_ciexyz_t tmp; libcolour_cielchuv_t tmp2; double L, C, h; if (to->white.X != from->white.X || to->white.Y != from->white.Y || to->white.Z != from->white.Z) { tmp.model = LIBCOLOUR_CIEXYZ; tmp2.model = LIBCOLOUR_CIELCHUV; tmp2.white = to->white; to_ciexyz((const libcolour_colour_t*)from, &tmp); to_cielchuv((const libcolour_colour_t*)&tmp, &tmp2); L = tmp2.L, C = tmp2.C, h = tmp2.h; } else { L = from->L, C = from->C, h = from->h; } to->L = L; to->u = C * cos(h); to->v = C * sin(h); } static void other_to_cieluv(const libcolour_colour_t* restrict from, libcolour_cieluv_t* restrict to) { libcolour_ciexyz_t tmp; tmp.model = LIBCOLOUR_CIEXYZ; to_ciexyz(from, &tmp); ciexyz_to_cieluv(&tmp, to); } static void to_cieluv(const libcolour_colour_t* restrict from, libcolour_cieluv_t* restrict to) { switch (from->model) { case LIBCOLOUR_CIEXYZ: ciexyz_to_cieluv(&from->ciexyz, to); return; case LIBCOLOUR_CIELCHUV: cielchuv_to_cieluv(&from->cielchuv, to); return; case LIBCOLOUR_CIELUV: if (to->white.X == from->cieluv.white.X && to->white.Y == from->cieluv.white.Y && to->white.Z == from->cieluv.white.Z) { *to = from->cieluv; return; } /* fall through */ default: other_to_cieluv(from, to); return; } } static void cieluv_to_cielchuv(const libcolour_cieluv_t* restrict from, libcolour_cielchuv_t* restrict to) { libcolour_cieluv_t tmp; double L, u, v; if (to->white.X != from->white.X || to->white.Y != from->white.Y || to->white.Z != from->white.Z) { tmp.model = LIBCOLOUR_CIELUV; tmp.white = to->white; to_cieluv((const libcolour_colour_t*)from, &tmp); L = tmp.L, u = tmp.u, v = tmp.v; } else { L = from->L, u = from->u, v = from->v; } to->L = L; to->C = sqrt(u * u + v * v); to->h = atan2(v, u); } static void other_to_cielchuv(const libcolour_colour_t* restrict from, libcolour_cielchuv_t* restrict to) { libcolour_cieluv_t tmp; tmp.model = LIBCOLOUR_CIELUV; tmp.white = to->white; to_cieluv(from, &tmp); cieluv_to_cielchuv(&tmp, to); } static void to_cielchuv(const libcolour_colour_t* restrict from, libcolour_cielchuv_t* restrict to) { switch (from->model) { case LIBCOLOUR_CIELUV: cieluv_to_cielchuv(&from->cieluv, to); return; case LIBCOLOUR_CIELCHUV: if (to->white.X == from->cielchuv.white.X && to->white.Y == from->cielchuv.white.Y && to->white.Z == from->cielchuv.white.Z) { *to = from->cielchuv; return; } /* fall through */ default: other_to_cielchuv(from, to); return; } } static void to_yiq(const libcolour_colour_t* restrict from, libcolour_yiq_t* restrict to) { /* TODO to_yiq: add direct conversion from and to CIE XYZ and use it as the default intermediary */ double r, g, b; libcolour_colour_t tmp = *from; switch (from->model) { case LIBCOLOUR_YIQ: *to = from->yiq; return; default: tmp.model = LIBCOLOUR_SRGB; tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); /* fall through */ case LIBCOLOUR_SRGB: if (tmp.srgb.with_gamma) { tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); } r = tmp.srgb.R; g = tmp.srgb.G; b = tmp.srgb.B; to->Y = r * 0.29893602129377540 + g * 0.5870430744511212 + b * 0.11402090425510336; to->I = r * 0.59594574307079930 - g * 0.2743886357457892 - b * 0.32155710732501010; to->Q = r * 0.21149734030682846 - g * 0.5229106903029739 + b * 0.31141334999614540; return; } } static void to_ydbdr(const libcolour_colour_t* restrict from, libcolour_ydbdr_t* restrict to) { /* TODO to_ydbdr: add direct conversion from and to CIE XYZ and use it as the default intermediary */ double r, g, b; libcolour_colour_t tmp = *from; switch (from->model) { case LIBCOLOUR_YDBDR: *to = from->ydbdr; return; case LIBCOLOUR_YUV: to->Y = tmp.yuv.Y; to->Db = tmp.yuv.U * 3.069; to->Dr = -tmp.yuv.V * 2.169; return; default: tmp.model = LIBCOLOUR_SRGB; tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); /* fall through */ case LIBCOLOUR_SRGB: if (tmp.srgb.with_gamma) { tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); } r = tmp.srgb.R; g = tmp.srgb.G; b = tmp.srgb.B; to->Y = r * 0.299 + g * 0.587 + b * 0.114; to->Db = -r * 0.450 - g * 0.883 + b * 1.333; to->Dr = -r * 1.333 + g * 1.116 + b * 0.217; return; } } static void to_yuv(const libcolour_colour_t* restrict from, libcolour_yuv_t* restrict to) { /* TODO to_yuv: add direct conversion from and to CIE XYZ and use it as the default intermediary */ /* TODO to_yub: add direct conversion from and to sRGB */ libcolour_colour_t tmp = *from; switch (from->model) { case LIBCOLOUR_YUV: *to = from->yuv; return; default: tmp.model = LIBCOLOUR_YDBDR; to_ydbdr(from, &tmp.ydbdr); /* fall through */ case LIBCOLOUR_YDBDR: to->Y = tmp.ydbdr.Y; to->U = tmp.ydbdr.Db / 3.069; to->V = -tmp.ydbdr.Dr / 2.169; return; } } static void to_ypbpr(const libcolour_colour_t* restrict from, libcolour_ypbpr_t* restrict to) { libcolour_colour_t tmp = *from; switch (from->model) { case LIBCOLOUR_YPBPR: *to = from->ypbpr; return; default: tmp.model = LIBCOLOUR_SRGB; tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); /* fall through */ case LIBCOLOUR_SRGB: if (tmp.srgb.with_gamma) { tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); } to->Y = tmp.srgb.R * 0.2126 + tmp.srgb.G * 0.7152 + tmp.srgb.B * 0.0722; to->Pb = tmp.srgb.B - to->Y; to->Pr = tmp.srgb.R - to->Y; return; } } static void to_ycgco(const libcolour_colour_t* restrict from, libcolour_ycgco_t* restrict to) { /* TODO to_ycgco: add direct conversion from and to CIE XYZ and use it as the default intermediary */ libcolour_colour_t tmp = *from; switch (from->model) { case LIBCOLOUR_YCGCO: *to = from->ycgco; return; default: tmp.model = LIBCOLOUR_SRGB; tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); /* fall through */ case LIBCOLOUR_SRGB: if (tmp.srgb.with_gamma) { tmp.srgb.with_gamma = 0; to_srgb(from, &tmp.srgb); } to->Y = tmp.srgb.R / 4 + tmp.srgb.G / 2 + tmp.srgb.B / 4; to->Cg = -tmp.srgb.R / 4 + tmp.srgb.G / 2 - tmp.srgb.B / 4; to->Co = tmp.srgb.R / 2 - tmp.srgb.B / 2; return; } } static void to_cie1960ucs(const libcolour_colour_t* restrict from, libcolour_cie1960ucs_t* restrict to) { libcolour_colour_t tmp = *from; double u, v, w, y; switch (from->model) { case LIBCOLOUR_CIE1960UCS: *to = from->cie1960ucs; return; case LIBCOLOUR_CIEUVW: u = from->cieuvw.U, v = from->cieuvw.V, w = from->cieuvw.W; y = (w + 17) / 25; y *= y * y; w *= 13; to->u = u / w + from->cieuvw.u0; to->v = v / w + from->cieuvw.v0; to->Y = y; return; default: tmp.model = LIBCOLOUR_CIEXYZ; to_ciexyz(from, &tmp.ciexyz); /* fall through */ case LIBCOLOUR_CIEXYZ: w = tmp.ciexyz.X + 15 * tmp.ciexyz.Y + 3 * tmp.ciexyz.Z; to->u = 4 * tmp.ciexyz.X / w; to->v = 6 * tmp.ciexyz.Y / w; to->Y = tmp.ciexyz.Y; return; } } static void to_cieuvw(const libcolour_colour_t* restrict from, libcolour_cieuvw_t* restrict to) { double U, V, W, w; libcolour_colour_t tmp = *from; switch (from->model) { case LIBCOLOUR_CIEUVW: U = from->cieuvw.U, V = from->cieuvw.V, W = from->cieuvw.W; w = W * 13; U += w * (from->cieuvw.u0 - to->u0); V += w * (from->cieuvw.v0 - to->v0); to->U = U, to->V = V, to->W = W; return; default: tmp.model = LIBCOLOUR_CIE1960UCS; to_cie1960ucs(from, &tmp.cie1960ucs); /* fall through */ case LIBCOLOUR_CIE1960UCS: U = tmp.cie1960ucs.u, V = tmp.cie1960ucs.v, W = tmp.cie1960ucs.Y; W = 25 * cbrt(W) - 17; w = W * 13; to->U = w * (U - to->u0); to->V = w * (V - to->v0); to->W = W; return; } } int libcolour_convert(const libcolour_colour_t* restrict from, libcolour_colour_t* restrict to) { #define X(C, T) 1 + if (from->model < 0 || from->model > LIBCOLOUR_LIST_MODELS 0) { #undef X errno = EINVAL; return -1; } switch (to->model) { case LIBCOLOUR_RGB: to_rgb(from, &to->rgb); break; case LIBCOLOUR_SRGB: to_srgb(from, &to->srgb); break; case LIBCOLOUR_CIEXYY: to_ciexyy(from, &to->ciexyy); break; case LIBCOLOUR_CIEXYZ: to_ciexyz(from, &to->ciexyz); break; case LIBCOLOUR_CIELAB: to_cielab(from, &to->cielab); break; case LIBCOLOUR_CIELUV: to_cieluv(from, &to->cieluv); break; case LIBCOLOUR_CIELCHUV: to_cielchuv(from, &to->cielchuv); break; case LIBCOLOUR_YIQ: to_yiq(from, &to->yiq); break; case LIBCOLOUR_YDBDR: to_ydbdr(from, &to->ydbdr); break; case LIBCOLOUR_YUV: to_yuv(from, &to->yuv); break; case LIBCOLOUR_YPBPR: to_ypbpr(from, &to->ypbpr); break; case LIBCOLOUR_YCGCO: to_ycgco(from, &to->ycgco); break; case LIBCOLOUR_CIE1960UCS: to_cie1960ucs(from, &to->cie1960ucs); break; case LIBCOLOUR_CIEUVW: to_cieuvw(from, &to->cieuvw); break; default: errno = EINVAL; return -1; } return 0; } double libcolour_srgb_encode(double x) { return x <= 0.0031306684425217108 ? 12.92 * x : 1.055 * pow(x, 1 / 2.4) - 0.055; } double libcolour_srgb_decode(double x) { return x <= 0.040448236277380506 ? x / 12.92 : pow((x + 0.055) / 1.055, 2.4); } int libcolour_delta_e(const libcolour_colour_t* a, const libcolour_colour_t* b, double* 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 = sqrt(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, double (*Mp)[n][m]) { #define M (*Mp) size_t r1, r2, c; double d; double* R1; double* 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; double 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; 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; double Sr, Sg, Sb; double 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(double)), M[0][3] = Sr; memcpy(M[1], M[1] + 3, 3 * sizeof(double)), M[1][3] = Sg; memcpy(M[2], M[2] + 3, 3 * sizeof(double)), 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; double M[3][6]; double 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(double)), M[0][3] = 1, M[0][4] = 0, M[0][5] = 0; memcpy(M[1], cs->M[1], 3 * sizeof(double)), M[1][3] = 0, M[1][4] = 1, M[1][5] = 0; memcpy(M[2], cs->M[2], 3 * sizeof(double)), 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(double)); memcpy(cs->Minv[1], M[1] + 3, 3 * sizeof(double)); memcpy(cs->Minv[2], M[2] + 3, 3 * sizeof(double)); return 0; } #define invert(Minv, M, n) invert_(n, &(Minv), &(M)) static int invert_(size_t n, double (*Minvp)[n][n], double (*Mp)[n][n]) { #define Minv (*Minvp) #define M (*Mp) #define J (*Jp) double J[n][2 * n] = alloca(sizeof(double[n][2 * n])); size_t i; for (i = 0; i < n; i++) { memcpy(J[i], M[i], n * sizeof(double)); memset(J[i] + n, 0, n * sizeof(double)); 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(double)); return 0; #undef J #undef M #undef Minv } static double transfer_function_l_star(double t) { t = (t > 216. / 24389.) ? cbrt(t) : t * 841. / 108. + 4. / 29.; return t * 1.16 - 0.16; } static double invtransfer_function_l_star(double t) { t = (t + 0.16) / 1.16; return (t > 6. / 29.) ? t * t * t : (t - 4. / 29.) * 108 / 841; } 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: case LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_HLG: /* TODO http://www.itu.int/dms_pubrec/itu-r/rec/bt/R-REC-BT.2100-0-201607-I!!PDF-E.pdf */ 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 = XVALUE, .y = YVALUE, .Y = 1} switch (space) { case LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MEASUREMENTS: if (get_matrices(cs)) return -1; return 0; case LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MATRIX: if (invert(cs->Minv, cs->M, 3) || get_primaries(cs)) return -1; return 0; case LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_INV_MATRIX: if (invert(cs->M, cs->Minv, 3) || get_primaries(cs)) return -1; 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 = 0.055; cs->gamma = 2.4; cs->slope = 12.92; cs->transition = 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 = 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 = 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 = 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 = 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 = 2.2; break; case LIBCOLOUR_RGB_COLOUR_SPACE_CIE_RGB: cs->red = XYY(0.7350, 0.2650); cs->green = XYY(0.2640, 0.7170); cs->blue = XYY(0.1670, 0.0090); cs->white = LIBCOLOUR_ILLUMINANT_E; cs->encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE; cs->gamma = 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 = 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 = 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 = 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 = 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 = 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 = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 = 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 = 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 / 0.45; cs->offset = 0.1115721957735072; cs->slope = 4.0; cs->transition = 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 / 0.45; cs->offset = 0.09929682680944; cs->slope = 4.5; cs->transition = 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 = 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)\ case C:\ if (buf)\ memcpy((char*)buf + sizeof(int), colour, sizeof(T));\ return sizeof(int) + sizeof(T); LIBCOLOUR_LIST_MODELS #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 = *(int*)buf; if (ver != MARSHAL_VERSION) { errno = EINVAL; return 0; } model = *(enum libcolour_model*)((char*)buf + sizeof(int)); switch (model) { #define X(C, T)\ case C:\ if (colour)\ memcpy(colour, (char*)buf + sizeof(int), sizeof(T));\ r = sizeof(int) + sizeof(T);\ break; LIBCOLOUR_LIST_MODELS #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; }