Flatten geometry to 2D before writing GeoPackage

Sites with an elevation produce 3D point geometry ([lon, lat, elev]).
GeoServer's GeoPackage reader then fails computing bounds:
  MismatchedDimensionException: Argument "WGS 84" has 3 dimensions, while 2 was expected.

Call force_2d() on the GeoDataFrame before to_file(GPKG). Elevation is
still carried as an attribute; only the geometry dimensionality changes.
The GeoJSON product on GCS is unaffected.

Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>

Author: jirhiker

orchestration/assets/products.py

@@ -172,6 +172,11 @@ def _geojson_to_geopackage(geojson_path: Path, layer_name: str, out_dir: Path) -
     if gdf.crs is None:
         gdf = gdf.set_crs("EPSG:4326")
 
+    # Sites with an elevation get 3D point geometry in the GeoJSON. GeoServer's
+    # GeoPackage reader rejects a 3D CRS ("WGS 84 has 3 dimensions") when
+    # computing bounds, so flatten to 2D — elevation remains an attribute.
+    gdf["geometry"] = gdf.geometry.force_2d()
+
     gpkg_path = out_dir / f"{layer_name}.gpkg"
     gdf.to_file(gpkg_path, driver="GPKG", layer=layer_name)
     return gpkg_path