geos 11.1.1

Rust bindings for GEOS C API
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

use crate::error::Error;
use crate::{ConstGeometry, Geom, Geometry as GGeometry};
use geo_types::Geometry;
use wkt;
use wkt::TryFromWkt;

use std::convert::TryFrom;

fn to_geo<T: Geom>(other: &T) -> Result<Geometry<f64>, Error> {
    // This is a first draft, it's very inefficient, we use wkt as a pivot format to
    // translate the geometry.
    // We should at least use wkb, or even better implement a direct translation
    let wkt_str = other.to_wkt()?;
    geo_types::Geometry::try_from_wkt_str(&wkt_str)
        .map_err(|e| Error::ConversionError(format!("impossible to read wkt: {e}")))
}

impl TryFrom<GGeometry> for Geometry<f64> {
    type Error = Error;

    fn try_from(other: GGeometry) -> Result<Self, Self::Error> {
        to_geo(&other)
    }
}

impl TryFrom<&GGeometry> for Geometry<f64> {
    type Error = Error;

    fn try_from(other: &GGeometry) -> Result<Self, Self::Error> {
        to_geo(other)
    }
}

impl TryFrom<ConstGeometry<'_>> for Geometry<f64> {
    type Error = Error;

    fn try_from(other: ConstGeometry<'_>) -> Result<Self, Self::Error> {
        to_geo(&other)
    }
}

impl TryFrom<&ConstGeometry<'_>> for Geometry<f64> {
    type Error = Error;

    fn try_from(other: &ConstGeometry<'_>) -> Result<Self, Self::Error> {
        to_geo(other)
    }
}

#[cfg(test)]
mod test {
    use crate::Geometry as GGeometry;
    use geo_types::{Coord, Geometry, LineString, MultiPoint, MultiPolygon, Point, Polygon};
    use std::convert::TryInto;

    fn coords(tuples: Vec<(f64, f64)>) -> Vec<Coord<f64>> {
        tuples.into_iter().map(Coord::from).collect()
    }

    #[test]
    fn geom_to_geo_polygon() {
        let poly = "MULTIPOLYGON(((0 0, 0 1, 1 1, 1 0, 0 0)))";
        let poly = GGeometry::new_from_wkt(poly).unwrap();

        let geo_polygon: Geometry<f64> = (&poly).try_into().unwrap();

        let exterior = LineString(coords(vec![
            (0., 0.),
            (0., 1.),
            (1., 1.),
            (1., 0.),
            (0., 0.),
        ]));
        let expected_poly = MultiPolygon(vec![Polygon::new(exterior, vec![])]);
        let expected: Geometry<_> = expected_poly.into();
        assert_eq!(expected, geo_polygon);
        // This check is to enforce that `TryFrom` is implemented for both reference and value.
        assert_eq!(expected, poly.try_into().unwrap());
    }

    #[test]
    fn geom_to_geo_multipoint() {
        let mp = "MULTIPOINT (0 0, 1 1)";
        let mp = GGeometry::new_from_wkt(mp).unwrap();

        let geo_multipoint: Geometry<f64> = (&mp).try_into().unwrap();

        let expected_multipoint = MultiPoint(vec![
            Point(Coord::from((0., 0.))),
            Point(Coord::from((1., 1.))),
        ]);
        let expected: Geometry<_> = expected_multipoint.into();
        assert_eq!(expected, geo_multipoint);
        // This check is to enforce that `TryFrom` is implemented for both reference and value.
        assert_eq!(expected, mp.try_into().unwrap());
    }
}