Oklab color space

The Oklab color space is a uniform color space for device independent color designed to improve perceptual uniformity, hue and lightness prediction, color blending, and usability while ensuring numerical stability and ease of implementation. Introduced by Björn Ottosson in December 2020, Oklab and its cylindrical counterpart, Oklch, have been included in the CSS Color Level 4 and Level 5 drafts for device-independent web colors since December 2021. They are supported by recent versions of major web browsers and allow the specification of wide-gamut P3 colors.
Oklab's model is fitted with improved color appearance data: CAM16 data for lightness and chroma, and IPT data for hue. The new fit addresses issues such as unexpected hue and lightness changes in blue colors present in the CIELAB color space, simplifying the creation of color schemes and smoother color gradients.
As Ottosson explained, he chose the name Oklab because the model does an OK (adequate) job and is based on the three color-space coordinates L, a, and b.
Oklab uses the same spatial structure as CIELAB, representing color using three components:
Like CIELCh, Oklch represents colors using:
Neutral greys, pure black and the reference white are achromatic, that is, C = 0 {\displaystyle C=0} , a = 0 {\displaystyle a=0} , b = 0 {\displaystyle b=0} , and h is undefined. Assigning any real value to their hue component has no effect on conversions between color spaces.
The perceptual color difference in Oklab is calculated as the Euclidean distance between the (L, a, b) coordinates.
Like CIELCh, the Cartesian coordinates a and b are converted to the polar coordinates C and h as follows:
C = a 2 + b 2 h = atan2 ⁡ ( b , a ) {\displaystyle {\begin{aligned}C&={\sqrt {a^{2}+b^{2}}}\\h&=\operatorname {atan2} (b,a)\end{aligned}}}
And the polar coordinates are converted to the Cartesian coordinates as follows:
a = C cos ⁡ ( h ) b = C sin ⁡ ( h ) {\displaystyle {\begin{aligned}a&=C\cos(h)\\b&=C\sin(h)\end{aligned}}}
Converting from CIE XYZ with a Standard Illuminant D65 involves:
Given:
M 1 = [ 0.8189330101 − 0.3618667424 − 0.1288597137 0.0329845436 − 0.9293118715 − 0.0361456387 0.0482003018 − 0.2643662691 − 0.6338517070 ] M 2 = [ 0.2104542553 − 0.7936177850 − 0.0040720468 1.9779984951 − 2.4285922050 − 0.4505937099 0.0259040371 − 0.7827717662 − 0.8086757660 ] {\displaystyle {\begin{aligned}\mathbf {M} _{1}&={\begin{bmatrix}0.8189330101&{\phantom {-}}0.3618667424&-0.1288597137\\0.0329845436&{\phantom {-}}0.9293118715&{\phantom {-}}0.0361456387\\0.0482003018&{\phantom {-}}0.2643662691&{\phantom {-}}0.6338517070\end{bmatrix}}\\\mathbf {M} _{2}&={\begin{bmatrix}0.2104542553&{\phantom {-}}0.7936177850&-0.0040720468\\1.9779984951&-2.4285922050&{\phantom {-}}0.4505937099\\0.0259040371&{\phantom {-}}0.7827717662&-0.8086757660\end{bmatrix}}\end{aligned}}}
Converting from sRGB requires first converting from sRGB to CIE XYZ with a Standard Illuminant D65. As the last step of this conversion is a linear map from linear RGB to CIE XYZ, the reference implementation directly employs the multiplied matrix representing the composition of the two linear maps:
Converting to CIE XYZ and sRGB simply involves applying the respective inverse functions in reverse order:
[ l ′ m ′ s ′ ] = M 2 − 1 [ L a b ] [ l m s ] = [ ( l ′ ) 3 ( m ′ ) 3 ( s ′ ) 3 ] [ X Y Z ] = M 1 − 1 [ l m s ] {\displaystyle {\begin{aligned}{\begin{bmatrix}l'\\m'\\s'\end{bmatrix}}&=\mathbf {M} _{2}^{-1}{\begin{bmatrix}L\\a\\b\end{bmatrix}}\\{\begin{bmatrix}l\\m\\s\end{bmatrix}}&={\begin{bmatrix}\left(l'\right)^{3}\\\left(m'\right)^{3}\\\left(s'\right)^{3}\end{bmatrix}}\\{\begin{bmatrix}X\\Y\\Z\end{bmatrix}}&=\mathbf {M} _{1}^{-1}{\begin{bmatrix}l\\m\\s\end{bmatrix}}\end{aligned}}}