Implementação do CIEDE2000 em SQL

Esta página apresenta uma implementação de referência da fórmula de diferença de cor CIEDE2000 em SQL. Se desejar obter uma correspondência exacta com implementações de terceiros até 10 casas decimais, poderá ter de fazer algumas alterações ao código-fonte, em particular comentando e descomentando algumas linhas, que podem ser aplicadas automaticamente através da ligação abaixo.

A função ΔE2000 em SQL

Consideremos a mais comum e académica (Sharma, 2005) das duas formulações.

-- This function written in SQL 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.

DELIMITER //

-- Delete any function of the same name that already exists
DROP FUNCTION IF EXISTS ciede_2000 //

-- 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.
CREATE FUNCTION ciede_2000(l_1 DOUBLE, a_1 DOUBLE, b_1 DOUBLE, l_2 DOUBLE, a_2 DOUBLE, b_2 DOUBLE)
RETURNS DOUBLE
DETERMINISTIC
NO SQL
BEGIN
	-- Working in SQL/PSM 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...
	DECLARE k_l, k_c, k_h DOUBLE DEFAULT 1.0;
	DECLARE n, c_1, c_2, h_1, h_2, h_m, h_d, r_t, p, t, l, c, h DOUBLE;

	SET n = (SQRT(a_1 * a_1 + b_1 * b_1) + SQRT(a_2 * a_2 + b_2 * b_2)) * 0.5;
	SET 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.
	SET n = 1.0 + 0.5 * (1.0 - SQRT(n / (n + 6103515625.0)));
	-- Application of the chroma correction factor.
	SET c_1 = SQRT(a_1 * a_1 * n * n + b_1 * b_1);
	SET 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.
	SET h_1 = COALESCE(ATAN2(b_1, a_1 * n), 0);
	SET h_2 = COALESCE(ATAN2(b_2, a_2 * n), 0);
	IF h_1 < 0 THEN SET h_1 = h_1 + 2 * PI(); END IF;
	IF h_2 < 0 THEN SET h_2 = h_2 + 2 * PI(); END IF;
	SET n = ABS(h_2 - h_1);
	-- Cross-implementation consistent rounding.
	IF PI() - 1E-14 < n AND n < PI() + 1E-14 THEN SET n = PI(); END IF;
	-- 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.
	SET h_m = (h_1 + h_2) * 0.5;
	SET h_d = (h_2 - h_1) * 0.5;
	IF PI() < n THEN
		SET h_d = h_d + 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.
		SET h_m = h_m + PI();
		-- SET h_m = h_m + CASE WHEN h_m < PI() THEN PI() ELSE -PI() END;
	END IF;
	SET p = 36.0 * h_m - 55.0 * PI();
	SET n = (c_1 + c_2) * 0.5;
	SET 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.
	SET r_t = -2.0 * SQRT(n / (n + 6103515625.0)) * SIN(PI() / 3.0 * EXP(p * p / (-25.0 * PI() * PI())));
	SET n = (l_1 + l_2) * 0.5;
	SET n = (n - 50.0) * (n - 50.0);
	-- Lightness.
	SET 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.
	SET t = 1.0	+ 0.24 * SIN(2.0 * h_m + PI() * 0.5)
			+ 0.32 * SIN(3.0 * h_m + 8.0 * PI() / 15.0)
			- 0.17 * SIN(h_m + PI() / 3.0)
			- 0.20 * SIN(4.0 * h_m + 3.0 * PI() / 20.0);
	SET n = c_1 + c_2;
	-- Hue.
	SET h = 2.0 * SQRT(c_1 * c_2) * SIN(h_d) / (k_h * (1.0 + 0.0075 * n * t));
	-- Chroma.
	SET 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);
END //

DELIMITER ;

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

-- L1 = 52.3   a1 = 21.9   b1 = 2.7
-- L2 = 53.8   a2 = 28.0   b2 = -3.1
-- CIE ΔE00 = 5.0119430211 (Bruce Lindbloom, Netflix’s VMAF, ...)
-- CIE ΔE00 = 5.0119254601 (Gaurav Sharma, OpenJDK, ...)
-- Deviation between implementations ≈ 1.8e-5

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

Precisão e fiabilidade do código fonte

A diferença entre as formulações de Sharma e Lindbloom nunca excede ±0,0003 no ΔE2000 final, o que corresponde à diferença habitual medida entre duas implementações de 32 bits e é imperceptível ao olho humano. Nossas implementações de 64 bits, todas consistentes entre si, garantem pelo menos 10 casas decimais corretas, de modo que a escolha de uma formulação em vez de outra depende principalmente da interoperabilidade desejada. A formulação que aparece por defeito nesta página é a mais utilizada (a sua microvantagem reside na sua ancoragem comunitária e na sua maior leveza do que a do seu análogo quando vectorizado).

✎ Se encontrar um comentário no código-fonte que não corresponda a outra língua, informe o autor do sítio, que estudará a sua sugestão e a incorporará no código-fonte.

Como é que se convertem cores RGB em L*a*b*?

Vá para a página AWK, C, Dart, Java, JavaScript, Kotlin, Lua, PHP, Python, Ruby ou Rust onde esse conversor (utilizando o iluminante D65) já está implementado para além da função de comparação de cores.

Intervalos de valores no CIELAB e interpretação do ΔE2000

No espaço de cor CIELAB, o componente L* representa a luminosidade e normalmente varia de 0 (preto) a 100 (branco). Os componentes a* e b* representam os eixos de cor: a* vai do verde ao vermelho, enquanto b* vai do azul ao amarelo. Na prática, os valores de a* e b* costumam estar entre -128 e +127, embora possam ultrapassar ligeiramente esses limites dependendo da conversão de cor.

Parâmetros k_l, k_c e k_h

Os parâmetros k_l, k_c e k_h são fatores de ponderação aplicados aos termos de luminosidade (ΔL*), croma (ΔC*) e matiz (ΔH*) na fórmula CIEDE2000. O seu valor padrão é 1, o que corresponde às condições de observação padrão recomendadas pela Comissão Internacional de Iluminação. Na prática, estes coeficientes são ajustados para refletir condições específicas: por exemplo, k_l = 2 é por vezes usado para dar mais peso às diferenças de luminosidade (comum em impressão), enquanto k_c ou k_h podem ser reduzidos para aumentar a tolerância a variações de saturação ou matiz conforme as exigências do controlo de qualidade. Dependendo do contexto, estes coeficientes variam tipicamente entre 0,5 e 2.

ΔE2000 (CIEDE2000) mede a diferença perceptível entre duas cores: 0 significa cores idênticas, e valores maiores (até cerca de 185 em casos extremos) indicam uma diferença mais significativa. Por exemplo, um valor ΔE2000 em torno de 5 indica cores próximas, enquanto em torno de 15 indica cores claramente diferentes.

Exemplo de utilização em SQL

-- Example usage of the ΔE*00 function in SQL

-- Color 1: l1 = 94.1   a1 = 30.7   b1 = 2.9
-- Color 2: l2 = 92.2   a2 = 26.4   b2 = -2.2

SELECT ciede_2000(l1, a1, b1, l2, a2, b2) AS delta_e;

-- .................................................. This shows a ΔE2000 of 3.8819773139
-- As explained in the comments, compliance with Gaurav Sharma would display 3.8819904826

Resultados dos testes

O driver escrito na linguagem C99, com 250 testes estáticos precisos, provou que esta função SQL é interoperável com a função CIEDE2000 disponível noutras linguagens de programação.

CIEDE2000 Verification Summary :
  First Verified Line : 57.7,61,95.1,50.49,-53,-91.7,72.48542739412765
             Duration : 44.06 s
            Successes : 10000000
               Errors : 0
      Average Delta E : 62.9405
    Average Deviation : 4.2545206732635951e-15
    Maximum Deviation : 1.1368683772161603e-13

Ficheiros para descarregar

Pode utilizar livremente estes ficheiros disponibilizados pelo Michel, mesmo para fins comerciais.

Comunidade

Se quiser deixar a sua opinião sobre este código fonte SQL ou sobre o CIEDE2000 em geral, o livro de visitas já contém 1 mensagens em português, e 9 mensagens no total, por isso diga-nos o que pensa.