Changeset 15368 for trunk/ppStack/src/ppStackMatch.c

Timestamp:

Oct 23, 2007, 1:52:49 PM (19 years ago)

Author:

Paul Price

Message:

Using fake image instead of individual fake sources. Was losing most of the suitable stamp stars because they were too close to other stars (or artifacts). In order to get around this, need to generate all sources that might be near a footprint; it's easier to just do everything on the image.

File:

• trunk/ppStack/src/ppStackMatch.c (modified) (1 diff)

trunk/ppStack/src/ppStackMatch.c

@@ -58 +58 @@
     }
 
+    // Generate a fake image to match to
+    float maxMag = -INFINITY;           // Maximum magnitude of sources
+    for (int i = 0; i < sources->n; i++) {
+        pmSource *source = sources->data[i]; // Source of interest
+        if (source->psfMag > maxMag) {
+            maxMag = source->psfMag;
+        }
+    }
+    pmReadout *fake = pmReadoutFakeFromSources(input->image->numCols, input->image->numRows,
+                                               sources, target, powf(10.0, -0.4 * maxMag));
+
+
     // Do the image matching
-    if (!pmSubtractionMatch(output, input, NULL, footprint, regionSize, spacing, threshold, sources,
+    if (!pmSubtractionMatch(output, input, fake, footprint, regionSize, spacing, threshold, sources,
                             stampsName, target, type, size, order, widths, orders, inner, ringsOrder,
                             binning, optimum, optWidths, optOrder, optThresh, iter, rej, maskBad,
                             maskBlank, badFrac)) {
         psError(PS_ERR_UNKNOWN, false, "Unable to match images.");
+        psFree(fake);
         psFree(optWidths);
         return false;
     }
+    psFree(fake);
     psFree(optWidths);