fast image operations
Diffstat (limited to 'src/drawing/poly.rs')
| -rw-r--r-- | src/drawing/poly.rs | 74 |
1 files changed, 74 insertions, 0 deletions
diff --git a/src/drawing/poly.rs b/src/drawing/poly.rs new file mode 100644 index 0000000..408aa47 --- /dev/null +++ b/src/drawing/poly.rs @@ -0,0 +1,74 @@ +//! draw polygons +use std::{ + cmp::{max, min}, + ops::{Deref, DerefMut}, +}; + +use crate::Image; + +impl<T: Deref<Target = [u8]> + DerefMut<Target = [u8]>, const CHANNELS: usize> Image<T, CHANNELS> { + /// Draws a filled polygon from a slice of points. Please close your poly. (first == last) + /// Borrowed from [imageproc](https://docs.rs/imageproc/latest/src/imageproc/drawing/polygon.rs.html#31), modified for less allocations. + /// ``` + /// # use fimg::Image; + /// let mut i = Image::alloc(10, 10); + /// i.points(&[(1, 8), (3, 1), (8, 1), (6, 6), (8, 8), (1, 8)], [255]); + /// # assert_eq!(i.buffer(), b"\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\xff\xff\xff\xff\xff\xff\x00\x00\x00\xff\xff\xff\xff\xff\xff\xff\x00\x00\x00\xff\xff\xff\xff\xff\xff\xff\x00\x00\x00\xff\xff\xff\xff\xff\xff\x00\x00\x00\xff\xff\xff\xff\xff\xff\xff\x00\x00\x00\xff\xff\xff\xff\xff\xff\x00\x00\x00\x00\xff\xff\xff\xff\xff\xff\xff\x00\x00\x00\xff\xff\xff\xff\xff\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00"); + /// ``` + pub fn points(&mut self, poly: &[(i32, i32)], c: [u8; CHANNELS]) { + if poly.len() <= 1 { + return; + } + let (mut y_max, mut y_min) = poly[..poly.len() - 1] + .iter() + .fold((i32::MIN, i32::MAX), |(max, min), &(_, y)| { + (y.max(max), y.min(min)) + }); + y_min = max(0, min(y_min, self.height() as i32 - 1)); + y_max = max(0, min(y_max, self.height() as i32 - 1)); + let mut intersections = vec![]; + for y in y_min..=y_max { + for [p0, p1] in poly.array_windows::<2>() { + if p0.1 <= y && p1.1 >= y || p1.1 <= y && p0.1 >= y { + if p0.1 == p1.1 { + intersections.push(p0.0); + intersections.push(p1.0); + } else if p0.1 == y || p1.1 == y { + if p1.1 > y { + intersections.push(p0.0); + } + if p0.1 > y { + intersections.push(p1.0); + } + } else { + let fraction = (y - p0.1) as f32 / (p1.1 - p0.1) as f32; + let inter = fraction.mul_add((p1.0 - p0.0) as f32, p0.0 as f32); + intersections.push(inter.round() as i32); + } + } + } + intersections.sort_unstable(); + // SAFETY: must. + unsafe { crate::assert_unchecked!(intersections.len() % 2 == 0) }; + for &[x, y_] in intersections.array_chunks::<2>() { + let mut from = min(x, self.width() as i32); + let mut to = min(y_, self.width() as i32 - 1); + if from < self.width() as i32 && to >= 0 { + // check bounds + from = max(0, from); + to = max(0, to); + + for x in from..=to { + // SAFETY: bounds are checked + unsafe { self.set_pixel(x as u32, y as u32, c) }; + } + } + } + intersections.clear(); + } + + for &[(x1, y1), (x2, y2)] in poly.array_windows::<2>() { + self.line((x1, y1), (x2, y2), c); + } + } +} |