etiles-core 0.0.1-alpha.9

Core primitives and operations for processing 3D Tiles 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

use ecoord::AxisAlignedBoundingCube;
use eproj::Coordinate3;
use nalgebra::{Point3, Vector3};

#[derive(Debug, Clone, Copy)]
pub struct BoundingCube {
    center: Point3<f64>,
    width: f64,
}

impl BoundingCube {
    pub fn new(center: Point3<f64>, width: f64) -> Self {
        Self { center, width }
    }

    pub fn from_axis_aligned_bounding_cube(bounding_cube: &AxisAlignedBoundingCube) -> Self {
        Self {
            center: bounding_cube.center(),
            width: bounding_cube.edge_length(),
        }
    }

    pub fn get_lower_bound(&self) -> Point3<f64> {
        self.center - Vector3::new(self.width, self.width, self.width)
    }

    pub fn get_upper_bound(&self) -> Point3<f64> {
        self.center + Vector3::new(self.width, self.width, self.width)
    }

    pub fn center_vector(&self) -> Vector3<f64> {
        self.center.coords
    }

    pub fn get_octant(&self, x_half: bool, y_half: bool, z_half: bool) -> BoundingCube {
        let octant_width = self.width / 2.0;
        let x_sign = if x_half { 1.0 } else { -1.0 };
        let y_sign = if y_half { 1.0 } else { -1.0 };
        let z_sign = if z_half { 1.0 } else { -1.0 };

        let octant_center = self.center
            + Vector3::new(
                octant_width * x_sign,
                octant_width * y_sign,
                octant_width * z_sign,
            );

        Self::new(octant_center, octant_width)
    }

    pub fn x_axis(&self) -> Vector3<f64> {
        Vector3::new(self.width / 2.0, 0.0, 0.0)
    }

    pub fn y_axis(&self) -> Vector3<f64> {
        Vector3::new(0.0, self.width / 2.0, 0.0)
    }

    pub fn z_axis(&self) -> Vector3<f64> {
        Vector3::new(0.0, 0.0, self.width / 2.0)
    }

    pub fn bounding_array(&self) -> [f64; 12] {
        let center_vec = self.center_vector();
        let x_axis_vec = self.x_axis();
        let y_axis_vec = self.y_axis();
        let z_axis_vec = self.z_axis();

        [
            center_vec.x,
            center_vec.y,
            center_vec.z,
            x_axis_vec.x,
            x_axis_vec.y,
            x_axis_vec.z,
            y_axis_vec.x,
            y_axis_vec.y,
            y_axis_vec.z,
            z_axis_vec.x,
            z_axis_vec.y,
            z_axis_vec.z,
        ]
    }
}

#[derive(Debug, Clone, Copy)]
pub struct BoundingRegion {
    south_west_min_height: Coordinate3,
    north_east_max_height: Coordinate3,
}

impl BoundingRegion {
    pub fn new(south_west_min_height: Coordinate3, north_east_max_height: Coordinate3) -> Self {
        Self {
            south_west_min_height,
            north_east_max_height,
        }
    }

    pub fn as_array(&self) -> [f64; 6] {
        let south_west_min_height_radian = self.south_west_min_height.to_radians();
        let north_east_max_height_radian = self.north_east_max_height.to_radians();

        [
            south_west_min_height_radian.x(),
            south_west_min_height_radian.y(),
            north_east_max_height_radian.x(),
            north_east_max_height_radian.y(),
            south_west_min_height_radian.z(),
            north_east_max_height_radian.z(),
        ]
    }
}