Changeset 14513 for trunk/psModules/src/imcombine/pmSubtraction.c

Timestamp:

Aug 15, 2007, 1:17:56 PM (19 years ago)

Author:

Paul Price

Message:

Adding function to set stamps based on a list of x,y positions (as opposed to find a set of stamps automatically). Adding another function to generate target stamps so that we can convolve to the target PSF.

File:

• trunk/psModules/src/imcombine/pmSubtraction.c (modified) (4 diffs)

trunk/psModules/src/imcombine/pmSubtraction.c

@@ -4 +4 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.36 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-08-14 02:14:30 $
+ *  @version $Revision: 1.37 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-08-15 23:17:56 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -542 +542 @@
                 stamp->status = PM_SUBTRACTION_STAMP_REJECTED;
                 psTrace("psModules.imcombine", 3, "Rejecting stamp %d (%d,%d) because of bad equation\n",
-                        i, stamp->x, stamp->y);
+                        i, (int)(stamp->x + 0.5), (int)(stamp->y + 0.5));
             } else {
                 stamp->status = PM_SUBTRACTION_STAMP_USED;
@@ -718 +718 @@
             deviations->data.F32[i] = sqrt(stats->sampleMean / 2.0);
             psTrace("psModules.imcombine", 1, "Deviation for stamp %d (%d,%d): %f\n",
-                    i, stamp->x, stamp->y, deviations->data.F32[i]);
+                    i, (int)(stamp->x + 0.5), (int)(stamp->y + 0.5), deviations->data.F32[i]);
             totalSquareDev += PS_SQR(deviations->data.F32[i]);
             numStamps++;
@@ -736 +736 @@
         if (stamp->status == PM_SUBTRACTION_STAMP_USED && deviations->data.F32[i] > limit) {
             // Mask out the stamp in the image so you it's not found again
-            psTrace("psModules.imcombine", 3, "Rejecting stamp %d (%d,%d)\n", i, stamp->x, stamp->y);
+            psTrace("psModules.imcombine", 3, "Rejecting stamp %d (%d,%d)\n", i,
+                    (int)(stamp->x + 0.5), (int)(stamp->y + 0.5));
             numRejected++;
             for (int y = stamp->y - footprint; y <= stamp->y + footprint; y++) {