sigwx-parser 0.1.0__tar.gz

sigwx_parser-0.1.0/PKG-INFO

@@ -0,0 +1,11 @@
+Metadata-Version: 2.4
+Name: sigwx-parser
+Version: 0.1.0
+Summary: A parser to convert WAFS IWXXM SIGWX XML files into GeoJSON FeatureCollections.
+Author-email: Miguel Ebersbach <ebersbachmiguel@gmail.com>
+Project-URL: Homepage, https://github.com/miguel-pub/sigwx-parser
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+Requires-Dist: lxml
+Requires-Dist: shapely
+Requires-Dist: antimeridian

sigwx_parser-0.1.0/pyproject.toml

@@ -0,0 +1,21 @@
+[build-system]
+requires = ["setuptools>=61.0"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "sigwx-parser"
+version = "0.1.0"
+description = "A parser to convert WAFS IWXXM SIGWX XML files into GeoJSON FeatureCollections."
+readme = "README.md"
+authors = [
+    { name="Miguel Ebersbach", email="ebersbachmiguel@gmail.com" },
+]
+requires-python = ">=3.10"
+dependencies = [
+    "lxml",
+    "shapely",
+    "antimeridian"
+]
+
+[project.urls]
+"Homepage" = "https://github.com/miguel-pub/sigwx-parser"

sigwx_parser-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

sigwx_parser-0.1.0/src/sigwx_parser/__init__.py

@@ -0,0 +1,2 @@
+from .parser import SigwxParser
+__all__ = ["SigwxParser"]

sigwx_parser-0.1.0/src/sigwx_parser/geometry.py

@@ -0,0 +1,67 @@
+import antimeridian
+from shapely.geometry import shape, mapping
+from shapely.geometry.polygon import orient
+
+def process_polygons(feature_type, polygon_patch, feature_id, upper_fl=None, lower_fl=None):
+    exterior = polygon_patch.find(".//exterior//posList").text.strip().split()
+    interiors = [int.text.strip().split() for int in polygon_patch.findall(".//interior//posList")]
+    exterior = [[float(y), float(x)] for x, y in zip(exterior[::2], exterior[1::2])]
+    interiors = [[[float(y), float(x)] for x, y in zip(int[::2], int[1::2])] for int in interiors]
+
+    feature = {
+        "type": "Feature",
+        "geometry": {
+            "type": "Polygon",
+            "coordinates": [exterior, *interiors],
+        },
+        "properties": {
+            "feature_id": feature_id,
+            "phenomenon": feature_type,
+            "upper_fl": upper_fl,
+            "lower_fl": lower_fl,
+        }
+    }
+    raw_shape = shape(feature["geometry"])
+    properly_wound_shape = orient(raw_shape, sign=1.0)
+    feature["geometry"] = mapping(properly_wound_shape)
+
+    fixed_feature = antimeridian.fix_geojson(feature)
+    return fixed_feature
+
+def process_points(feature_type, point, feature_id, upper_fl=None, lower_fl=None):
+    coords = point.find(".//pos").text.strip().split()
+    coords[1], coords[0] = float(coords[0]), float(coords[1])  # Swap to [lat, lon]
+    feature = {
+        "type": "Feature",
+        "geometry": {
+            "type": "Point",
+            "coordinates": coords,
+        },
+        "properties": {
+            "feature_id": feature_id,
+            "phenomenon": feature_type,
+            "upper_fl": upper_fl,
+            "lower_fl": lower_fl,
+        }
+    }
+    return feature
+
+def process_lines(feature_type, line, feature_id, upper_fl=None, lower_fl=None):
+    coords = line.find(".//posList").text.strip().split()
+    coords = [[float(y), float(x)] for x, y in zip(coords[::2], coords[1::2])]
+    feature = {
+        "type": "Feature",
+        "geometry": {
+            "type": "LineString",
+            "coordinates": coords,
+        },
+        "properties": {
+            "feature_id": feature_id,
+            "phenomenon": feature_type,
+            "upper_fl": upper_fl,
+            "lower_fl": lower_fl,
+        }
+    }
+    fixed_feature = antimeridian.fix_geojson(feature)
+    return fixed_feature
+

sigwx_parser-0.1.0/src/sigwx_parser/parser.py

@@ -0,0 +1,95 @@
+import json
+from pathlib import Path
+from geometry import process_points, process_lines, process_polygons
+from lxml import etree
+
+class SigwxParser:
+    def __init__(self, file_path: str, xml_content: str = None):
+        self.file_path = file_path
+        self.xml_content = xml_content
+        self.root = None
+
+        if file_path:
+            path = Path(file_path)
+            if not path.is_file():
+                raise FileNotFoundError(f"Could not find a file at: {path.absolute()}")
+            raw_text = path.read_text(encoding='utf-8')
+        elif xml_content:
+            raw_text = xml_content
+        else:
+            raise ValueError("Either file_path or xml_content must be provided.")
+        self.root = self._clean_xml_string(raw_text)
+
+    @staticmethod
+    def _clean_xml_string(raw_text: str):
+        clean_text = raw_text.replace('\uFEFF', '').strip()
+        parser = etree.XMLParser(recover=True)
+        root = etree.fromstring(clean_text.encode('utf-8'), parser=parser)
+
+        for elem in root.iter():
+            elem.tag = etree.QName(elem).localname
+
+            clean_attribs = {etree.QName(k).localname: v for k, v in elem.attrib.items()}
+            elem.attrib.clear()
+            elem.attrib.update(clean_attribs)
+
+        etree.cleanup_namespaces(root)
+
+        return root
+    @staticmethod
+    def _safe_text(element, default="N/A"):
+        if element is not None and element.text is not None:
+            return element.text.strip()
+        return default
+    def parse(self):
+        forecasts = self.root.xpath("//WAFSSignificantWeatherForecast")
+        collections = []
+        for f in forecasts:
+            valid_time = f.find(".//beginPosition").text.strip()
+            issue_time = f.find(".//timePosition").text.strip()
+            feature_collection = {
+                "type": "FeatureCollection",
+                "features": []
+            }
+            met_features = f.xpath(".//MeteorologicalFeature")
+
+            for feature in met_features:
+                feature_type = feature.find(".//phenomenon").get("href").strip().split("/")[-1] or "N/A"
+                feature_id = feature.get("id").strip() if feature.get("id") else "N/A"
+                upper_fl = self._safe_text(feature.find(".//upperElevation"))
+                lower_fl = self._safe_text(feature.find(".//lowerElevation"))
+
+
+                if feature.find(".//PolygonPatch") is not None:
+                    processed_polygon = process_polygons(feature_type=feature_type,
+                                                         polygon_patch=feature.find(".//PolygonPatch"),
+                                                         feature_id=feature_id,
+                                                         upper_fl=upper_fl,
+                                                         lower_fl=lower_fl)
+
+                    feature_collection.get("features").append(processed_polygon)
+
+                elif feature.find(".//Point") is not None:
+                    processed_point = process_points(feature_type=feature_type,
+                                                     point=feature.find(".//Point"),
+                                                     feature_id=feature_id,
+                                                     upper_fl=upper_fl,
+                                                     lower_fl=lower_fl)
+                    feature_collection.get("features").append(processed_point)
+                else:
+                    processed_line = process_lines(feature_type=feature_type,
+                                                   line=feature.find(".//CubicSpline"),
+                                                   feature_id=feature_id,
+                                                   upper_fl=upper_fl,
+                                                   lower_fl=lower_fl)
+                    feature_collection.get("features").append(processed_line)
+
+            collections.append({
+                "valid_time": valid_time,
+                "issue_time": issue_time,
+                "feature_collection": feature_collection})
+        return collections
+
+
+
+

sigwx_parser-0.1.0/src/sigwx_parser.egg-info/PKG-INFO

@@ -0,0 +1,11 @@
+Metadata-Version: 2.4
+Name: sigwx-parser
+Version: 0.1.0
+Summary: A parser to convert WAFS IWXXM SIGWX XML files into GeoJSON FeatureCollections.
+Author-email: Miguel Ebersbach <ebersbachmiguel@gmail.com>
+Project-URL: Homepage, https://github.com/miguel-pub/sigwx-parser
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+Requires-Dist: lxml
+Requires-Dist: shapely
+Requires-Dist: antimeridian

sigwx_parser-0.1.0/src/sigwx_parser.egg-info/SOURCES.txt

@@ -0,0 +1,10 @@
+README.md
+pyproject.toml
+src/sigwx_parser/__init__.py
+src/sigwx_parser/geometry.py
+src/sigwx_parser/parser.py
+src/sigwx_parser.egg-info/PKG-INFO
+src/sigwx_parser.egg-info/SOURCES.txt
+src/sigwx_parser.egg-info/dependency_links.txt
+src/sigwx_parser.egg-info/requires.txt
+src/sigwx_parser.egg-info/top_level.txt

sigwx_parser-0.1.0/src/sigwx_parser.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

sigwx_parser-0.1.0/src/sigwx_parser.egg-info/requires.txt

@@ -0,0 +1,3 @@
+lxml
+shapely
+antimeridian

sigwx_parser-0.1.0/src/sigwx_parser.egg-info/top_level.txt

@@ -0,0 +1 @@
+sigwx_parser