Implement geo_traits for geojson types.

Cargo.toml

@@ -12,11 +12,14 @@ edition = "2018"
 
 [features]
 default = ["geo-types"]
+geo-traits = ["dep:geo-traits", "dep:bytemuck"]
 
 [dependencies]
 serde = { version="~1.0", features = ["derive"] }
 serde_json = "~1.0"
 geo-types = { version = "0.7.13", features = ["serde"], optional = true }
+geo-traits = { version = "0.1.0", optional = true }
+bytemuck = { version = "1", features = ["derive"], optional = true }
 thiserror = "1.0.20"
 log = "0.4.17"
 

src/geo_traits_impl/coord.rs

@@ -0,0 +1,47 @@
+use super::PointType;
+use geo_traits::Dimensions;
+
+impl geo_traits::CoordTrait for PointType {
+    type T = f64;
+
+    fn dim(&self) -> Dimensions {
+        match self.0.len() {
+            0 | 1 => panic!("Position must have at least 2 dimensions"),
+            2 => Dimensions::Xy,
+            3 => Dimensions::Xyz,
+            _ => Dimensions::Unknown(self.0.len()),
+        }
+    }
+
+    fn x(&self) -> Self::T {
+        self.0[0]
+    }
+
+    fn y(&self) -> Self::T {
+        self.0[1]
+    }
+
+    fn nth_unchecked(&self, n: usize) -> Self::T {
+        self.0[n]
+    }
+}
+
+impl geo_traits::CoordTrait for &PointType {
+    type T = f64;
+
+    fn dim(&self) -> Dimensions {
+        PointType::dim(self)
+    }
+
+    fn x(&self) -> Self::T {
+        PointType::x(self)
+    }
+
+    fn y(&self) -> Self::T {
+        PointType::y(self)
+    }
+
+    fn nth_unchecked(&self, n: usize) -> Self::T {
+        PointType::nth_unchecked(self, n)
+    }
+}