1use crate::{ToAlpha, XyzD65, round_dp};
2use core::fmt;
3
4#[derive(Debug, Clone, Copy, PartialEq)]
11pub struct Oklab {
12 pub lightness: f64,
13 pub a: f64,
14 pub b: f64,
15 pub alpha: f32,
16}
17
18impl Oklab {
19 pub fn new(lightness: f64, a: f64, b: f64, alpha: f32) -> Self {
20 Self {
21 lightness: lightness.clamp(0.0, 100.0),
22 a: a.clamp(-128.0, 127.0),
23 b: b.clamp(-128.0, 127.0),
24 alpha: alpha.clamp(0.0, 100.0),
25 }
26 }
27}
28
29impl ToAlpha for Oklab {
30 fn to_alpha(&self) -> f32 {
31 self.alpha
32 }
33}
34
35impl fmt::Display for Oklab {
36 fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
37 let Self { lightness, a, b, alpha } = self;
38 write!(f, "oklab({} {} {}", round_dp(*lightness, 5), round_dp(*a, 3), round_dp(*b, 3))?;
39 if *alpha < 100.0 {
40 write!(f, " / {}", round_dp(*alpha as f64, 2))?;
41 }
42 write!(f, ")")
43 }
44}
45
46impl From<XyzD65> for Oklab {
47 fn from(value: XyzD65) -> Self {
48 let XyzD65 { x, y, z, alpha } = value;
49 let x = x / 100.0;
50 let y = y / 100.0;
51 let z = z / 100.0;
52 let l = x * 0.819_022_437_996_703 + y * 0.3619062600528904 + z * (-0.1288737815209879);
53 let m = x * 0.0329836539323885 + y * 0.9292868615863434 + z * 0.0361446663506424;
54 let s = x * 0.0481771893596242 + y * 0.2642395317527308 + z * 0.6335478284694309;
55 let l = l.cbrt();
56 let m = m.cbrt();
57 let s = s.cbrt();
58 let lightness = l * 0.210_454_268_309_314 + m * 0.7936177747023054 + s * (-0.0040720430116193);
59 let a = l * 1.9779985324311684 + m * (-2.428_592_242_048_58) + s * 0.450_593_709_617_411;
60 let b = l * 0.0259040424655478 + m * 0.7827717124575296 + s * (-0.8086757549230774);
61 Oklab::new(lightness, a, b, alpha)
62 }
63}
64
65impl From<Oklab> for XyzD65 {
66 fn from(value: Oklab) -> Self {
67 let Oklab { lightness, a, b, alpha } = value;
68 let l = lightness * 1.0000000000000000 + a * 0.3963377773761749 + b * 0.2158037573099136;
69 let m = lightness * 1.0000000000000000 + a * (-0.1055613458156586) + b * (-0.0638541728258133);
70 let l = l.powi(3);
71 let s = lightness * 1.0000000000000000 + a * (-0.0894841775298119) + b * (-1.2914855480194092);
72 let m = m.powi(3);
73 let s = s.powi(3);
74 let x = l * 1.2268798758459243 + m * (-0.5578149944602171) + s * 0.2813910456659647;
75 let y = l * (-0.0405757452148008) + m * 1.112_286_803_280_317 + s * (-0.0717110580655164);
76 let z = l * (-0.0763729366746601) + m * (-0.4214933324022432) + s * 1.5869240198367816;
77 XyzD65::new(x * 100.0, y * 100.0, z * 100.0, alpha)
78 }
79}