geometry-kernel 1.0.3

A Rust-first geometry kernel for deterministic planar geometry operations
Documentation 
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

use crate::buffer::{self, BufferOptions};
use crate::canonicalize::canonicalize_polygon;
use crate::error::Result;
use crate::kernel::{GeometryKernel, KernelBackend};
use crate::overlay;
use crate::precision::PrecisionModel;
use crate::predicates::{polygon_area, segment_intersection, SegmentIntersection};
use crate::types::{Coord, LineString, MultiPolygon, Polygon};

#[derive(Debug, Clone, Copy, PartialEq)]
pub struct PureRustKernel {
    precision: PrecisionModel,
}

impl PureRustKernel {
    pub const fn new(precision: PrecisionModel) -> Self {
        Self { precision }
    }
}

impl Default for PureRustKernel {
    fn default() -> Self {
        Self::new(PrecisionModel::floating())
    }
}

impl GeometryKernel for PureRustKernel {
    fn backend(&self) -> KernelBackend {
        KernelBackend::PureRust
    }

    fn precision(&self) -> PrecisionModel {
        self.precision
    }

    fn canonicalize_polygon(&self, polygon: &Polygon) -> Polygon {
        canonicalize_polygon(polygon, self.precision)
    }

    fn polygon_area(&self, polygon: &Polygon) -> Result<f64> {
        Ok(polygon_area(polygon))
    }

    fn largest_polygon(&self, multi_polygon: &MultiPolygon) -> Result<Option<Polygon>> {
        Ok(multi_polygon
            .polygons
            .iter()
            .max_by(|left, right| {
                polygon_area(left)
                    .partial_cmp(&polygon_area(right))
                    .unwrap_or(std::cmp::Ordering::Equal)
            })
            .cloned())
    }

    fn buffer_polygon(
        &self,
        polygon: &Polygon,
        distance: f64,
        options: BufferOptions,
    ) -> Result<MultiPolygon> {
        buffer::buffer_polygon(polygon, distance, options, self.precision)
    }

    fn line_buffer(
        &self,
        line: &LineString,
        distance: f64,
        options: BufferOptions,
    ) -> Result<MultiPolygon> {
        buffer::line_buffer(line, distance, options, self.precision)
    }

    fn intersection(&self, subject: &MultiPolygon, clip: &MultiPolygon) -> Result<MultiPolygon> {
        overlay::intersection(subject, clip, self.precision)
    }

    fn difference(&self, subject: &MultiPolygon, clip: &MultiPolygon) -> Result<MultiPolygon> {
        overlay::difference(subject, clip, self.precision)
    }

    fn line_polygon_intersections(
        &self,
        line: &LineString,
        polygon: &Polygon,
    ) -> Result<Vec<Coord>> {
        let mut points = Vec::new();

        for (line_start, line_end) in line.segments() {
            for ring in std::iter::once(&polygon.exterior).chain(polygon.holes.iter()) {
                for (ring_start, ring_end) in ring.segments() {
                    match segment_intersection(
                        line_start,
                        line_end,
                        ring_start,
                        ring_end,
                        self.precision,
                    ) {
                        Some(SegmentIntersection::Point(point)) => {
                            push_unique(&mut points, point, self.precision);
                        }
                        Some(SegmentIntersection::Overlap(start, end)) => {
                            push_unique(&mut points, start, self.precision);
                            push_unique(&mut points, end, self.precision);
                        }
                        None => {}
                    }
                }
            }
        }

        Ok(points)
    }
}

fn push_unique(points: &mut Vec<Coord>, point: Coord, precision: PrecisionModel) {
    let point = precision.snap_coord(point);
    if !points
        .iter()
        .any(|existing| precision.same_coord(*existing, point))
    {
        points.push(point);
    }
}