cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 1 | /* |
Cristy | 7ce65e7 | 2015-12-12 18:03:16 -0500 | [diff] [blame] | 2 | Copyright 1999-2016 ImageMagick Studio LLC, a non-profit organization |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 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 | */ |
Cristy | 83bceaa | 2016-06-03 20:39:35 -0400 | [diff] [blame^] | 18 | #ifndef MAGICKCORE_GEM_PRIVATE_H |
| 19 | #define MAGICKCORE_GEM_PRIVATE_H |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 20 | |
cristy | 0b9c0d8 | 2015-05-24 22:27:57 +0000 | [diff] [blame] | 21 | #include "MagickCore/pixel-accessor.h" |
cristy | 6398ec7 | 2013-11-28 02:00:27 +0000 | [diff] [blame] | 22 | |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 23 | #if defined(__cplusplus) || defined(c_plusplus) |
| 24 | extern "C" { |
| 25 | #endif |
| 26 | |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 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) |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 32 | |
cristy | 8ea8122 | 2011-09-04 10:33:32 +0000 | [diff] [blame] | 33 | extern MagickPrivate double |
| 34 | GenerateDifferentialNoise(RandomInfo *,const Quantum,const NoiseType, |
cristy | 9ed1f81 | 2011-10-08 02:00:08 +0000 | [diff] [blame] | 35 | const double); |
cristy | 8ea8122 | 2011-09-04 10:33:32 +0000 | [diff] [blame] | 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 | |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 42 | extern MagickPrivate void |
cristy | 722fc0c | 2012-08-04 23:15:43 +0000 | [diff] [blame] | 43 | ConvertHCLToRGB(const double,const double,const double,double *,double *, |
| 44 | double *), |
cristy | 9e2436a | 2013-05-02 20:35:59 +0000 | [diff] [blame] | 45 | ConvertHCLpToRGB(const double,const double,const double,double *,double *, |
| 46 | double *), |
cristy | 0a39a5c | 2012-06-27 12:51:45 +0000 | [diff] [blame] | 47 | ConvertHSBToRGB(const double,const double,const double,double *,double *, |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 48 | double *), |
cristy | af64eb2 | 2013-05-02 14:07:10 +0000 | [diff] [blame] | 49 | ConvertHSIToRGB(const double,const double,const double,double *,double *, |
| 50 | double *), |
cristy | 246c313 | 2013-05-02 16:35:53 +0000 | [diff] [blame] | 51 | ConvertHSVToRGB(const double,const double,const double,double *,double *, |
| 52 | double *), |
cristy | 0a39a5c | 2012-06-27 12:51:45 +0000 | [diff] [blame] | 53 | ConvertHWBToRGB(const double,const double,const double,double *,double *, |
cristy | 3094b7f | 2011-10-01 23:18:02 +0000 | [diff] [blame] | 54 | double *), |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 55 | ConvertLCHabToRGB(const double,const double,const double,double *,double *, |
| 56 | double *), |
| 57 | ConvertLCHuvToRGB(const double,const double,const double,double *,double *, |
cristy | b285095 | 2013-04-04 19:00:52 +0000 | [diff] [blame] | 58 | double *), |
cristy | 722fc0c | 2012-08-04 23:15:43 +0000 | [diff] [blame] | 59 | ConvertRGBToHCL(const double,const double,const double,double *,double *, |
| 60 | double *), |
cristy | 9e2436a | 2013-05-02 20:35:59 +0000 | [diff] [blame] | 61 | ConvertRGBToHCLp(const double,const double,const double,double *,double *, |
| 62 | double *), |
cristy | 0a39a5c | 2012-06-27 12:51:45 +0000 | [diff] [blame] | 63 | ConvertRGBToHSB(const double,const double,const double,double *,double *, |
cristy | 3094b7f | 2011-10-01 23:18:02 +0000 | [diff] [blame] | 64 | double *), |
cristy | af64eb2 | 2013-05-02 14:07:10 +0000 | [diff] [blame] | 65 | ConvertRGBToHSI(const double,const double,const double,double *,double *, |
| 66 | double *), |
cristy | 246c313 | 2013-05-02 16:35:53 +0000 | [diff] [blame] | 67 | ConvertRGBToHSV(const double,const double,const double,double *,double *, |
| 68 | double *), |
cristy | 0a39a5c | 2012-06-27 12:51:45 +0000 | [diff] [blame] | 69 | ConvertRGBToHWB(const double,const double,const double,double *,double *, |
cristy | b285095 | 2013-04-04 19:00:52 +0000 | [diff] [blame] | 70 | double *), |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 71 | ConvertRGBToLCHab(const double,const double,const double,double *,double *, |
| 72 | double *), |
| 73 | ConvertRGBToLCHuv(const double,const double,const double,double *,double *, |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 74 | double *); |
| 75 | |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 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); |
cristy | 35605e9 | 2013-05-06 11:33:57 +0000 | [diff] [blame] | 87 | y=(L+16.0)/116.0; |
| 88 | x=y+a/500.0; |
| 89 | z=y-b/200.0; |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 90 | if ((x*x*x) > CIEEpsilon) |
| 91 | x=(x*x*x); |
| 92 | else |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 93 | x=(116.0*x-16.0)/CIEK; |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 94 | if ((y*y*y) > CIEEpsilon) |
| 95 | y=(y*y*y); |
| 96 | else |
cristy | 35605e9 | 2013-05-06 11:33:57 +0000 | [diff] [blame] | 97 | y=L/CIEK; |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 98 | if ((z*z*z) > CIEEpsilon) |
| 99 | z=(z*z*z); |
| 100 | else |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 101 | z=(116.0*z-16.0)/CIEK; |
| 102 | *X=D65X*x; |
| 103 | *Y=D65Y*y; |
| 104 | *Z=D65Z*z; |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 105 | } |
| 106 | |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 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); |
cristy | a871691 | 2013-05-06 12:04:10 +0000 | [diff] [blame] | 113 | if (L > (CIEK*CIEEpsilon)) |
| 114 | *Y=(double) pow((L+16.0)/116.0,3.0); |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 115 | else |
cristy | a871691 | 2013-05-06 12:04:10 +0000 | [diff] [blame] | 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))+ |
cristy | 887d513 | 2013-05-06 20:05:35 +0000 | [diff] [blame] | 118 | 5.0*(*Y))/((((52.0*L/(u+13.0*L*(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z))))-1.0)/ |
cristy | a871691 | 2013-05-06 12:04:10 +0000 | [diff] [blame] | 119 | 3.0)-(-1.0/3.0)); |
cristy | 887d513 | 2013-05-06 20:05:35 +0000 | [diff] [blame] | 120 | *Z=(*X*(((52.0*L/(u+13.0*L*(4.0*D65X/(D65X+15.0*D65Y+3.0*D65Z))))-1.0)/3.0))- |
cristy | a871691 | 2013-05-06 12:04:10 +0000 | [diff] [blame] | 121 | 5.0*(*Y); |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 122 | } |
| 123 | |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 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 | |
cristy | 8eccf57 | 2013-05-03 17:30:54 +0000 | [diff] [blame] | 132 | /* |
| 133 | Convert RGB to XYZ colorspace. |
| 134 | */ |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 135 | assert(X != (double *) NULL); |
| 136 | assert(Y != (double *) NULL); |
| 137 | assert(Z != (double *) NULL); |
cristy | 58f7537 | 2013-05-05 17:34:43 +0000 | [diff] [blame] | 138 | r=QuantumScale*DecodePixelGamma(red); |
| 139 | g=QuantumScale*DecodePixelGamma(green); |
| 140 | b=QuantumScale*DecodePixelGamma(blue); |
cristy | d65a94d | 2015-01-28 22:20:48 +0000 | [diff] [blame] | 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; |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 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; |
cristy | c32cb02 | 2013-05-06 13:32:31 +0000 | [diff] [blame] | 153 | |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 154 | assert(L != (double *) NULL); |
| 155 | assert(a != (double *) NULL); |
| 156 | assert(b != (double *) NULL); |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 157 | if ((X/D65X) > CIEEpsilon) |
| 158 | x=pow(X/D65X,1.0/3.0); |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 159 | else |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 160 | x=(CIEK*X/D65X+16.0)/116.0; |
| 161 | if ((Y/D65Y) > CIEEpsilon) |
| 162 | y=pow(Y/D65Y,1.0/3.0); |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 163 | else |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 164 | y=(CIEK*Y/D65Y+16.0)/116.0; |
| 165 | if ((Z/D65Z) > CIEEpsilon) |
| 166 | z=pow(Z/D65Z,1.0/3.0); |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 167 | else |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 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; |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 172 | } |
| 173 | |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 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); |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 183 | if ((Y/D65Y) > CIEEpsilon) |
| 184 | *L=(double) (116.0*pow(Y/D65Y,1.0/3.0)-16.0); |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 185 | else |
cristy | 15ec6a6 | 2013-04-20 01:11:51 +0000 | [diff] [blame] | 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; |
cristy | df42b17 | 2013-04-05 13:35:37 +0000 | [diff] [blame] | 193 | } |
| 194 | |
cristy | 48348e7 | 2013-05-07 01:21:28 +0000 | [diff] [blame] | 195 | static inline void ConvertXYZToRGB(const double X,const double Y,const double Z, |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 196 | double *red,double *green,double *blue) |
| 197 | { |
| 198 | double |
| 199 | b, |
| 200 | g, |
| 201 | r; |
| 202 | |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 203 | assert(red != (double *) NULL); |
| 204 | assert(green != (double *) NULL); |
| 205 | assert(blue != (double *) NULL); |
cristy | d65a94d | 2015-01-28 22:20:48 +0000 | [diff] [blame] | 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; |
cristy | 58f7537 | 2013-05-05 17:34:43 +0000 | [diff] [blame] | 209 | *red=EncodePixelGamma(QuantumRange*r); |
| 210 | *green=EncodePixelGamma(QuantumRange*g); |
| 211 | *blue=EncodePixelGamma(QuantumRange*b); |
cristy | d803bd0 | 2013-04-03 19:10:51 +0000 | [diff] [blame] | 212 | } |
| 213 | |
cristy | d1dd6e4 | 2011-09-04 01:46:08 +0000 | [diff] [blame] | 214 | #if defined(__cplusplus) || defined(c_plusplus) |
| 215 | } |
| 216 | #endif |
| 217 | |
| 218 | #endif |