feat: try to create a WCS object

iris_pipeline/resample/resample.py

@@ -5,6 +5,9 @@
 # try first importing everything from JWST
 from jwst import datamodels
 
+from astropy.io import fits
+from astropy.wcs import WCS
+
 from . import gwcs_drizzle
 from . import resample_utils
 from jwst.model_blender import blendmeta
@@ -47,6 +50,10 @@ def __init__(self, input_models, output=None, **pars):
             filename for output
         """
         self.input_models = input_models
+        for i in input_models:
+            with fits.open("/home/azonca/p/issues/202002_mosaic/example_field/iris_sim_gc_filterKN3_fix.fits") as f:
+                wcs = WCS(f[0].header)
+            i.meta.wcs = wcs
         self.drizpars = pars
         if output is None:
             output = input_models.meta.resample.output