color_gradient/interpolation/
mod.rs

1use std::{collections::BTreeMap, ops::Range};
2
3type Resolution = u16;
4
5/// A interpolator that interpolates values in range.
6#[derive(Clone, Debug)]
7#[cfg_attr(feature = "serde", derive(serde::Serialize, serde::Deserialize))]
8pub struct Interpolator {
9    keys: BTreeMap<Resolution, f32>,
10    lhs: f32,
11    rhs: f32,
12}
13
14impl Default for Interpolator {
15    fn default() -> Self {
16        Self { keys: Default::default(), lhs: 0.0, rhs: 1.0 }
17    }
18}
19
20impl Interpolator {
21    /// Create a new interpolator with the given lower and upper bounds.
22    ///
23    /// # Examples
24    ///
25    /// ```
26    /// # use color_gradient::Interpolator;
27    /// let interpolator = Interpolator::new(0.0, 1.0);
28    /// assert_eq!(interpolator.head(), 0.0);
29    /// assert_eq!(interpolator.tail(), 1.0);
30    /// ```
31    pub fn new(lower: f32, upper: f32) -> Self {
32        Self { keys: Default::default(), lhs: lower, rhs: upper }
33    }
34    /// Get the left most value of the interpolator, if not set, the lower bound is returned.
35    ///
36    /// # Examples
37    ///
38    /// ```
39    /// # use color_gradient::Interpolator;
40    /// let mut interpolator = Interpolator::default();
41    /// assert_eq!(interpolator.head(), 0.0);
42    /// interpolator.insert(0, 0.5);
43    /// assert_eq!(interpolator.head(), 0.5);
44    /// ```
45    pub fn head(&self) -> f32 {
46        self.keys.get(&Resolution::MIN).copied().unwrap_or(self.lhs)
47    }
48    /// Get the right most value of the interpolator, if not set, the upper bound is returned.
49    ///
50    /// # Examples
51    ///
52    /// ```
53    /// # use color_gradient::Interpolator;
54    /// let mut interpolator = Interpolator::default();
55    /// assert_eq!(interpolator.tail(), 1.0);
56    /// interpolator.insert(65535, 0.5);
57    /// assert_eq!(interpolator.tail(), 0.5);
58    /// ```
59    pub fn tail(&self) -> f32 {
60        self.keys.get(&Resolution::MAX).copied().unwrap_or(self.rhs)
61    }
62    /// Insert a new key value pair into the interpolator, overwriting any existing value.
63    ///
64    /// # Examples
65    ///
66    /// ```
67    /// # use color_gradient::{Interpolator};
68    /// let mut gradient = Interpolator::default();
69    /// gradient.insert(0, 0.5);
70    /// gradient.insert(65535, 0.5);
71    /// assert_eq!(gradient.head(), 0.5);
72    /// assert_eq!(gradient.tail(), 0.5);
73    /// ```
74    pub fn insert(&mut self, key: Resolution, value: f32) {
75        self.keys.insert(key, value);
76    }
77    /// Remove a key value pair from the interpolator.
78    ///
79    /// # Examples
80    ///
81    /// ```
82    /// # use color_gradient::{Interpolator};
83    /// let mut gradient = Interpolator::default();
84    /// gradient.insert(0, 0.5);
85    /// assert_eq!(gradient.head(), 0.5);
86    /// gradient.remove(0);
87    /// assert_eq!(gradient.head(), 0.0);
88    /// ```
89    pub fn remove(&mut self, key: Resolution) {
90        self.keys.remove(&key);
91    }
92    /// Clear all key value pairs from the interpolator.
93    ///
94    /// # Examples
95    ///
96    /// ```
97    /// # use color_gradient::{Interpolator};
98    /// ```
99    pub fn clear(&mut self) {
100        self.keys.clear();
101    }
102}
103
104impl Interpolator {
105    pub(crate) fn get_ratio(range: &Range<f32>, value: f32) -> u16 {
106        if value <= range.start {
107            0
108        }
109        else if value >= range.end {
110            65535
111        }
112        else {
113            let ratio = (value - range.start) / (range.end - range.start);
114            (ratio * 65535.0) as u16
115        }
116    }
117
118    /// Get zero-order interpolation, that is, the first number greater than ratio
119    pub fn get_step(&self, key: Resolution) -> f32 {
120        if key == Resolution::MIN {
121            self.head()
122        }
123        else if key == Resolution::MAX {
124            self.tail()
125        }
126        else {
127            self.keys.range(..=key).next_back().map(|(_, v)| *v).unwrap_or(self.lhs)
128        }
129    }
130    /// Get first-order linear interpolation
131    pub fn get_linear(&self, key: Resolution) -> f32 {
132        if key == Resolution::MIN {
133            self.head()
134        }
135        else if key == Resolution::MAX {
136            self.tail()
137        }
138        else {
139            let s1 = (&Resolution::MIN, &self.head());
140            let s2 = (&Resolution::MAX, &self.tail());
141            let (k1, v1) = self.keys.range(..=key).next_back().unwrap_or(s1);
142            let (k2, v2) = self.keys.range(key..).next().unwrap_or(s2);
143            if k1 == k2 { *v1 } else { v1 + (v2 - v1) * (key - k1) as f32 / (k2 - k1) as f32 }
144        }
145    }
146    /// Get bezier interpolation
147    pub fn get_bezier(&self, key: Resolution) -> f32 {
148        if key == Resolution::MIN {
149            self.head()
150        }
151        else if key == Resolution::MAX {
152            self.tail()
153        }
154        else {
155            // bezier 1D
156            let s1 = (&Resolution::MIN, &self.head());
157            let s2 = (&Resolution::MAX, &self.tail());
158            let (k1, v1) = self.keys.range(..=key).next_back().unwrap_or(s1);
159            let (k2, v2) = self.keys.range(key..).next().unwrap_or(s2);
160            let t = (key - k1) as f32 / (k2 - k1) as f32;
161            let v = (1.0 - t) * v1 + t * v2;
162            v
163        }
164    }
165}