2 Copyright 1999-2017 ImageMagick Studio LLC, a non-profit organization
3 dedicated to making software imaging solutions freely available.
5 You may not use this file except in compliance with the License.
6 obtain a copy of the License at
8 https://www.imagemagick.org/script/license.php
10 Unless required by applicable law or agreed to in writing, software
11 distributed under the License is distributed on an "AS IS" BASIS,
12 WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 See the License for the specific language governing permissions and
14 limitations under the License.
16 MagickCore private graphic gems methods.
18 #ifndef MAGICKCORE_GEM_PRIVATE_H
19 #define MAGICKCORE_GEM_PRIVATE_H
21 #include "MagickCore/pixel-accessor.h"
23 #if defined(__cplusplus) || defined(c_plusplus)
30 #define CIEEpsilon (216.0/24389.0)
31 #define CIEK (24389.0/27.0)
33 extern MagickPrivate double
34 GenerateDifferentialNoise(RandomInfo *,const Quantum,const NoiseType,
37 extern MagickPrivate size_t
38 GetOptimalKernelWidth(const double,const double),
39 GetOptimalKernelWidth1D(const double,const double),
40 GetOptimalKernelWidth2D(const double,const double);
42 extern MagickPrivate void
43 ConvertHCLToRGB(const double,const double,const double,double *,double *,
45 ConvertHCLpToRGB(const double,const double,const double,double *,double *,
47 ConvertHSBToRGB(const double,const double,const double,double *,double *,
49 ConvertHSIToRGB(const double,const double,const double,double *,double *,
51 ConvertHSVToRGB(const double,const double,const double,double *,double *,
53 ConvertHWBToRGB(const double,const double,const double,double *,double *,
55 ConvertLCHabToRGB(const double,const double,const double,double *,double *,
57 ConvertLCHuvToRGB(const double,const double,const double,double *,double *,
59 ConvertRGBToHCL(const double,const double,const double,double *,double *,
61 ConvertRGBToHCLp(const double,const double,const double,double *,double *,
63 ConvertRGBToHSB(const double,const double,const double,double *,double *,
65 ConvertRGBToHSI(const double,const double,const double,double *,double *,
67 ConvertRGBToHSV(const double,const double,const double,double *,double *,
69 ConvertRGBToHWB(const double,const double,const double,double *,double *,
71 ConvertRGBToLCHab(const double,const double,const double,double *,double *,
73 ConvertRGBToLCHuv(const double,const double,const double,double *,double *,
76 static inline void ConvertLabToXYZ(const double L,const double a,const double b,
77 double *X,double *Y,double *Z)
84 assert(X != (double *) NULL);
85 assert(Y != (double *) NULL);
86 assert(Z != (double *) NULL);
90 if ((x*x*x) > CIEEpsilon)
93 x=(116.0*x-16.0)/CIEK;
94 if ((y*y*y) > CIEEpsilon)
98 if ((z*z*z) > CIEEpsilon)
101 z=(116.0*z-16.0)/CIEK;
107 static inline void ConvertLuvToXYZ(const double L,const double u,const double v,
108 double *X,double *Y,double *Z)
110 assert(X != (double *) NULL);
111 assert(Y != (double *) NULL);
112 assert(Z != (double *) NULL);
113 if (L > (CIEK*CIEEpsilon))
114 *Y=(double) pow((L+16.0)/116.0,3.0);
117 *X=((*Y*((39.0*L/(v+13.0*L*(9.0*D65Y/(D65X+15.0*D65Y+3.0*D65Z))))-5.0))+
118 5.0*(*Y))/((((52.0*L/(u+13.0*L*(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z))))-1.0)/
120 *Z=(*X*(((52.0*L/(u+13.0*L*(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z))))-1.0)/3.0))-
124 static inline void ConvertRGBToXYZ(const double red,const double green,
125 const double blue,double *X,double *Y,double *Z)
133 Convert RGB to XYZ colorspace.
135 assert(X != (double *) NULL);
136 assert(Y != (double *) NULL);
137 assert(Z != (double *) NULL);
138 r=QuantumScale*DecodePixelGamma(red);
139 g=QuantumScale*DecodePixelGamma(green);
140 b=QuantumScale*DecodePixelGamma(blue);
141 *X=0.4124564*r+0.3575761*g+0.1804375*b;
142 *Y=0.2126729*r+0.7151522*g+0.0721750*b;
143 *Z=0.0193339*r+0.1191920*g+0.9503041*b;
146 static inline void ConvertXYZToLab(const double X,const double Y,const double Z,
147 double *L,double *a,double *b)
154 assert(L != (double *) NULL);
155 assert(a != (double *) NULL);
156 assert(b != (double *) NULL);
157 if ((X/D65X) > CIEEpsilon)
158 x=pow(X/D65X,1.0/3.0);
160 x=(CIEK*X/D65X+16.0)/116.0;
161 if ((Y/D65Y) > CIEEpsilon)
162 y=pow(Y/D65Y,1.0/3.0);
164 y=(CIEK*Y/D65Y+16.0)/116.0;
165 if ((Z/D65Z) > CIEEpsilon)
166 z=pow(Z/D65Z,1.0/3.0);
168 z=(CIEK*Z/D65Z+16.0)/116.0;
169 *L=((116.0*y)-16.0)/100.0;
170 *a=(500.0*(x-y))/255.0+0.5;
171 *b=(200.0*(y-z))/255.0+0.5;
174 static inline void ConvertXYZToLuv(const double X,const double Y,const double Z,
175 double *L,double *u,double *v)
180 assert(L != (double *) NULL);
181 assert(u != (double *) NULL);
182 assert(v != (double *) NULL);
183 if ((Y/D65Y) > CIEEpsilon)
184 *L=(double) (116.0*pow(Y/D65Y,1.0/3.0)-16.0);
187 alpha=PerceptibleReciprocal(X+15.0*Y+3.0*Z);
188 *u=13.0*(*L)*((4.0*alpha*X)-(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z)));
189 *v=13.0*(*L)*((9.0*alpha*Y)-(9.0*D65Y/(D65X+15.0*D65Y+3.0*D65Z)));
195 static inline void ConvertXYZToRGB(const double X,const double Y,const double Z,
196 double *red,double *green,double *blue)
203 assert(red != (double *) NULL);
204 assert(green != (double *) NULL);
205 assert(blue != (double *) NULL);
206 r=3.2404542*X-1.5371385*Y-0.4985314*Z;
207 g=(-0.9692660)*X+1.8760108*Y+0.0415560*Z;
208 b=0.0556434*X-0.2040259*Y+1.0572252*Z;
209 *red=EncodePixelGamma(QuantumRange*r);
210 *green=EncodePixelGamma(QuantumRange*g);
211 *blue=EncodePixelGamma(QuantumRange*b);
214 #if defined(__cplusplus) || defined(c_plusplus)