Fix division throwing error in some cases

GispoCoding · Dec 3, 2024 · fd0cb2f · fd0cb2f
1 parent 060c2a0
commit fd0cb2f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/eis_toolkit/vector_processing/idw_interpolation.py b/eis_toolkit/vector_processing/idw_interpolation.py
@@ -67,7 +67,7 @@ def _idw_core(x, y, z, xi, yi: np.ndarray, power: Number, search_radius: Optiona
             over += z[n] / dist
             under += 1.0 / dist
 
-    interpolated_values = over / under
+    interpolated_values = np.divide(over, under, out=np.full_like(over, np.nan), where=under != 0)
     return interpolated_values