chromashift/
oklab.rs

1use crate::{XyzD65, round_dp};
2use core::fmt;
3
4/// A more adequate expression of LAB, in the CIE colour space.
5/// The components are:
6/// - L - a number between 0.0 and 100.0
7/// - A - a number between -128.0 and +127.0
8/// - B - a number between -128.0 and +127.0
9/// - Alpha - a number between 0.0 and 100.0
10#[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 fmt::Display for Oklab {
30	fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
31		let Self { lightness, a, b, alpha } = self;
32		write!(f, "oklab({} {} {}", round_dp(*lightness, 5), round_dp(*a, 3), round_dp(*b, 3))?;
33		if *alpha < 100.0 {
34			write!(f, " / {}", round_dp(*alpha as f64, 2))?;
35		}
36		write!(f, ")")
37	}
38}
39
40impl From<XyzD65> for Oklab {
41	fn from(value: XyzD65) -> Self {
42		let XyzD65 { x, y, z, alpha } = value;
43		let x = x / 100.0;
44		let y = y / 100.0;
45		let z = z / 100.0;
46		let l = x * 0.819_022_437_996_703 + y * 0.3619062600528904 + z * (-0.1288737815209879);
47		let m = x * 0.0329836539323885 + y * 0.9292868615863434 + z * 0.0361446663506424;
48		let s = x * 0.0481771893596242 + y * 0.2642395317527308 + z * 0.6335478284694309;
49		let l = l.cbrt();
50		let m = m.cbrt();
51		let s = s.cbrt();
52		let lightness = l * 0.210_454_268_309_314 + m * 0.7936177747023054 + s * (-0.0040720430116193);
53		let a = l * 1.9779985324311684 + m * (-2.428_592_242_048_58) + s * 0.450_593_709_617_411;
54		let b = l * 0.0259040424655478 + m * 0.7827717124575296 + s * (-0.8086757549230774);
55		Oklab::new(lightness, a, b, alpha)
56	}
57}
58
59impl From<Oklab> for XyzD65 {
60	fn from(value: Oklab) -> Self {
61		let Oklab { lightness, a, b, alpha } = value;
62		let l = lightness * 1.0000000000000000 + a * 0.3963377773761749 + b * 0.2158037573099136;
63		let m = lightness * 1.0000000000000000 + a * (-0.1055613458156586) + b * (-0.0638541728258133);
64		let l = l.powi(3);
65		let s = lightness * 1.0000000000000000 + a * (-0.0894841775298119) + b * (-1.2914855480194092);
66		let m = m.powi(3);
67		let s = s.powi(3);
68		let x = l * 1.2268798758459243 + m * (-0.5578149944602171) + s * 0.2813910456659647;
69		let y = l * (-0.0405757452148008) + m * 1.112_286_803_280_317 + s * (-0.0717110580655164);
70		let z = l * (-0.0763729366746601) + m * (-0.4214933324022432) + s * 1.5869240198367816;
71		XyzD65::new(x * 100.0, y * 100.0, z * 100.0, alpha)
72	}
73}
chromashift/oklab.rs

chromashift/
oklab.rs
