Changeset 26305 for branches/pap/ppStack/src/ppStackMatch.c

Timestamp:

Dec 2, 2009, 3:49:31 AM (16 years ago)

Author:

Paul Price

Message:

Dump sources before generating fake image.

File:

• branches/pap/ppStack/src/ppStackMatch.c (modified) (4 diffs)

branches/pap/ppStack/src/ppStackMatch.c

@@ -18 +18 @@
 #define COVAR_FRAC 0.01                 // Truncation fraction for covariance matrix
 
-//#define TESTING                         // Enable debugging output
+#define TESTING                         // Enable debugging output
 
 #ifdef TESTING
@@ -64 +64 @@
     return;
 }
+
+#ifdef TESTING
+static void dumpSources(const char *root, // Root name for source file
+                        psArray *sources// Sources to dump
+    )
+{
+    static int number = 0;              // Index of sources
+    psString name = NULL;               // Output name
+    psStringAppend(&name, "%s_%d.sources", root, number);
+    FILE file = fopen(name, "w");       // File to write
+    psFree(name);
+    for (int i = 0; i < sources->n; i++) {
+        pmSource *source = sources->data[i]; // Source of interest
+        if (!source) {
+            continue;
+        }
+        float x, y;                     // Coordinates of source
+        coordsFromSource(&x, &y, source);
+        float mag = source->psfMag;     // Magnitude of source
+        fprintf(file, "%f %f %f\n", x, y, mag);
+    }
+    fclose(file);
+    return;
+}
+#endif
+
 
 static psArray *stackSourcesFilter(psArray *sources, // Source list to filter
@@ -343 +369 @@
             psArray *stampSources = stackSourcesFilter(options->sourceLists->data[index],
                                                        footprint); // Filtered list of sources
-
+#ifdef TESTING
+            dumpSources("fakeSources", stampSources);
+#endif
             if (!pmReadoutFakeFromSources(fake, readout->image->numCols, readout->image->numRows,
                                           stampSources, SOURCE_MASK, NULL, NULL, options->psf,
@@ -356 +384 @@
             fake->mask = psImageCopy(NULL, readout->mask, PS_TYPE_IMAGE_MASK);
 
+#if 0
             // Add the background into the target image
+            // XXX Is this desirable?  There should be zero background on the input!
             psImage *bgImage = stackBackgroundModel(readout, config); // Image of background
             psBinaryOp(fake->image, fake->image, "+", bgImage);
             psFree(bgImage);
+#endif
 
 #ifdef TESTING