def _run(self, img):
imgid = self.get_imgid(img)
_logger.debug('correcting flat in %s', imgid)
_logger.debug('flat mean is %f', self.flat_stats)
data = self.datamodel.get_data(img)
data = array.correct_flatfield(data, self.flatdata, dtype=self.dtype)
# FIXME
img[0].data = data
return img
def run(self, img):
imgid = self.get_imgid(img)
_logger.debug('correcting flat in %s', imgid)
_logger.debug('flat mean is %f', self.flat_stats)
data = self.datamodel.get_data(img)
data = array.correct_flatfield(data, self.flatdata, dtype=self.dtype)
# FIXME
img[0].data = data
hdr = img['primary'].header
hdr['NUM-FF'] = self.calibid
hdr['history'] = 'Flat-field correction with {}'.format(self.calibid)
hdr['history'] = 'Flat-field correction time {}'.format(datetime.datetime.utcnow().isoformat())
hdr['history'] = 'Flat-field correction mean {}'.format(self.flat_stats)
return img
def correct_superflat(self, frame, fitted, step=0, save=True):
frame.flat_corrected = name_skyflat_proc(frame.label, step)
if save:
shutil.copyfile(frame.resized_base, frame.flat_corrected)
else:
os.rename(frame.resized_base, frame.flat_corrected)
self.logger.info("Step %d, SF: apply superflat to frame %s", step,
frame.flat_corrected)
with fits.open(frame.flat_corrected, mode='update') as hdulist:
data = hdulist['primary'].data
datar = data[frame.valid_region]
data[frame.valid_region] = narray.correct_flatfield(datar, fitted)
frame.lastname = frame.flat_corrected
def correct_superflat(self, frame, fitted, step=0, save=True):
frame.flat_corrected = name_skyflat_proc(frame.label, step)
if save:
shutil.copyfile(frame.resized_base, frame.flat_corrected)
else:
os.rename(frame.resized_base, frame.flat_corrected)
self.logger.info("Step %d, SF: apply superflat to frame %s",
step, frame.flat_corrected)
with fits.open(frame.flat_corrected, mode='update') as hdulist:
data = hdulist['primary'].data
datar = data[frame.valid_region]
data[frame.valid_region] = narray.correct_flatfield(datar, fitted)
frame.lastname = frame.flat_corrected
def correct_superflat(self, frame, fitted, step=0, save=True):
frame.flat_corrected = name_skyflat_proc(frame.baselabel, step)
if save:
shutil.copyfile(frame.resized_base, frame.flat_corrected)
else:
os.rename(frame.resized_base, frame.flat_corrected)
_logger.info("Step %d, SF: apply superflat to frame %s", step,
frame.flat_corrected)
with fits.open(frame.flat_corrected, mode='update') as hdulist:
data = hdulist['primary'].data
datar = data[frame.valid_region]
data[frame.valid_region] = correct_flatfield(datar, fitted)
frame.lastname = frame.flat_corrected
# FIXME: plotting
try:
self.figure_image(data[frame.valid_region], frame)
except ValueError:
_logger.warning('Problem plotting %s', frame.lastname)