chromashift/
lab.rs

1use crate::{ToAlpha, XyzD50, round_dp};
2use core::fmt;
3
4const D50X: f64 = 96.4220;
5const D50Y: f64 = 100.0;
6const D50Z: f64 = 82.5210;
7
8/// An CIE defined colour space representing L - perceptual lightness, and two axes A & B.
9/// The components are:
10/// - L - a number between 0.0 and 100.0
11/// - A - a number between -125.0 and +125.0
12/// - B - a number between -125.0 and +125.0
13/// - Alpha - a number between 0.0 and 100.0
14#[derive(Debug, Clone, Copy, PartialEq)]
15pub struct Lab {
16	pub lightness: f64,
17	pub a: f64,
18	pub b: f64,
19	pub alpha: f32,
20}
21
22impl Lab {
23	pub fn new(lightness: f64, a: f64, b: f64, alpha: f32) -> Self {
24		Self {
25			lightness: lightness.clamp(0.0, 100.0),
26			a: a.clamp(-125.0, 125.0),
27			b: b.clamp(-125.0, 125.0),
28			alpha: alpha.clamp(0.0, 100.0),
29		}
30	}
31}
32
33impl ToAlpha for Lab {
34	fn to_alpha(&self) -> f32 {
35		self.alpha
36	}
37}
38
39impl fmt::Display for Lab {
40	fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
41		let Self { lightness, a, b, alpha } = self;
42		write!(f, "lab({} {} {}", round_dp(*lightness, 2), round_dp(*a, 3), round_dp(*b, 3))?;
43		if *alpha < 100.0 {
44			write!(f, " / {}", round_dp(*alpha as f64, 2))?;
45		}
46		write!(f, ")")
47	}
48}
49
50impl From<XyzD50> for Lab {
51	fn from(value: XyzD50) -> Self {
52		let XyzD50 { x, y, z, alpha } = value;
53		let x = x / D50X;
54		let y = y / D50Y;
55		let z = z / D50Z;
56		let epsilon = 216.0 / 24389.0; // 6^3/29^3
57		let kappa = 24389.0 / 27.0; // 29^3/3^3
58		let fx = if x > epsilon { x.cbrt() } else { (kappa * x + 16.0) / 116.0 };
59		let fy = if y > epsilon { y.cbrt() } else { (kappa * y + 16.0) / 116.0 };
60		let fz = if z > epsilon { z.cbrt() } else { (kappa * z + 16.0) / 116.0 };
61		let lightness = 116.0 * fy - 16.0;
62		let a = 500.0 * (fx - fy);
63		let b = 200.0 * (fy - fz);
64		Lab::new(lightness, a, b, alpha)
65	}
66}
67
68impl From<Lab> for XyzD50 {
69	fn from(value: Lab) -> Self {
70		let Lab { lightness, a, b, alpha } = value;
71		let epsilon = 216.0 / 24389.0; // 6^3/29^3
72		let kappa = 24389.0 / 27.0; // 29^3/3^3
73		let fy = (lightness + 16.0) / 116.0;
74		let fx = a / 500.0 + fy;
75		let fz = fy - b / 200.0;
76		let x = if fx.powi(3) > epsilon { fx.powi(3) } else { (116.0 * fx - 16.0) / kappa };
77		let y = if lightness > kappa * epsilon { ((lightness + 16.0) / 116.0).powi(3) } else { lightness / kappa };
78		let z = if fz.powi(3) > epsilon { fz.powi(3) } else { (116.0 * fz - 16.0) / kappa };
79		XyzD50::new(x * D50X, y * D50Y, z * D50Z, alpha)
80	}
81}
chromashift/lab.rs

chromashift/
lab.rs
