h3o 0.9.4

A Rust implementation of the H3 geospatial indexing system.
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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169

use super::*;
use geo::{line_string, polygon};

#[test]
fn test_simple() {
    // ┏━┓
    // ┗━┛
    #[rustfmt::skip]
        let rings = vec![
            line_string![(x: 1., y: 1.), (x: 1., y: 3.), (x: 3., y: 3.), (x: 3., y: 1.), (x: 1., y: 1.)],
        ];
    // Expect a single polygon without hole.
    #[rustfmt::skip]
        let expected = MultiPolygon(vec![
            polygon![(x: 1., y: 1.), (x: 1., y: 3.), (x: 3., y: 3.), (x: 3., y: 1.), (x: 1., y: 1.)],
        ]);

    let res = MultiPolygon::from(RingHierarchy::new(rings));

    assert_eq!(res, expected);
}

#[test]
fn test_with_one_hole() {
    // ┏━━━┓
    // ┃┏━┓┃
    // ┃┗━┛┃
    // ┗━━━┛
    #[rustfmt::skip]
        let rings = vec![
            line_string![(x: 2., y: 2.), (x: 2., y: 4.), (x: 4., y: 4.), (x: 4., y: 2.), (x: 2., y: 2.)],
            line_string![(x: 1., y: 1.), (x: 1., y: 5.), (x: 5., y: 5.), (x: 5., y: 1.), (x: 1., y: 1.)],
        ];
    // Expect a single polygon with a hole.
    #[rustfmt::skip]
        let expected = MultiPolygon(vec![polygon![
            exterior:   [(x: 1., y: 1.), (x: 1., y: 5.), (x: 5., y: 5.), (x: 5., y: 1.), (x: 1., y: 1.)],
            interiors: [[(x: 2., y: 2.), (x: 2., y: 4.), (x: 4., y: 4.), (x: 4., y: 2.), (x: 2., y: 2.)]]
        ]]);

    let res = MultiPolygon::from(RingHierarchy::new(rings));

    assert_eq!(res, expected);
}

#[test]
fn test_with_multiple_holes() {
    // ┏━━━━━━━┓
    // ┃┏━┓ ┏━┓┃
    // ┃┗━┛ ┗━┛┃
    // ┗━━━━━━━┛
    #[rustfmt::skip]
        let rings = vec![
            line_string![(x: 3., y: 4.), (x: 3., y: 6.), (x: 4., y: 6.), (x: 4., y: 4.), (x: 3., y: 4.)],
            line_string![(x: 1., y: 3.), (x: 1., y: 7.), (x: 9., y: 7.), (x: 9., y: 3.), (x: 1., y: 3.)],
            line_string![(x: 6., y: 4.), (x: 6., y: 6.), (x: 7., y: 6.), (x: 7., y: 4.), (x: 6., y: 4.)],
        ];
    // Expect a single polygon with two holes.
    #[rustfmt::skip]
        let expected = MultiPolygon(vec![polygon![
            exterior: [(x: 1., y: 3.), (x: 1., y: 7.), (x: 9., y: 7.), (x: 9., y: 3.), (x: 1., y: 3.)],
            interiors: [
                [(x: 3., y: 4.), (x: 3., y: 6.), (x: 4., y: 6.), (x: 4., y: 4.), (x: 3., y: 4.)],
                [(x: 6., y: 4.), (x: 6., y: 6.), (x: 7., y: 6.), (x: 7., y: 4.), (x: 6., y: 4.)],
            ]
        ]]);
    let res = MultiPolygon::from(RingHierarchy::new(rings));

    assert_eq!(res, expected);
}

#[test]
fn test_with_multiple_outers() {
    // ┏━┓ ┏━┓
    // ┗━┛ ┗━┛
    #[rustfmt::skip]
        let rings = vec![
            line_string![(x: 1., y: 2.), (x: 1., y: 4.), (x: 2., y: 4.), (x: 2., y: 2.), (x: 1., y: 2.)],
            line_string![(x: 4., y: 2.), (x: 4., y: 4.), (x: 5., y: 4.), (x: 5., y: 2.), (x: 4., y: 2.)],
        ];
    // Expect a two polygons without holes.
    #[rustfmt::skip]
        let expected = MultiPolygon(vec![
            polygon![(x: 1., y: 2.), (x: 1., y: 4.), (x: 2., y: 4.), (x: 2., y: 2.), (x: 1., y: 2.)],
            polygon![(x: 4., y: 2.), (x: 4., y: 4.), (x: 5., y: 4.), (x: 5., y: 2.), (x: 4., y: 2.)],
        ]);
    let res = MultiPolygon::from(RingHierarchy::new(rings));

    assert_eq!(res, expected);
}

#[test]
fn test_nested() {
    // ┏━━━━━┓
    // ┃┏━━━┓┃
    // ┃┃┏━┓┃┃
    // ┃┃┗━┛┃┃
    // ┃┗━━━┛┃
    // ┗━━━━━┛
    #[rustfmt::skip]
        let rings = vec![
            line_string![(x: 3., y: 3.), (x: 3., y: 5.), (x: 5., y: 5.), (x: 5., y: 3.), (x: 3., y: 3.)],
            line_string![(x: 2., y: 2.), (x: 2., y: 6.), (x: 6., y: 6.), (x: 6., y: 2.), (x: 2., y: 2.)],
            line_string![(x: 1., y: 1.), (x: 1., y: 7.), (x: 7., y: 7.), (x: 7., y: 1.), (x: 1., y: 1.)],
        ];
    // Expect a polygon with a hole (outermost) + a polygon without hole
    // (innermost one)
    #[rustfmt::skip]
        let expected = MultiPolygon(vec![
            polygon![
                exterior:   [(x: 1., y: 1.), (x: 1., y: 7.), (x: 7., y: 7.), (x: 7., y: 1.), (x: 1., y: 1.)],
                interiors: [[(x: 2., y: 2.), (x: 2., y: 6.), (x: 6., y: 6.), (x: 6., y: 2.), (x: 2., y: 2.)]]
            ],
            polygon![(x: 3., y: 3.), (x: 3., y: 5.), (x: 5., y: 5.), (x: 5., y: 3.), (x: 3., y: 3.)],
        ]);

    let res = MultiPolygon::from(RingHierarchy::new(rings));

    assert_eq!(res, expected);
}

#[test]
fn test_gamut() {
    // ┏━━━━━━━┓
    // ┃┏━━━━━┓┃  ┏━━━━━━━┓
    // ┃┃┏━━━┓┃┃  ┃┏━┓ ┏━┓┃
    // ┃┃┃┏━┓┃┃┃  ┃┗━┛ ┗━┛┃
    // ┃┃┃┗━┛┃┃┃  ┗━━━━━━━┛
    // ┃┃┗━━━┛┃┃
    // ┃┗━━━━━┛┃     ┏━┓
    // ┗━━━━━━━┛     ┗━┛
    #[rustfmt::skip]
        let rings = vec![
            // Nested.
            line_string![(x: 1., y: 1.), (x: 1., y: 9.), (x: 9., y: 9.), (x: 9., y: 1.), (x: 1., y: 1.)],
            line_string![(x: 2., y: 2.), (x: 2., y: 8.), (x: 8., y: 8.), (x: 8., y: 2.), (x: 2., y: 2.)],
            line_string![(x: 3., y: 3.), (x: 3., y: 7.), (x: 7., y: 7.), (x: 7., y: 3.), (x: 3., y: 3.)],
            line_string![(x: 4., y: 4.), (x: 4., y: 6.), (x: 6., y: 6.), (x: 6., y: 4.), (x: 4., y: 4.)],
            // Multiple holes.
            line_string![(x: 13., y: 14.), (x: 13., y: 16.), (x: 14., y: 16.), (x: 14., y: 14.), (x: 13., y: 14.)],
            line_string![(x: 11., y: 13.), (x: 11., y: 17.), (x: 19., y: 17.), (x: 19., y: 13.), (x: 11., y: 13.)],
            line_string![(x: 16., y: 14.), (x: 16., y: 16.), (x: 17., y: 16.), (x: 17., y: 14.), (x: 16., y: 14.)],
            // Simple
            line_string![(x: 19., y: 19.), (x: 19., y: 21.), (x: 21., y: 21.), (x: 21., y: 19.), (x: 19., y: 19.)],
        ];
    #[rustfmt::skip]
        let expected = MultiPolygon(vec![
            polygon![
                exterior:   [(x: 1., y: 1.), (x: 1., y: 9.), (x: 9., y: 9.), (x: 9., y: 1.), (x: 1., y: 1.)],
                interiors: [[(x: 2., y: 2.), (x: 2., y: 8.), (x: 8., y: 8.), (x: 8., y: 2.), (x: 2., y: 2.)]]
            ],
            polygon![
                exterior: [(x: 11., y: 13.), (x: 11., y: 17.), (x: 19., y: 17.), (x: 19., y: 13.), (x: 11., y: 13.)],
                interiors: [
                    [(x: 13., y: 14.), (x: 13., y: 16.), (x: 14., y: 16.), (x: 14., y: 14.), (x: 13., y: 14.)],
                    [(x: 16., y: 14.), (x: 16., y: 16.), (x: 17., y: 16.), (x: 17., y: 14.), (x: 16., y: 14.)],
                ]
            ],
            polygon![(x: 19., y: 19.), (x: 19., y: 21.), (x: 21., y: 21.), (x: 21., y: 19.), (x: 19., y: 19.)],
            polygon![
                exterior:   [(x: 3., y: 3.), (x: 3., y: 7.), (x: 7., y: 7.), (x: 7., y: 3.), (x: 3., y: 3.)],
                interiors: [[(x: 4., y: 4.), (x: 4., y: 6.), (x: 6., y: 6.), (x: 6., y: 4.), (x: 4., y: 4.)]]
            ],
        ]);

    let res = MultiPolygon::from(RingHierarchy::new(rings));

    assert_eq!(res, expected);
}