Changeset 35455 for trunk/psModules/src/imcombine/pmStack.c

Timestamp:

Apr 30, 2013, 6:10:41 PM (13 years ago)

Author:

watersc1

Message:

Changes to stacking code:

• New input FWHM clipping algorithm using mixture models
• New recipes/reductions that further reduce FWHM clipping
• pmSubtractionRejectStamps now fits flux/chi^{2 values in log-log space to prevent outliers biasing the fit (and causing images to be rejected due to a small number of bad stamps)}
• pmPSF now correctly chooses between useReff forms as appropriate
• pmPSF change fixes circularization issues in pmReadoutFakeFromSources
• pmPSFEnvelope somewhat more clear about circularization
• Changed debugging in pmSubtractionEquation to add information.

I do not see any conflicts with the ppSub code.

File:

• trunk/psModules/src/imcombine/pmStack.c (modified) (2 diffs)

trunk/psModules/src/imcombine/pmStack.c

@@ -41 +41 @@
 /* #define TEST_X 5745                       // x coordinate to examine */
 /* #define TEST_Y 5331                       // y coordinate to examine */
-#define TEST_X 25
-#define TEST_Y 25
+// #define TEST_X 972
+// #define TEST_Y 3213
+#define TEST_X 3289
+#define TEST_Y 4810
 #define TEST_RADIUS 0.5                 // Radius to examine
 # endif
@@ -752 +754 @@
     if (PS_SQR(x - TEST_X) + PS_SQR(y - TEST_Y) <= PS_SQR(TEST_RADIUS)) {
         for (int i = 0; i < numGood; i++) {
-            fprintf(stderr,"Input %d, pixel %d,%d (%" PRIu16 "): %f %f (%f) %f %f %d %x %x -> %x %x\n",
-                    i, x, y, pixelSources->data.U16[i], pixelData->data.F32[i], pixelVariances->data.F32[i],
-                    addVariance->data.F32[i], pixelWeights->data.F32[i], pixelExps->data.F32[i],
-                    pixelSuspects->data.U8[i], badMaskBits, suspectMaskBits, *badMask, *goodMask);
+            fprintf(stderr,"Input %d, pixel %d,%d (%" PRIu16 "): %f %f (%g) %g %f %d %x %x -> %x %x\n",
+                    i, x, y, pixelSources->data.U16[i],
+                    pixelData->data.F32[i], pixelVariances->data.F32[i],
+                    addVariance->data.F32[i],
+                    pixelWeights->data.F32[i], pixelExps->data.F32[i],
+                    pixelSuspects->data.U8[i],
+                    badMaskBits, suspectMaskBits,
+                    *badMask, *goodMask);
         }
     }