zerometry 0.3.0

Make some zerocopy operations on serialized geo data
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

use insta::assert_compact_debug_snapshot;

use crate::{Coord, RelationBetweenShapes, Zerometry, Zoint};

#[test]
fn test_mono_multipolygon_contains_points() {
    // This test highlights a bug I found on cellulite where a point was considered contained in an h3 cell even though it clearly is not.
    let buffer = [6.0197316417968105, 49.63676497357687];
    let point = Zoint::new(unsafe { Coord::from_slice(&buffer) });
    println!("point: {}", print_geojson(&point.into()));

    let wrong_multipolygon = geo_types::MultiPolygon::new(vec![geo_types::Polygon::new(
        geo_types::LineString::from(vec![
            geo_types::Point::new(-6.436337296790293, 55.37739041554851),
            geo_types::Point::new(-4.889760342933786, 51.22372845966178),
            geo_types::Point::new(1.188509553443464, 49.47027919866874),
            geo_types::Point::new(3.6300086390995316, 50.610463312569514),
            geo_types::Point::new(6.259687055981991, 51.96477015603749),
            geo_types::Point::new(5.52364654929031, 55.70676846515227),
            geo_types::Point::new(-0.9315871635106105, 57.689497374592854),
            geo_types::Point::new(-6.436337296790293, 55.37739041554851),
        ]),
        Vec::new(),
    )]);
    let mut wrong_buffer = Vec::new();
    Zerometry::write_from_geometry(
        &mut wrong_buffer,
        &geo_types::Geometry::MultiPolygon(wrong_multipolygon.clone()),
    )
    .unwrap();
    let wrong_multipolygon = unsafe { Zerometry::from_bytes(&wrong_buffer).unwrap() };
    let wrong = wrong_multipolygon.all_relation(&point);
    assert_compact_debug_snapshot!(wrong, @"OutputRelation { contains: Some(false), strict_contains: Some(false), contained: Some(false), strict_contained: Some(false), intersect: Some(false), disjoint: Some(true) }");

    let right_multipolygon = geo_types::MultiPolygon::new(vec![geo_types::Polygon::new(
        geo_types::LineString::from(vec![
            geo_types::Point::new(7.509948481928903, 43.78660935394501),
            geo_types::Point::new(12.677317810359865, 46.406957457982266),
            geo_types::Point::new(12.345747364400054, 50.55427508726938),
            geo_types::Point::new(6.259687055981991, 51.96477015603749),
            geo_types::Point::new(3.630008639099554, 50.610463312569486),
            geo_types::Point::new(1.188509553443464, 49.47027919866874),
            geo_types::Point::new(2.026568965384611, 45.18424868970644),
            geo_types::Point::new(7.509948481928903, 43.78660935394501),
        ]),
        Vec::new(),
    )]);

    let mut right_buffer = Vec::new();
    Zerometry::write_from_geometry(
        &mut right_buffer,
        &geo_types::Geometry::MultiPolygon(right_multipolygon.clone()),
    )
    .unwrap();
    let right_multipolygon = unsafe { Zerometry::from_bytes(&right_buffer).unwrap() };
    println!("right_multipolygon: {}", print_geojson(&right_multipolygon));

    let right = right_multipolygon.all_relation(&point);
    assert_compact_debug_snapshot!(right, @"OutputRelation { contains: Some(true), strict_contains: Some(true), contained: Some(false), strict_contained: Some(false), intersect: Some(false), disjoint: Some(false) }");
}

fn print_geojson(geometry: &Zerometry) -> String {
    geojson::GeoJson::Geometry(geojson::Geometry::new(geojson::Value::from(
        &geometry.to_geo(),
    )))
    .to_string_pretty()
    .unwrap()
}