geo-repair 0.13.0

Fix invalid GIS geometries in parallel using rewritten GEOS/JTS repair module
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

use super::prep::PreparedLines;
use crate::core::MakeValidError;
use geo::Coord;
use rustc_hash::FxHashMap;
use spade::handles::FixedVertexHandle;
use spade::{ConstrainedDelaunayTriangulation, Triangulation};

pub(crate) fn build(
    prepared: &PreparedLines,
) -> Result<ConstrainedDelaunayTriangulation<Coord<f64>>, MakeValidError> {
    let mut cdt = ConstrainedDelaunayTriangulation::<Coord<f64>>::new();
    let mut pos_to_handle: FxHashMap<(u64, u64), FixedVertexHandle> = FxHashMap::default();

    for &line in &prepared.lines {
        let start_key = (line.start.x.to_bits(), line.start.y.to_bits());
        let start = *pos_to_handle
            .entry(start_key)
            .or_insert_with(|| cdt.insert(line.start).expect("vertex insertion failed"));

        let end_key = (line.end.x.to_bits(), line.end.y.to_bits());
        let end = *pos_to_handle
            .entry(end_key)
            .or_insert_with(|| cdt.insert(line.end).expect("vertex insertion failed"));

        if start != end && cdt.can_add_constraint(start, end) {
            cdt.add_constraint(start, end);
        }
    }

    Ok(cdt)
}