Реализация CIEDE2000 на PHP

На этой странице представлена эталонная реализация формулы цветовой разницы CIEDE2000 на PHP. Если вы хотите получить точное совпадение со сторонними реализациями до 10 знаков после запятой, вам может потребоваться внести некоторые изменения в исходный код, в частности, закомментировать и разкомментировать несколько строк, которые могут быть применены автоматически по ссылке ниже.

Функция ΔE2000 в PHP

Рассмотрим более распространенную и академическую (Sharma, 2005) из двух формулировок.

<?php

// This function written in PHP is not affiliated with the CIE (International Commission on Illumination),
// and is released into the public domain. It is provided "as is" without any warranty, express or implied.

// The classic CIE ΔE2000 implementation, which operates on two L*a*b* colors, and returns their difference.
// "l" ranges from 0 to 100, while "a" and "b" are unbounded and commonly clamped to the range of -128 to 127.
function ciede_2000($l_1, $a_1, $b_1, $l_2, $a_2, $b_2) {
	// Working in PHP with the CIEDE2000 color-difference formula.
	// k_l, k_c, k_h are parametric factors to be adjusted according to
	// different viewing parameters such as textures, backgrounds...
	$k_l = $k_c = $k_h = 1.0;
	$n = (sqrt($a_1 * $a_1 + $b_1 * $b_1) + sqrt($a_2 * $a_2 + $b_2 * $b_2)) * 0.5;
	$n = $n * $n * $n * $n * $n * $n * $n;
	// A factor involving chroma raised to the power of 7 designed to make
	// the influence of chroma on the total color difference more accurate.
	$n = 1.0 + 0.5 * (1.0 - sqrt($n / ($n + 6103515625.0)));
	// Application of the chroma correction factor.
	$c_1 = sqrt($a_1 * $a_1 * $n * $n + $b_1 * $b_1);
	$c_2 = sqrt($a_2 * $a_2 * $n * $n + $b_2 * $b_2);
	// atan2 is preferred over atan because it accurately computes the angle of
	// a point (x, y) in all quadrants, handling the signs of both coordinates.
	$h_1 = atan2($b_1, $a_1 * $n);
	$h_2 = atan2($b_2, $a_2 * $n);
	$h_1 += 2.0 * M_PI * ($h_1 < 0.0);
	$h_2 += 2.0 * M_PI * ($h_2 < 0.0);
	$n = abs($h_2 - $h_1);
	// Cross-implementation consistent rounding.
	if (M_PI - 1E-14 < $n && $n < M_PI + 1E-14)
		$n = M_PI;
	// When the hue angles lie in different quadrants, the straightforward
	// average can produce a mean that incorrectly suggests a hue angle in
	// the wrong quadrant, the next lines handle this issue.
	$h_m = ($h_1 + $h_2) * 0.5;
	$h_d = ($h_2 - $h_1) * 0.5;
	if (M_PI < $n) {
		$h_d += M_PI;
		// 📜 Sharma’s formulation doesn’t use the next line, but the one after it,
		// and these two variants differ by ±0.0003 on the final color differences.
		$h_m += M_PI;
		// $h_m += $h_m < M_PI ? M_PI : -M_PI;
	}
	$p = 36.0 * $h_m - 55.0 * M_PI;
	$n = ($c_1 + $c_2) * 0.5;
	$n = $n * $n * $n * $n * $n * $n * $n;
	// The hue rotation correction term is designed to account for the
	// non-linear behavior of hue differences in the blue region.
	$r_t = -2.0 * sqrt($n / ($n + 6103515625.0))
			   * sin(M_PI / 3.0 * exp($p * $p / (-25.0 * M_PI * M_PI)));
	$n = ($l_1 + $l_2) * 0.5;
	$n = ($n - 50.0) * ($n - 50.0);
	// Lightness.
	$l = ($l_2 - $l_1) / ($k_l * (1.0 + 0.015 * $n / sqrt(20.0 + $n)));
	// These coefficients adjust the impact of different harmonic
	// components on the hue difference calculation.
	$t = 1.0	+ 0.24 * sin(2.0 * $h_m + M_PI * 0.5)
			+ 0.32 * sin(3.0 * $h_m + 8.0 * M_PI / 15.0)
			- 0.17 * sin($h_m + M_PI / 3.0)
			- 0.20 * sin(4.0 * $h_m + 3.0 * M_PI / 20.0);
	$n = $c_1 + $c_2;
	// Hue.
	$h = 2.0 * sqrt($c_1 * $c_2) * sin($h_d) / ($k_h * (1.0 + 0.0075 * $n * $t));
	// Chroma.
	$c = ($c_2 - $c_1) / ($k_c * (1.0 + 0.0225 * $n));
	// Returning the square root ensures that dE00 accurately reflects the
	// geometric distance in color space, which can range from 0 to around 185.
	return sqrt($l * $l + $h * $h + $c * $c + $c * $h * $r_t);
}

// GitHub Project : https://github.com/michel-leonard/ciede2000-color-matching
//   Online Tests : https://michel-leonard.github.io/ciede2000-color-matching

// L1 = 39.9   a1 = 45.8   b1 = -2.0
// L2 = 40.7   a2 = 51.5   b2 = 2.4
// CIE ΔE00 = 2.9464564049 (Bruce Lindbloom, Netflix’s VMAF, ...)
// CIE ΔE00 = 2.9464695839 (Gaurav Sharma, OpenJDK, ...)
// Deviation between implementations ≈ 1.3e-5

// See the source code comments for easy switching between these two widely used ΔE*00 implementation variants.

Точность и надежность исходного кода

Разница между формулировками Шармы и Линдблума никогда не превышает ±0,0003 по финальному ΔE2000, что соответствует обычной разнице между двумя 32-битными реализациями и незаметно для невооружённого глаза. Наши 64-битные реализации, все согласованные между собой, гарантируют как минимум 10 правильных десятичных знаков, поэтому выбор одной формулы вместо другой в основном зависит от желаемой совместимости. Наиболее часто используется формула, которая по умолчанию отображается на этой странице (ее микродостоинство в том, что она основана на сообществе и при векторизации легче своего аналога).

✎ Если вы обнаружили в исходном коде комментарий, который не соответствует другому языку, пожалуйста, сообщите об этом автору сайта, который изучит ваше предложение и включит его в исходный код.

Как преобразовать цвета RGB в L*a*b*?

Для преобразования необходимо использовать промежуточное цветовое пространство XYZ, и если вам нужна помощь, исходный код приведен в нижней части этой страницы (используется точка белого D65, формализованная в 1964 году).

Диапазоны значений в CIELAB и интерпретация ΔE2000

В цветовом пространстве CIELAB компонент L* обозначает светлоту и обычно изменяется от 0 (черный) до 100 (белый). Компоненты a* и b* описывают цветовые оси: a* идет от зеленого к красному, а b* — от синего к желтому. На практике значения a* и b* чаще всего находятся в диапазоне от -128 до +127, хотя в зависимости от конверсии они могут немного выходить за эти пределы.

Параметры k_l, k_c и k_h

Параметры k_l, k_c и k_h являются весовыми коэффициентами, применяемыми к компонентам яркости (ΔL*), хроматичности (ΔC*) и оттенка (ΔH*) в формуле CIEDE2000. Их значение по умолчанию — 1, что соответствует стандартным условиям наблюдения, рекомендованным Международная комиссия по освещению. На практике эти коэффициенты настраиваются для отражения конкретных условий: например, k_l = 2 иногда используется для придания большего веса различиям яркости (часто в печати), тогда как k_c или k_h могут быть уменьшены для повышения допустимости вариаций насыщенности или оттенка в зависимости от требований контроля качества. В зависимости от контекста эти коэффициенты обычно варьируются от 0,5 до 2.

ΔE2000 (CIEDE2000) измеряет воспринимаемую разницу между двумя цветами: 0 означает идентичные цвета, а более высокие значения (до примерно 185 в крайних случаях) показывают большую разницу. Например, значение ΔE2000 около 5 означает близкие цвета, а около 15 — явно разные.

Пример использования в PHP

// Define color 1 in L*a*b* space
$l1 = 57.0;
$a1 = 17.0;
$b1 = -3.6;

// Define color 2 in L*a*b* space
$l2 = 56.6;
$a2 = 21.0;
$b2 = 4.6;

// Calculate the ΔE*00 color difference
$deltaE = ciede_2000($l1, $a1, $b1, $l2, $a2, $b2);
echo $deltaE;

// .................................................. This shows a ΔE2000 of 5.9908949846
// As explained in the comments, compliance with Gaurav Sharma would display 5.9909105284

Результаты испытаний

Драйвер, написанный на языке C99, с помощью 250 точных статических тестов доказал, что эта функция PHP совместима с функцией CIEDE2000, доступной на других языках программирования.

CIEDE2000 Verification Summary :
  First Verified Line : 34,30,-17,69.8,60,-33,36.861058302008
             Duration : 48.97 s
            Successes : 10000000
               Errors : 0
      Average Delta E : 62.9515
    Average Deviation : 4.5546872117219198e-13
    Maximum Deviation : 5.0590642786119133e-12

Файлы для загрузки

Не стесняйтесь использовать эти файлы, предоставленные Мишелем, даже в коммерческих целях.

Автор веб-сайта и исходного кода колориметрии

ktm620enduro

Сообщество

Если вы хотите оставить своё мнение об этом исходном коде PHP или о CIEDE2000 в целом, то в гостевой книге уже есть 1 сообщения на русском языке и 9 сообщений в целом, так что, пожалуйста, дайте нам знать, что вы думаете по этому поводу.