18#ifndef MAGICKCORE_GEM_PRIVATE_H
19#define MAGICKCORE_GEM_PRIVATE_H
21#include "magick/pixel-accessor.h"
22#include "magick/visual-effects.h"
24#if defined(__cplusplus) || defined(c_plusplus)
31#define CIEEpsilon (216.0/24389.0)
32#define CIEK (24389.0/27.0)
34static inline void ConvertLabToXYZ(
const double L,
const double a,
const double b,
35 double *X,
double *Y,
double *Z)
42 assert(X != (
double *) NULL);
43 assert(Y != (
double *) NULL);
44 assert(Z != (
double *) NULL);
48 if ((x*x*x) > CIEEpsilon)
51 x=(116.0*x-16.0)/CIEK;
52 if (L > (CIEK*CIEEpsilon))
56 if ((z*z*z) > CIEEpsilon)
59 z=(116.0*z-16.0)/CIEK;
65static inline void ConvertXYZToLuv(
const double X,
const double Y,
const double Z,
66 double *L,
double *u,
double *v)
71 assert(L != (
double *) NULL);
72 assert(u != (
double *) NULL);
73 assert(v != (
double *) NULL);
74 if ((Y/D65Y) > CIEEpsilon)
75 *L=(double) (116.0*pow(Y/D65Y,1.0/3.0)-16.0);
78 alpha=PerceptibleReciprocal(X+15.0*Y+3.0*Z);
79 *u=13.0*(*L)*((4.0*alpha*X)-(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z)));
80 *v=13.0*(*L)*((9.0*alpha*Y)-(9.0*D65Y/(D65X+15.0*D65Y+3.0*D65Z)));
86static inline void ConvertRGBToXYZ(
const Quantum red,
const Quantum green,
87 const Quantum blue,
double *X,
double *Y,
double *Z)
94 assert(X != (
double *) NULL);
95 assert(Y != (
double *) NULL);
96 assert(Z != (
double *) NULL);
97 r=QuantumScale*DecodePixelGamma((MagickRealType) red);
98 g=QuantumScale*DecodePixelGamma((MagickRealType) green);
99 b=QuantumScale*DecodePixelGamma((MagickRealType) blue);
100 *X=(0.41239079926595*r)+(0.35758433938387*g)+(0.18048078840183*b);
101 *Y=(0.21263900587151*r)+(0.71516867876775*g)+(0.072192315360733*b);
102 *Z=(0.019330818715591*r)+(0.11919477979462*g)+(0.95053215224966*b);
105static inline void ConvertXYZToLab(
const double X,
const double Y,
const double Z,
106 double *L,
double *a,
double *b)
113 assert(L != (
double *) NULL);
114 assert(a != (
double *) NULL);
115 assert(b != (
double *) NULL);
116 if ((X/D65X) > CIEEpsilon)
117 x=pow(X/D65X,1.0/3.0);
119 x=(CIEK*X/D65X+16.0)/116.0;
120 if ((Y/D65Y) > CIEEpsilon)
121 y=pow(Y/D65Y,1.0/3.0);
123 y=(CIEK*Y/D65Y+16.0)/116.0;
124 if ((Z/D65Z) > CIEEpsilon)
125 z=pow(Z/D65Z,1.0/3.0);
127 z=(CIEK*Z/D65Z+16.0)/116.0;
128 *L=((116.0*y)-16.0)/100.0;
129 *a=(500.0*(x-y))/255.0+0.5;
130 *b=(200.0*(y-z))/255.0+0.5;
133static inline void ConvertLuvToXYZ(
const double L,
const double u,
const double v,
134 double *X,
double *Y,
double *Z)
139 assert(X != (
double *) NULL);
140 assert(Y != (
double *) NULL);
141 assert(Z != (
double *) NULL);
142 if (L > (CIEK*CIEEpsilon))
143 *Y=(double) pow((L+16.0)/116.0,3.0);
146 gamma=PerceptibleReciprocal((((52.0*L*PerceptibleReciprocal(u+13.0*L*
147 (4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z))))-1.0)/3.0)-(-1.0/3.0));
148 *X=gamma*((*Y*((39.0*L*PerceptibleReciprocal(v+13.0*L*(9.0*D65Y/
149 (D65X+15.0*D65Y+3.0*D65Z))))-5.0))+5.0*(*Y));
150 *Z=(*X*(((52.0*L*PerceptibleReciprocal(u+13.0*L*(4.0*D65X/
151 (D65X+15.0*D65Y+3.0*D65Z))))-1.0)/3.0))-5.0*(*Y);
154static inline void ConvertXYZToRGB(
const double X,
const double Y,
const double Z,
155 Quantum *red,Quantum *green,Quantum *blue)
162 assert(red != (Quantum *) NULL);
163 assert(green != (Quantum *) NULL);
164 assert(blue != (Quantum *) NULL);
165 r=(3.240969941904521*X)+(-1.537383177570093*Y)+(-0.498610760293*Z);
166 g=(-0.96924363628087*X)+(1.87596750150772*Y)+(0.041555057407175*Z);
167 b=(0.055630079696993*X)+(-0.20397695888897*Y)+(1.056971514242878*Z);
168 *red=ClampToQuantum((MagickRealType) EncodePixelGamma((MagickRealType)
170 *green=ClampToQuantum((MagickRealType) EncodePixelGamma((MagickRealType)
172 *blue=ClampToQuantum((MagickRealType) EncodePixelGamma((MagickRealType)
176#if defined(__cplusplus) || defined(c_plusplus)