granicus.if.org Git - imagemagick/blob - MagickCore/gem-private.h

   1 /*
   2   Copyright 1999-2013 ImageMagick Studio LLC, a non-profit organization
   3   dedicated to making software imaging solutions freely available.
   4
   5   You may not use this file except in compliance with the License.
   6   obtain a copy of the License at
   7
   8     http://www.imagemagick.org/script/license.php
   9
  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.
  15
  16   MagickCore private graphic gems methods.
  17 */
  18 #ifndef _MAGICKCORE_GEM_PRIVATE_H
  19 #define _MAGICKCORE_GEM_PRIVATE_H
  20
  21 #if defined(__cplusplus) || defined(c_plusplus)
  22 extern "C" {
  23 #endif
  24
  25 #include "MagickCore/pixel-private.h"
  26
  27 #define D65X  0.950456
  28 #define D65Y  1.0
  29 #define D65Z  1.088754
  30 #define CIEEpsilon  (216.0/24389.0)
  31 #define CIEK  (24389.0/27.0)
  32
  33 extern MagickPrivate double
  34   GenerateDifferentialNoise(RandomInfo *,const Quantum,const NoiseType,
  35     const double);
  36
  37 extern MagickPrivate size_t
  38   GetOptimalKernelWidth(const double,const double),
  39   GetOptimalKernelWidth1D(const double,const double),
  40   GetOptimalKernelWidth2D(const double,const double);
  41
  42 extern MagickPrivate void
  43   ConvertHCLToRGB(const double,const double,const double,double *,double *,
  44     double *),
  45   ConvertHCLpToRGB(const double,const double,const double,double *,double *,
  46     double *),
  47   ConvertHSBToRGB(const double,const double,const double,double *,double *,
  48     double *),
  49   ConvertHSIToRGB(const double,const double,const double,double *,double *,
  50     double *),
  51   ConvertHSVToRGB(const double,const double,const double,double *,double *,
  52     double *),
  53   ConvertHWBToRGB(const double,const double,const double,double *,double *,
  54     double *),
  55   ConvertLCHabToRGB(const double,const double,const double,double *,double *,
  56     double *),
  57   ConvertLCHuvToRGB(const double,const double,const double,double *,double *,
  58     double *),
  59   ConvertRGBToHCL(const double,const double,const double,double *,double *,
  60     double *),
  61   ConvertRGBToHCLp(const double,const double,const double,double *,double *,
  62     double *),
  63   ConvertRGBToHSB(const double,const double,const double,double *,double *,
  64     double *),
  65   ConvertRGBToHSI(const double,const double,const double,double *,double *,
  66     double *),
  67   ConvertRGBToHSV(const double,const double,const double,double *,double *,
  68     double *),
  69   ConvertRGBToHWB(const double,const double,const double,double *,double *,
  70     double *),
  71   ConvertRGBToLCHab(const double,const double,const double,double *,double *,
  72     double *),
  73   ConvertRGBToLCHuv(const double,const double,const double,double *,double *,
  74     double *);
  75
  76 static inline void ConvertLabToXYZ(const double L,const double a,const double b,
  77   double *X,double *Y,double *Z)
  78 {
  79   double
  80     x,
  81     y,
  82     z;
  83
  84   assert(X != (double *) NULL);
  85   assert(Y != (double *) NULL);
  86   assert(Z != (double *) NULL);
  87   y=(100.0*L+16.0)/116.0;
  88   x=y+255.0*(a-0.5)/500.0;
  89   z=y-255.0*(b-0.5)/200.0;
  90   if ((x*x*x) > CIEEpsilon)
  91     x=(x*x*x);
  92   else
  93     x=(116.0*x-16.0)/CIEK;
  94   if ((y*y*y) > CIEEpsilon)
  95     y=(y*y*y);
  96   else
  97     y=(100.0*L)/CIEK;
  98   if ((z*z*z) > CIEEpsilon)
  99     z=(z*z*z);
 100   else
 101     z=(116.0*z-16.0)/CIEK;
 102   *X=D65X*x;
 103   *Y=D65Y*y;
 104   *Z=D65Z*z;
 105 }
 106
 107 static inline void ConvertLuvToXYZ(const double L,const double u,const double v,
 108   double *X,double *Y,double *Z)
 109 {
 110   assert(X != (double *) NULL);
 111   assert(Y != (double *) NULL);
 112   assert(Z != (double *) NULL);
 113   if ((100.0*L) > (CIEK*CIEEpsilon))
 114     *Y=(double) pow(((100.0*L)+16.0)/116.0,3.0);
 115   else
 116     *Y=(100.0*L)/CIEK;
 117   *X=((*Y*((39.0*(100.0*L)/((262.0*v-140.0)+13.0*(100.0*L)*(9.0*D65Y/
 118     (D65X+15.0*D65Y+3.0*D65Z))))-5.0))+5.0*(*Y))/((((52.0f*(100.0*L)/
 119     ((354.0*u-134.0)+13.0*(100.0*L)*(4.0*D65X/(D65X+15.0*D65Y+3.0*
 120     D65Z))))-1.0f)/3.0)-(-1.0f/3.0));
 121   *Z=(*X*(((52.0f*(100.0*L)/((354.0*u-134.0)+13.0*(100.0*L)*(4.0*D65X/
 122     (D65X+15.0*D65Y+3.0*D65Z))))-1.0f)/3.0))-5.0*(*Y);
 123 }
 124
 125 static inline void ConvertRGBToXYZ(const double red,const double green,
 126   const double blue,double *X,double *Y,double *Z)
 127 {
 128   double
 129     b,
 130     g,
 131     r;
 132
 133   assert(X != (double *) NULL);
 134   assert(Y != (double *) NULL);
 135   assert(Z != (double *) NULL);
 136   r=QuantumScale*red;
 137   g=QuantumScale*green;
 138   b=QuantumScale*blue;
 139   *X=0.41239558896741421610*r+0.35758343076371481710*g+0.18049264738170157350*b;
 140   *Y=0.21258623078559555160*r+0.71517030370341084990*g+0.07220049864333622685*b;
 141   *Z=0.01929721549174694484*r+0.11918386458084853180*g+0.95049712513157976600*b;
 142 }
 143
 144 static inline void ConvertXYZToLab(const double X,const double Y,const double Z,
 145   double *L,double *a,double *b)
 146 {
 147   double
 148     x,
 149     y,
 150     z;
 151   assert(L != (double *) NULL);
 152   assert(a != (double *) NULL);
 153   assert(b != (double *) NULL);
 154   if ((X/D65X) > CIEEpsilon)
 155     x=pow(X/D65X,1.0/3.0);
 156   else
 157     x=(CIEK*X/D65X+16.0)/116.0;
 158   if ((Y/D65Y) > CIEEpsilon)
 159     y=pow(Y/D65Y,1.0/3.0);
 160   else
 161     y=(CIEK*Y/D65Y+16.0)/116.0;
 162   if ((Z/D65Z) > CIEEpsilon)
 163     z=pow(Z/D65Z,1.0/3.0);
 164   else
 165     z=(CIEK*Z/D65Z+16.0)/116.0;
 166   *L=((116.0*y)-16.0)/100.0;
 167   *a=(500.0*(x-y))/255.0+0.5;
 168   *b=(200.0*(y-z))/255.0+0.5;
 169 }
 170
 171 static inline void ConvertXYZToLuv(const double X,const double Y,const double Z,
 172   double *L,double *u,double *v)
 173 {
 174   double
 175     alpha;
 176
 177   assert(L != (double *) NULL);
 178   assert(u != (double *) NULL);
 179   assert(v != (double *) NULL);
 180   if ((Y/D65Y) > CIEEpsilon)
 181     *L=(double) (116.0*pow(Y/D65Y,1.0/3.0)-16.0);
 182   else
 183     *L=CIEK*(Y/D65Y);
 184   alpha=PerceptibleReciprocal(X+15.0*Y+3.0*Z);
 185   *u=13.0*(*L)*((4.0*alpha*X)-(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z)));
 186   *v=13.0*(*L)*((9.0*alpha*Y)-(9.0*D65Y/(D65X+15.0*D65Y+3.0*D65Z)));
 187   *L/=100.0;
 188   *u=(*u+134.0)/354.0;
 189   *v=(*v+140.0)/262.0;
 190 }
 191
 192 static inline void ConvertXYZToRGB(const double x,const double y,const double z,
 193   double *red,double *green,double *blue)
 194 {
 195   double
 196     b,
 197     g,
 198     r;
 199
 200   /*
 201     Convert XYZ to RGB colorspace.
 202   */
 203   assert(red != (double *) NULL);
 204   assert(green != (double *) NULL);
 205   assert(blue != (double *) NULL);
 206   r=3.2406*x-1.5372*y-0.4986*z;
 207   g=(-0.9689*x+1.8758*y+0.0415*z);
 208   b=0.0557*x-0.2040*y+1.0570*z;
 209   *red=QuantumRange*r;
 210   *green=QuantumRange*g;
 211   *blue=QuantumRange*b;
 212 }
 213
 214
 215 #if defined(__cplusplus) || defined(c_plusplus)
 216 }
 217 #endif
 218
 219 #endif