rawshift_image/data/
icc.rs

1//! ICC profile constants.
2//!
3//! Reference values for ICC color management, including the Profile
4//! Connection Space (PCS) white point and standard color space primaries.
5
6/// ICC Profile Connection Space (PCS) white point: D50.
7///
8/// XYZ tristimulus values (Y=1 normalized) per ICC specification.
9pub const PCS_WHITEPOINT: [f64; 3] = [0.9505, 1.0000, 1.0890];
10
11/// sRGB color space primaries in CIE 1931 xy chromaticity.
12///
13/// Format: [[Rx, Ry], [Gx, Gy], [Bx, By]]
14pub const SRGB_PRIMARIES: [[f64; 2]; 3] = [
15    [0.6400, 0.3300], // Red
16    [0.3000, 0.6000], // Green
17    [0.1500, 0.0600], // Blue
18];
19
20/// sRGB to XYZ (D65) matrix (3x3, row-major).
21pub const SRGB_TO_XYZ_D65: [f64; 9] = [
22    0.4124564, 0.3575761, 0.1804375, 0.2126729, 0.7151522, 0.0721750, 0.0193339, 0.1191920,
23    0.9503041,
24];
25
26/// XYZ (D65) to sRGB matrix (3x3, row-major).
27pub const XYZ_D65_TO_SRGB: [f64; 9] = [
28    3.2404542, -1.5371385, -0.4985314, -0.9692660, 1.8760108, 0.0415560, 0.0556434, -0.2040259,
29    1.0572252,
30];
31
32/// Adobe RGB (1998) primaries in CIE 1931 xy chromaticity.
33pub const ADOBE_RGB_PRIMARIES: [[f64; 2]; 3] = [
34    [0.6400, 0.3300], // Red
35    [0.2100, 0.7100], // Green
36    [0.1500, 0.0600], // Blue
37];
38
39/// Bradford chromatic adaptation matrix: D65 to D50.
40///
41/// Used to adapt color data from D65 (common camera reference)
42/// to D50 (ICC PCS).
43pub const BRADFORD_D65_TO_D50: [f64; 9] = [
44    1.0478112, 0.0228866, -0.0501270, 0.0295424, 0.9904844, -0.0170491, -0.0092345, 0.0150436,
45    0.7521316,
46];
47
48/// Bradford chromatic adaptation matrix: D50 to D65.
49pub const BRADFORD_D50_TO_D65: [f64; 9] = [
50    0.9555766, -0.0230393, 0.0631636, -0.0282895, 1.0099416, 0.0210077, 0.0122982, -0.0204830,
51    1.3299098,
52];
53
54#[cfg(test)]
55mod tests {
56    use super::*;
57
58    #[test]
59    fn test_srgb_matrix_identity_white() {
60        // sRGB matrix applied to D65 white point should give ~[1,1,1] in sRGB
61        let xyz = [0.9505, 1.0000, 1.0890]; // D50, not D65 — just verify structure
62        let r =
63            XYZ_D65_TO_SRGB[0] * xyz[0] + XYZ_D65_TO_SRGB[1] * xyz[1] + XYZ_D65_TO_SRGB[2] * xyz[2];
64        // Not exactly 1.0 since xyz here is D50, but should be a reasonable value
65        assert!(r > 0.0 && r < 2.0);
66    }
67
68    #[test]
69    fn test_bradford_roundtrip() {
70        // D65->D50->D65 should approximate identity
71        let input = [0.5, 0.6, 0.7];
72        // Forward: D65 -> D50
73        let mid = [
74            BRADFORD_D65_TO_D50[0] * input[0]
75                + BRADFORD_D65_TO_D50[1] * input[1]
76                + BRADFORD_D65_TO_D50[2] * input[2],
77            BRADFORD_D65_TO_D50[3] * input[0]
78                + BRADFORD_D65_TO_D50[4] * input[1]
79                + BRADFORD_D65_TO_D50[5] * input[2],
80            BRADFORD_D65_TO_D50[6] * input[0]
81                + BRADFORD_D65_TO_D50[7] * input[1]
82                + BRADFORD_D65_TO_D50[8] * input[2],
83        ];
84        // Inverse: D50 -> D65
85        let output = [
86            BRADFORD_D50_TO_D65[0] * mid[0]
87                + BRADFORD_D50_TO_D65[1] * mid[1]
88                + BRADFORD_D50_TO_D65[2] * mid[2],
89            BRADFORD_D50_TO_D65[3] * mid[0]
90                + BRADFORD_D50_TO_D65[4] * mid[1]
91                + BRADFORD_D50_TO_D65[5] * mid[2],
92            BRADFORD_D50_TO_D65[6] * mid[0]
93                + BRADFORD_D50_TO_D65[7] * mid[1]
94                + BRADFORD_D50_TO_D65[8] * mid[2],
95        ];
96        for i in 0..3 {
97            assert!(
98                (output[i] - input[i]).abs() < 1e-4,
99                "roundtrip failed at index {i}: {:.6} != {:.6}",
100                output[i],
101                input[i]
102            );
103        }
104    }
105}
rawshift_image/data/icc.rs

rawshift_image/data/
icc.rs
