Delete null island points from the hull calculation

catenarytransit · Apr 1, 2024 · dcf0c4a · dcf0c4a
1 parent cfef04b
commit dcf0c4a
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/maple/gtfs_handlers/hull_from_gtfs.rs b/src/maple/gtfs_handlers/hull_from_gtfs.rs
@@ -12,6 +12,10 @@ pub fn hull_from_gtfs(gtfs: &gtfs_structures::Gtfs) -> Option<Polygon> {
                 .map(|(id, points)| {
                     points
                         .iter()
+                        .filter(|point| {
+                            let is_null_island = f64::abs(0. - point.latitude) < 0.000001 && f64::abs(0. - point.latitude) < 0.000001;
+                            !is_null_island
+                        })
                         .map(|point| Point::new(point.longitude, point.latitude))
                 })
                 .flatten()
@@ -27,6 +31,10 @@ pub fn hull_from_gtfs(gtfs: &gtfs_structures::Gtfs) -> Option<Polygon> {
                         .stops
                         .iter()
                         .filter(|(_, stop)| stop.longitude.is_some() && stop.latitude.is_some())
+                        .filter(|(_, stop)| {
+                            let is_null_island = f64::abs(0. - stop.latitude.unwrap()) < 0.000001 && f64::abs(0. - stop.latitude.unwrap()) < 0.000001;
+                            !is_null_island
+                        })
                         .map(|(_, stop)| {
                             Point::new(stop.longitude.unwrap(), stop.latitude.unwrap())
                         })