polytope-python 2.1.8__tar.gz → 2.1.10__tar.gz

{polytope_python-2.1.8/polytope_python.egg-info → polytope_python-2.1.10}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: polytope-python
-Version: 2.1.8
+Version: 2.1.10
 Summary: Polytope datacube feature extraction library
 Author-email: "European Centre for Medium-Range Weather Forecasts (ECMWF)" <software.support@ecmwf.int>
 Maintainer-email: James Hawkes <James.Hawkes@ecmwf.int>, Mathilde Leuridan <Mathilde.Leuridan@ecmwf.int>

{polytope_python-2.1.8 → polytope_python-2.1.10}/polytope_feature/_version.py

@@ -1,2 +1,2 @@
 # Do not change! Do not track in version control!
-__version__ = "2.1.8"
+__version__ = "2.1.10"

{polytope_python-2.1.8 → polytope_python-2.1.10}/polytope_feature/engine/quadtree_slicer.py

@@ -19,6 +19,7 @@ class QuadTreeSlicer(Engine):
         # TODO: maybe we create the quadtree as soon as we have an unstructured slicer type and return it
         # to the slicer somehow?
         quad_tree = QuadTree()
+        # NOTE: the points here are assumed to be lat/lon implicitly
         points = [tuple(point) for point in points]
         quad_tree.build_point_tree(points)
         self.points = points
@@ -27,17 +28,29 @@ class QuadTreeSlicer(Engine):
     def extract_single(self, datacube, polytope):
         # extract a single polygon
         # if need to find nearest points, then take alternative slicing method using quadtree to find nearest point
+        axes = polytope.axes()
+        assert len(axes) == 2
+        assert "latitude" in axes and "longitude" in axes
+        revert_axes = not (list(axes) == ["latitude", "longitude"])
         if use_rust:
             if len(datacube.nearest_search) == 0:
-                polytope_points = [tuple(point) for point in polytope.points]
+                if revert_axes:
+                    polytope_points = [tuple(reversed(point)) for point in polytope.points]
+                else:
+                    polytope_points = [tuple(point) for point in polytope.points]
                 polygon_points = self.quad_tree.query_polygon(self.points, 0, polytope_points)
             else:
                 k = datacube.nearest_search[tuple(polytope.axes())][1]
-                nn_points = [tuple(pt) for pt in datacube.nearest_search[tuple(polytope.axes())][0]]
+                if revert_axes:
+                    nn_points = [tuple(reversed(pt)) for pt in datacube.nearest_search[tuple(polytope.axes())][0]]
+                else:
+                    nn_points = [tuple(pt) for pt in datacube.nearest_search[tuple(polytope.axes())][0]]
                 polygon_points = []
                 for nn_pt in nn_points:
                     polygon_points.extend(self.quad_tree.k_nearest_neighbor(nn_pt, k, self.points))
         else:
+            if revert_axes:
+                polytope.points = [tuple(reversed(point)) for point in polytope.points]
             polygon_points = self.quad_tree.query_polygon(polytope)
         return polygon_points
 

{polytope_python-2.1.8 → polytope_python-2.1.10}/polytope_feature/options.py

@@ -1,12 +1,11 @@
 import argparse
+import logging
 from abc import ABC
 from typing import Dict, List, Literal, Optional, Tuple, Union
 
 from conflator import ConfigModel, Conflator
 from pydantic import ConfigDict
 
-from polytope_feature.datacube.switching_grid_helper import lookup_grid_config
-
 
 class TransformationConfig(ConfigModel):
     model_config = ConfigDict(extra="forbid")
@@ -111,9 +110,12 @@ class PolytopeOptions(ABC):
                 replaced = replace_grid_config_in_options(config_options, pre_path)
                 if replaced:
                     axis_config = config_options.axis_config
-            except Exception:
-                # Fail silently and continue with original config
-                pass
+            except Exception as e:
+                logging.warning(
+                    "Dynamic grid replacement failed for georef '%s': %s. Using static grid config.",
+                    pre_path.get("georef", "unknown"),
+                    e,
+                )
         return (
             axis_config,
             compressed_axes_config,
@@ -151,6 +153,8 @@ def gridspec_to_grid_config(gridspec, md5hash):
 
 
 def replace_grid_config_in_options(options, req):
+    from polytope_feature.datacube.switching_grid_helper import lookup_grid_config
+
     gridspec, md5hash = lookup_grid_config(req)
     grid_config = gridspec_to_grid_config(gridspec, md5hash)
     if grid_config is not None:

{polytope_python-2.1.8 → polytope_python-2.1.10/polytope_python.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: polytope-python
-Version: 2.1.8
+Version: 2.1.10
 Summary: Polytope datacube feature extraction library
 Author-email: "European Centre for Medium-Range Weather Forecasts (ECMWF)" <software.support@ecmwf.int>
 Maintainer-email: James Hawkes <James.Hawkes@ecmwf.int>, Mathilde Leuridan <Mathilde.Leuridan@ecmwf.int>