Merge pull request #23 from svank/interpolation-points

Align interpolation points to data points
punch-mission · Apr 6, 2023 · 6fdcba1 · 6fdcba1
2 parents 9ba9a28 + 61d5123
commit 6fdcba1
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/regularizepsf/fitter.py b/regularizepsf/fitter.py
@@ -346,12 +346,12 @@ def generator():
                 interpolator = RectBivariateSpline(np.arange(image.shape[0]), 
                                                    np.arange(image.shape[1]), 
                                                    image)
-                image = interpolator(np.linspace(0, 
-                                                 image.shape[0], 
-                                                 image.shape[0] * interpolation_scale),
-                                     np.linspace(0, 
-                                                 image.shape[1], 
-                                                 image.shape[1] * interpolation_scale))
+                image = interpolator(np.linspace(0,
+                                                 image.shape[0] - 1,
+                                                 1 + (image.shape[0] - 1) * interpolation_scale),
+                                     np.linspace(0,
+                                                 image.shape[1] - 1,
+                                                 1 + (image.shape[1] - 1) * interpolation_scale))
 
             # find stars using SEP
             background = sep.Background(image)