Changeset 19268

Timestamp:

Aug 28, 2008, 2:48:22 PM (18 years ago)

Author:

Paul Price

Message:

Moving the inspect/reject debugging mask output to here, so that the
entire image can be written instead of just chunks (as was done by
having it in ppStackReadout). Also ensure masked regions really look
masked, via pmReadoutMaskApply.

File:

• trunk/ppStack/src/ppStackLoop.c (modified) (3 diffs)

trunk/ppStack/src/ppStackLoop.c

@@ -12 +12 @@
 #include "ppStack.h"
 
-//#define TESTING
+#define TESTING
 
 #define WCS_TOLERANCE 0.001             // Tolerance for WCS
@@ -619 +619 @@
             psLogMsg("ppStack", PS_LOG_INFO, "%ld total pixels to inspect from image %d", inspect->n, i);
 
+#ifdef TESTING
+            {
+                psImage *pix = psPixelsToMask(NULL, inspect, psRegionSet(0, outRO->image->numCols - 1,
+                                                                         0, outRO->image->numRows - 1), 0xff);
+                psString name = NULL;           // Name of image
+                psStringAppend(&name, "inspect_%03d.fits", i);
+                psFits *fits = psFitsOpen(name, "w");
+                psFree(name);
+                psFitsWriteImage(fits, NULL, pix, 0, NULL);
+                psFree(pix);
+                psFitsClose(fits);
+            }
+#endif
+
             psPixels *reject = pmStackReject(inspect, NULL, threshold, subRegions->data[i],
                                              subKernels->data[i]); // Pixels to reject
+
+#ifdef TESTING
+            {
+                psImage *pix = psPixelsToMask(NULL, reject, psRegionSet(0, outRO->image->numCols - 1,
+                                                                        0, outRO->image->numRows - 1), 0xff);
+                psString name = NULL;           // Name of image
+                psStringAppend(&name, "reject_%03d.fits", i);
+                psFits *fits = psFitsOpen(name, "w");
+                psFree(name);
+                psFitsWriteImage(fits, NULL, pix, 0, NULL);
+                psFree(pix);
+                psFitsClose(fits);
+            }
+#endif
+
             psFree(inspect);
             if (!reject) {
@@ -727 +756 @@
         }
         psFree(readouts);
+
+        // Ensure masked regions really look masked
+        if (!pmReadoutMaskApply(outRO, maskBad)) {
+            psWarning("Unable to apply mask");
+        }
 
         // Close up