From 22db4d7f25f7b1e4cf90c540bf3fe4288621bac6 Mon Sep 17 00:00:00 2001 From: Julien Guy Date: Thu, 30 Jan 2025 15:37:14 -0800 Subject: [PATCH] cfinder as argument of subtract_peramp_overscan --- py/desispec/ccdcalib.py | 9 +++++++-- 1 file changed, 7 insertions(+), 2 deletions(-) diff --git a/py/desispec/ccdcalib.py b/py/desispec/ccdcalib.py index 03fe64e0e..384ed6d4b 100755 --- a/py/desispec/ccdcalib.py +++ b/py/desispec/ccdcalib.py @@ -280,7 +280,7 @@ def compute_bias_file(rawfiles, outfile, camera, explistfile=None, image=fitsfile[camera].data.astype("float64") - subtract_peramp_overscan(image, image_header) + subtract_peramp_overscan(image, image_header, cfinder) if shape is None : shape=image.shape @@ -747,9 +747,14 @@ def compare_bias(rawfile, biasfile1, biasfile2, ny=8, nx=40): image, hdr = fitsio.read(rawfile, ext=cam1, header=True) + primary_hdr = None + with fitsio.FITS(rawfile) as fx: + primary_hdr = io.raw.read_raw_primary_header(fx) + cfinder = CalibFinder([primary_hdr, hdr],fallback_on_dark_not_found=True) + #- subtract constant per-amp overscan region image = image.astype(float) - subtract_peramp_overscan(image, hdr) + subtract_peramp_overscan(image, hdr, cfinder) #- calculate differences per-amp readout_mode = get_readout_mode(hdr)