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
use types::{CoordinateType, MultiPolygon, Polygon};
use algorithm::winding_order::{Winding, WindingOrder};
pub trait Orient<T> {
fn orient(&self, orientation: Direction) -> Self;
}
impl<T> Orient<T> for Polygon<T>
where
T: CoordinateType,
{
fn orient(&self, direction: Direction) -> Polygon<T> {
orient(self, direction)
}
}
impl<T> Orient<T> for MultiPolygon<T>
where
T: CoordinateType,
{
fn orient(&self, direction: Direction) -> MultiPolygon<T> {
MultiPolygon(self.0.iter().map(|poly| poly.orient(direction)).collect())
}
}
#[derive(Copy, Clone, Debug)]
pub enum Direction {
Default,
Reversed,
}
fn orient<T>(poly: &Polygon<T>, direction: Direction) -> Polygon<T>
where
T: CoordinateType,
{
let interiors = poly.interiors.iter().map(|l| l.clone_to_winding_order(
match direction {
Direction::Default => WindingOrder::Clockwise,
Direction::Reversed => WindingOrder::CounterClockwise,
})).collect();
let ext_ring = poly.exterior.clone_to_winding_order(match direction {
Direction::Default => WindingOrder::CounterClockwise,
Direction::Reversed => WindingOrder::Clockwise,
});
Polygon::new(ext_ring, interiors)
}
#[cfg(test)]
mod test {
use types::{LineString, Point, Polygon};
use super::*;
#[test]
fn test_polygon_orientation() {
let points_ext_raw = vec![(1.0, 0.0), (0.0, 1.0), (1.0, 2.0), (2.0, 1.0), (1.0, 0.0)];
let points_int_raw = vec![(1.0, 0.5), (1.5, 1.0), (1.0, 1.5), (0.5, 1.0), (1.0, 0.5)];
let points_ext = points_ext_raw
.iter()
.map(|e| Point::new(e.0, e.1))
.collect::<Vec<_>>();
let points_int = points_int_raw
.iter()
.map(|e| Point::new(e.0, e.1))
.collect::<Vec<_>>();
let poly1 = Polygon::new(LineString(points_ext), vec![LineString(points_int)]);
let oriented_ext = vec![(1.0, 0.0), (2.0, 1.0), (1.0, 2.0), (0.0, 1.0), (1.0, 0.0)];
let oriented_ext_ls = LineString(
oriented_ext
.iter()
.map(|e| Point::new(e.0, e.1))
.collect::<Vec<_>>(),
);
let oriented_int_raw = vec![(1.0, 0.5), (0.5, 1.0), (1.0, 1.5), (1.5, 1.0), (1.0, 0.5)];
let oriented_int_ls = LineString(
oriented_int_raw
.iter()
.map(|e| Point::new(e.0, e.1))
.collect::<Vec<_>>(),
);
let oriented = orient(&poly1, Direction::Default);
assert_eq!(oriented.exterior.0, oriented_ext_ls.0);
assert_eq!(oriented.interiors[0].0, oriented_int_ls.0);
}
}