Changeset 16602 for branches/pap_branch_080207/ppStack/src/ppStackReadout.c

Timestamp:

Feb 22, 2008, 9:20:13 AM (18 years ago)

Author:

Paul Price

Message:

Fixing debugging outputs so that the full range is output (incremental
reads were overwriting each other).

File:

• branches/pap_branch_080207/ppStack/src/ppStackReadout.c (modified) (6 diffs)

branches/pap_branch_080207/ppStack/src/ppStackReadout.c

@@ -10 +10 @@
 #include "ppStack.h"
 
-
 #define WCS_TOLERANCE 0.001             // Tolerance for WCS
 
 //#define REJECTION_FILES                 // Write rejection mask?
 //#define INSPECTION_FILES                // Write inspection mask?
+
+static int sectionNum = 0;              // Section number; for debugging outputs
+
 
 bool ppStackReadout(const pmConfig *config, pmReadout *outRO, const psArray *readouts,
@@ -26 +28 @@
     assert(readouts->n == regions->n);
     assert(regions->n == kernels->n);
+
 
     // Get the recipe values
@@ -113 +116 @@
 #ifdef INSPECTION_FILES
     {
-        psFits *fits = psFitsOpen("combined.fits", "w");
+        psString name = NULL;           // Name of image
+        psStringAppend(&name, "combined_%03d.fits", sectionNum);
+        psFits *fits = psFitsOpen(name, "w");
+        psFree(name);
         psFitsWriteImage(fits, NULL, outRO->image, 0, NULL);
         psFitsClose(fits);
@@ -121 +127 @@
         pmStackData *data = stack->data[i]; // Data for this image
         psImage *inspected = psPixelsToMask(NULL, data->pixels,
-                                            psRegionSet(0, data->readout->image->numCols,
-                                                        0, data->readout->image->numRows),
+                                            psRegionSet(0, data->readout->image->numCols - 1,
+                                                        0, data->readout->image->numRows - 1),
                                             maskBlank);
         psString name = NULL;           // Name of image
-        psStringAppend(&name, "inspect%03d.fits", i);
+        psStringAppend(&name, "inspect_%03d_%03d.fits", sectionNum, i);
         psFits *fits = psFitsOpen(name, "w");
         psFree(name);
@@ -156 +162 @@
         }
         psImage *rejected = psPixelsToMask(NULL, data->pixels,
-                                           psRegionSet(0, data->readout->image->numCols,
-                                                       0, data->readout->image->numRows),
+                                           psRegionSet(0, data->readout->image->numCols - 1,
+                                                       0, data->readout->image->numRows - 1),
                                            maskBlank);
         psString name = NULL;           // Name of image
-        psStringAppend(&name, "reject%03d.fits", i);
+        psStringAppend(&name, "reject_%03d_%03d.fits", sectionNum, i);
         psFits *fits = psFitsOpen(name, "w");
         psFree(name);
@@ -194 +200 @@
     psFree(stack);
 
+    sectionNum++;
+
     return true;
 }