aboutsummaryrefslogtreecommitdiffstats
path: root/src
diff options
context:
space:
mode:
Diffstat (limited to '')
-rw-r--r--src/convert.c1068
-rw-r--r--src/libcolour.c717
-rw-r--r--src/libcolour.h301
-rw-r--r--src/test.c657
4 files changed, 0 insertions, 2743 deletions
diff --git a/src/convert.c b/src/convert.c
deleted file mode 100644
index 6d9be16..0000000
--- a/src/convert.c
+++ /dev/null
@@ -1,1068 +0,0 @@
-/**
- * Copyright © 2016 Mattias Andrée <maandree@kth.se>
- *
- * 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 <http://www.gnu.org/licenses/>.
- */
-#include "libcolour.h"
-
-#include <errno.h>
-#include <math.h>
-#include <stdio.h>
-#include <stdlib.h>
-#include <string.h>
-
-
-#define WASDIV0(x) (isinf(x) || isnan(x))
-#define PI (3.14159265358979323846)
-#define PI2 (2 * PI)
-
-
-
-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;
-}
-
-#define REGULAR(s, t) ((t) <= (s)->transition ? (s)->slope * (t) : (1 + (s)->offset) * pow((t), 1 / (s)->gamma) - (s)->offset)
-#define INVREGULAR(s, t) ((t) <= (s)->transitioninv ? (t) / (s)->slope : pow(((t) + (s)->offset) / (1 + (s)->offset), (s)->gamma))
-
-static void rgb_encode(libcolour_rgb_t* restrict colour, const libcolour_rgb_t* restrict space)
-{
- double r_sign = 1, g_sign = 1, b_sign = 1;
- switch (space->encoding_type) {
- case LIBCOLOUR_ENCODING_TYPE_LINEAR:
- break;
- case LIBCOLOUR_ENCODING_TYPE_SIMPLE:
- case LIBCOLOUR_ENCODING_TYPE_REGULAR:
- if (colour->R < 0) r_sign = -1, colour->R = -colour->R;
- if (colour->G < 0) g_sign = -1, colour->G = -colour->G;
- if (colour->B < 0) b_sign = -1, colour->B = -colour->B;
- if (space->encoding_type == LIBCOLOUR_ENCODING_TYPE_SIMPLE) {
- colour->R = pow(colour->R, 1 / space->gamma);
- colour->G = pow(colour->G, 1 / space->gamma);
- colour->B = pow(colour->B, 1 / space->gamma);
- } else {
- colour->R = REGULAR(space, colour->R);
- colour->G = REGULAR(space, colour->G);
- colour->B = REGULAR(space, colour->B);
- }
- colour->R *= r_sign;
- colour->G *= g_sign;
- colour->B *= b_sign;
- break;
- case LIBCOLOUR_ENCODING_TYPE_CUSTOM:
- colour->R = (space->to_encoded_red)(colour->R);
- colour->G = (space->to_encoded_green)(colour->G);
- colour->B = (space->to_encoded_blue)(colour->B);
- break;
- default:
- fprintf(stderr, "libcolour: invalid encoding type\n");
- abort();
- }
-}
-
-static void rgb_decode(libcolour_rgb_t* restrict colour, const libcolour_rgb_t* restrict space)
-{
- double r_sign = 1, g_sign = 1, b_sign = 1;
- switch (space->encoding_type) {
- case LIBCOLOUR_ENCODING_TYPE_LINEAR:
- break;
- case LIBCOLOUR_ENCODING_TYPE_SIMPLE:
- case LIBCOLOUR_ENCODING_TYPE_REGULAR:
- if (colour->R < 0) r_sign = -1, colour->R = -colour->R;
- if (colour->G < 0) g_sign = -1, colour->G = -colour->G;
- if (colour->B < 0) b_sign = -1, colour->B = -colour->B;
- if (space->encoding_type == LIBCOLOUR_ENCODING_TYPE_SIMPLE) {
- colour->R = pow(colour->R, space->gamma);
- colour->G = pow(colour->G, space->gamma);
- colour->B = pow(colour->B, space->gamma);
- } else {
- colour->R = INVREGULAR(space, colour->R);
- colour->G = INVREGULAR(space, colour->G);
- colour->B = INVREGULAR(space, colour->B);
- }
- colour->R *= r_sign;
- colour->G *= g_sign;
- colour->B *= b_sign;
- break;
- case LIBCOLOUR_ENCODING_TYPE_CUSTOM:
- colour->R = (space->to_decoded_red)(colour->R);
- colour->G = (space->to_decoded_green)(colour->G);
- colour->B = (space->to_decoded_blue)(colour->B);
- break;
- default:
- fprintf(stderr, "libcolour: invalid encoding type\n");
- abort();
- }
-}
-
-#undef REGULAR
-#undef INVREGULAR
-
-static int rgb_same_gamma(const libcolour_rgb_t* restrict a, const libcolour_rgb_t* restrict b)
-{
- if (a->encoding_type != b->encoding_type)
- return 0;
- switch (a->encoding_type) {
- case LIBCOLOUR_ENCODING_TYPE_SIMPLE:
- return a->gamma == b->gamma;
- case LIBCOLOUR_ENCODING_TYPE_REGULAR:
- return a->gamma == b->gamma &&
- a->offset == b->offset &&
- a->slope == b->slope &&
- a->transition == b->transition;
- case LIBCOLOUR_ENCODING_TYPE_CUSTOM:
- return a->to_encoded_red == b->to_encoded_red &&
- a->to_encoded_green == b->to_encoded_green &&
- a->to_encoded_blue == b->to_encoded_blue &&
- a->to_decoded_red == b->to_decoded_red &&
- a->to_decoded_green == b->to_decoded_green &&
- a->to_decoded_blue == b->to_decoded_blue;
- default:
- return 1;
- }
-}
-
-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->R = from->rgb.R;
- to->G = from->rgb.G;
- to->B = from->rgb.B;
- if (have_gamma && with_gamma && !rgb_same_gamma(&from->rgb, to))
- rgb_decode(to, &from->rgb), have_gamma = 0;
- 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) {
- rgb_encode(to, to);
- } else {
- rgb_decode(to, &from->rgb);
- }
- }
-}
-
-
-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;
-#define MULTIPLY(CX, CY, CZ) ((CX) * X + (CY) * Y + (CZ) * Z)
- to->R = MULTIPLY(3.240446254647737500675930277794, -1.537134761820080575134284117667, -0.498530193022728718155178739835);
- to->G = MULTIPLY(-0.969266606244679751469561779231, 1.876011959788370209167851498933, 0.041556042214430065351304932619);
- to->B = MULTIPLY(0.055643503564352832235773149705, -0.204026179735960239147729566866, 1.057226567722703292062647051353);
-#undef MULTIPLY
-}
-
-static inline double srgb_encode(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t <= 0.0031306684425217108 ? 12.92 * t : 1.055 * pow(t, 1 / 2.4) - 0.055;
- return t * sign;
-}
-
-static inline double srgb_decode(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t <= 0.0031306684425217108 * 12.92 ? t / 12.92 : pow((t + 0.055) / 1.055, 2.4);
- return t * sign;
-}
-
-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 = srgb_encode(from->R);
- to->G = srgb_encode(from->G);
- to->B = srgb_encode(from->B);
- } else {
- to->R = srgb_decode(from->R);
- to->G = srgb_decode(from->G);
- to->B = 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.95629483232089407263032398986979387700557708740234 + Q * 0.62102512544472865396727456754888407886028289794922;
- to->G = Y - I * 0.27212147408397735492968649850809015333652496337891 - Q * 0.64738095351761570928061928498209454119205474853516;
- to->B = Y - I * 1.10698990856712820018969978264067322015762329101562 + Q * 1.70461497549882934343656870623817667365074157714844;
- /* ⎛1 1000 cos 33° / 877 1000 sin 33° / 877 ⎞
- ⎜1 9500 sin 33° / 24067 - 299000 * cos 33° / 514799 -9500 cos 33° / 24067 - 299000 * sin 33° / 514799⎟
- ⎝1 -250 sin 33° / 123 250 cos 33° / 123 ⎠
- */
-}
-
-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 yuv_to_srgb(const libcolour_yuv_t* restrict from, libcolour_srgb_t* restrict to)
-{
- double Y = from->Y, U = from->U, V = from->V;
- to->R = Y + U * 0.00028328010485821202317155420580263580632163211703 + V * 1.14070449590558520291949662350816652178764343261719;
- to->G = Y - U * 0.39630886669497211727275498560629785060882568359375 - V * 0.58107364288228224857846271333983168005943298339844;
- to->B = Y + U * 2.03990003507541306504435851820744574069976806640625 + V * 0.00017179031692307700847528739718228507626918144524;
-}
-
-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_YUV:
- yuv_to_srgb(&from->yuv, 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;
- }
-#define MULTIPLIY(CR, CG, CB) ((CR) * R + (CG) * G + (CB) * B)
- to->X = MULTIPLIY(0.412457445582367576708548995157, 0.357575865245515878143578447634, 0.180437247826399665973085006954);
- to->Y = MULTIPLIY(0.212673370378408277403536885686, 0.715151730491031756287156895269, 0.072174899130559869164791564344);
- to->Z = MULTIPLIY(0.019333942761673460208893260415, 0.119191955081838593666354597644, 0.950302838552371742508739771438);
-#undef MULTIPLIY
-}
-
-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);
- to->Z = cielab_finv(to->Z) * 1.08883;
-}
-
-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 yiq_to_ciexyz(const libcolour_yiq_t* restrict from, libcolour_ciexyz_t* restrict to)
-{
-#define MULTIPLIY(CY, CI, CQ) ((CY) * from->Y + (CI) * from->I + (CQ) * from->Q)
- to->X = MULTIPLIY(0.95047055865428309306963683411595411598682403564453, 0.09738463974178063198294807989441324025392532348633, 0.33223466706854809515903070860076695680618286132812);
- to->Y = MULTIPLIY(0.99999999999999988897769753748434595763683319091797, -0.07112658301916767455974621725545148365199565887451, -0.20786968876398304040264974901219829916954040527344);
- to->Z = MULTIPLIY(1.08882873639588373393394249433185905218124389648438, -1.06592139332461721679123911599162966012954711914062, 1.55474471255181900808395312196807935833930969238281);
-#undef MULTIPLIY
-}
-
-static void ydbdr_to_ciexyz(const libcolour_ydbdr_t* restrict from, libcolour_ciexyz_t* restrict to)
-{
-#define MULTIPLIY(CY, CDB, CDR) ((CY) * from->Y + (CDB) * from->Db + (CDR) * from->Dr)
- to->X = MULTIPLIY(0.95047055865428309306963683411595411598682403564453, 0.07379612357298576119646327242662664502859115600586, -0.12113653724874726136384595065464964136481285095215);
- to->Y = MULTIPLIY(0.99999999999999988897769753748434595763683319091797, -0.04435684145428285540813106990754022262990474700928, 0.07973534004202506575431641522300196811556816101074);
- to->Z = MULTIPLIY(1.08882873639588373393394249433185905218124389648438, 0.61625657933494271123464613992837257683277130126953, 0.02168821359337728266192257819966471288353204727173);
-#undef MULTIPLIY
-}
-
-static void yuv_to_ciexyz(const libcolour_yuv_t* restrict from, libcolour_ciexyz_t* restrict to)
-{
-#define MULTIPLIY(CY, CU, CV) ((CY) * from->Y + (CU) * from->U + (CV) * from->V)
- to->X = MULTIPLIY(0.95047055865428309306963683411595411598682403564453, 0.22648030324549334180161963558930438011884689331055, 0.26274514929253273143316960158699657768011093139648);
- to->Y = MULTIPLIY(0.99999999999999988897769753748434595763683319091797, -0.13613114642319409930415474718756740912795066833496, -0.17294595255115238763288232348713790997862815856934);
- to->Z = MULTIPLIY(1.08882873639588373393394249433185905218124389648438, 1.89129144197893928058817891724174842238426208496094, -0.04704173528403532422714761196402832865715026855469);
-#undef MULTIPLIY
-}
-
-static void ypbpr_to_ciexyz(const libcolour_ypbpr_t* restrict from, libcolour_ciexyz_t* restrict to)
-{
-#define MULTIPLIY(CY, CPB, CPR) ((CY) * from->Y + (CPB) * from->Pb + (CPR) * from->Pr)
- to->X = MULTIPLIY(0.95047055865428298204733437160030007362365722656250, 0.14433968424876231217979238863335922360420227050781, 0.30616461986760712399302519770571961998939514160156);
- to->Y = MULTIPLIY(0.99999999999999977795539507496869191527366638183594, -0.00002022802471486617736928792510298080742359161377, 0.00008771894888734421691367515450110659003257751465);
- to->Z = MULTIPLIY(1.08882873639588373393394249433185905218124389648438, 0.93827031735982591165168287261622026562690734863281, -0.01609699914324668607035206946420657914131879806519);
-#undef MULTIPLIY
-}
-
-static void ycgco_to_ciexyz(const libcolour_ycgco_t* restrict from, libcolour_ciexyz_t* restrict to)
-{
-#define MULTIPLIY(CY, CCG, CCO) ((CY) * from->Y + (CCG) * from->Cg + (CCO) * from->Co)
- to->X = MULTIPLIY(0.95047055865428309306963683411595411598682403564453, -0.23531882816325136453805555447615915909409523010254, 0.23202019775596791073546398820326430723071098327637);
- to->Y = MULTIPLIY(0.99999999999999988897769753748434595763683319091797, 0.43030346098206362359661625305307097733020782470703, 0.14049847124784842211653312915586866438388824462891);
- to->Z = MULTIPLIY(1.08882873639588373393394249433185905218124389648438, -0.85044482623220662986796014592982828617095947265625, -0.93096889579069830311652822274481877684593200683594);
-#undef MULTIPLIY
-}
-
-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);
- break;
- case LIBCOLOUR_SRGB:
- srgb_to_ciexyz(&from->srgb, to);
- break;
- case LIBCOLOUR_CIEXYY:
- ciexyy_to_ciexyz(&from->ciexyy, to);
- break;
- case LIBCOLOUR_CIEXYZ:
- *to = from->ciexyz;
- break;
- case LIBCOLOUR_CIELAB:
- cielab_to_ciexyz(&from->cielab, to);
- break;
- case LIBCOLOUR_CIELUV:
- cieluv_to_ciexyz(&from->cieluv, to);
- break;
- case LIBCOLOUR_CIELCHUV:
- cielchuv_to_ciexyz(&from->cielchuv, to);
- break;
- case LIBCOLOUR_YIQ:
- yiq_to_ciexyz(&from->yiq, to);
- break;
- case LIBCOLOUR_YDBDR:
- ydbdr_to_ciexyz(&from->ydbdr, to);
- break;
- case LIBCOLOUR_YUV:
- yuv_to_ciexyz(&from->yuv, to);
- break;
- case LIBCOLOUR_YPBPR:
- ypbpr_to_ciexyz(&from->ypbpr, to);
- break;
- case LIBCOLOUR_YCGCO:
- ycgco_to_ciexyz(&from->ycgco, to);
- break;
- case LIBCOLOUR_CIE1960UCS:
- cie1960ucs_to_ciexyz(&from->cie1960ucs, to);
- break;
- case LIBCOLOUR_CIEUVW:
- cieuvw_to_ciexyz(&from->cieuvw, to);
- break;
- default:
- other_to_ciexyz(from, to);
- break;
- }
-}
-
-
-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;
- Z /= 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;
- tmp2.one_revolution = PI2;
- 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 * PI2 / from->one_revolution;
- }
- 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);
- break;
- case LIBCOLOUR_CIELCHUV:
- cielchuv_to_cieluv(&from->cielchuv, to);
- break;
- 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;
- break;
- }
- /* fall through */
- default:
- other_to_cieluv(from, to);
- break;
- }
-}
-
-
-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) / PI2 * to->one_revolution;
- if (!isnan(to->h) && !isinf(to->h) && (to->h < 0))
- to->h += to->one_revolution;
-}
-
-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)
-{
- double one_revolution;
- switch (from->model) {
- case LIBCOLOUR_CIELUV:
- cieluv_to_cielchuv(&from->cieluv, to);
- break;
- 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) {
- if (to->one_revolution == from->cielchuv.one_revolution) {
- *to = from->cielchuv;
- } else {
- one_revolution = to->one_revolution;
- *to = from->cielchuv;
- to->one_revolution = one_revolution;
- to->h = to->h / from->cielchuv.one_revolution * one_revolution;
- }
- break;
- }
- /* fall through */
- default:
- other_to_cielchuv(from, to);
- break;
- }
-}
-
-
-static void to_yiq(const libcolour_colour_t* restrict from, libcolour_yiq_t* restrict to)
-{
- double x, y, z;
- libcolour_colour_t tmp = *from;
- switch (from->model) {
- case LIBCOLOUR_YIQ:
- *to = from->yiq;
- break;
- case LIBCOLOUR_SRGB:
- if (tmp.srgb.with_gamma) {
- tmp.srgb.with_gamma = 0;
- to_srgb(from, &tmp.srgb);
- }
- x = tmp.srgb.R;
- y = tmp.srgb.G;
- z = tmp.srgb.B;
- to->Y = x * 0.299 + y * 0.587 + z * 0.114;
- to->I = x * 0.59571613491277464191853141528554260730743408203125 /* (0.877 cos 33°)(1 - 0.299) - (0.492 sin 33°)(-0.299) */
- - y * 0.27445283783925644716106262421817518770694732666016 /* (0.877 cos 33°)(-0.587) - (0.492 sin 33°)(-0.587) */
- - z * 0.32126329707351808373516632855171337723731994628906; /* (0.877 cos 33°)(-0.114) - (0.492 sin 33°)(1 - 0.114) */
- to->Q = x * 0.21145640212011795888713550084503367543220520019531 /* (0.877 sin 33°)(1 - 0.299) + (0.492 cos 33°)(-0.299) */
- - y * 0.52259104529161115593183239980135113000869750976562 /* (0.877 sin 33°)(-0.587) + (0.492 cos 33°)(-0.587) */
- + z * 0.31113464317149330806699936147197149693965911865234; /* (0.877 sin 33°)(-0.114) + (0.492 cos 33°)(1 - 0.114) */
- break;
- default:
- tmp.model = LIBCOLOUR_CIEXYZ;
- to_ciexyz(from, &tmp.ciexyz);
- /* fall through */
- case LIBCOLOUR_CIEXYZ:
- x = tmp.ciexyz.X;
- y = tmp.ciexyz.Y;
- z = tmp.ciexyz.Z;
-#define MULTIPLIY(CX, CY, CZ) ((CX) * x + (CY) * y + (CZ) * z)
- to->Y = MULTIPLIY(0.40627729168038273499519164033699780702590942382812, 0.61835674212166968910509012857801280915737152099609, -0.00414330221353725880462093300593551248311996459961);
- to->I = MULTIPLIY(2.17852787350219845308174626552499830722808837890625, -1.36502666214454104753883711964590474963188171386719, -0.64803574634025240541745915834326297044754028320312);
- to->Q = MULTIPLIY(1.20905577682138853923277110880007967352867126464844, -1.36890364998339797431015085749095305800437927246094, 0.20180559439597040016778350945969577878713607788086);
-#undef MULTIPLIY
- break;
- }
-}
-
-
-static void to_ydbdr(const libcolour_colour_t* restrict from, libcolour_ydbdr_t* restrict to)
-{
- double x, y, z;
- libcolour_colour_t tmp = *from;
- switch (from->model) {
- case LIBCOLOUR_YDBDR:
- *to = from->ydbdr;
- break;
- case LIBCOLOUR_SRGB:
- if (tmp.srgb.with_gamma) {
- tmp.srgb.with_gamma = 0;
- to_srgb(from, &tmp.srgb);
- }
- x = tmp.srgb.R;
- y = tmp.srgb.G;
- z = tmp.srgb.B;
- to->Y = x * 0.299 + y * 0.587 + z * 0.114;
- to->Db = -x * 0.450 - y * 0.883 + z * 1.333;
- to->Dr = -x * 1.333 + y * 1.116 + z * 0.217;
- break;
- case LIBCOLOUR_YUV:
- to->Y = tmp.yuv.Y;
- to->Db = tmp.yuv.U * 3.069;
- to->Dr = -tmp.yuv.V * 2.169;
- break;
- default:
- tmp.model = LIBCOLOUR_CIEXYZ;
- to_ciexyz(from, &tmp.ciexyz);
- /* fall through */
- case LIBCOLOUR_CIEXYZ:
- x = tmp.ciexyz.X;
- y = tmp.ciexyz.Y;
- z = tmp.ciexyz.Z;
-#define MULTIPLY(CX, CY, CZ) ((CX) * x + (CY) * y + (CZ) * z)
- to->Y = MULTIPLY(0.40627729168038273499519164033699780702590942382812, 0.61835674212166968910509012857801280915737152099609, -0.00414330221353725880462093300593551248311996459961);
- to->Db = MULTIPLY(-0.52816561102614745237815441214479506015777587890625, -1.23677481526212962315014465275453403592109680175781, 1.59692761635924940222253098909277468919754028320312);
- to->Dr = MULTIPLY(-5.38914174974103143966885909321717917919158935546875, 4.09835630362728497999569299281574785709381103515625, 0.94033545560642795013706063400604762136936187744141);
-#undef MULTIPLY
- break;
- }
-}
-
-
-static void to_yuv(const libcolour_colour_t* restrict from, libcolour_yuv_t* restrict to)
-{
- libcolour_colour_t tmp = *from;
- switch (from->model) {
- case LIBCOLOUR_YUV:
- *to = from->yuv;
- break;
- case LIBCOLOUR_YDBDR:
- to->Y = from->ydbdr.Y;
- to->U = from->ydbdr.Db / 3.069;
- to->V = from->ydbdr.Dr / -2.169;
- break;
- case LIBCOLOUR_SRGB:
- if (tmp.srgb.with_gamma) {
- tmp.srgb.with_gamma = 0;
- to_srgb(from, &tmp.srgb);
- }
-#define MULTIPLY(CR, CG, CB) ((CR) * tmp.srgb.R + (CG) * tmp.srgb.G + (CB) * tmp.srgb.B)
- to->Y = MULTIPLY(0.29899999999999998800959133404830936342477798461914, 0.58699999999999996624922005139524117112159729003906, 0.11400000000000000410782519111307919956743717193604);
- to->U = MULTIPLY(-0.14662756598240470062854967636667424812912940979004, -0.28771586836102963635752871596196200698614120483398, 0.43434343434343436474165400795754976570606231689453);
- to->V = MULTIPLY(0.61456892577224520035628074765554629266262054443359, -0.51452282157676354490405401520547457039356231689453, -0.10004610419548178035231700278018251992762088775635);
-#undef MULTIPLY
- break;
- default:
- tmp.model = LIBCOLOUR_CIEXYZ;
- to_ciexyz(from, &tmp.ciexyz);
- /* fall through */
- case LIBCOLOUR_CIEXYZ:
-#define MULTIPLY(CX, CY, CZ) ((CX) * tmp.ciexyz.X + (CY) * tmp.ciexyz.Y + (CZ) * tmp.ciexyz.Z)
- to->Y = MULTIPLY(0.40627729168038273499519164033699780702590942382812, 0.61835674212166968910509012857801280915737152099609, -0.00414330221353725880462093300593551248311996459961);
- to->U = MULTIPLY(-0.17209697328971895746718701047939248383045196533203, -0.40298951295605395239718404809536878019571304321289, 0.52034135430408912093014350830344483256340026855469);
- to->V = MULTIPLY(2.48462044709130047692724474472925066947937011718750, -1.88951420176453876997868519538315013051033020019531, -0.43353409663735725798616726933687459677457809448242);
-#undef MULTIPLY
- break;
- }
-}
-
-
-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;
- break;
- 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;
- break;
- default:
- tmp.model = LIBCOLOUR_CIEXYZ;
- to_ciexyz(from, &tmp.ciexyz);
- /* fall through */
- case LIBCOLOUR_CIEXYZ:
-#define MULTIPLY(CX, CY, CZ) ((CX) * tmp.ciexyz.X + (CY) * tmp.ciexyz.Y + (CZ) * tmp.ciexyz.Z)
- to->Y = MULTIPLY(-0.00028314209073960778378920011277841695118695497513, 1.00019821310075673892470149439759552478790283203125, 0.00006512054470741990286342115723527967929840087891);
- to->Pb = MULTIPLY(0.05592664565509243568275365987574332393705844879150, -1.20422439283671711685030913940863683819770812988281, 1.05716144717799576113748116767965257167816162109375);
- to->Pr = MULTIPLY(3.24072939673847715269516811531502753496170043945312, -2.53733297492083753610359053709544241428375244140625, -0.49859531356743613805804216099204495549201965332031);
-#undef MULTIPLY
- break;
- }
-}
-
-
-static void to_ycgco(const libcolour_colour_t* restrict from, libcolour_ycgco_t* restrict to)
-{
- libcolour_colour_t tmp = *from;
- switch (from->model) {
- case LIBCOLOUR_YCGCO:
- *to = from->ycgco;
- break;
- 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;
- break;
- default:
- tmp.model = LIBCOLOUR_CIEXYZ;
- to_ciexyz(from, &tmp.ciexyz);
- /* fall through */
- case LIBCOLOUR_CIEXYZ:
-#define MULTIPLY(CX, CY, CZ) ((CX) * tmp.ciexyz.X + (CY) * tmp.ciexyz.Y + (CZ) * tmp.ciexyz.Z)
- to->Y = MULTIPLY(0.33938913643068269188063368346774950623512268066406, 0.50271574450517486631895280879689380526542663574219, 0.16045211478220866574417868832824751734733581542969);
- to->Cg = MULTIPLY(-1.30865574267536244335019546269904822111129760742188, 1.37329621528319534284889869013568386435508728027344, -0.11889607256777862120955546743061859160661697387695);
- to->Co = MULTIPLY(1.59240137554169236544510113162687048316001892089844, -0.66655429104206020962664069884340278804302215576172, -0.77787838037271606062006412685150280594825744628906);
-#undef MULTIPLY
- break;
- }
-}
-
-
-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;
- break;
- 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;
- break;
- }
-}
-
-
-
-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;
-}
-
diff --git a/src/libcolour.c b/src/libcolour.c
deleted file mode 100644
index 867e280..0000000
--- a/src/libcolour.c
+++ /dev/null
@@ -1,717 +0,0 @@
-/**
- * Copyright © 2016 Mattias Andrée <maandree@kth.se>
- *
- * 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 <http://www.gnu.org/licenses/>.
- */
-#include "libcolour.h"
-
-#include <alloca.h>
-#include <errno.h>
-#include <math.h>
-#include <string.h>
-
-
-#define MARSHAL_VERSION 0
-
-
-
-double libcolour_srgb_encode(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t <= 0.0031306684425217108 ? 12.92 * t : 1.055 * pow(t, 1 / 2.4) - 0.055;
- return t * sign;
-}
-
-double libcolour_srgb_decode(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t <= 0.0031306684425217108 * 12.92 ? t / 12.92 : pow((t + 0.055) / 1.055, 2.4);
- return t * sign;
-}
-
-
-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;
- if (colour->cielchuv.one_revolution == 0.)
- colour->cielchuv.one_revolution = 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;
- 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)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t > 216. / 24389. ? 1.16 * cbrt(t) - 0.16 : t * 24389. / 2700.;
- return t * sign;
-}
-
-static double invtransfer_function_l_star(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t > 0.08 ? (((1000000. * t + 480000.) * t + 76800.) * t + 4096.) / 1560896.0 : t * 2700. / 24389.;
- return t * sign;
-}
-
-static double transfer_function_oetf_hlg(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = 12 * t <= 1 ? sqrt(3 * t) : 0.17883277 * log(t - 0.02372241) + 1.004293468902569985701234145381;
- return t * sign;
-}
-
-static double invtransfer_function_oetf_hlg(double t)
-{
- double sign = 1;
- if (t < 0) {
- t = -t;
- sign = -1;
- }
- t = t <= 0.5 ? t * t / 3 : exp(t - 1.004293468902569985701234145381) / 0.17883277 + 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 = XVALUE, .y = 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 = 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.2740, 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;
-}
-
diff --git a/src/libcolour.h b/src/libcolour.h
deleted file mode 100644
index 6387151..0000000
--- a/src/libcolour.h
+++ /dev/null
@@ -1,301 +0,0 @@
-/**
- * Copyright © 2016 Mattias Andrée <maandree@kth.se>
- *
- * 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 <http://www.gnu.org/licenses/>.
- */
-#ifndef LIBCOLOUR_H
-#define LIBCOLOUR_H
-
-
-#include <stddef.h>
-
-
-
-#define LIBCOLOUR_ILLUMINANT__(x__, y__) (libcolour_ciexyy_t){.model = LIBCOLOUR_CIEXYY, .x = x__, .y = y__, .Y = 1}
-
-#define LIBCOLOUR_ILLUMINANT_A LIBCOLOUR_ILLUMINANT__(0.447573514098910552050369915378, 0.407439444306660847328060981454)
-#define LIBCOLOUR_ILLUMINANT_B LIBCOLOUR_ILLUMINANT__(0.348407693041403399014654951316, 0.351617234807268863594487129376)
-#define LIBCOLOUR_ILLUMINANT_C LIBCOLOUR_ILLUMINANT__(0.310058473730255412803558101587, 0.316149707523236456196968902077)
-#define LIBCOLOUR_ILLUMINANT_D50 LIBCOLOUR_ILLUMINANT__(0.345668037029273123028616510055, 0.358496838937619077825047497754)
-#define LIBCOLOUR_ILLUMINANT_D55 LIBCOLOUR_ILLUMINANT__(0.332424102468830251488896010414, 0.347428039087666229445261478759)
-#define LIBCOLOUR_ILLUMINANT_D65 LIBCOLOUR_ILLUMINANT__(0.312726871026564878786047074755, 0.329023206641284038376227272238)
-#define LIBCOLOUR_ILLUMINANT_D75 LIBCOLOUR_ILLUMINANT__(0.299022300412497055166483050925, 0.314852737888341893679466920730)
-#define LIBCOLOUR_ILLUMINANT_E LIBCOLOUR_ILLUMINANT__(1. / 3, 1. / 3)
-#define LIBCOLOUR_ILLUMINANT_F1 LIBCOLOUR_ILLUMINANT__(0.313062433035651010992950205036, 0.337106477918307445573731229160)
-#define LIBCOLOUR_ILLUMINANT_F2 LIBCOLOUR_ILLUMINANT__(0.372068154452825539113547392844, 0.375122558203110079144693145281)
-#define LIBCOLOUR_ILLUMINANT_F3 LIBCOLOUR_ILLUMINANT__(0.409090035308107391465171076561, 0.394117134255365986206243178458)
-#define LIBCOLOUR_ILLUMINANT_F4 LIBCOLOUR_ILLUMINANT__(0.440181095827666568620628595454, 0.403090691158138336724903183494)
-#define LIBCOLOUR_ILLUMINANT_F5 LIBCOLOUR_ILLUMINANT__(0.313756583095696484075887155996, 0.345160794752101929283583103825)
-#define LIBCOLOUR_ILLUMINANT_F6 LIBCOLOUR_ILLUMINANT__(0.377882361062687466279896852939, 0.388192885537868959122675960316)
-#define LIBCOLOUR_ILLUMINANT_F7 LIBCOLOUR_ILLUMINANT__(0.312852472915475354753311876266, 0.329174178033567632617462095368)
-#define LIBCOLOUR_ILLUMINANT_F8 LIBCOLOUR_ILLUMINANT__(0.345805753550315952971061506105, 0.358617583214377477762724311106)
-#define LIBCOLOUR_ILLUMINANT_F9 LIBCOLOUR_ILLUMINANT__(0.374105245592801061160770359493, 0.372672400924498159469067104510)
-#define LIBCOLOUR_ILLUMINANT_F10 LIBCOLOUR_ILLUMINANT__(0.346086913993929323751785886998, 0.358751605952200347537939251197)
-#define LIBCOLOUR_ILLUMINANT_F11 LIBCOLOUR_ILLUMINANT__(0.380537485483030235577928124258, 0.376915309293930078649026427229)
-#define LIBCOLOUR_ILLUMINANT_F12 LIBCOLOUR_ILLUMINANT__(0.437023901312296902954557253906, 0.404214327891585678553809657387)
-
-#define LIBCOLOUR_ILLUMINANT_A_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.451173939693730152722395132514, 0.405936604212625562482230634487)
-#define LIBCOLOUR_ILLUMINANT_B_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.349819801494100579564161535018, 0.352687989927865819250740742063)
-#define LIBCOLOUR_ILLUMINANT_C_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.310388663270034004248998371622, 0.319050711366790695766582075521)
-#define LIBCOLOUR_ILLUMINANT_D50_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.347729429961154856698613002663, 0.359522508516545380441442603114)
-#define LIBCOLOUR_ILLUMINANT_D55_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.334116336430253457745465084372, 0.348766090975953568786849245953)
-#define LIBCOLOUR_ILLUMINANT_D65_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.313823646938709621689866935412, 0.330998985489933561510156323493)
-#define LIBCOLOUR_ILLUMINANT_D75_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.299679971345752860223399238748, 0.317403239854836705102769656150)
-#define LIBCOLOUR_ILLUMINANT_E_10DEG_OBS LIBCOLOUR_ILLUMINANT__(1. / 3, 1. / 3)
-#define LIBCOLOUR_ILLUMINANT_F1_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.318098801070991199502202562144, 0.335489451474129951602520804954)
-#define LIBCOLOUR_ILLUMINANT_F2_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.379274832262508854174853922814, 0.367227934400669309145115448700)
-#define LIBCOLOUR_ILLUMINANT_F3_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.417644682102624287267644831445, 0.383124504918675723441623404142)
-#define LIBCOLOUR_ILLUMINANT_F4_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.449247699162001246087072559021, 0.390605475879083674506375700730)
-#define LIBCOLOUR_ILLUMINANT_F5_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.319739939104951298443069163113, 0.342367055369128092667807550242)
-#define LIBCOLOUR_ILLUMINANT_F6_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.386626908526034762658696308790, 0.378372201588893453116924092683)
-#define LIBCOLOUR_ILLUMINANT_F7_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.315645637312390425766039925293, 0.329508145132134222521358424274)
-#define LIBCOLOUR_ILLUMINANT_F8_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.348965563721531868424108324689, 0.359317299140994528272585739614)
-#define LIBCOLOUR_ILLUMINANT_F9_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.378258900384649654480284652891, 0.370371375730762564248976786985)
-#define LIBCOLOUR_ILLUMINANT_F10_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.350893389986753234666139178444, 0.354302210111646531665030579461)
-#define LIBCOLOUR_ILLUMINANT_F11_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.385435391037903751776383387551, 0.371094786781121399599214782938)
-#define LIBCOLOUR_ILLUMINANT_F12_10DEG_OBS LIBCOLOUR_ILLUMINANT__(0.442654456042513022584472537346, 0.397060737666593277506166259627)
-
-
-#define LIBCOLOUR_LIST_MODELS\
- X(LIBCOLOUR_RGB, libcolour_rgb_t)\
- X(LIBCOLOUR_SRGB, libcolour_srgb_t)\
- X(LIBCOLOUR_CIEXYY, libcolour_ciexyy_t)\
- X(LIBCOLOUR_CIEXYZ, libcolour_ciexyz_t)\
- X(LIBCOLOUR_CIELAB, libcolour_cielab_t)\
- X(LIBCOLOUR_CIELUV, libcolour_cieluv_t)\
- X(LIBCOLOUR_CIELCHUV, libcolour_cielchuv_t)\
- X(LIBCOLOUR_YIQ, libcolour_yiq_t)\
- X(LIBCOLOUR_YDBDR, libcolour_ydbdr_t)\
- X(LIBCOLOUR_YUV, libcolour_yuv_t)\
- X(LIBCOLOUR_YPBPR, libcolour_ypbpr_t)\
- X(LIBCOLOUR_YCGCO, libcolour_ycgco_t)\
- X(LIBCOLOUR_CIE1960UCS, libcolour_cie1960ucs_t)\
- X(LIBCOLOUR_CIEUVW, libcolour_cieuvw_t)
-
-
-typedef enum libcolour_model {
-#define X(C, T) C,
- LIBCOLOUR_LIST_MODELS
-#undef X
-} libcolour_model_t;
-
-
-typedef enum libcolour_encoding_type {
- LIBCOLOUR_ENCODING_TYPE_LINEAR,
- LIBCOLOUR_ENCODING_TYPE_SIMPLE,
- LIBCOLOUR_ENCODING_TYPE_REGULAR,
- LIBCOLOUR_ENCODING_TYPE_CUSTOM
-} libcolour_encoding_type;
-
-
-typedef enum libcolour_rgb_colour_space {
- LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MEASUREMENTS,
- LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MATRIX,
- LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_INV_MATRIX,
- LIBCOLOUR_RGB_COLOUR_SPACE_SRGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_ADOBE_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_APPLE_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_BEST_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_BETA_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_BRUCE_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_CIE_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_COLORMATCH_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_D65,
- LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_THEATER,
- LIBCOLOUR_RGB_COLOUR_SPACE_DON_RGB_4,
- LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2,
- LIBCOLOUR_RGB_COLOUR_SPACE_EKTA_SPACE_PS5,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_625_LINE,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_525_LINE,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_709,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2020,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_PQ,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_PQ,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_PQ,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_HLG,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_HLG,
- LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_HLG,
- LIBCOLOUR_RGB_COLOUR_SPACE_LIGHTROOM_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_NTSC_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_PAL_SECAM_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_PROPHOTO_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_SGI_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_240M_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_C_RGB,
- LIBCOLOUR_RGB_COLOUR_SPACE_WIDE_GAMUT_RGB
-} libcolour_rgb_colour_space_t;
-
-
-typedef struct libcolour_srgb {
- enum libcolour_model model;
- double R;
- double G;
- double B;
- int with_gamma;
-} libcolour_srgb_t;
-
-typedef struct libcolour_ciexyy {
- enum libcolour_model model;
- double x;
- double y;
- double Y;
-} libcolour_ciexyy_t;
-
-typedef struct libcolour_ciexyz {
- enum libcolour_model model;
- double X;
- double Y;
- double Z;
-} libcolour_ciexyz_t;
-
-typedef struct libcolour_cielab {
- enum libcolour_model model;
- double L;
- double a;
- double b;
-} libcolour_cielab_t;
-
-typedef struct libcolour_yiq {
- enum libcolour_model model;
- double Y;
- double I;
- double Q;
-} libcolour_yiq_t;
-
-typedef struct libcolour_ydbdr {
- enum libcolour_model model;
- double Y;
- double Db;
- double Dr;
-} libcolour_ydbdr_t;
-
-typedef struct libcolour_yuv {
- enum libcolour_model model;
- double Y;
- double U;
- double V;
-} libcolour_yuv_t;
-
-typedef struct libcolour_ypbpr {
- enum libcolour_model model;
- double Y;
- double Pb;
- double Pr;
-} libcolour_ypbpr_t;
-
-typedef struct libcolour_ycgco {
- enum libcolour_model model;
- double Y;
- double Cg;
- double Co;
-} libcolour_ycgco_t;
-
-typedef struct libcolour_cie1960ucs {
- enum libcolour_model model;
- double u;
- double v;
- double Y;
-} libcolour_cie1960ucs_t;
-
-typedef struct libcolour_cieuvw {
- enum libcolour_model model;
- double U;
- double V;
- double W;
- double u0;
- double v0;
-} libcolour_cieuvw_t;
-
-typedef struct libcolour_cieluv {
- enum libcolour_model model;
- double L;
- double u;
- double v;
- struct libcolour_ciexyz white;
-} libcolour_cieluv_t;
-
-typedef struct libcolour_cielchuv {
- enum libcolour_model model;
- double L;
- double C;
- double h;
- struct libcolour_ciexyz white;
- double one_revolution;
-} libcolour_cielchuv_t;
-
-typedef struct libcolour_rgb {
- enum libcolour_model model;
- double R;
- double G;
- double B;
- int with_gamma;
- enum libcolour_encoding_type encoding_type;
- double gamma;
- double offset;
- double slope;
- double transition;
- double transitioninv;
- double (*to_encoded_red)(double);
- double (*to_decoded_red)(double);
- double (*to_encoded_green)(double);
- double (*to_decoded_green)(double);
- double (*to_encoded_blue)(double);
- double (*to_decoded_blue)(double);
- struct libcolour_ciexyy red;
- struct libcolour_ciexyy green;
- struct libcolour_ciexyy blue;
- struct libcolour_ciexyy white;
- double white_r;
- double white_g;
- double white_b;
- double M[3][3];
- double Minv[3][3];
- enum libcolour_rgb_colour_space colour_space;
-} libcolour_rgb_t;
-
-
-typedef union libcolour_colour {
- enum libcolour_model model;
- struct libcolour_rgb rgb;
- struct libcolour_srgb srgb;
- struct libcolour_ciexyy ciexyy;
- struct libcolour_ciexyz ciexyz;
- struct libcolour_cielab cielab;
- struct libcolour_cieluv cieluv;
- struct libcolour_cielchuv cielchuv;
- struct libcolour_yiq yiq;
- struct libcolour_ydbdr ydbdr;
- struct libcolour_yuv yuv;
- struct libcolour_ypbpr ypbpr;
- struct libcolour_ycgco ycgco;
- struct libcolour_cie1960ucs cie1960ucs;
- struct libcolour_cieuvw cieuvw;
-} libcolour_colour_t;
-
-
-
-int libcolour_convert(const libcolour_colour_t* restrict from, libcolour_colour_t* restrict to);
-double libcolour_srgb_encode(double);
-double libcolour_srgb_decode(double);
-int libcolour_delta_e(const libcolour_colour_t*, const libcolour_colour_t*, double*);
-int libcolour_proper(libcolour_colour_t*);
-int libcolour_get_rgb_colour_space(libcolour_rgb_t*, libcolour_rgb_colour_space_t);
-size_t libcolour_marshal(const libcolour_colour_t*, void*);
-size_t libcolour_unmarshal(libcolour_colour_t*, const void*);
-
-
-
-#endif
-
diff --git a/src/test.c b/src/test.c
deleted file mode 100644
index ad7acfc..0000000
--- a/src/test.c
+++ /dev/null
@@ -1,657 +0,0 @@
-/**
- * Copyright © 2016 Mattias Andrée <maandree@kth.se>
- *
- * 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 <http://www.gnu.org/licenses/>.
- */
-#include "libcolour.h"
-
-#include <math.h>
-#include <stdio.h>
-#include <string.h>
-
-
-#ifndef SKIP_2CONVERT
-static int test_2convert(libcolour_colour_t* c1, libcolour_colour_t* c2, libcolour_colour_t* c3)
-{
- double ch1, ch2, ch3;
- if (libcolour_convert(c1, c2))
- return -1;
- if (libcolour_convert(c2, c3))
- return -1;
- if ((c1->model == LIBCOLOUR_CIELCHUV) &&
- !isnan(c1->cielchuv.h) && !isnan(c3->cielchuv.h) &&
- !isinf(c1->cielchuv.h) && !isinf(c3->cielchuv.h) &&
- (c3->cielchuv.h / c3->cielchuv.one_revolution > 0.99999))
- c3->cielchuv.h = 0;
- ch1 = c1->srgb.R - c3->srgb.R;
- ch2 = c1->srgb.G - c3->srgb.G;
- ch3 = c1->srgb.B - c3->srgb.B;
- ch1 *= ch1;
- ch2 *= ch2;
- ch3 *= ch3;
- switch (c1->model) {
- case LIBCOLOUR_CIEXYY:
- if (c1->ciexyy.Y == 0)
- ch1 = ch2 = 0;
- if (c1->ciexyy.y == 0)
- ch1 = ch2 = ch3 = 0;
- break;
- case LIBCOLOUR_CIEXYZ:
- if (c1->ciexyz.Y == 0)
- ch1 = ch3 = 0;
- break;
- case LIBCOLOUR_CIE1960UCS:
- if (c1->cie1960ucs.v == 0)
- ch1 = ch2 = ch3 = 0;
- if (c1->cie1960ucs.Y == 0)
- ch1 = ch2 = 0;
- break;
- case LIBCOLOUR_CIEUVW:
- if (c1->cieuvw.W == 0)
- ch1 = ch2 = 0;
- break;
- case LIBCOLOUR_CIELCHUV:
- if (c1->cielchuv.C == 0)
- ch1 = ch3 = 0;
- break;
- case LIBCOLOUR_CIELAB:
- if (c1->cielchuv.L == 0)
- ch2 = ch3 = 0;
- break;
- default:
- break;
- }
- if (ch1 > 0.000001 ||
- ch2 > 0.000001 ||
- ch3 > 0.000001)
- return 0;
- return 1;
-}
-
-
-static int test_2convert_11(libcolour_colour_t* c1, libcolour_model_t model)
-{
- libcolour_colour_t c2, c3;
- int r1, r2;
- c3 = *c1;
- c2.model = model;
- switch (model) {
- case LIBCOLOUR_RGB:
- c2.rgb.with_gamma = 0;
- if (libcolour_get_rgb_colour_space(&c2.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB) < 0)
- return -1;
- return test_2convert(c1, &c2, &c3);
- case LIBCOLOUR_SRGB:
- c2.srgb.with_gamma = 0;
- if (r1 = test_2convert(c1, &c2, &c3), r1 < 0)
- return -1;
- c2.srgb.with_gamma = 1;
- if (r2 = test_2convert(c1, &c2, &c3), r2 < 0)
- return -1;
- return r1 & r2;
- case LIBCOLOUR_CIELCHUV:
- c2.cielchuv.one_revolution = 360.;
- /* fall though */
- case LIBCOLOUR_CIELUV:
- c2.cieluv.white.model = LIBCOLOUR_CIEXYZ;
- c2.cieluv.white.X = 1.0294;
- c2.cieluv.white.Y = 1;
- c2.cieluv.white.Z = 0.9118;
- return test_2convert(c1, &c2, &c3);
- case LIBCOLOUR_CIEUVW:
- c2.cieuvw.u0 = 0.37;
- c2.cieuvw.v0 = 0.30;
- return test_2convert(c1, &c2, &c3);
- default:
- return test_2convert(c1, &c2, &c3);
- }
-}
-
-
-static int test_2convert_1n(libcolour_model_t model, const char* model_name, double ch1, double ch2, double ch3)
-{
- libcolour_colour_t c1;
- int run, r, rc = 1;
-
- c1.model = model;
- switch (model) {
- case LIBCOLOUR_CIELAB:
- c1.cielab.L = ch1 * 100, c1.cielab.a = ch2 * 100, c1.cielab.b = ch3 * 100;
- break;
- case LIBCOLOUR_CIELCHUV:
- if (ch3 > 0.9999)
- return 1;
- /* fall though */
- default:
- c1.srgb.R = ch1, c1.srgb.G = ch2, c1.srgb.B = ch3;
- break;
- }
-
- for (run = 0;; run++) {
- switch (model) {
- case LIBCOLOUR_RGB:
- if (run == 2)
- return rc;
- c1.rgb.with_gamma = run;
- if (libcolour_get_rgb_colour_space(&c1.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB) < 0)
- return -1;
- break;
- case LIBCOLOUR_SRGB:
- if (run == 2)
- return rc;
- c1.srgb.with_gamma = run;
- break;
- case LIBCOLOUR_CIELCHUV:
- c1.cielchuv.one_revolution = run < 2 ? 360. : 2 * 3.14159265358979323846;
- c1.cielchuv.h *= c1.cielchuv.one_revolution;
- /* fall through */
- case LIBCOLOUR_CIELUV:
- c1.cieluv.white.model = LIBCOLOUR_CIEXYZ;
- if (run == 0) {
- c1.cieluv.white.X = 1.0294;
- c1.cieluv.white.Y = 1;
- c1.cieluv.white.Z = 0.9118;
- } else if (run == 1) {
- c1.cieluv.white.X = 1.03;
- c1.cieluv.white.Y = 0.8;
- c1.cieluv.white.Z = 0.92;
- } else if (run == 2 && model == LIBCOLOUR_CIELCHUV) {
- c1.cieluv.white.X = 1.0294;
- c1.cieluv.white.Y = 1;
- c1.cieluv.white.Z = 0.9118;
- } else if (run == 3 && model == LIBCOLOUR_CIELCHUV) {
- c1.cieluv.white.X = 1.03;
- c1.cieluv.white.Y = 0.8;
- c1.cieluv.white.Z = 0.92;
- } else {
- return rc;
- }
- break;
- case LIBCOLOUR_CIEUVW:
- if (run == 0) {
- c1.cieuvw.u0 = 0.37;
- c1.cieuvw.v0 = 0.30;
- } else if (run == 1) {
- c1.cieuvw.u0 = 0.47;
- c1.cieuvw.v0 = 0.31;
- } else {
- return rc;
- }
- break;
- default:
- if (run == 1)
- return rc;
- break;
- }
-
-#define X(ENUM, TYPE)\
- r = test_2convert_11(&c1, ENUM);\
- if (r < 0)\
- return -1;\
- if (!r)\
- printf("%s -> %s -> %s failed at run %i with (%lf, %lf, %lf)\n",\
- model_name, #ENUM, model_name, run, ch1, ch2, ch3), rc = 0;
- LIBCOLOUR_LIST_MODELS;
-#undef X
- }
-}
-
-
-static int test_2convert_nm(double ch1, double ch2, double ch3)
-{
- int r, rc = 1;
-
-#define X(ENUM, TYPE)\
- r = test_2convert_1n(ENUM, #ENUM, ch1, ch2, ch3);\
- if (r < 0)\
- return -1;\
- if (!r)\
- rc = 0;
- LIBCOLOUR_LIST_MODELS;
-#undef X
-
- return rc;
-}
-
-
-static int test_2convert_nm_all(void)
-{
-#define N 11
- int r, rc = 1;
- double ch1, ch2, ch3;
- int i, j, k;
-
- for (i = 0; i < N; i++) {
- for (j = 0; j < N; j++) {
- for (k = 0; k < N; k++) {
- ch1 = ((i + N / 2) % N) / (N - 1);
- ch2 = ((j + N / 2) % N) / (N - 1);
- ch3 = ((k + N / 2) % N) / (N - 1);
-
- if (ch1 > 0.999) ch1 = 1;
- if (ch2 > 0.999) ch2 = 1;
- if (ch3 > 0.999) ch3 = 1;
-
- r = test_2convert_nm(ch1, ch2, ch3);
- if (r < 0)
- return -1;
- if (!r)
- rc = 0;
- }
- }
- }
-
- return rc;
-#undef N
-}
-#endif
-
-
-static inline int ftest(double a, double b, double l)
-{
- return a - l <= b && b <= a + l;
-}
-
-
-#define test_rgb(COLOUR_SPACE, ...) test_rgb_(COLOUR_SPACE, #COLOUR_SPACE, __VA_ARGS__)
-static int test_rgb_(enum libcolour_rgb_colour_space colour_space, const char* name,
- double r, double g, double b, double x, double y, double z,
- double xx, double yy, double zz)
-{
- libcolour_colour_t c1, c2;
-
- c1.model = LIBCOLOUR_RGB;
- c2.model = LIBCOLOUR_CIEXYZ;
- c1.rgb.with_gamma = 0;
- c1.rgb.R = r, c1.rgb.G = g, c1.rgb.B = b;
- if (libcolour_get_rgb_colour_space(&c1.rgb, colour_space) ||
- libcolour_convert(&c1, &c2)) {
- printf("%s failed without gamma\n", name);
- return 0;
- } else if (!ftest(c2.ciexyz.X, x, 0.0001) ||
- !ftest(c2.ciexyz.Y, y, 0.0001) ||
- !ftest(c2.ciexyz.Z, z, 0.0001)) {
- printf("%s failed without gamma\n", name);
- return 0;
- }
-
- c1.model = LIBCOLOUR_RGB;
- c2.model = LIBCOLOUR_CIEXYZ;
- c1.rgb.with_gamma = 1;
- c1.rgb.R = r, c1.rgb.G = g, c1.rgb.B = b;
- if (libcolour_get_rgb_colour_space(&c1.rgb, colour_space) ||
- libcolour_convert(&c1, &c2)) {
- printf("%s failed with gamma\n", name);
- return 0;
- } else if (!ftest(c2.ciexyz.X, xx, 0.0001) ||
- !ftest(c2.ciexyz.Y, yy, 0.0001) ||
- !ftest(c2.ciexyz.Z, zz, 0.0001)) {
- printf("%s failed with gamma\n", name);
- return 0;
- }
-
- if (libcolour_convert(&c2, &c1) ||
- !ftest(c1.rgb.R, r, 0.0000001) ||
- !ftest(c1.rgb.G, g, 0.0000001) ||
- !ftest(c1.rgb.B, b, 0.0000001)) {
- printf("%s failed, encode\n", name);
- return 0;
- }
-
- return 1;
-}
-
-
-static int test_1convert(libcolour_colour_t *c1, libcolour_model_t model, const char* name1,
- const char* name2, double ch1, double ch2, double ch3, double d1,
- double d2, double d3)
-{
- libcolour_colour_t c2, c3;
- int rc = 1;
- c3 = *c1;
- c2.model = model;
- switch (model) {
- case LIBCOLOUR_CIELCHUV:
- c2.cielchuv.one_revolution = 360;
- /* fall-through */
- case LIBCOLOUR_CIELUV:
- c2.cieluv.white.X = 0.9504705586542831;
- c2.cieluv.white.Y = 1;
- c2.cieluv.white.Z = 1.0888287363958837;
- break;
- default:
- break;
- }
- if (libcolour_convert(c1, &c2))
- return -1;
- if (!ftest(c2.srgb.R, ch1, d1) ||
- !ftest(c2.srgb.G, ch2, d2) ||
- !ftest(c2.srgb.B, ch3, d3)) {
- printf("%s -> %s failed\n", name1, name2), rc = 0;
- rc = 0;
- }
- if (libcolour_convert(&c2, &c3))
- return -1;
- if (!ftest(c3.srgb.R, c1->srgb.R, 0.00000001) ||
- !ftest(c3.srgb.G, c1->srgb.G, 0.00000001) ||
- !ftest(c3.srgb.B, c1->srgb.B, 0.00000001)) {
- printf("%s -> %s failed\n", name2, name1), rc = 0;
- rc = 0;
- }
- return rc;
-}
-
-
-/**
- * Test libcolour
- *
- * @return 0: All tests passed
- * 1: At least one test fail
- * 2: An error occurred
- */
-int main(int argc, char* argv[])
-{
- int r, rc = 0;
- libcolour_colour_t c1, c2, c3, c4;
- double t1, t2;
- size_t n;
- char buf[sizeof(int) + sizeof(libcolour_rgb_t)];
-
-#ifndef SKIP_2CONVERT
- r = test_2convert_nm_all();
- if (r < 0)
- goto fail;
- if (!r)
- rc = 1;
-#else
- (void) r;
-#endif
-
- c1.model = LIBCOLOUR_SRGB;
- c1.srgb.R = 0.02, c1.srgb.G = 0.5, c1.srgb.B = 0.9;
- c1.srgb.with_gamma = 0;
- c2.model = LIBCOLOUR_RGB;
- c2.rgb.with_gamma = 0;
- if (libcolour_get_rgb_colour_space(&c2.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB)) {
- printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
- goto colour_spaces_done;
- }
- if (c2.rgb.white_r != 1 || c2.rgb.white_g != 1 || c2.rgb.white_b != 1 ||
- c2.rgb.red.model != LIBCOLOUR_CIEXYY || c2.rgb.green.model != LIBCOLOUR_CIEXYY ||
- c2.rgb.blue.model != LIBCOLOUR_CIEXYY || c2.rgb.white.model != LIBCOLOUR_CIEXYY ||
- c2.rgb.colour_space != LIBCOLOUR_RGB_COLOUR_SPACE_SRGB || c2.rgb.gamma != 2.4 ||
- c2.rgb.encoding_type != LIBCOLOUR_ENCODING_TYPE_REGULAR || c2.rgb.slope != 12.92
- || c2.rgb.offset != 0.055 || c2.rgb.transitioninv != c2.rgb.transition * c2.rgb.slope ||
- !ftest(c2.rgb.transition, 0.00313067, 0.00000001)) {
- printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
- goto colour_spaces_done;
- }
- t1 = c2.rgb.transition * c2.rgb.slope;
- t2 = (1 + c2.rgb.offset) * pow(c2.rgb.transition, 1 / c2.rgb.gamma) - c2.rgb.offset;
- if (!ftest(t1, t2, 0.0000001) || c2.rgb.white.Y != 1 ||
- c2.rgb.red.x != 0.64 || c2.rgb.red.y != 0.33 || !ftest(c2.rgb.red.Y, 0.21265, 0.00005) ||
- c2.rgb.green.x != 0.30 || c2.rgb.green.y != 0.60 || !ftest(c2.rgb.green.Y, 0.71515, 0.00005) ||
- c2.rgb.blue.x != 0.15 || c2.rgb.blue.y != 0.06 || !ftest(c2.rgb.blue.Y, 0.07215, 0.00005) ||
- !ftest(c2.rgb.white.x, 0.312727, 0.000001) || !ftest(c2.rgb.white.y, 0.329023, 0.000001) ||
- !ftest(c2.rgb.red.Y + c2.rgb.green.Y + c2.rgb.blue.Y, 1, 0.00000001)) {
- printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
- goto colour_spaces_done;
- }
- if (libcolour_convert(&c1, &c2) ||
- !ftest(c1.srgb.R, c2.rgb.R, 0.0000001) ||
- !ftest(c1.srgb.G, c2.rgb.G, 0.0000001) ||
- !ftest(c1.srgb.B, c2.rgb.B, 0.0000001)) {
- printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
- goto colour_spaces_done;
- }
- if (libcolour_convert(&c2, &c1) ||
- !ftest(c1.srgb.R, c2.rgb.R, 0.0000001) ||
- !ftest(c1.srgb.G, c2.rgb.G, 0.0000001) ||
- !ftest(c1.srgb.B, c2.rgb.B, 0.0000001)) {
- printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
- goto colour_spaces_done;
- }
-
- c1.model = c2.model = LIBCOLOUR_CIELAB;
- c1.cielab.L = 30, c2.cielab.L = 80;
- c1.cielab.a = 40, c2.cielab.a = -10;
- c1.cielab.b = -50, c2.cielab.b = -40;
- if (libcolour_delta_e(&c1, &c2, &t1) || !ftest(t1, 71.4142842854285, 0.0000001))
- printf("libcolour_delta_e failed\n"), rc = 0;
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_ADOBE_RGB, 0.3, 0.2, 0.9,
- 0.379497, 0.282430, 0.914245, 0.195429, 0.098925, 0.790020);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_APPLE_RGB, 0.3, 0.2, 0.9,
- 0.364212, 0.282789, 0.866008, 0.221570, 0.134028, 0.773782);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_BEST_RGB, 0.3, 0.2, 0.9,
- 0.345007, 0.246779, 0.736029, 0.151406, 0.064655, 0.647212);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_BETA_RGB, 0.3, 0.2, 0.9,
- 0.342837, 0.253386, 0.714198, 0.146437, 0.066824, 0.623382);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_BRUCE_RGB, 0.3, 0.2, 0.9,
- 0.368857, 0.276910, 0.915272, 0.191184, 0.096699, 0.791486);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_CIE_RGB, 0.3, 0.2, 0.9,
- 0.389293, 0.225188, 0.892857, 0.202678, 0.044608, 0.785312);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_COLORMATCH_RGB, 0.3, 0.2, 0.9,
- 0.337557, 0.274378, 0.651989, 0.186858, 0.123210, 0.581381);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_DON_RGB_4, 0.3, 0.2, 0.9,
- 0.344990, 0.251411, 0.727872, 0.150504, 0.066348, 0.638058);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB, 0.3, 0.2, 0.9,
- 0.353021, 0.286400, 0.695202, 0.196733, 0.134157, 0.630279);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2, 0.3, 0.2, 0.9,
- 0.353021, 0.286400, 0.695202, 0.149594, 0.097238, 0.579927);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_EKTA_SPACE_PS5, 0.3, 0.2, 0.9,
- 0.320377, 0.229160, 0.713291, 0.127134, 0.043253, 0.622392);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_NTSC_RGB, 0.3, 0.2, 0.9,
- 0.397081, 0.310031, 1.017821, 0.235650, 0.156675, 0.906707);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_PAL_SECAM_RGB, 0.3, 0.2, 0.9,
- 0.357972, 0.272130, 0.877151, 0.208580, 0.120322, 0.769056);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_PROPHOTO_RGB, 0.3, 0.2, 0.9,
- 0.294559, 0.228864, 0.742689, 0.124735, 0.072340, 0.682655);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_C_RGB, 0.3, 0.2, 0.9,
- 0.363595, 0.281822, 0.890350, 0.216773, 0.131311, 0.783348);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_WIDE_GAMUT_RGB, 0.3, 0.2, 0.9,
- 0.367485, 0.237631, 0.706442, 0.170318, 0.052665, 0.614915);
-
- rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_LIGHTROOM_RGB, 0.3, 0.2, 0.9,
- 0.294559, 0.228864, 0.742689, 0.294559, 0.228864, 0.742689);
-
- /*
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_D65
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_THEATER
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_625_LINE
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_525_LINE
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_709
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2020
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_PQ
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_PQ
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_PQ
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_HLG
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_HLG
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_HLG
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_SGI_RGB
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_240M_RGB
-
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MEASUREMENTS
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MATRIX
- TODO test LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_INV_MATRIX
- */
-
- /* TODO test transfer functions more rigorously */
-
- c1.model = c2.model = LIBCOLOUR_RGB;
- c3.model = c4.model = LIBCOLOUR_SRGB;
- c1.rgb.with_gamma = c3.srgb.with_gamma = 0;
- c2.rgb.with_gamma = c4.srgb.with_gamma = 1;
- if (libcolour_get_rgb_colour_space(&c1.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB))
- goto fail;
- if (libcolour_get_rgb_colour_space(&c2.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB))
- goto fail;
- for (t1 = 0; t1 <= 1; t1 += 0.001) {
- c1.rgb.R = c3.srgb.R = t1 - 1;
- c1.rgb.G = c3.srgb.G = t1;
- c1.rgb.B = c3.srgb.B = t1 + 1;
- if (libcolour_convert(&c1, &c2))
- goto fail;
- if (libcolour_convert(&c3, &c4))
- goto fail;
- if (!ftest(c2.rgb.R, c4.srgb.R, 0.00000000001) ||
- !ftest(c2.rgb.G, c4.srgb.G, 0.00000000001) ||
- !ftest(c2.rgb.B, c4.srgb.B, 0.00000000001) ||
- !ftest(libcolour_srgb_encode(c1.rgb.R), c2.srgb.R, 0.00000000001) ||
- !ftest(libcolour_srgb_encode(c1.rgb.G), c2.srgb.G, 0.00000000001) ||
- !ftest(libcolour_srgb_encode(c1.rgb.B), c2.srgb.B, 0.00000000001)) {
- printf("libcolour_srgb_encode failed\n"), rc = 0;
- goto colour_spaces_done;
- }
- if (libcolour_convert(&c2, &c1))
- goto fail;
- if (libcolour_convert(&c4, &c3))
- goto fail;
- if (!ftest(c1.rgb.R, c3.srgb.R, 0.00000000001) ||
- !ftest(c1.rgb.G, c3.srgb.G, 0.00000000001) ||
- !ftest(c1.rgb.B, c3.srgb.B, 0.00000000001) ||
- !ftest(libcolour_srgb_decode(c2.rgb.R), c1.srgb.R, 0.00000000001) ||
- !ftest(libcolour_srgb_decode(c2.rgb.G), c1.srgb.G, 0.00000000001) ||
- !ftest(libcolour_srgb_decode(c2.rgb.B), c1.srgb.B, 0.00000000001)) {
- printf("libcolour_srgb_decode failed\n"), rc = 0;
- goto colour_spaces_done;
- }
- }
-
- c1.rgb.encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
- c1.rgb.gamma = 2.2;
- c1.rgb.with_gamma = 0;
- c1.rgb.R = 0.25;
- c1.rgb.G = 0.5;
- c1.rgb.B = 0.75;
- c4 = c3 = c2 = c1;
- c2.rgb.with_gamma = c3.srgb.with_gamma = 1;
- c3.rgb.gamma = c4.rgb.gamma = 1.8;
- if (libcolour_convert(&c1, &c2))
- goto fail;
- if (libcolour_convert(&c2, &c3))
- goto fail;
- if (libcolour_convert(&c3, &c4))
- goto fail;
- if (!ftest(c1.rgb.R, c4.srgb.R, 0.0000001) ||
- !ftest(c1.rgb.G, c4.srgb.G, 0.0000001) ||
- !ftest(c1.rgb.B, c4.srgb.B, 0.0000001)) {
- printf("libcolour_convert failed to convert between two transfer functions\n"), rc = 0;
- goto colour_spaces_done;
- }
-
- if (libcolour_convert(&c1, &c4))
- goto fail;
- if (!ftest(c1.rgb.R, c4.srgb.R, 0.0000001) ||
- !ftest(c1.rgb.G, c4.srgb.G, 0.0000001) ||
- !ftest(c1.rgb.B, c4.srgb.B, 0.0000001)) {
- printf("libcolour_convert failed to convert when two different transfer functions are not applied\n"), rc = 0;
- goto colour_spaces_done;
- }
-
-#define TEST(TO, CH1, CH2, CH3, D1, D2, D3)\
- r = test_1convert(&c1, TO, "LIBCOLOUR_SRGB", #TO, CH1, CH2, CH3, D1, D2, D3);\
- if (r < 0)\
- goto fail;\
- if (!r)\
- rc = 1;
-
- c1.model = LIBCOLOUR_SRGB;
- c1.srgb.with_gamma = 0;
- c1.srgb.R = 0.3, c1.srgb.G = 0.2, c1.srgb.B = 0.7;
- TEST(LIBCOLOUR_CIEXYZ, 0.321558, 0.257355, 0.694851, 0.000001, 00000001, 0.000001);
- TEST(LIBCOLOUR_CIEXYY, 0.252447, 0.202043, 0.257355, 0.000001, 0.000001, 0.000001);
- TEST(LIBCOLOUR_CIELAB, 57.7851, 30.3599, -44.9740, 0.0001, 0.0001, 0.0001);
- TEST(LIBCOLOUR_CIELUV, 57.7851, 5.5723, -74.1571, 0.0001, 0.0001, 0.0001);
- TEST(LIBCOLOUR_CIELCHUV, 57.7851, 74.3662, 274.2972, 0.0001, 0.0001, 0.0001);
- /*
- * TODO convert test from sRGB to:
- * LChab 57.7851, 54.2622, 304.0215
- * YDbDr
- * YPbPr
- * YCgCo
- * CIE 1960 UCS
- * CIEUVW
- */
- c1.srgb.with_gamma = 0;
- c1.srgb.R = 32 / 255., c1.srgb.G = 65 / 255., c1.srgb.B = 32 / 255.;
- TEST(LIBCOLOUR_YIQ, 0.201, -0.035, -0.068, 0.001, 0.001, 0.001);
- TEST(LIBCOLOUR_YUV, 0.201, -0.037, -0.067, 0.001, 0.001, 0.001);
- /* XXX YCbCr(60, 118, 116) */
- /* XXX HSV(120, 0.508, 0.255) */
- /* XXX HSI(120, 0.256, 0.169) */
- /* XXX HSL(120, 0.340, 0.190) */
-
- if (rc)
- goto colour_spaces_done;
-
-#undef TEST
- colour_spaces_done:
-
- memset(&c1, 0, sizeof(c1));
- memset(&c2, 0, sizeof(c2));
- c1.model = LIBCOLOUR_RGB;
- if (libcolour_get_rgb_colour_space(&c1.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2))
- goto fail;
- c3 = c1;
- c1.rgb.to_encoded_red = NULL;
- c1.rgb.to_encoded_green = NULL;
- c1.rgb.to_encoded_blue = NULL;
- c1.rgb.to_decoded_red = NULL;
- c1.rgb.to_decoded_green = NULL;
- c1.rgb.to_decoded_blue = NULL;
- if (libcolour_marshal(&c1, NULL) > sizeof(buf)) {
- printf("libcolour_marshal failed\n"), rc = 0;
- goto marshal_done;
- }
- n = libcolour_marshal(&c1, buf);
- if (n > sizeof(buf)) {
- printf("libcolour_marshal failed\n"), rc = 0;
- goto marshal_done;
- }
- if (libcolour_unmarshal(NULL, buf) != n) {
- printf("libcolour_unmarshal failed\n"), rc = 0;
- goto marshal_done;
- }
- if (libcolour_unmarshal(&c2, buf) != n) {
- printf("libcolour_unmarshal failed\n"), rc = 0;
- goto marshal_done;
- }
- if (memcmp(&c2, &c3, sizeof(c2))) {
- printf("libcolour_(un)marshal failed\n"), rc = 0;
- goto marshal_done;
- }
- marshal_done:
-
- return rc;
- fail:
- perror(*argv);
- return 2;
- (void) argc;
-}