Fix depr warning in fringing variation level

LSSTDESC · Mar 19, 2024 · 7237452 · 7237452
1 parent 2a0f1e4
commit 7237452
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 7 deletions.
diff --git a/data/fringing_data/skyline_var.fits b/data/fringing_data/skyline_var.fits
diff --git a/data/fringing_data/skyline_var.pkl b/data/fringing_data/skyline_var.pkl
diff --git a/imsim/sky_model.py b/imsim/sky_model.py
@@ -1,11 +1,11 @@
-
+from astropy.io import fits
 import copy
 import warnings
 import numpy as np
 import galsim
 import pickle
 from galsim.config import InputLoader, RegisterInputType, RegisterValueType
-from scipy.interpolate import RegularGridInterpolator
+from scipy.interpolate import RegularGridInterpolator, RectBivariateSpline
 import os
 from .meta_data import data_dir
 
@@ -198,11 +198,17 @@ def fringe_variation_level(self):
 
         if self.spatial_vary:
             # Load 2d interpolator for OH spatial variation
-            filename = os.path.join(data_dir, 'fringing_data',
-                                        'skyline_var.pkl')
-            with open(filename, 'rb') as f:
-                interp = pickle.load(f)
-
+            filename = os.path.join(
+                data_dir,
+                'fringing_data',
+                'skyline_var.fits'
+            )
+            hdu = fits.open(filename)[0]
+            z = hdu.data
+            nx, ny = z.shape
+            x = np.linspace(hdu.header['XMIN'], hdu.header['XMAX'], nx)
+            y = np.linspace(hdu.header['YMIN'], hdu.header['YMAX'], ny)
+            interp = RectBivariateSpline(x, y, z)
             dx, dy = self.boresight.project(self.true_center)
             # calculated OH flux level wrst the center of focal plane.
             level = interp(dx.deg,dy.deg)/interp(0,0)