From 4e9733b0df65cfb11948b2c71ad2b332a55d360b Mon Sep 17 00:00:00 2001 From: Thibault Barnouin Date: Fri, 14 Mar 2025 14:30:30 +0100 Subject: [PATCH] Save the raw total flux image as PrimaryHDU rebase display on main --- package/FOC_reduction.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/package/FOC_reduction.py b/package/FOC_reduction.py index 1d4443d..c2da366 100755 --- a/package/FOC_reduction.py +++ b/package/FOC_reduction.py @@ -46,8 +46,8 @@ def main(target=None, proposal_id=None, infiles=None, output_dir="./data", crop= display_bkg = True # Data binning - pxsize = 4 - pxscale = "px" # pixel, arcsec or full + pxsize = 0.05 + pxscale = "arcsec" # pixel, arcsec or full rebin_operation = "sum" # sum or average # Alignement @@ -60,8 +60,8 @@ def main(target=None, proposal_id=None, infiles=None, output_dir="./data", crop= # Smoothing smoothing_function = "combine" # gaussian_after, weighted_gaussian_after, gaussian, weighted_gaussian or combine - smoothing_FWHM = 1.5 # If None, no smoothing is done - smoothing_scale = "px" # pixel or arcsec + smoothing_FWHM = 0.075 # If None, no smoothing is done + smoothing_scale = "arcsec" # pixel or arcsec # Rotation rotate_data = False @@ -274,7 +274,7 @@ def main(target=None, proposal_id=None, infiles=None, output_dir="./data", crop= flux_data=np.array([flux_data, flux_error, flux_mask]), flux_head=flux_head, ) - outfiles.append("/".join([data_folder, Stokes_hdul[0].header["FILENAME"] + ".fits"])) + outfiles.append("/".join([data_folder, Stokes_hdul["I_STOKES"].header["FILENAME"] + ".fits"])) # Step 5: # crop to desired region of interest (roi)