From a3697b9d1e9b4af0c71c772907c1e7d6cf71954f Mon Sep 17 00:00:00 2001
From: CSTEZCAN <33690601+CSTEZCAN@users.noreply.github.com>
Date: Sun, 9 Mar 2025 16:19:04 +0000
Subject: [PATCH] Update image_handling.py

threw this error when using with phase_hologram_github_logo_generation_and_reconstruction.py on a jupyter notebook;
"NotImplementedError: interp2d has been removed in SciPy 1.14.0."

updated.
---
 diffractsim/util/image_handling.py | 22 +++++++++++++---------
 1 file changed, 13 insertions(+), 9 deletions(-)

diff --git a/diffractsim/util/image_handling.py b/diffractsim/util/image_handling.py
index a6b40aa..26d1479 100644
--- a/diffractsim/util/image_handling.py
+++ b/diffractsim/util/image_handling.py
@@ -55,16 +55,20 @@ def convert_graymap_image_to_hsvmap_image(img):
 
 
 
-def resize_array(img_array, new_shape):
+from scipy.interpolate import RegularGridInterpolator
 
+def resize_array(img_array, new_shape):
     Ny, Nx = img_array.shape
+    x = np.linspace(0, 1, Nx)
+    y = np.linspace(0, 1, Ny)
+    
+    interpolator = RegularGridInterpolator((y, x), img_array, method="cubic", bounds_error=False, fill_value=0)
+    
+    new_x = np.linspace(0, 1, new_shape[1])
+    new_y = np.linspace(0, 1, new_shape[0])
+    new_grid = np.meshgrid(new_y, new_x, indexing='ij')
+    
+    resize_img_array = interpolator(np.array([new_grid[0].flatten(), new_grid[1].flatten()]).T).reshape(new_shape)
     
-    from scipy.interpolate import interp2d
-    fun = interp2d(
-        np.linspace(0, 1, Nx),
-        np.linspace(0, 1, Ny),
-        img_array,
-        kind="cubic",
-    )
-    resize_img_array = fun(np.linspace(0, 1, new_shape[1]), np.linspace(0, 1, new_shape[0]))
     return resize_img_array
+