egml-core 0.0.2-alpha.3

Core primitives and operations for processing GML 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
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

use crate::Error;
use crate::model::geometry::primitives::{
    AbstractSurface, AsAbstractSurface, AsAbstractSurfaceMut, SurfacePatchArrayProperty,
    TriangulatedSurface,
};
use crate::model::geometry::{DirectPosition, Envelope};
use nalgebra::Isometry3;

/// A 2-D geometry composed of one or more surface patches.
///
/// Corresponds to `gml:Surface` in [OGC 07-036 ยง10.5.10](https://docs.ogc.org/is/07-036/07-036.pdf).  Patches are stored as
/// a [`SurfacePatchArrayProperty`] and may be of mixed kinds (polygons, triangles).
#[derive(Debug, Clone, PartialEq)]
pub struct Surface {
    pub(crate) abstract_surface: AbstractSurface,
    patches: SurfacePatchArrayProperty,
}

impl Surface {
    /// Creates a new `Surface` from its abstract base and patch array.
    pub fn new(abstract_surface: AbstractSurface, patches: SurfacePatchArrayProperty) -> Self {
        Surface {
            abstract_surface,
            patches,
        }
    }

    pub(crate) fn into_patches(self) -> SurfacePatchArrayProperty {
        self.patches
    }

    /// Decomposes this surface into a [`TriangulatedSurface`] by triangulating each patch.
    ///
    /// # Errors
    ///
    /// Returns [`Error::TriangulationFailed`] if any patch cannot be triangulated.
    pub fn triangulate(&self) -> Result<TriangulatedSurface, Error> {
        let a = self.patches.patches.iter().map(|patch| patch.triangulate());
        let b = a.collect::<Result<Vec<TriangulatedSurface>, Error>>()?;
        TriangulatedSurface::from_triangulated_surfaces(b)
    }

    /// Returns the union of the bounding boxes of all patches.
    pub fn compute_envelope(&self) -> Envelope {
        self.patches.compute_envelope()
    }

    pub fn apply_transform(&mut self, m: &Isometry3<f64>) {
        self.patches.apply_transform(m)
    }

    pub fn points(&self) -> Vec<&DirectPosition> {
        todo!("needs to be implemented")
    }
}

/// Object-safe read accessor for [`Surface`] data.
pub trait AsSurface: AsAbstractSurface {
    /// Returns a reference to the underlying [`Surface`].
    fn surface(&self) -> &Surface;

    /// Returns the patch array of this surface.
    fn patches(&self) -> &SurfacePatchArrayProperty {
        &self.surface().patches
    }

    fn patches_len(&self) -> usize {
        self.patches().patches_len()
    }
}

/// Mutable companion to [`AsSurface`].
pub trait AsSurfaceMut: AsSurface + AsAbstractSurfaceMut {
    /// Returns a mutable reference to the underlying [`Surface`].
    fn surface_mut(&mut self) -> &mut Surface;

    fn patches_mut(&mut self) -> &mut SurfacePatchArrayProperty {
        &mut self.surface_mut().patches
    }
}

impl AsSurface for Surface {
    fn surface(&self) -> &Surface {
        self
    }
}

impl AsSurfaceMut for Surface {
    fn surface_mut(&mut self) -> &mut Surface {
        self
    }
}

#[doc(hidden)]
#[macro_export]
macro_rules! impl_surface_traits {
    ($type:ty) => {
        $crate::impl_abstract_surface_traits!($type);

        impl $crate::model::geometry::primitives::AsAbstractSurface for $type {
            fn abstract_surface(&self) -> &$crate::model::geometry::primitives::AbstractSurface {
                use $crate::model::geometry::primitives::AsSurface;
                &self.surface().abstract_surface
            }
        }

        impl $crate::model::geometry::primitives::AsAbstractSurfaceMut for $type {
            fn abstract_surface_mut(
                &mut self,
            ) -> &mut $crate::model::geometry::primitives::AbstractSurface {
                use $crate::model::geometry::primitives::AsSurfaceMut;
                &mut self.surface_mut().abstract_surface
            }
        }
    };
}

impl_surface_traits!(Surface);