Implémentation de CIEDE2000 en VBA

Cette page présente une implémentation de référence de la formule de différence de couleur CIEDE2000 en VBA. Si vous souhaitez obtenir une correspondance exacte avec des implémentations tierces jusqu’à 10 décimales, vous pourriez avoir à apporter des modifications au code source, notamment en commentant et décommentant quelques lignes, ce qui peut être appliqué automatiquement via le lien suivant.

La fonction ΔE2000 en VBA

Considérons la plus courante et académique (Sharma, 2005) des deux formulations.

' Cette fonction écrite en VBA est placée dans le domaine public et
' n’est pas affiliée à la CIE (Commission Internationale de l’Éclairage).

' Bonne pratique pour partager un code source.
Option Explicit

' L’implémentation CIEDE2000 classique qui accepte deux couleurs L*a*b* et renvoie leur différence.
' La composante "L" varie de 0 à 100. "a" et "b", non bornées, sont souvent projetées entre -128 et 127.
Public Function ciede_2000(l_1 As Double, a_1 As Double, b_1 As Double, l_2 As Double, a_2 As Double, b_2 As Double) As Double
	' Traite la formule de différence de couleurs CIEDE2000 en VBA.
	' k_l, k_c et k_h sont des facteurs paramétriques qu’on ajuste
	' selon des recommandations propres au secteur industriel.
	Const M_PI = 3.14159265358979323846264338328, k_l = 1.0, k_c = 1.0, k_h = 1.0
	Dim n As Double, c_1 As Double, c_2 As Double, h_1 As Double, h_2 As Double
	Dim h_m As Double, h_d As Double, p As Double, r_t As Double, l As Double
	Dim t As Double, h As Double, c As Double
	n = (Sqr(a_1 * a_1 + b_1 * b_1) + Sqr(a_2 * a_2 + b_2 * b_2)) * 0.5
	n = n * n * n * n * n * n * n
	' Un facteur impliquant la chroma moyenne à la puissance 7,
	' conçu pour modéliser plus précisément son influence.
	n = 1.0 + 0.5 * (1.0 - Sqr(n / (n + 6103515625.0)))
	' Application du facteur de correction de la chroma pour compenser sa non-linéarité.
	c_1 = Sqr(a_1 * a_1 * n * n + b_1 * b_1)
	c_2 = Sqr(a_2 * a_2 * n * n + b_2 * b_2)
	' Utilisation de 14 lignes pour simuler atan2, car VBA n’a pas cette fonction.
	If 0.0 < a_1 Then
		h_1 = Atn(b_1 / (a_1 * n)) - (b_1 < 0.0) * 2.0 * M_PI
	ElseIf a_1 < 0.0 Then
		h_1 = Atn(b_1 / (a_1 * n)) + M_PI
	Else
		h_1 = M_PI + ((0.0 < b_1) - (b_1 < 0.0)) * 0.5 * M_PI
	End If
	If 0.0 < a_2 Then
		h_2 = Atn(b_2 / (a_2 * n)) - (b_2 < 0.0) * 2.0 * M_PI
	ElseIf a_2 < 0.0 Then
		h_2 = Atn(b_2 / (a_2 * n)) + M_PI
	Else
		h_2 = M_PI + ((0.0 < b_2) - (b_2 < 0.0)) * 0.5 * M_PI
	End If
	' Le calcul avec atan2 adapté à CIEDE2000 est terminé.
	n = Abs(h_2 - h_1)
	' Prévient le branchement de dépendre du RoundingMode du langage de programmation.
	If M_PI - 1E-14 < n And n < M_PI + 1E-14 Then n = M_PI
	' Lorsque les angles de teinte sont dans différents quadrants,
	' la moyenne arithmétique simple peut donner un angle incorrect,
	' les lignes suivantes prennent en compte cette correction angulaire.
	h_m = (h_1 + h_2) * 0.5
	h_d = (h_2 - h_1) * 0.5
	If M_PI < n Then
		h_d = h_d + M_PI
		' 📜 La formulation de Sharma n’utilise pas la ligne suivante, mais plutôt celle d’après.
		' Note : ces deux variantes ne diffèrent que de ±0,0003 sur la différence de couleur finale.
		h_m = h_m + M_PI
		' If h_m < M_PI Then h_m = h_m + M_PI Else h_m = h_m - M_PI
	End If
	p = 36.0 * h_m - 55.0 * M_PI
	n = (c_1 + c_2) * 0.5
	n = n * n * n * n * n * n * n
	' Le terme de correction de la rotation de teinte ajuste le comportement
	' de l’algorithme, d’autant plus si la comparaison porte sur des teintes bleues.
	r_t = -2.0 * Sqr(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)
	' Luminosité.
	l = (l_2 - l_1) / (k_l * (1.0 + 0.015 * n / Sqr(20.0 + n)))
	' Ces coefficients modulent l’influence des composantes
	' harmoniques dans le calcul de la différence de teinte.
	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.2 * Sin(4.0 * h_m + 3.0 * M_PI / 20.0)
	n = c_1 + c_2
	' Teinte.
	h = 2.0 * Sqr(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))
	' Retourner la racine carrée assure que dE00 représente une distance
	' géométrique (comprise entre 0 et environ 185) dans l’espace CIELAB.
	ciede_2000 = Sqr(l * l + h * h + c * c + c * h * r_t)
End Function

'    Projet GitHub : https://github.com/michel-leonard/ciede2000-color-matching
'   Tests en ligne : https://michel-leonard.github.io/ciede2000-color-matching

' L1 = 85.3   a1 = 21.6   b1 = 3.9
' L2 = 87.7   a2 = 27.3   b2 = -4.2
' CIE ΔE00 = 6.1457265560 (Bruce Lindbloom, Netflix’s VMAF, ...)
' CIE ΔE00 = 6.1457093985 (Gaurav Sharma, OpenJDK, ...)
' Écart entre les implémentations ≈ 1.7e-5

' Voir les commentaires du code source pour passer d’une de ces variantes d’implémentation de ΔE*00 à l’autre.

Précision et fiabilité du code source

La différence entre les formulations de Sharma et Lindbloom ne dépasse jamais ±0,0003 sur le ΔE2000 final, ce qui correspond à l’écart habituel mesuré entre deux implémentations 32 bits et est imperceptible à l’œil humain. Nos implémentations 64 bits, toutes cohérentes entre elles, garantissent au moins 10 décimales correctes, de sorte que le choix d’une formulation plutôt qu’une autre dépend surtout de l’interopérabilité souhaitée. La formulation qui apparait par défaut sur cette page est la plus couramment utilisée (son micro-avantage réside dans son ancrage communautaire et dans sa légèreté supérieure à celle de son analogue lorsque vectorisée).

✎ Si vous trouvez dans le code source un commentaire qui ne correspond pas à une autre langue, veuillez en informer l’auteur du site, qui étudiera votre suggestion et l’intégrera dans le code source.

Comment convertir les couleurs RGB en L*a*b* ?

Rendez-vous sur la page AWK, C, Dart, Java, JavaScript, Kotlin, Lua, PHP, Python, Ruby ou Rust où un tel convertisseur (utilisant l’illuminant D65) est déjà implémenté en plus de la fonction de comparaison de couleurs.

Plages de valeurs dans CIELAB et interprétation du ΔE2000

Dans l’espace colorimétrique CIELAB, la composante L* représente la luminosité et varie de 0 (noir) à 100 (blanc). Les composantes a* et b* décrivent les axes de couleur : a* s’étend du vert au rouge, tandis que b* va du bleu au jaune. Dans la pratique, les valeurs de a* et b* se situent dans la plage -128 à +127, même si elles peuvent légèrement dépasser ces limites selon les conversions colorimétriques.

Le ΔE2000 (CIEDE2000) quantifie la différence perceptuelle entre deux couleurs : 0 signifie deux couleurs identiques, et des valeurs plus élevées (jusqu’à environ 185 dans les cas extrêmes) indiquent une différence plus marquée. Par exemple, une valeur ΔE2000 autour de 5 correspond à des couleurs proches, tandis qu’une valeur autour de 15 correspond à des couleurs clairement distinctes.

Exemple d’utilisation en VBA

' Compute the Delta E (CIEDE2000) color difference between two L*a*b* colors in VBA

' Color 1: L1 = 96.4   a1 = 45.4   b1 = -3.0
' Color 2: L2 = 95.9   a2 = 50.4   b2 = 3.7

Dim deltaE As Double
deltaE = ciede_2000(l1, a1, b1, l2, a2, b2)
Debug.Print deltaE

' .................................................. This shows a ΔE2000 of 3.7852068385
' As explained in the comments, compliance with Gaurav Sharma would display 3.7852203702

Résultats des tests

Le driver écrit en langage C99, doté de 250 tests statiques précis, a prouvé que cette fonction VBA est interopérable avec la fonction CIEDE2000 mise à disposition dans les autres langages de programmation.

CIEDE2000 Verification Summary :
  First Verified Line : 31.8,56.39,27.1,61,85,111.57,37.83481076432601
             Duration : 43.31 s
            Successes : 10000000
               Errors : 0
      Average Delta E : 62.9544
    Average Deviation : 8.8259009656255217e-15
    Maximum Deviation : 2.9842794901924208e-13

Fichiers à télécharger

Utilisez librement ces fichiers mis à disposition par Michel, même à des fins commerciales.

Communauté

Si vous souhaitez laisser votre avis sur ce code source VBA ou sur CIEDE2000 en général, le livre d’or contient déjà 1 messages en français, et 9 messages au total, alors laissez-nous savoir ce que vous pensez.