Don't use Delaunay with fewer than 4 points

metno · Jan 16, 2025 · 5a6d517 · 5a6d517
1 parent 96465b8
commit 5a6d517
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/bris/outputs/verif.py b/bris/outputs/verif.py
@@ -80,12 +80,12 @@ def __init__(
         self.opoints = gridpp.Points(self.obs_lats, self.obs_lons, self.obs_altitudes)
         self.opoints_tuple = np.column_stack((self.obs_lats, self.obs_lons))
 
-
-        if not self._is_gridded_input:
+        self.triangulation = self.ipoints_tuple
+        if not self._is_gridded_input and len(self.ipoints_tuple[0]) > 3:
+            # This speeds up interpolation from irregular points to observation points
+            # but Delaunay needs enough points for this to work
             self.triangulation = Delaunay(self.ipoints_tuple)
-        else:
-            self.triangulation = None
-
+
         # The intermediate will only store the final output locations
         intermediate_pm = PredictMetadata(
             [variable],