地下埋設物を対象とした平面直角座標系への対応

nusamai-citygml/src/geometry.rs

@@ -84,6 +84,7 @@ pub struct SurfaceSpan {
 #[derive(Default)]
 pub(crate) struct GeometryCollector {
     pub vertices: indexmap::IndexSet<[u64; 3], ahash::RandomState>,
+    pub epsg: Option<EpsgCode>,
     pub multipolygon: MultiPolygon<'static, u32>,
     pub multilinestring: MultiLineString<'static, u32>,
     pub multipoint: MultiPoint<'static, u32>,
@@ -132,8 +133,14 @@ impl GeometryCollector {
             ]);
         }
 
+        let epsg = if let Some(epsg) = self.epsg {
+            epsg
+        } else {
+            EPSG_JGD2011_GEOGRAPHIC_3D
+        };
+
         GeometryStore {
-            epsg: EPSG_JGD2011_GEOGRAPHIC_3D,
+            epsg,
             vertices,
             multipolygon: self.multipolygon,
             multilinestring: self.multilinestring,

nusamai-citygml/src/parser.rs

@@ -1,5 +1,6 @@
 use std::{io::BufRead, mem, str};
 
+use nusamai_projection::crs::EpsgCode;
 use quick_xml::{
     events::{BytesStart, Event},
     name::{Namespace, ResolveResult::Bound},
@@ -514,19 +515,29 @@
         lod: u8,
     ) -> Result<(), ParseError> {
         let mut surface_id = None;
+        let mut epsg = None;
         loop {
             match self.reader.read_event_into(&mut self.state.buf1) {
                 Ok(Event::Start(start)) => {
                     let (nsres, localname) = self.reader.resolve_element(start.name());
                     let poly_begin = self.state.geometry_collector.multipolygon.len();
 
-                    // surface id
                     for attr in start.attributes().flatten() {
                         let (nsres, localname) = self.reader.resolve_attribute(attr.key);
+                        // surface id
                         if nsres == Bound(GML31_NS) && localname.as_ref() == b"id" {
                             let id = String::from_utf8_lossy(attr.value.as_ref()).to_string();
                             surface_id = Some(self.state.context.id_to_integer_id(id));
-                            break;
+                        }
+                        // epsg
+                        if localname.as_ref() == b"srsName" {
+                            let srsname = String::from_utf8_lossy(attr.value.as_ref()).to_string();
+                            let epsg_head = "http://www.opengis.net/def/crs/EPSG/0/";
+                            if srsname.starts_with(epsg_head) {
+                                if let Some(stripped) = srsname.strip_prefix(epsg_head) {
+                                    epsg = stripped.parse::<EpsgCode>().ok();
+                                }
+                            }
                         }
                     }
 
@@ -595,6 +606,10 @@
                                     end: poly_end as u32,
                                 });
                         }
+
+                        if let Some(epsg) = epsg {
+                            self.state.geometry_collector.epsg = Some(epsg);
+                        }
                     }
                 }
                 Ok(Event::End(_)) => break,