OpenGeodeWeb-Back 6.16.1__tar.gz → 6.17.0__tar.gz

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: OpenGeodeWeb-Back
-Version: 6.16.1
+Version: 6.17.0
 Summary: OpenGeodeWeb-Back is an open source framework that proposes handy python functions and wrappers for the OpenGeode ecosystem
 Author-email: Geode-solutions <team-web@geode-solutions.com>
 Project-URL: Homepage, https://github.com/Geode-solutions/OpenGeodeWeb-Back
@@ -13,6 +13,7 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: asgiref==3.*
 Requires-Dist: blinker==1.*
+Requires-Dist: certifi==2026.7.22
 Requires-Dist: click==8.*
 Requires-Dist: defusedxml==0.7.1
 Requires-Dist: flask[async]==3.*
@@ -28,6 +29,7 @@ Requires-Dist: opengeode-geosciencesio==5.13.1
 Requires-Dist: opengeode-inspector==6.12.4
 Requires-Dist: opengeode-inspector-geosciences==1.2.4
 Requires-Dist: opengeode-io==7.11.2
+Requires-Dist: pillow==13.0.0
 Requires-Dist: types-defusedxml==0.7.0.20260504
 Requires-Dist: werkzeug==3.1.8
 Requires-Dist: opengeodeweb-microservice==1.*,>=1.3.2

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "OpenGeodeWeb-Back"
-version = "6.16.1"
+version = "6.17.0"
 dynamic = ["dependencies"]
 authors = [{ name = "Geode-solutions", email = "team-web@geode-solutions.com" }]
 description = "OpenGeodeWeb-Back is an open source framework that proposes handy python functions and wrappers for the OpenGeode ecosystem"
@@ -33,4 +33,6 @@ where = ["src"]
 "opengeodeweb_back.routes.schemas" = ["*.json"]
 "opengeodeweb_back.routes.models.schemas" = ["*.json"]
 "opengeodeweb_back.routes.create.schemas" = ["*.json"]
+"opengeodeweb_back.routes.crs.schemas" = ["*.json"]
+"opengeodeweb_back.routes.basemap.schemas" = ["*.json"]
 "opengeodeweb_back" = ["py.typed"]

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/requirements.txt

@@ -8,6 +8,8 @@ asgiref==3.*
     # via flask
 blinker==1.*
     # via flask
+certifi==2026.7.22
+    # via -r requirements.in
 click==8.*
     # via flask
 defusedxml==0.7.1
@@ -62,6 +64,8 @@ opengeode-io==7.11.2
     #   -r requirements.in
     #   geode-viewables
     #   opengeode-geosciencesio
+pillow==13.0.0
+    # via -r requirements.in
 types-defusedxml==0.7.0.20260504
     # via -r requirements.in
 werkzeug==3.1.8

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/src/OpenGeodeWeb_Back.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: OpenGeodeWeb-Back
-Version: 6.16.1
+Version: 6.17.0
 Summary: OpenGeodeWeb-Back is an open source framework that proposes handy python functions and wrappers for the OpenGeode ecosystem
 Author-email: Geode-solutions <team-web@geode-solutions.com>
 Project-URL: Homepage, https://github.com/Geode-solutions/OpenGeodeWeb-Back
@@ -13,6 +13,7 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: asgiref==3.*
 Requires-Dist: blinker==1.*
+Requires-Dist: certifi==2026.7.22
 Requires-Dist: click==8.*
 Requires-Dist: defusedxml==0.7.1
 Requires-Dist: flask[async]==3.*
@@ -28,6 +29,7 @@ Requires-Dist: opengeode-geosciencesio==5.13.1
 Requires-Dist: opengeode-inspector==6.12.4
 Requires-Dist: opengeode-inspector-geosciences==1.2.4
 Requires-Dist: opengeode-io==7.11.2
+Requires-Dist: pillow==13.0.0
 Requires-Dist: types-defusedxml==0.7.0.20260504
 Requires-Dist: werkzeug==3.1.8
 Requires-Dist: opengeodeweb-microservice==1.*,>=1.3.2

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/src/OpenGeodeWeb_Back.egg-info/SOURCES.txt

@@ -16,7 +16,16 @@ src/opengeodeweb_back/py.typed
 src/opengeodeweb_back/test_utils.py
 src/opengeodeweb_back/typed_route.py
 src/opengeodeweb_back/utils_functions.py
+src/opengeodeweb_back/basemap/__init__.py
+src/opengeodeweb_back/basemap/basemap_mesh.py
+src/opengeodeweb_back/basemap/basemap_providers.py
+src/opengeodeweb_back/basemap/basemap_tiles.py
+src/opengeodeweb_back/crs/__init__.py
+src/opengeodeweb_back/crs/crs_alignment.py
+src/opengeodeweb_back/crs/crs_assignment.py
+src/opengeodeweb_back/crs/crs_functions.py
 src/opengeodeweb_back/geode_objects/__init__.py
+src/opengeodeweb_back/geode_objects/crs_coordinates.py
 src/opengeodeweb_back/geode_objects/geode_brep.py
 src/opengeodeweb_back/geode_objects/geode_cross_section.py
 src/opengeodeweb_back/geode_objects/geode_edged_curve2d.py
@@ -53,15 +62,37 @@ src/opengeodeweb_back/geode_objects/geode_triangulated_surface3d.py
 src/opengeodeweb_back/geode_objects/geode_vertex_set.py
 src/opengeodeweb_back/routes/__init__.py
 src/opengeodeweb_back/routes/blueprint_routes.py
+src/opengeodeweb_back/routes/basemap/__init__.py
+src/opengeodeweb_back/routes/basemap/blueprint_basemap.py
+src/opengeodeweb_back/routes/basemap/schemas/__init__.py
+src/opengeodeweb_back/routes/basemap/schemas/build.json
+src/opengeodeweb_back/routes/basemap/schemas/build.py
+src/opengeodeweb_back/routes/basemap/schemas/delete.json
+src/opengeodeweb_back/routes/basemap/schemas/delete.py
 src/opengeodeweb_back/routes/create/__init__.py
 src/opengeodeweb_back/routes/create/blueprint_create.py
 src/opengeodeweb_back/routes/create/schemas/__init__.py
 src/opengeodeweb_back/routes/create/schemas/edged_curve.json
 src/opengeodeweb_back/routes/create/schemas/edged_curve.py
+src/opengeodeweb_back/routes/create/schemas/horizon_stack.json
+src/opengeodeweb_back/routes/create/schemas/horizon_stack.py
 src/opengeodeweb_back/routes/create/schemas/point_set.json
 src/opengeodeweb_back/routes/create/schemas/point_set.py
 src/opengeodeweb_back/routes/create/schemas/polygonal_surface.json
 src/opengeodeweb_back/routes/create/schemas/polygonal_surface.py
+src/opengeodeweb_back/routes/crs/__init__.py
+src/opengeodeweb_back/routes/crs/blueprint_crs.py
+src/opengeodeweb_back/routes/crs/schemas/__init__.py
+src/opengeodeweb_back/routes/crs/schemas/active.json
+src/opengeodeweb_back/routes/crs/schemas/active.py
+src/opengeodeweb_back/routes/crs/schemas/align.json
+src/opengeodeweb_back/routes/crs/schemas/align.py
+src/opengeodeweb_back/routes/crs/schemas/assign_geographic.json
+src/opengeodeweb_back/routes/crs/schemas/assign_geographic.py
+src/opengeodeweb_back/routes/crs/schemas/convert_geographic.json
+src/opengeodeweb_back/routes/crs/schemas/convert_geographic.py
+src/opengeodeweb_back/routes/crs/schemas/geographic_list.json
+src/opengeodeweb_back/routes/crs/schemas/geographic_list.py
 src/opengeodeweb_back/routes/schemas/__init__.py
 src/opengeodeweb_back/routes/schemas/allowed_files.json
 src/opengeodeweb_back/routes/schemas/allowed_files.py
@@ -85,8 +116,8 @@ src/opengeodeweb_back/routes/schemas/geode_object_inheritance.json
 src/opengeodeweb_back/routes/schemas/geode_object_inheritance.py
 src/opengeodeweb_back/routes/schemas/geode_objects_and_output_extensions.json
 src/opengeodeweb_back/routes/schemas/geode_objects_and_output_extensions.py
-src/opengeodeweb_back/routes/schemas/geographic_coordinate_systems.json
-src/opengeodeweb_back/routes/schemas/geographic_coordinate_systems.py
+src/opengeodeweb_back/routes/schemas/horizon_stack_info.json
+src/opengeodeweb_back/routes/schemas/horizon_stack_info.py
 src/opengeodeweb_back/routes/schemas/import_project.json
 src/opengeodeweb_back/routes/schemas/import_project.py
 src/opengeodeweb_back/routes/schemas/kill.json
@@ -133,7 +164,9 @@ src/opengeodeweb_back/routes/schemas/vertex_attribute_names.json
 src/opengeodeweb_back/routes/schemas/vertex_attribute_names.py
 src/opengeodeweb_back/routes/schemas/vertex_attribute_range.json
 src/opengeodeweb_back/routes/schemas/vertex_attribute_range.py
+tests/test_basemap.py
 tests/test_create_routes.py
+tests/test_crs_routes.py
 tests/test_events.py
 tests/test_geode_functions.py
 tests/test_routes.py

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/src/OpenGeodeWeb_Back.egg-info/requires.txt

@@ -1,5 +1,6 @@
 asgiref==3.*
 blinker==1.*
+certifi==2026.7.22
 click==8.*
 defusedxml==0.7.1
 flask[async]==3.*
@@ -15,6 +16,7 @@ opengeode-geosciencesio==5.13.1
 opengeode-inspector==6.12.4
 opengeode-inspector-geosciences==1.2.4
 opengeode-io==7.11.2
+pillow==13.0.0
 types-defusedxml==0.7.0.20260504
 werkzeug==3.1.8
 opengeodeweb-microservice==1.*,>=1.3.2

{opengeodeweb_back-6.16.1 → opengeodeweb_back-6.17.0}/src/opengeodeweb_back/app.py

@@ -19,7 +19,9 @@ from werkzeug.exceptions import HTTPException
 
 from opengeodeweb_back import app_config, utils_functions
 from opengeodeweb_back.routes import blueprint_routes
+from opengeodeweb_back.routes.basemap import blueprint_basemap
 from opengeodeweb_back.routes.create import blueprint_create
+from opengeodeweb_back.routes.crs import blueprint_crs
 
 if TYPE_CHECKING:
     from collections.abc import Generator
@@ -171,6 +173,16 @@ def register_ogw_back_blueprints(app: flask.Flask) -> None:
         url_prefix="/opengeodeweb_back/create",
         name="opengeodeweb_create",
     )
+    app.register_blueprint(
+        blueprint_crs.routes,
+        url_prefix="/opengeodeweb_back/crs",
+        name="opengeodeweb_crs",
+    )
+    app.register_blueprint(
+        blueprint_basemap.routes,
+        url_prefix="/opengeodeweb_back/basemap",
+        name="opengeodeweb_basemap",
+    )
 
 
 def run_server(app: Flask) -> None:

opengeodeweb_back-6.17.0/src/opengeodeweb_back/basemap/basemap_mesh.py

@@ -0,0 +1,93 @@
+from __future__ import annotations
+
+# Standard library imports
+import math
+from typing import TYPE_CHECKING
+
+# Third party imports
+import opengeode as og
+import opengeode_geosciences as og_geosciences
+
+# Local application imports
+from opengeodeweb_back.basemap.basemap_tiles import HALF_WORLD
+from opengeodeweb_back.crs.crs_functions import convert_xy_corners
+from opengeodeweb_back.geode_objects.geode_triangulated_surface3d import (
+    GeodeTriangulatedSurface3D,
+)
+
+if TYPE_CHECKING:
+    from opengeodeweb_back.basemap.basemap_tiles import TileRange
+
+GRID_SIZE = 32
+MESH_NAME = "basemap"
+UV_ATTRIBUTE = "basemap_uv"
+WEB_MERCATOR = og_geosciences.GeographicCoordinateSystemInfo(
+    "EPSG", "3857", "WGS 84 / Pseudo-Mercator"
+)
+OUTSIDE_MESSAGE = "Area outside the map projection limits"
+
+
+class OutsideProjectionError(Exception):
+    pass
+
+
+def bounds_to_3857(
+    bounds: tuple[float, float, float, float],
+    crs: og_geosciences.GeographicCoordinateSystemInfo,
+) -> tuple[float, float, float, float]:
+    """Bounding box in Web Mercator of the four corners of bounds, given in crs."""
+    try:
+        corners = convert_xy_corners(bounds, crs, WEB_MERCATOR)
+    except RuntimeError as error:  # PROJ rejects coordinates outside either CRS domain
+        raise OutsideProjectionError(OUTSIDE_MESSAGE) from error
+    xs = [x for x, _ in corners]
+    ys = [y for _, y in corners]
+    projected = (min(xs), min(ys), max(xs), max(ys))
+    if any(not abs(value) <= HALF_WORLD for value in projected):  # also rejects NaN
+        raise OutsideProjectionError(OUTSIDE_MESSAGE)
+    return projected
+
+
+def build_mesh(
+    tiles: TileRange,
+    crs: og_geosciences.GeographicCoordinateSystemInfo,
+    z_value: float,
+) -> GeodeTriangulatedSurface3D:
+    """Grid covering the tiles, with texture coordinates, converted from Web Mercator to crs."""
+    x0, y0, x1, y1 = tiles.bounds_3857()
+    mesh = GeodeTriangulatedSurface3D()
+    builder = mesh.builder()
+    builder.set_name(MESH_NAME)
+    attributes = mesh.triangulated_surface.vertex_attribute_manager()
+    values = og.AttributeValuesArrayDouble2()
+    values.default_value = [0.0, 0.0]
+    uv_id = attributes.create_attribute_variable_arraydouble2(
+        UV_ATTRIBUTE, values, og.AttributeProperties()
+    )
+    uv = attributes.find_attribute_variable_arraydouble2(uv_id)
+    last = GRID_SIZE - 1
+    for row in range(GRID_SIZE):
+        for column in range(GRID_SIZE):
+            u, v = column / last, row / last
+            point = og.Point3D([x0 + u * (x1 - x0), y0 + v * (y1 - y0), z_value])
+            uv.set_value(builder.create_point(point), [u, v])
+    for row in range(last):
+        for column in range(last):
+            corner = row * GRID_SIZE + column
+            builder.create_triangle([corner, corner + 1, corner + GRID_SIZE + 1])
+            builder.create_triangle([corner, corner + GRID_SIZE + 1, corner + GRID_SIZE])
+    builder.compute_polygon_adjacencies()
+    mesh.assign_crs(WEB_MERCATOR.authority_code(), WEB_MERCATOR)
+    # Assigning then converting to the same CRS name is rejected by OpenGeode
+    if crs.authority_code() == WEB_MERCATOR.authority_code():
+        return mesh
+    try:
+        mesh.convert_crs(crs.authority_code(), crs)
+    except RuntimeError as error:  # PROJ rejects tiles outside the project CRS domain
+        raise OutsideProjectionError(OUTSIDE_MESSAGE) from error
+    surface = mesh.triangulated_surface
+    for vertex in range(surface.nb_vertices()):
+        point = surface.point(vertex)
+        if not all(math.isfinite(point.value(axis)) for axis in range(3)):
+            raise OutsideProjectionError(OUTSIDE_MESSAGE)
+    return mesh

opengeodeweb_back-6.17.0/src/opengeodeweb_back/basemap/basemap_providers.py

@@ -0,0 +1,100 @@
+from __future__ import annotations
+
+# Standard library imports
+import os
+from dataclasses import dataclass
+from typing import TYPE_CHECKING
+
+if TYPE_CHECKING:
+    from collections.abc import Container, Iterable
+
+USER_AGENT = "Vease (+https://geode-solutions.com)"
+# Every source is cut in 256 px Web Mercator tiles: they share one tile grid
+TILE_SIZE = 256
+
+
+@dataclass(frozen=True)
+class Source:
+    url_template: str
+    max_zoom: int
+    attribution: str
+
+
+_GEOPLATEFORME = (
+    "https://data.geopf.fr/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER={layer}"
+    "&STYLE=normal&TILEMATRIXSET=PM&TILEMATRIX={{z}}&TILEROW={{y}}&TILECOL={{x}}&FORMAT={format}"
+)
+_IGN = "© IGN"
+
+SOURCES: dict[str, Source] = {
+    "eox_s2cloudless_2016": Source(
+        "https://tiles.maps.eox.at/wmts/1.0.0/s2cloudless_3857/default/g/{z}/{y}/{x}.jpg",
+        17,
+        "Sentinel-2 cloudless 2016 by EOX IT Services GmbH "
+        "(Contains modified Copernicus Sentinel data 2016)",
+    ),
+    "ign_orthophotos": Source(
+        _GEOPLATEFORME.format(layer="ORTHOIMAGERY.ORTHOPHOTOS", format="image/jpeg"), 19, _IGN
+    ),
+    # Volunteer servers: to be replaced before production, see the basemap design
+    "osm": Source(
+        "https://tile.openstreetmap.org/{z}/{x}/{y}.png", 19, "© OpenStreetMap contributors"
+    ),
+    "ign_plan": Source(
+        _GEOPLATEFORME.format(layer="GEOGRAPHICALGRIDSYSTEMS.PLANIGNV2", format="image/png"),
+        19,
+        _IGN,
+    ),
+    "ign_roads": Source(
+        _GEOPLATEFORME.format(layer="TRANSPORTNETWORKS.ROADS", format="image/png"), 18, _IGN
+    ),
+    "eox_overlay": Source(
+        "https://tiles.maps.eox.at/wmts/1.0.0/overlay_3857/default/g/{z}/{y}/{x}.png",
+        14,
+        "© OpenStreetMap contributors, rendering © EOX",
+    ),
+}
+
+# The layers of each map, from the bottom up. A layer draws the first of its sources serving a
+# tile: the next one is only fetched when it serves none
+MAPS: dict[str, tuple[tuple[str, ...], ...]] = {
+    # The IGN covers the world imagery or map where it has tiles (France)
+    "satellite": (("eox_s2cloudless_2016",), ("ign_orthophotos",)),
+    "plan": (("osm",), ("ign_plan",)),
+    # Transparent: both sources would draw the same roads twice, EOX outside France
+    "routes": (("ign_roads", "eox_overlay"),),
+}
+
+
+def url_template(source: str) -> str:
+    # Local servers for tests: one variable per source, else one for all where {source} is replaced
+    template = os.environ.get(f"BASEMAP_TILE_URL_{source.upper()}")
+    if template is not None:
+        return template
+    shared = os.environ.get("BASEMAP_TILE_URL")
+    if shared is not None:
+        return shared.replace("{source}", source)
+    return SOURCES[source].url_template
+
+
+def sorted_maps(maps: Iterable[str]) -> list[str]:
+    requested = set(maps)
+    return [name for name in MAPS if name in requested]
+
+
+def map_layers(maps: list[str], served: Container[str]) -> dict[str, list[str]]:
+    """The served sources of each sorted map, from the bottom up: maps without any are left out."""
+    layers = {
+        name: [source for layer in MAPS[name] for source in layer if source in served]
+        for name in maps
+    }
+    return {name: sources for name, sources in layers.items() if sources}
+
+
+def max_zoom(sources: Iterable[str]) -> int:
+    # The grid follows the finest source: coarser ones are overzoomed
+    return max(SOURCES[source].max_zoom for source in sources)
+
+
+def attribution(sources: Iterable[str]) -> str:
+    return " | ".join(dict.fromkeys(SOURCES[source].attribution for source in sources))