cazan-common 0.1.1

The common library to `cazan-cli` and `cazanw`
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

use crate::geometry::{distance, Point, Polygon};

pub fn rdp(polygon: &Polygon, epsilon: f64) -> Polygon {
    let mut dmax = 0.0;
    let mut index = 0;
    let end = polygon.len() - 1;

    for i in 1..end {
        let d = perpendicular_distance(polygon[i], polygon[0], polygon[end]);
        if d > dmax {
            index = i;
            dmax = d;
        }
    }

    if dmax > epsilon {
        let mut results = rdp(&Vec::from(&polygon[..=index]), epsilon);
        results.pop();
        results.extend(rdp(&Vec::from(&polygon[index..]), epsilon));
        results
    } else {
        vec![polygon[0], polygon[end]]
    }
}

fn perpendicular_distance(p: Point, a: Point, b: Point) -> f64 {
    let num = ((b.y as i32 - a.y as i32) * p.x as i32) - ((b.x as i32 - a.x as i32) * p.y as i32)
        + (b.x as i32 * a.y as i32)
        - (b.y as i32 * a.x as i32);
    let num = num.abs() as f64;
    let den = distance(a, b) as f64;
    num / den
}