|
1
|
|
|
<?php |
|
2
|
|
|
|
|
3
|
|
|
/* |
|
4
|
|
|
|
|
5
|
|
|
Color Space : v1.25 : 2009.04.20 |
|
6
|
|
|
———————————————————————————————— |
|
7
|
|
|
XYZ <-> Luv <-> LCHuv |
|
8
|
|
|
XYZ <-> Lab <-> LCHab |
|
9
|
|
|
RGB <-> XYZ <-> xyY |
|
10
|
|
|
RGB <-> HSL |
|
11
|
|
|
RGB <-> HSV <-> RYB |
|
12
|
|
|
RGB <-> CMY <-> CMYK |
|
13
|
|
|
RGB <-> HEX <-> STRING |
|
14
|
|
|
|
|
15
|
|
|
*/ |
|
16
|
|
|
|
|
17
|
|
|
class ColorSpace |
|
18
|
|
|
{ |
|
19
|
|
|
private $RYB_H = []; |
|
20
|
|
|
private $H_RYB = []; |
|
21
|
|
|
private $_XYZ_RGB = []; |
|
22
|
|
|
private $White; |
|
23
|
|
|
private $_xyY; |
|
24
|
|
|
|
|
25
|
|
|
public function __construct() |
|
26
|
|
|
{ |
|
27
|
|
|
$wheel = [ |
|
28
|
|
|
[0, 0], |
|
29
|
|
|
[15, 8], // Red |
|
30
|
|
|
[30, 17], |
|
31
|
|
|
[45, 26], // Orange |
|
32
|
|
|
[60, 34], |
|
33
|
|
|
[75, 41], // Yellow |
|
34
|
|
|
[90, 48], |
|
35
|
|
|
[105, 54], // Lime |
|
36
|
|
|
[120, 60], |
|
37
|
|
|
[135, 81], // Green |
|
38
|
|
|
[150, 103], |
|
39
|
|
|
[165, 123], // Teal |
|
40
|
|
|
[180, 138], |
|
41
|
|
|
[195, 155], // Cyan |
|
42
|
|
|
[210, 171], |
|
43
|
|
|
[225, 187], // Azure |
|
44
|
|
|
[240, 204], |
|
45
|
|
|
[255, 219], // Blue |
|
46
|
|
|
[270, 234], |
|
47
|
|
|
[285, 251], // Indigo |
|
48
|
|
|
[300, 267], |
|
49
|
|
|
[315, 282], // Purple |
|
50
|
|
|
[330, 298], |
|
51
|
|
|
[345, 329], // Pink |
|
52
|
|
|
[360, 0], |
|
53
|
|
|
]; |
|
54
|
|
|
|
|
55
|
|
|
$a; |
|
|
|
|
|
|
56
|
|
|
$b; |
|
|
|
|
|
|
57
|
|
|
$i; |
|
|
|
|
|
|
58
|
|
|
$H; |
|
|
|
|
|
|
59
|
|
|
$_H; |
|
|
|
|
|
|
60
|
|
|
$H_; |
|
|
|
|
|
|
61
|
|
|
|
|
62
|
|
|
for ($H = 0; $H < 360; $H++) { |
|
63
|
|
|
$H_ = false; |
|
64
|
|
|
$_H = false; |
|
65
|
|
|
|
|
66
|
|
|
for ($i = 0; $i < 24; $i++) { |
|
67
|
|
|
$a = $wheel[$i]; |
|
68
|
|
|
$b = $wheel[$i + 1]; |
|
69
|
|
|
|
|
70
|
|
|
if ($b && $b[1] < $a[1]) { |
|
|
|
|
|
|
71
|
|
|
$b[1] += 360; |
|
72
|
|
|
} |
|
73
|
|
|
|
|
74
|
|
View Code Duplication |
if (!$H_ && $a[0] <= $H && $b[0] > $H) { |
|
|
|
|
|
|
75
|
|
|
$this->H_RYB[$H] = (($a[1] + ($b[1] - $a[1]) * ($H - $a[0]) / ($b[0] - $a[0])) % 360); |
|
76
|
|
|
$H_ = true; |
|
77
|
|
|
} |
|
78
|
|
|
|
|
79
|
|
View Code Duplication |
if (!$_H && $a[1] <= $H && $b[1] > $H) { |
|
|
|
|
|
|
80
|
|
|
$this->RYB_H[$H] = (($a[0] + ($b[0] - $a[0]) * ($H - $a[1]) / ($b[1] - $a[1])) % 360); |
|
81
|
|
|
$_H = true; |
|
82
|
|
|
} |
|
83
|
|
|
|
|
84
|
|
|
if ($H_ == true && _H == true) { |
|
|
|
|
|
|
85
|
|
|
break; |
|
86
|
|
|
} |
|
87
|
|
|
} |
|
88
|
|
|
} |
|
89
|
|
|
$this->White = $this->illuminant('2', 'D65'); |
|
90
|
|
|
$this->profile('sRGB'); |
|
91
|
|
|
// echo "this->White = "; |
|
92
|
|
|
// print_r($this->White); |
|
93
|
|
|
} |
|
94
|
|
|
|
|
95
|
|
View Code Duplication |
public function RYB_HSV($o) |
|
|
|
|
|
|
96
|
|
|
{ |
|
97
|
|
|
$n = floor($o['H']); |
|
98
|
|
|
$x = $n > 0 ? $o['H'] % $n : 0; |
|
99
|
|
|
$a = $this->RYB_H[$n % 360]; |
|
100
|
|
|
$b = $this->RYB_H[ceil($o['H']) % 360]; |
|
101
|
|
|
|
|
102
|
|
|
return [ |
|
103
|
|
|
'H' => $a + ($b - $a) * $x, |
|
104
|
|
|
'S' => $o['S'], |
|
105
|
|
|
'V' => $o['V'], |
|
106
|
|
|
]; |
|
107
|
|
|
} |
|
108
|
|
|
|
|
109
|
|
View Code Duplication |
public function HSV_RYB($o) |
|
|
|
|
|
|
110
|
|
|
{ |
|
111
|
|
|
$n = floor($o['H']); |
|
112
|
|
|
$x = $n > 0 ? $o['H'] % $n : 0; |
|
113
|
|
|
$a = $this->H_RYB[$n % 360]; |
|
114
|
|
|
$b = $this->H_RYB[ceil($o['H']) % 360]; |
|
115
|
|
|
|
|
116
|
|
|
return [ |
|
117
|
|
|
'H' => $a + ($b - $a) * $x, |
|
118
|
|
|
'S' => $o['S'], |
|
119
|
|
|
'V' => $o['V'], |
|
120
|
|
|
]; |
|
121
|
|
|
} |
|
122
|
|
|
|
|
123
|
|
|
public function STRING_HEX($o) |
|
124
|
|
|
{ |
|
125
|
|
|
return 0 + ('0x'.$o); |
|
126
|
|
|
} |
|
127
|
|
|
|
|
128
|
|
|
public function HEX_STRING($o) |
|
129
|
|
|
{ |
|
130
|
|
|
$str = sprintf('%X', $o); |
|
131
|
|
|
$n = strlen($str); |
|
132
|
|
|
while ($n < 6) { |
|
133
|
|
|
$str = '0'.$str; |
|
134
|
|
|
$n++; |
|
135
|
|
|
} |
|
136
|
|
|
|
|
137
|
|
|
return $str; |
|
138
|
|
|
} |
|
139
|
|
|
|
|
140
|
|
|
public function HEX_RGB($o) |
|
141
|
|
|
{ |
|
142
|
|
|
return [ |
|
143
|
|
|
'R' => ($o >> 16), |
|
144
|
|
|
'G' => ($o >> 8) & 0xFF, |
|
145
|
|
|
'B' => $o & 0xFF, |
|
146
|
|
|
]; |
|
147
|
|
|
} |
|
148
|
|
|
|
|
149
|
|
|
public function HEX32_RGBA($o) |
|
150
|
|
|
{ |
|
151
|
|
|
return [ |
|
152
|
|
|
'R' => $o >> 16 & 0xFF, |
|
153
|
|
|
'G' => $o >> 9 & 0xFF, |
|
154
|
|
|
'B' => $o & 0xFF, |
|
155
|
|
|
'A' => $o >> 24, |
|
156
|
|
|
]; |
|
157
|
|
|
} |
|
158
|
|
|
|
|
159
|
|
|
public function RGBA_HEX32($o) |
|
160
|
|
|
{ |
|
161
|
|
|
return ($o['A'] << 24 | $o['R'] << 16 | $o['G'] << 8 | $o['B']) >> 0; |
|
162
|
|
|
} |
|
163
|
|
|
|
|
164
|
|
|
public function HEX32_rgbaa($o) |
|
165
|
|
|
{ |
|
166
|
|
|
return 'rgba('.($o >> 16 & 0xFF).','. |
|
167
|
|
|
($o >> 8 & 0xFF).','.($o & 0xFF).','. |
|
168
|
|
|
(($o >> 24) / 255).')'; |
|
169
|
|
|
} |
|
170
|
|
|
|
|
171
|
|
|
/* |
|
172
|
|
|
function RGB_STRING($o) |
|
173
|
|
|
{ |
|
174
|
|
|
return $this->HEX_STRING($this->RGB_HEX($o)); |
|
175
|
|
|
} |
|
176
|
|
|
function RGB_rgbaa($o) |
|
177
|
|
|
{ |
|
178
|
|
|
return 'rgba('. $o['R'] .','.$o['G'].','.$o['B'].','.$o['A'].')'; |
|
179
|
|
|
} |
|
180
|
|
|
function HEX_rgba($o) |
|
181
|
|
|
{ |
|
182
|
|
|
return $this->HEX32_rgbaa($o); |
|
183
|
|
|
} |
|
184
|
|
|
*/ |
|
185
|
|
|
public function RGB_HEX($o) |
|
186
|
|
|
{ |
|
187
|
|
|
return $o['R'] << 16 | $o['G'] << 8 | $o['B']; |
|
188
|
|
|
} |
|
189
|
|
|
|
|
190
|
|
|
public function RGB_CMY($o) |
|
191
|
|
|
{ |
|
192
|
|
|
return [ |
|
193
|
|
|
'C' => 1 - ($o['R'] / 255), |
|
194
|
|
|
'M' => 1 - ($o['G'] / 255), |
|
195
|
|
|
'Y' => 1 - ($o['B'] / 255), |
|
196
|
|
|
]; |
|
197
|
|
|
} |
|
198
|
|
|
|
|
199
|
|
|
public function RGB_HSL($o) |
|
200
|
|
|
{ |
|
201
|
|
|
$_R = $o['R'] / 255; |
|
202
|
|
|
$_G = $o['G'] / 255; |
|
203
|
|
|
$_B = $o['B'] / 255; |
|
204
|
|
|
$min = min($_R, $_G, $_B); |
|
205
|
|
|
$max = max($_R, $_G, $_B); |
|
206
|
|
|
$D = $max - $min; |
|
207
|
|
|
$L = ($max + $min) / 2; |
|
208
|
|
|
|
|
209
|
|
|
if ($D == 0) { |
|
210
|
|
|
$H = 0; |
|
211
|
|
|
$S = 0; |
|
212
|
|
|
} // No Chroma |
|
213
|
|
|
|
|
214
|
|
|
else { |
|
215
|
|
View Code Duplication |
if ($L < 0.5) { |
|
|
|
|
|
|
216
|
|
|
$S = $D / ($max + $min); |
|
217
|
|
|
} else { |
|
218
|
|
|
$S = $D / (2 - $max - $min); |
|
219
|
|
|
} |
|
220
|
|
|
|
|
221
|
|
|
$DR = ((($max - $_R) / 6) + ($D / 2)) / $D; |
|
222
|
|
|
$DG = ((($max - $_G) / 6) + ($D / 2)) / $D; |
|
223
|
|
|
$DB = ((($max - $_B) / 6) + ($D / 2)) / $D; |
|
224
|
|
|
|
|
225
|
|
View Code Duplication |
if ($_R == $max) { |
|
|
|
|
|
|
226
|
|
|
$H = $DB - $DG; |
|
227
|
|
|
} else { |
|
228
|
|
|
if ($_G == $max) { |
|
229
|
|
|
$H = (1 / 3) + $DR - $DB; |
|
230
|
|
|
} else { |
|
231
|
|
|
if ($_B == $max) { |
|
232
|
|
|
$H = (2 / 3) + $DG - $DR; |
|
233
|
|
|
} |
|
234
|
|
|
} |
|
235
|
|
|
} |
|
236
|
|
|
|
|
237
|
|
|
if ($H < 0) { |
|
238
|
|
|
$H += 1; |
|
|
|
|
|
|
239
|
|
|
} |
|
240
|
|
|
if ($H > 1) { |
|
241
|
|
|
$H -= 1; |
|
242
|
|
|
} |
|
243
|
|
|
} |
|
244
|
|
|
|
|
245
|
|
|
return [ |
|
246
|
|
|
'H' => $H * 360, |
|
247
|
|
|
'S' => $S * 100, |
|
248
|
|
|
'L' => $L * 100, |
|
249
|
|
|
]; |
|
250
|
|
|
} |
|
251
|
|
|
|
|
252
|
|
|
public function RGB_HSV($o) |
|
253
|
|
|
{ |
|
254
|
|
|
$_R = $o['R'] / 255; |
|
255
|
|
|
$_G = $o['G'] / 255; |
|
256
|
|
|
$_B = $o['B'] / 255; |
|
257
|
|
|
|
|
258
|
|
|
$min = min($_R, $_G, $_B); |
|
259
|
|
|
$max = max($_R, $_G, $_B); |
|
260
|
|
|
$D = $max - $min; |
|
261
|
|
|
$V = $max; |
|
262
|
|
|
|
|
263
|
|
|
if ($D == 0) { |
|
264
|
|
|
$H = 0; |
|
265
|
|
|
$S = 0; |
|
266
|
|
|
} // No chroma |
|
267
|
|
|
|
|
268
|
|
|
else { // Chromatic data |
|
269
|
|
|
|
|
270
|
|
|
$S = $D / $max; |
|
271
|
|
|
|
|
272
|
|
|
$DR = ((($max - $_R) / 6) + ($D / 2)) / $D; |
|
273
|
|
|
$DG = ((($max - $_G) / 6) + ($D / 2)) / $D; |
|
274
|
|
|
$DB = ((($max - $_B) / 6) + ($D / 2)) / $D; |
|
275
|
|
|
|
|
276
|
|
View Code Duplication |
if ($_R == $max) { |
|
|
|
|
|
|
277
|
|
|
$H = $DB - $DG; |
|
278
|
|
|
} else { |
|
279
|
|
|
if ($_G == $max) { |
|
280
|
|
|
$H = (1 / 3) + $DR - $DB; |
|
281
|
|
|
} else { |
|
282
|
|
|
if ($_B == $max) { |
|
283
|
|
|
$H = (2 / 3) + $DG - $DR; |
|
284
|
|
|
} |
|
285
|
|
|
} |
|
286
|
|
|
} |
|
287
|
|
|
|
|
288
|
|
|
if ($H < 0) { |
|
289
|
|
|
$H += 1; |
|
|
|
|
|
|
290
|
|
|
} |
|
291
|
|
|
if ($H > 1) { |
|
292
|
|
|
$H -= 1; |
|
293
|
|
|
} |
|
294
|
|
|
} |
|
295
|
|
|
|
|
296
|
|
|
return ['H' => $H * 360, 'S' => $S * 100, 'V' => $V * 100]; |
|
297
|
|
|
} |
|
298
|
|
|
|
|
299
|
|
|
public function RGB_XYZ($o) |
|
300
|
|
|
{ |
|
301
|
|
|
$M = $this->_RGB_XYZ; |
|
302
|
|
|
$z = []; |
|
303
|
|
|
|
|
304
|
|
|
$R = $o['R'] / 255; |
|
305
|
|
|
$G = $o['G'] / 255; |
|
306
|
|
|
$B = $o['B'] / 255; |
|
307
|
|
|
|
|
308
|
|
|
if ($this->_Space == 'sRGB') { |
|
309
|
|
|
$R = ($R > 0.04045) ? pow((($R + 0.055) / 1.055), 2.4) : $R / 12.92; |
|
310
|
|
|
$G = ($G > 0.04045) ? pow((($G + 0.055) / 1.055), 2.4) : $G / 12.92; |
|
311
|
|
|
$B = ($B > 0.04045) ? pow((($B + 0.055) / 1.055), 2.4) : $B / 12.92; |
|
312
|
|
|
} else { |
|
313
|
|
|
$R = pow($R, $this->_Gamma); |
|
314
|
|
|
$G = pow($G, $this->_Gamma); |
|
315
|
|
|
$B = pow($B, $this->_Gamma); |
|
316
|
|
|
} |
|
317
|
|
|
|
|
318
|
|
|
$z['X'] = $R * $M[0] + $G * $M[3] + $B * $M[6]; |
|
319
|
|
|
$z['Y'] = $R * $M[1] + $G * $M[4] + $B * $M[7]; |
|
320
|
|
|
$z['Z'] = $R * $M[2] + $G * $M[5] + $B * $M[8]; |
|
321
|
|
|
|
|
322
|
|
|
return $z; |
|
323
|
|
|
} |
|
324
|
|
|
|
|
325
|
|
|
public function CMY_RGB($o) |
|
326
|
|
|
{ |
|
327
|
|
|
return [ |
|
328
|
|
|
'R' => max(0, (1 - $o['C']) * 255), |
|
329
|
|
|
'G' => max(0, (1 - $o['M']) * 255), |
|
330
|
|
|
'B' => max(0, (1 - $o['Y']) * 255), |
|
331
|
|
|
]; |
|
332
|
|
|
} |
|
333
|
|
|
|
|
334
|
|
|
public function CMY_CMYK($o) |
|
335
|
|
|
{ |
|
336
|
|
|
$C = $o['C']; |
|
337
|
|
|
$M = $o['M']; |
|
338
|
|
|
$Y = $o['Y']; |
|
339
|
|
|
$K = min($Y, min($M, min($C, 1))); |
|
340
|
|
|
|
|
341
|
|
|
$C = round(($C - $K) / (1 - $K) * 100); |
|
342
|
|
|
$M = round(($M - $K) / (1 - $K) * 100); |
|
343
|
|
|
$Y = round(($Y - $K) / (1 - $K) * 100); |
|
344
|
|
|
$K = round($K * 100); |
|
345
|
|
|
|
|
346
|
|
|
return ['C' => $C, 'M' => $M, 'Y' => $Y, 'K' => $K]; |
|
347
|
|
|
} |
|
348
|
|
|
|
|
349
|
|
|
// CMYK = C: Cyan / M: Magenta / Y: Yellow / K: Key (black) |
|
350
|
|
|
|
|
351
|
|
|
public function CMYK_CMY($o) |
|
352
|
|
|
{ |
|
353
|
|
|
return [ |
|
354
|
|
|
'C' => ($o['C'] * (1 - $o['K']) + $o['K']), |
|
355
|
|
|
'M' => ($o['M'] * (1 - $o['K']) + $o['K']), |
|
356
|
|
|
'Y' => ($o['Y'] * (1 - $o['K']) + $o['K']), |
|
357
|
|
|
]; |
|
358
|
|
|
} |
|
359
|
|
|
|
|
360
|
|
|
// HSL (1978) = H: Hue / S: Saturation / L: Lightess |
|
361
|
|
|
|
|
362
|
|
View Code Duplication |
public function Hue_2_RGB($v1, $v2, $vH) |
|
|
|
|
|
|
363
|
|
|
{ |
|
364
|
|
|
if ($vH < 0) { |
|
365
|
|
|
$vH += 1; |
|
366
|
|
|
} |
|
367
|
|
|
if ($vH > 1) { |
|
368
|
|
|
$vH -= 1; |
|
369
|
|
|
} |
|
370
|
|
|
if ((6 * $vH) < 1) { |
|
371
|
|
|
return $v1 + ($v2 - $v1) * 6 * $vH; |
|
372
|
|
|
} |
|
373
|
|
|
if ((2 * $vH) < 1) { |
|
374
|
|
|
return $v2; |
|
375
|
|
|
} |
|
376
|
|
|
if ((3 * $vH) < 2) { |
|
377
|
|
|
return $v1 + ($v2 - $v1) * ((2 / 3) - $vH) * 6; |
|
378
|
|
|
} |
|
379
|
|
|
|
|
380
|
|
|
return $v1; |
|
381
|
|
|
} |
|
382
|
|
|
|
|
383
|
|
|
public function HSL_RGB($o) |
|
384
|
|
|
{ |
|
385
|
|
|
$H = $o['H'] / 360; |
|
386
|
|
|
$S = $o['S'] / 100; |
|
387
|
|
|
$L = $o['L'] / 100; |
|
388
|
|
|
$R; |
|
|
|
|
|
|
389
|
|
|
$G; |
|
|
|
|
|
|
390
|
|
|
$B; |
|
|
|
|
|
|
391
|
|
|
$_1; |
|
|
|
|
|
|
392
|
|
|
$_2; |
|
|
|
|
|
|
393
|
|
|
|
|
394
|
|
|
if ($S == 0) { // HSL from 0 to 1 |
|
395
|
|
|
|
|
396
|
|
|
$R = $L * 255; |
|
397
|
|
|
$G = $L * 255; |
|
398
|
|
|
$B = $L * 255; |
|
399
|
|
|
} else { |
|
400
|
|
|
if ($L < 0.5) { |
|
401
|
|
|
$_2 = $L * (1 + $S); |
|
402
|
|
|
} else { |
|
403
|
|
|
$_2 = ($L + $S) - ($S * $L); |
|
404
|
|
|
} |
|
405
|
|
|
|
|
406
|
|
|
$_1 = 2 * $L - $_2; |
|
407
|
|
|
|
|
408
|
|
|
$R = 255 * $this->Hue_2_RGB($_1, $_2, $H + (1 / 3)); |
|
409
|
|
|
$G = 255 * $this->Hue_2_RGB($_1, $_2, $H); |
|
410
|
|
|
$B = 255 * $this->Hue_2_RGB($_1, $_2, $H - (1 / 3)); |
|
411
|
|
|
} |
|
412
|
|
|
|
|
413
|
|
|
return ['R' => $R, 'G' => $G, 'B' => $B]; |
|
414
|
|
|
} |
|
415
|
|
|
|
|
416
|
|
|
// HSV (1978) = H: Hue / S: Saturation / V: Value |
|
417
|
|
|
|
|
418
|
|
|
public function HSV_RGB($o) |
|
419
|
|
|
{ |
|
420
|
|
|
$H = $o['H'] / 360; |
|
421
|
|
|
$S = $o['S'] / 100; |
|
422
|
|
|
$V = $o['V'] / 100; |
|
423
|
|
|
$R; |
|
|
|
|
|
|
424
|
|
|
$G; |
|
|
|
|
|
|
425
|
|
|
$B; |
|
|
|
|
|
|
426
|
|
|
|
|
427
|
|
|
if ($S == 0) { |
|
428
|
|
|
$R = $G = $B = round($V * 255); |
|
429
|
|
|
} else { |
|
430
|
|
|
if ($H >= 1) { |
|
431
|
|
|
$H = 0; |
|
432
|
|
|
} |
|
433
|
|
|
|
|
434
|
|
|
$H = 6 * $H; |
|
435
|
|
|
$D = $H - floor($H); |
|
436
|
|
|
$A = round(255 * $V * (1 - $S)); |
|
437
|
|
|
$B = round(255 * $V * (1 - ($S * $D))); |
|
438
|
|
|
$C = round(255 * $V * (1 - ($S * (1 - $D)))); |
|
439
|
|
|
$V = round(255 * $V); |
|
440
|
|
|
|
|
441
|
|
|
switch (floor($H)) { |
|
442
|
|
|
|
|
443
|
|
|
case 0: |
|
444
|
|
|
$R = $V; |
|
445
|
|
|
$G = $C; |
|
446
|
|
|
$B = $A; |
|
447
|
|
|
break; |
|
448
|
|
|
case 1: |
|
449
|
|
|
$R = $B; |
|
450
|
|
|
$G = $V; |
|
451
|
|
|
$B = $A; |
|
452
|
|
|
break; |
|
453
|
|
|
case 2: |
|
454
|
|
|
$R = $A; |
|
455
|
|
|
$G = $V; |
|
456
|
|
|
$B = $C; |
|
457
|
|
|
break; |
|
458
|
|
|
case 3: |
|
459
|
|
|
$R = $A; |
|
460
|
|
|
$G = $B; |
|
461
|
|
|
$B = $V; |
|
462
|
|
|
break; |
|
463
|
|
|
case 4: |
|
464
|
|
|
$R = $C; |
|
465
|
|
|
$G = $A; |
|
466
|
|
|
$B = $V; |
|
467
|
|
|
break; |
|
468
|
|
|
case 5: |
|
469
|
|
|
$R = $V; |
|
470
|
|
|
$G = $A; |
|
471
|
|
|
$B = $B; |
|
|
|
|
|
|
472
|
|
|
break; |
|
473
|
|
|
} |
|
474
|
|
|
} |
|
475
|
|
|
|
|
476
|
|
|
return ['R' => $R, 'G' => $G, 'B' => $B]; |
|
|
|
|
|
|
477
|
|
|
} |
|
478
|
|
|
|
|
479
|
|
|
// CIE (Commission International de L’Eclairage) |
|
480
|
|
|
|
|
481
|
|
|
// CIE-XYZ (1931) = Y: Luminescence / XZ: Spectral Weighting Curves (Spectral Locus) |
|
482
|
|
|
|
|
483
|
|
|
public function XYZ_RGB($o) |
|
484
|
|
|
{ |
|
485
|
|
|
$M = $this->_XYZ_RGB; |
|
486
|
|
|
$z = []; |
|
487
|
|
|
|
|
488
|
|
|
$z['R'] = $o['X'] * $M[0] + $o['Y'] * $M[3] + $o['Z'] * $M[6]; |
|
489
|
|
|
$z['G'] = $o['X'] * $M[1] + $o['Y'] * $M[4] + $o['Z'] * $M[7]; |
|
490
|
|
|
$z['B'] = $o['X'] * $M[2] + $o['Y'] * $M[5] + $o['Z'] * $M[8]; |
|
491
|
|
|
|
|
492
|
|
|
if ($this->_Space == 'sRGB') { |
|
493
|
|
|
$z['R'] = ($z['R'] > 0.0031308) ? (1.055 * pow($z['R'], 1 / 2.4)) - 0.055 : 12.92 * $z['R']; |
|
494
|
|
|
$z['G'] = ($z['G'] > 0.0031308) ? (1.055 * pow($z['G'], 1 / 2.4)) - 0.055 : 12.92 * $z['G']; |
|
495
|
|
|
$z['B'] = ($z['B'] > 0.0031308) ? (1.055 * pow($z['B'], 1 / 2.4)) - 0.055 : 12.92 * $z['B']; |
|
496
|
|
|
} else { |
|
497
|
|
|
$z['R'] = pow($z['R'], 1 / $this->_Gamma); |
|
498
|
|
|
$z['G'] = pow($z['G'], 1 / $this->_Gamma); |
|
499
|
|
|
$z['B'] = pow($z['B'], 1 / $this->_Gamma); |
|
500
|
|
|
} |
|
501
|
|
|
|
|
502
|
|
|
return ['R' => round($z['R'] * 255), 'G' => round($z['G'] * 255), 'B' => round($z['B'] * 255)]; |
|
503
|
|
|
} |
|
504
|
|
|
|
|
505
|
|
|
public function XYZ_xyY($o) |
|
506
|
|
|
{ |
|
507
|
|
|
$n = $o['X'] + $o['Y'] + $o['Z']; |
|
508
|
|
|
|
|
509
|
|
|
if ($n == 0) { |
|
510
|
|
|
return ['x' => 0, 'y' => 0, 'Y' => $o['Y']]; |
|
511
|
|
|
} |
|
512
|
|
|
|
|
513
|
|
|
return ['x' => $o['X'] / $n, 'y' => $o['Y'] / $n, 'Y' => $o['Y']]; |
|
514
|
|
|
} |
|
515
|
|
|
|
|
516
|
|
|
public function XYZ_HLab($o) |
|
517
|
|
|
{ |
|
518
|
|
|
$n = sqrt($o['Y']); |
|
519
|
|
|
|
|
520
|
|
|
return [ |
|
521
|
|
|
'L' => 10 * $n, |
|
522
|
|
|
'a' => 17.5 * (((1.02 * $o['X']) - $o['Y']) / $n), |
|
523
|
|
|
'b' => 7 * (($o['Y'] - (0.847 * $o['Z'])) / $n), |
|
524
|
|
|
]; |
|
525
|
|
|
} |
|
526
|
|
|
|
|
527
|
|
|
public function XYZ_Lab($o) |
|
528
|
|
|
{ |
|
529
|
|
|
$r = $this->White; |
|
530
|
|
|
|
|
531
|
|
|
function fu($n) |
|
|
|
|
|
|
532
|
|
|
{ |
|
533
|
|
|
if ($n > 0.008856) { |
|
534
|
|
|
return pow($n, 1 / 3); |
|
535
|
|
|
} else { |
|
536
|
|
|
return (7.787 * $n) + (16 / 116); |
|
537
|
|
|
} |
|
538
|
|
|
} |
|
539
|
|
|
|
|
540
|
|
|
$X = fu($o['X'] / $r['X']); |
|
541
|
|
|
$Y = fu($o['Y'] / $r['Y']); |
|
542
|
|
|
$Z = fu($o['Z'] / $r['Z']); |
|
543
|
|
|
|
|
544
|
|
|
return ['L' => (116 * $Y) - 16, 'a' => 500 * ($X - $Y), 'b' => 200 * ($Y - $Z)]; |
|
545
|
|
|
} |
|
546
|
|
|
|
|
547
|
|
|
public function XYZ_Luv($o) |
|
548
|
|
|
{ |
|
549
|
|
|
$r = $this->White; |
|
550
|
|
|
|
|
551
|
|
|
$U = (4 * $o['X']) / ($o['X'] + (15 * $o['Y']) + (3 * $o['Z'])); |
|
552
|
|
|
$V = (9 * $o['Y']) / ($o['X'] + (15 * $o['Y']) + (3 * $o['Z'])); |
|
553
|
|
|
|
|
554
|
|
|
if ($o['Y'] > 0.008856) { |
|
555
|
|
|
$o['Y'] = pow($o['Y'], 1 / 3); |
|
556
|
|
|
} else { |
|
557
|
|
|
$o['Y'] = (7.787 * $o['Y']) + (16 / 116); |
|
558
|
|
|
} |
|
559
|
|
|
|
|
560
|
|
|
$_L = (116 * $o['Y']) - 16; |
|
561
|
|
|
$_U = (4 * $r['X']) / ($r['X'] + (15 * $r['Y']) + (3 * $r['Z'])); |
|
562
|
|
|
$_V = (9 * $r['Y']) / ($r['X'] + (15 * $r['Y']) + (3 * $r['Z'])); |
|
563
|
|
|
|
|
564
|
|
|
return ['L' => $_L, 'u' => 13 * $_L * ($U - $_U), 'v' => 13 * $_L * ($V - $_V)]; |
|
565
|
|
|
} |
|
566
|
|
|
|
|
567
|
|
|
// CIE-xyY (1931) = Y: Luminescence / xy: Chromaticity Co-ordinates (Spectral Locus) |
|
568
|
|
|
|
|
569
|
|
|
public function xyY_XYZ($o) |
|
570
|
|
|
{ |
|
571
|
|
|
return [ |
|
572
|
|
|
'X' => ($o['x'] * $o['Y']) / $o['y'], |
|
573
|
|
|
'Y' => $o['Y'], |
|
574
|
|
|
'Z' => ((1 - $o['x'] - $o['y']) * $o['Y']) / $o['y'], |
|
575
|
|
|
]; |
|
576
|
|
|
} |
|
577
|
|
|
|
|
578
|
|
|
// Hunter-L*ab (1948) = L: Lightness / ab: Color-opponent Dimensions |
|
579
|
|
|
|
|
580
|
|
|
public function HLab_XYZ($o) |
|
581
|
|
|
{ |
|
582
|
|
|
$_Y = $o['L'] / 10; |
|
583
|
|
|
$_X = ($o['a'] / 17.5) * ($o['L'] / 10); |
|
584
|
|
|
$_Z = ($o['b'] / 7) * ($o['L'] / 10); |
|
585
|
|
|
|
|
586
|
|
|
$Y = pow($_Y, 2); |
|
587
|
|
|
$X = ($_X + $Y) / 1.02; |
|
588
|
|
|
$Z = -1 * ($_Z - $Y) / 0.847; |
|
589
|
|
|
|
|
590
|
|
|
return ['X' => $X, 'Y' => $Y, 'Z' => $Z]; |
|
591
|
|
|
} |
|
592
|
|
|
|
|
593
|
|
|
// CIE-L*ab (1976) = L: Luminescence / a: Red / Green / b: Blue / Yellow |
|
594
|
|
|
|
|
595
|
|
|
public function Lab_XYZ($o) |
|
596
|
|
|
{ |
|
597
|
|
|
$r = $this->White; |
|
598
|
|
|
|
|
599
|
|
|
$Y = ($o['L'] + 16) / 116; |
|
600
|
|
|
$X = $o['a'] / 500 + $Y; |
|
601
|
|
|
$Z = $Y - $o['b'] / 200; |
|
602
|
|
|
|
|
603
|
|
|
$Y = pow($Y, 3) > 0.008856 ? pow($Y, 3) : ($Y - 16 / 116) / 7.787; |
|
604
|
|
|
$X = pow($X, 3) > 0.008856 ? pow($X, 3) : ($X - 16 / 116) / 7.787; |
|
605
|
|
|
$Z = pow($Z, 3) > 0.008856 ? pow($Z, 3) : ($Z - 16 / 116) / 7.787; |
|
606
|
|
|
|
|
607
|
|
|
return ['X' => $r['X'] * $X, 'Y' => $r['Y'] * $Y, $Z => $r['Z'] * $Z]; |
|
608
|
|
|
} |
|
609
|
|
|
|
|
610
|
|
View Code Duplication |
public function Lab_LCHab($o) |
|
|
|
|
|
|
611
|
|
|
{ |
|
612
|
|
|
$H = atan2($o['b'], $o['a']) * (180 / PI); |
|
613
|
|
|
|
|
614
|
|
|
if ($H < 0) { |
|
615
|
|
|
$H += 360; |
|
616
|
|
|
} else { |
|
617
|
|
|
if ($H > 360) { |
|
618
|
|
|
$H -= 360; |
|
619
|
|
|
} |
|
620
|
|
|
} |
|
621
|
|
|
|
|
622
|
|
|
return ['L' => $o['L'], 'C' => sqrt($o['a'] * $o['a'] + $o['b'] * $o['b']), 'H' => $H]; |
|
623
|
|
|
} |
|
624
|
|
|
|
|
625
|
|
|
// CIE-L*uv (1976) = L: Luminescence / u: Saturation / v: Hue |
|
626
|
|
|
|
|
627
|
|
|
public function Luv_XYZ($o) |
|
628
|
|
|
{ |
|
629
|
|
|
$r = $White; |
|
|
|
|
|
|
630
|
|
|
|
|
631
|
|
|
$Y = ($o['L'] + 16) / 116; |
|
632
|
|
|
$Y = (pow($Y, 3) > 0.008856) ? pow($Y, 3) : (($Y - 16 / 116) / 7.787); |
|
633
|
|
|
|
|
634
|
|
|
$_U = (4 * $r['X']) / ($r['X'] + (15 * $r['Y']) + (3 * $r['Z'])); |
|
635
|
|
|
$_V = (9 * $r['Y']) / ($r['X'] + (15 * $r['Y']) + (3 * $r['Z'])); |
|
636
|
|
|
|
|
637
|
|
|
$U = $o['u'] / (13 * $o['L']) + $_U; |
|
638
|
|
|
$V = $o['v'] / (13 * $o['L']) + $_V; |
|
639
|
|
|
|
|
640
|
|
|
$X = -(9 * $Y * $U) / (($U - 4) * $V - $U * $V); |
|
641
|
|
|
$Z = (9 * $Y - (15 * $V * $Y) - ($V * $X)) / (3 * $V); |
|
642
|
|
|
|
|
643
|
|
|
return ['X' => $X, 'Y' => $Y, 'Z' => $Z]; |
|
644
|
|
|
} |
|
645
|
|
|
|
|
646
|
|
View Code Duplication |
public function Luv_LCHuv($o) |
|
|
|
|
|
|
647
|
|
|
{ |
|
648
|
|
|
$H = atan2($o['v'], $o['u']) * (180 / PI); |
|
649
|
|
|
|
|
650
|
|
|
if ($H < 0) { |
|
651
|
|
|
$H += 360; |
|
652
|
|
|
} else { |
|
653
|
|
|
if ($H > 360) { |
|
654
|
|
|
$H -= 360; |
|
655
|
|
|
} |
|
656
|
|
|
} |
|
657
|
|
|
|
|
658
|
|
|
return ['L' => $o['L'], 'C' => sqrt($o['u'] * $o['u'] + $o['v'] * $o['v']), 'H' => $H]; |
|
659
|
|
|
} |
|
660
|
|
|
|
|
661
|
|
|
// CIE-L*CH (1986) = L: Luminescece / C: Chromacity / H: Hue |
|
662
|
|
|
|
|
663
|
|
|
public function LCHab_Lab($o) |
|
664
|
|
|
{ |
|
665
|
|
|
$rad = $o['H'] * (PI / 180); |
|
666
|
|
|
|
|
667
|
|
|
return ['L' => $o['L'], 'a' => cos($rad) * $o['C'], 'b' => sin($rad) * $o['C']]; |
|
668
|
|
|
} |
|
669
|
|
|
|
|
670
|
|
|
public function LCHuv_Luv($o) |
|
671
|
|
|
{ |
|
672
|
|
|
$rad = $o['H'] * (PI / 180); |
|
673
|
|
|
|
|
674
|
|
|
return ['L' => $o['L'], 'u' => cos($rad) * $o['C'], 'v' => sin($rad) * $o['C']]; |
|
675
|
|
|
} |
|
676
|
|
|
|
|
677
|
|
|
public function adapt($o, $type) |
|
678
|
|
|
{ |
|
679
|
|
|
$r = [ // Adaption methods |
|
680
|
|
|
'XYZ scaling' => [ |
|
681
|
|
|
'A' => [[1, 0, 0], [0, 1, 0], [0, 0, 1]], |
|
682
|
|
|
'Z' => [[1, 0, 0], [0, 1, 0], [0, 0, 1]], |
|
683
|
|
|
], |
|
684
|
|
|
'Von Kries' => [ |
|
685
|
|
|
'A' => [[0.400240, -0.226300, 0], [0.707600, 1.165320, 0], [-0.080810, 0.045700, 0.918220]], |
|
686
|
|
|
'Z' => [[1.859936, 0.361191, 0], [-1.129382, 0.638812, 0], [0.219897, -0.000006, 1.089064]], |
|
687
|
|
|
], |
|
688
|
|
|
'Bradford' => [ |
|
689
|
|
|
'A' => [ |
|
690
|
|
|
[0.895100, 0.26640000, -0.16139900], |
|
691
|
|
|
[-0.75019900, 1.71350, 0.0367000], |
|
692
|
|
|
[0.03889900, -0.0685000, 1.02960000], |
|
693
|
|
|
], |
|
694
|
|
|
'Z' => [ |
|
695
|
|
|
[0.986993, -0.14705399, 0.15996299], |
|
696
|
|
|
[0.43230499, 0.51836, 0.0492912], |
|
697
|
|
|
[-0.00852866, 0.0400428, 0.96848699], |
|
698
|
|
|
], |
|
699
|
|
|
], |
|
700
|
|
|
]; |
|
701
|
|
|
|
|
702
|
|
|
$WS = $this->_xyY; |
|
703
|
|
|
$WD = $this->White; // White Point Source + Destination |
|
704
|
|
|
// echo 'WS = '; |
|
705
|
|
|
// print_r($WS); |
|
706
|
|
|
// echo 'WD = '; |
|
707
|
|
|
// print_r($WD); |
|
708
|
|
|
|
|
709
|
|
|
$A = $r[$type]['A']; |
|
710
|
|
|
$Z = $r[$type]['Z']; // Load Matrices |
|
711
|
|
|
|
|
712
|
|
|
$CRD = $this->multiply($A, [[$WD['X']], [$WD['Y']], [$WD['Z']]]); // Convert to cone responce domain |
|
713
|
|
|
$CRS = $this->multiply($A, [[$WS['X']], [$WS['Y']], [$WS['Z']]]); |
|
714
|
|
|
//print_r($CRD); |
|
715
|
|
|
//print_r($CRS); |
|
716
|
|
|
|
|
717
|
|
|
$M = [ |
|
718
|
|
|
[$CRD[0][0] / $CRS[0][0], 0, 0], |
|
719
|
|
|
[0, $CRD[1][0] / $CRS[1][0], 0], |
|
720
|
|
|
[0, 0, $CRD[2][0] / $CRS[2][0]], |
|
721
|
|
|
]; // Scale Vectors |
|
722
|
|
|
|
|
723
|
|
|
$z = $this->multiply($Z, |
|
724
|
|
|
$this->multiply($M, $this->multiply($A, [[$o['X']], [$o['Y']], [$o['Z']]]))); // Back to XYZ |
|
725
|
|
|
|
|
726
|
|
|
return ['X' => $z[0][0], 'Y' => $z[1][0], 'Z' => $z[2][0]]; |
|
727
|
|
|
} |
|
728
|
|
|
|
|
729
|
|
|
public function f($o) |
|
730
|
|
|
{ |
|
731
|
|
|
$x = $this->xyY_XYZ($o); |
|
732
|
|
|
// echo 'o = '; |
|
733
|
|
|
// print_r($o); |
|
734
|
|
|
// echo 'x = '; |
|
735
|
|
|
// print_r($x); |
|
736
|
|
|
return $this->adapt($x, 'Bradford'); |
|
737
|
|
|
} |
|
738
|
|
|
|
|
739
|
|
|
public function illuminant($observer, $type) |
|
740
|
|
|
{ |
|
741
|
|
|
$o = $this->_illuminant[$type]; |
|
742
|
|
|
|
|
743
|
|
|
$o = ($observer == 2) ? ['x' => $o[0], 'y' => $o[1], 'Y' => 1] : ['x' => $o[2], 'y' => $o[3], 'Y' => 1]; |
|
744
|
|
|
|
|
745
|
|
|
//print_r($o); |
|
746
|
|
|
|
|
747
|
|
|
return $this->xyY_XYZ($o); |
|
748
|
|
|
} |
|
749
|
|
|
|
|
750
|
|
|
private $_Space; |
|
751
|
|
|
private $_Gamma; |
|
752
|
|
|
private $_White; |
|
753
|
|
|
private $_Matrix; |
|
754
|
|
|
|
|
755
|
|
|
public function profile($i) |
|
756
|
|
|
{ |
|
757
|
|
|
$m = $this->_profile[$i]; |
|
758
|
|
|
//print_r($m); |
|
759
|
|
|
|
|
760
|
|
|
$this->_Space = $i; |
|
761
|
|
|
$this->_Gamma = $m[0]; |
|
762
|
|
|
$this->_White = $m[1]; |
|
763
|
|
|
$this->_Matrix = $m; |
|
764
|
|
|
|
|
765
|
|
|
// Input Illuminant |
|
766
|
|
|
|
|
767
|
|
|
$this->_xyY = $this->illuminant('2', $m[1]); |
|
768
|
|
|
//print_r($this->_xyY); |
|
769
|
|
|
|
|
770
|
|
|
$R = $this->f(['x' => $m[2], 'y' => $m[3], 'Y' => $m[4]]); |
|
771
|
|
|
$G = $this->f(['x' => $m[5], 'y' => $m[6], 'Y' => $m[7]]); |
|
772
|
|
|
$B = $this->f(['x' => $m[8], 'y' => $m[9], 'Y' => $m[10]]); |
|
773
|
|
|
// print_r($R); |
|
774
|
|
|
// print_r($G); |
|
775
|
|
|
// print_r($B); |
|
776
|
|
|
|
|
777
|
|
|
$this->_RGB_XYZ = [$R['X'], $R['Y'], $R['Z'], $G['X'], $G['Y'], $G['Z'], $B['X'], $B['Y'], $B['Z']]; |
|
778
|
|
|
// print_r($this->_RGB_XYZ); |
|
779
|
|
|
|
|
780
|
|
|
$this->_XYZ_RGB = $this->inverse($this->_RGB_XYZ); |
|
781
|
|
|
} |
|
782
|
|
|
|
|
783
|
|
|
private $_RGB_XYZ; |
|
784
|
|
|
|
|
785
|
|
|
private $_profile = [ // [ Gamma, Illuminant, Matrix ] |
|
786
|
|
|
'Adobe (1998)' => [2.2, 'D65', 0.64, 0.33, 0.297361, 0.21, 0.71, 0.627355, 0.15, 0.06, 0.075285], |
|
787
|
|
|
// Adobe |
|
788
|
|
|
'Apple RGB' => [1.8, 'D65', 0.625, 0.34, 0.244634, 0.28, 0.595, 0.672034, 0.155, 0.07, 0.083332], |
|
789
|
|
|
// Apple, a.k.a. SGI |
|
790
|
|
|
'BestRGB' => [2.2, 'D50', 0.7347, 0.2653, 0.228457, 0.215, 0.775, 0.737352, 0.13, 0.035, 0.034191], |
|
791
|
|
|
// Don Hutcheson |
|
792
|
|
|
'Beta RGB' => [2.2, 'D50', 0.6888, 0.3112, 0.303273, 0.1986, 0.7551, 0.663786, 0.1265, 0.0352, 0.032941], |
|
793
|
|
|
// Bruce Lindbloom |
|
794
|
|
|
'Bruce RGB' => [2.2, 'D65', 0.64, 0.33, 0.240995, 0.28, 0.65, 0.683554, 0.15, 0.06, 0.075452], |
|
795
|
|
|
// Bruce Fraser |
|
796
|
|
|
'CIE RGB' => [2.2, 'E', 0.735, 0.265, 0.176204, 0.274, 0.717, 0.812985, 0.167, 0.009, 0.010811], |
|
797
|
|
|
// CIE |
|
798
|
|
|
'ColorMatch' => [1.8, 'D50', 0.63, 0.34, 0.274884, 0.295, 0.605, 0.658132, 0.15, 0.075, 0.066985], |
|
799
|
|
|
// Radius |
|
800
|
|
|
'DonRGB4' => [2.2, 'D50', 0.696, 0.3, 0.27835, 0.215, 0.765, 0.68797, 0.13, 0.035, 0.03368], |
|
801
|
|
|
'eciRGB' => [1.8, 'D50', 0.67, 0.33, 0.32025, 0.21, 0.71, 0.602071, 0.14, 0.08, 0.077679], |
|
802
|
|
|
// European Colour Initiative |
|
803
|
|
|
'Ekta Space PS5' => [2.2, 'D50', 0.695, 0.305, 0.260629, 0.26, 0.7, 0.734946, 0.11, 0.005, 0.004425], |
|
804
|
|
|
// Joseph Holmes |
|
805
|
|
|
'Generic RGB' => [1.8, 'D65', 0.6295, 0.3407, 0.232546, 0.2949, 0.6055, 0.672501, 0.1551, 0.0762, 0.094952], |
|
806
|
|
|
'HDTV (HD-CIF)' => [1.95, 'D65', 0.64, 0.33, 0.212673, 0.3, 0.6, 0.715152, 0.15, 0.06, 0.072175], |
|
807
|
|
|
// a.k.a. ITU-R BT.701 |
|
808
|
|
|
'NTSC' => [2.2, 'C', 0.67, 0.33, 0.298839, 0.21, 0.71, 0.586811, 0.14, 0.08, 0.11435], |
|
809
|
|
|
// National Television System Committee (NTSC), a.k.a. Y'I'Q' |
|
810
|
|
|
'PAL / SECAM' => [2.2, 'D65', 0.64, 0.33, 0.222021, 0.29, 0.6, 0.706645, 0.15, 0.06, 0.071334], |
|
811
|
|
|
// European Broadcasting Union (EBU), a.k.a. Y'U'V' |
|
812
|
|
|
'ProPhoto' => [1.8, 'D50', 0.7347, 0.2653, 0.28804, 0.1596, 0.8404, 0.711874, 0.0366, 0.0001, 0.000086], |
|
813
|
|
|
// Kodak, a.k.a. ROMM RGB |
|
814
|
|
|
'SGI' => [1.47, 'D65', 0.625, 0.34, 0.244651, 0.28, 0.595, 0.672030, 0.155, 0.07, 0.083319], |
|
815
|
|
|
'SMPTE-240M' => [1.92, 'D65', 0.63, 0.34, 0.212413, 0.31, 0.595, 0.701044, 0.155, 0.07, 0.086543], |
|
816
|
|
|
'SMPTE-C' => [2.2, 'D65', 0.63, 0.34, 0.212395, 0.31, 0.595, 0.701049, 0.155, 0.07, 0.086556], |
|
817
|
|
|
// Society of Motion Picture and Television Engineers (SMPTE) |
|
818
|
|
|
'sRGB' => [2.2, 'D65', 0.64, 0.33, 0.212656, 0.3, 0.6, 0.715158, 0.15, 0.06, 0.072186], |
|
819
|
|
|
// Microsoft & Hewlett - Packard |
|
820
|
|
|
'Wide Gamut' => [2.2, 'D50', 0.7347, 0.2653, 0.258187, 0.1152, 0.8264, 0.724938, 0.1566, 0.0177, 0.016875], |
|
821
|
|
|
// Adobe |
|
822
|
|
|
]; |
|
823
|
|
|
|
|
824
|
|
|
private $_illuminant = [ // [ x2°, y2°, x10°, y10°, CCT (Kelvin) ] |
|
825
|
|
|
|
|
826
|
|
|
'A' => [0.44757, 0.40745, 0.45117, 0.40594, 2856], // Incandescent tungsten |
|
827
|
|
|
'B' => [0.34842, 0.35161, 0.3498, 0.3527, 4874], // Obsolete, direct sunlight at noon |
|
828
|
|
|
'C' => [0.31006, 0.31616, 0.31039, 0.31905, 6774], // Obsolete, north sky daylight |
|
829
|
|
|
'D50' => [0.34567, 0.35850, 0.34773, 0.35952, 5003], // ICC Profile PCS. Horizon light. |
|
830
|
|
|
'D55' => [0.33242, 0.34743, 0.33411, 0.34877, 5503], // Compromise between incandescent and daylight |
|
831
|
|
|
'D65' => [0.31271, 0.32902, 0.31382, 0.33100, 6504], // Daylight, sRGB color space |
|
832
|
|
|
'D75' => [0.29902, 0.31485, 0.29968, 0.31740, 7504], // North sky day light |
|
833
|
|
|
'E' => [0.33333, 0.33333, 0.33333, 0.33333, 5454], // Equal energy |
|
834
|
|
|
'F1' => [0.31310, 0.33727, 0.31811, 0.33559, 6430], // Daylight Fluorescent |
|
835
|
|
|
'F2' => [0.37208, 0.37529, 0.37925, 0.36733, 4230], // Cool White Fluorescent |
|
836
|
|
|
'F3' => [0.40910, 0.39430, 0.41761, 0.38324, 3450], // White Fluorescent |
|
837
|
|
|
'F4' => [0.44018, 0.40329, 0.44920, 0.39074, 2940], // Warm White Fluorescent |
|
838
|
|
|
'F5' => [0.31379, 0.34531, 0.31975, 0.34246, 6350], // Daylight Fluorescent |
|
839
|
|
|
'F6' => [0.37790, 0.38835, 0.38660, 0.37847, 4150], // Lite White Fluorescent |
|
840
|
|
|
'F7' => [0.31292, 0.32933, 0.31569, 0.32960, 6500], // D65 simulator, day light simulator |
|
841
|
|
|
'F8' => [0.34588, 0.35875, 0.34902, 0.35939, 5000], // D50 simulator, Sylvania F40 Design |
|
842
|
|
|
'F9' => [0.37417, 0.37281, 0.37829, 0.37045, 4150], // Cool White Deluxe Fluorescent |
|
843
|
|
|
'F10' => [0.34609, 0.35986, 0.35090, 0.35444, 5000], // Philips TL85, Ultralume 50 |
|
844
|
|
|
'F11' => [0.38052, 0.37713, 0.38541, 0.37123, 4000], // Philips TL84, Ultralume 40 |
|
845
|
|
|
'F12' => [0.43695, 0.40441, 0.44256, 0.39717, 3000], |
|
846
|
|
|
]; // Philips TL83, Ultralume 30 |
|
847
|
|
|
|
|
848
|
|
|
public function multiply($m1, $m2) |
|
849
|
|
|
{ |
|
850
|
|
|
$ni = count($m1); |
|
851
|
|
|
$ki = $ni; |
|
852
|
|
|
$i; |
|
|
|
|
|
|
853
|
|
|
$nj; |
|
|
|
|
|
|
854
|
|
|
$kj = count($m2[0]); |
|
855
|
|
|
$j; |
|
|
|
|
|
|
856
|
|
|
$cols = count($m1[0]); |
|
857
|
|
|
$M = []; |
|
858
|
|
|
$sum; |
|
|
|
|
|
|
859
|
|
|
$nc; |
|
|
|
|
|
|
860
|
|
|
$c; |
|
|
|
|
|
|
861
|
|
|
do { |
|
862
|
|
|
$i = $ki - $ni; |
|
863
|
|
|
$M[$i] = []; |
|
864
|
|
|
$nj = $kj; |
|
865
|
|
|
do { |
|
866
|
|
|
$j = $kj - $nj; |
|
867
|
|
|
$sum = 0; |
|
868
|
|
|
$nc = $cols; |
|
869
|
|
|
do { |
|
870
|
|
|
$c = $cols - $nc; |
|
871
|
|
|
$sum += $m1[$i][$c] * $m2[$c][$j]; |
|
872
|
|
|
} while (($nc -= 1)); |
|
873
|
|
|
$M[$i][$j] = $sum; |
|
874
|
|
|
} while (($nj -= 1)); |
|
875
|
|
|
} while (($ni -= 1)); |
|
876
|
|
|
|
|
877
|
|
|
return $M; |
|
878
|
|
|
} |
|
879
|
|
|
|
|
880
|
|
|
public function determinant($m) |
|
881
|
|
|
{ // 3x3 |
|
882
|
|
|
|
|
883
|
|
|
// print_r($m); |
|
884
|
|
|
return $m[0] * ($m[4] * $m[8] - $m[5] * $m[7]) - |
|
885
|
|
|
$m[1] * ($m[3] * $m[8] - $m[5] * $m[6]) + |
|
886
|
|
|
$m[2] * ($m[3] * $m[7] - $m[4] * $m[6]); |
|
887
|
|
|
} |
|
888
|
|
|
|
|
889
|
|
|
public function inverse($m) |
|
890
|
|
|
{ // 3x3 |
|
891
|
|
|
|
|
892
|
|
|
$d = 1.0 / $this->determinant($m); |
|
893
|
|
|
|
|
894
|
|
|
return [ |
|
895
|
|
|
$d * ($m[4] * $m[8] - $m[5] * $m[7]), |
|
896
|
|
|
$d * (-1 * ($m[1] * $m[8] - $m[2] * $m[7])), |
|
897
|
|
|
$d * ($m[1] * $m[5] - $m[2] * $m[4]), |
|
898
|
|
|
$d * (-1 * ($m[3] * $m[8] - $m[5] * $m[6])), |
|
899
|
|
|
$d * ($m[0] * $m[8] - $m[2] * $m[6]), |
|
900
|
|
|
$d * (-1 * ($m[0] * $m[5] - $m[2] * $m[3])), |
|
901
|
|
|
$d * ($m[3] * $m[7] - $m[4] * $m[6]), |
|
902
|
|
|
$d * (-1 * ($m[0] * $m[7] - $m[1] * $m[6])), |
|
903
|
|
|
$d * ($m[0] * $m[4] - $m[1] * $m[3]), |
|
904
|
|
|
]; |
|
905
|
|
|
} |
|
906
|
|
|
} |
|
907
|
|
|
|
This error can happen if you refactor code and forget to move the variable initialization.
Let’s take a look at a simple example:
The above code is perfectly fine. Now imagine that we re-order the statements:
In that case,
$xwould be read before it is initialized. This was a very basic example, however the principle is the same for the found issue.