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

   1 /*
   2   Copyright 1999-2014 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 #include "MagickCore/pixel-private.h"
  22
  23 #if defined(__cplusplus) || defined(c_plusplus)
  24 extern "C" {
  25 #endif
  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=(L+16.0)/116.0;
  88   x=y+a/500.0;
  89   z=y-b/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=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 (L > (CIEK*CIEEpsilon))
 114     *Y=(double) pow((L+16.0)/116.0,3.0);
 115   else
 116     *Y=L/CIEK;
 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)/
 119     3.0)-(-1.0/3.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))-
 121     5.0*(*Y);
 122 }
 123
 124 static inline void ConvertRGBToXYZ(const double red,const double green,
 125   const double blue,double *X,double *Y,double *Z)
 126 {
 127   double
 128     b,
 129     g,
 130     r;
 131
 132   /*
 133     Convert RGB to XYZ colorspace.
 134   */
 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.41239558896741421610*r+0.35758343076371481710*g+0.18049264738170157350*b;
 142   *Y=0.21258623078559555160*r+0.71517030370341084990*g+0.07220049864333622685*b;
 143   *Z=0.01929721549174694484*r+0.11918386458084853180*g+0.95049712513157976600*b;
 144 }
 145
 146 static inline void ConvertXYZToLab(const double X,const double Y,const double Z,
 147   double *L,double *a,double *b)
 148 {
 149   double
 150     x,
 151     y,
 152     z;
 153
 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);
 159   else
 160     x=(CIEK*X/D65X+16.0)/116.0;
 161   if ((Y/D65Y) > CIEEpsilon)
 162     y=pow(Y/D65Y,1.0/3.0);
 163   else
 164     y=(CIEK*Y/D65Y+16.0)/116.0;
 165   if ((Z/D65Z) > CIEEpsilon)
 166     z=pow(Z/D65Z,1.0/3.0);
 167   else
 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;
 172 }
 173
 174 static inline void ConvertXYZToLuv(const double X,const double Y,const double Z,
 175   double *L,double *u,double *v)
 176 {
 177   double
 178     alpha;
 179
 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);
 185   else
 186     *L=CIEK*(Y/D65Y);
 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)));
 190   *L/=100.0;
 191   *u=(*u+134.0)/354.0;
 192   *v=(*v+140.0)/262.0;
 193 }
 194
 195 static inline void ConvertXYZToRGB(const double X,const double Y,const double Z,
 196   double *red,double *green,double *blue)
 197 {
 198   double
 199     b,
 200     g,
 201     r;
 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=EncodePixelGamma(QuantumRange*r);
 210   *green=EncodePixelGamma(QuantumRange*g);
 211   *blue=EncodePixelGamma(QuantumRange*b);
 212 }
 213
 214 #if defined(__cplusplus) || defined(c_plusplus)
 215 }
 216 #endif
 217
 218 #endif