1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160
use std::f64;
use crate::{Point, Result, Transform, Vector};
/// Minimum and maximum bounds in three dimensions.
#[derive(Clone, Copy, Debug, PartialEq)]
pub struct Bounds {
/// The minimum values.
pub min: Vector<f64>,
/// The maximum values.
pub max: Vector<f64>,
}
impl Bounds {
/// Grows the bounds to encompass this point in xyz space.
///
/// # Examples
///
/// ```
/// # use las::{Bounds, Point};
/// let point = Point { x: 1., y: 2., z: 3., ..Default::default() };
/// let mut bounds = Bounds { ..Default::default() };
/// bounds.grow(&point);
/// assert_eq!(1., bounds.min.x);
/// ```
pub fn grow(&mut self, point: &Point) {
if point.x < self.min.x {
self.min.x = point.x;
}
if point.y < self.min.y {
self.min.y = point.y;
}
if point.z < self.min.z {
self.min.z = point.z;
}
if point.x > self.max.x {
self.max.x = point.x;
}
if point.y > self.max.y {
self.max.y = point.y;
}
if point.z > self.max.z {
self.max.z = point.z;
}
}
/// Transform the bounds to be compatible with the chosen transform. Otherwise, points may lay outside of the bounding box due to floating-point issues.
///
/// # Example
///
/// ```
/// use las::{Bounds, Transform, Vector};
///
/// let bounds = Bounds {
/// min: Vector {
/// x: -2.7868618965148926,
/// y: -0.9322229027748108,
/// z: -5.8063459396362305,
/// },
/// max: Vector {
/// x: 0.6091402173042297,
/// y: 1.5428568124771118,
/// z: -0.09441471844911575,
/// },
/// };
///
/// // Currently, the default scale is 0.001.
/// let new_bounds = bounds.adapt(&Default::default()).unwrap();
/// assert_eq!(new_bounds.max.z, -0.094);
/// ```
pub fn adapt(&self, transform: &Vector<Transform>) -> Result<Self> {
fn inner_convert(value: f64, transform: &Transform) -> Result<f64> {
// During saving, an instance with +-inf is saved. We must consider for this corner case.
if value.is_infinite() {
return Ok(value);
}
Ok(transform.direct(transform.inverse(value)?))
}
Ok(Self {
min: Vector {
x: inner_convert(self.min.x, &transform.x)?,
y: inner_convert(self.min.y, &transform.y)?,
z: inner_convert(self.min.z, &transform.z)?,
},
max: Vector {
x: inner_convert(self.max.x, &transform.x)?,
y: inner_convert(self.max.y, &transform.y)?,
z: inner_convert(self.max.z, &transform.z)?,
},
})
}
}
impl Default for Bounds {
fn default() -> Bounds {
Bounds {
min: Vector {
x: f64::INFINITY,
y: f64::INFINITY,
z: f64::INFINITY,
},
max: Vector {
x: f64::NEG_INFINITY,
y: f64::NEG_INFINITY,
z: f64::NEG_INFINITY,
},
}
}
}
#[cfg(test)]
mod tests {
use super::*;
use crate::Point;
#[test]
fn grow() {
let mut bounds = Bounds {
..Default::default()
};
bounds.grow(&Point {
x: 1.,
y: 2.,
z: 3.,
..Default::default()
});
assert_eq!(1., bounds.min.x);
assert_eq!(1., bounds.max.x);
assert_eq!(2., bounds.min.y);
assert_eq!(2., bounds.max.y);
assert_eq!(3., bounds.min.z);
assert_eq!(3., bounds.max.z);
bounds.grow(&Point {
x: 0.,
y: 1.,
z: 2.,
..Default::default()
});
assert_eq!(0., bounds.min.x);
assert_eq!(1., bounds.max.x);
assert_eq!(1., bounds.min.y);
assert_eq!(2., bounds.max.y);
assert_eq!(2., bounds.min.z);
assert_eq!(3., bounds.max.z);
bounds.grow(&Point {
x: 2.,
y: 3.,
z: 4.,
..Default::default()
});
assert_eq!(0., bounds.min.x);
assert_eq!(2., bounds.max.x);
assert_eq!(1., bounds.min.y);
assert_eq!(3., bounds.max.y);
assert_eq!(2., bounds.min.z);
assert_eq!(4., bounds.max.z);
}
}