isr: mark amp as bad if entirely saturated

If the amp is completely saturated, overscan subtraction may fail.
lsst · Feb 16, 2017 · 479ed1d · 479ed1d
1 parent 105b842
commit 479ed1d
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/python/lsst/obs/subaru/isr.py b/python/lsst/obs/subaru/isr.py
@@ -247,6 +247,11 @@ def runDataRef(self, sensorRef):
                 self.saturationDetection(ccdExposure, amp)
             if self.config.doSuspect and not badAmp:
                 self.suspectDetection(ccdExposure, amp)
+            ampMask = afwImage.MaskU(ccdExposure.getMaskedImage().getMask(), amp.getRawDataBBox(),
+                                     afwImage.PARENT)
+            maskVal = ampMask.getPlaneBitMask([self.config.saturatedMaskName, self.config.suspectMaskName])
+            if numpy.all(ampMask.getArray() & maskVal > 0):
+                badAmp = True
             if self.config.doOverscan and not badAmp:
                 ampImage = afwImage.MaskedImageF(ccdExposure.getMaskedImage(), amp.getRawDataBBox(),
                                                  afwImage.PARENT)