def _run(self, img): imgid = self.get_imgid(img) _logger.debug('correcting sky in %s', imgid) _logger.debug('sky mean is %f', self.calib_stats) data = self.datamodel.get_data(img) data = array.correct_sky(data, self.skydata, dtype=self.dtype) # FIXME img[0].data = data return img
def run(self, img): imgid = self.get_imgid(img) if self.datamodel.do_sky_correction(img): _logger.debug('correcting sky in %s', imgid) _logger.debug('sky mean is %f', self.calib_stats) data = self.datamodel.get_data(img) data = array.correct_sky(data, self.skydata, dtype=self.dtype) # FIXME img[0].data = data hdr = img['primary'].header hdr['NUM-SK'] = self.calibid hdr['history'] = 'Sky subtraction with {}'.format(self.calibid) hdr['history'] = 'Sky subtraction time {}'.format(datetime.datetime.utcnow().isoformat()) hdr['history'] = 'Sky subtraction mean {}'.format(self.calib_stats) else: _logger.debug('skip sky correction in %s', imgid) return img