Changeset 15664

Timestamp:

Nov 20, 2007, 3:04:54 PM (18 years ago)

Author:

Paul Price

Message:

Mask pixel if flat-field is not finite.

File:

• trunk/psModules/src/detrend/pmFlatField.c (modified) (1 diff)

trunk/psModules/src/detrend/pmFlatField.c

@@ -77 +77 @@
         for (int i = 0; i < inImage->numCols; i++) { \
             ps##TYPE flatValue = flatImage->data.TYPE[j + yOffset][i + xOffset]; \
-            if (flatValue <= 0.0 || (flatMask && flatMask->data.U8[j + yOffset][i + xOffset])) { \
+            if (!isfinite(flatValue) || flatValue <= 0.0 ||
+                (flatMask && flatMask->data.U8[j + yOffset][i + xOffset])) { \
                 if (inMask) { \
                     inMask->data.PS_TYPE_MASK_DATA[j][i] |= badFlat; \