Changeset 34247 for branches/eam_branches/ipp-20120627/psphot/src/psphotStackMatchPSFs.c

Timestamp:

Jul 31, 2012, 11:50:18 AM (14 years ago)

Author:

eugene

Message:

merge changes from trunk

Location:

branches/eam_branches/ipp-20120627/psphot

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/psphotStackMatchPSFs.c (modified) (2 diffs)

branches/eam_branches/ipp-20120627/psphot/src/psphotStackMatchPSFs.c

@@ -87 +87 @@
     }
 
-    // copy the header data from Src to Out
-    // pmHDU *hduSrc = pmHDUFromReadout(readoutSrc);
-    // psAssert (hduSrc, "input missing hdu?");
-
-
     // set NAN pixels to 'SAT'
     psImageMaskType maskSat = pmConfigMaskGet("SAT", config);
@@ -116 +111 @@
     // save the output fwhm values in the readout->analysis.  we may have / will have multiple output PSF sizes,
     // so we save this in a vector.  if the vector is not yet defined, create it
+    // Skip this if the readout deconvolution fraction was over the limit.
     // NOTE: fwhmValues as defined here has 1 + nMatched PSF : 0 == unmatched
     psVector *fwhmValues = psVectorAllocEmpty(10, PS_TYPE_F32);
     psVectorAppend(fwhmValues, NAN); // XXX this corresponds to the unmatched image set
-    for (int i = 0; i < options->targetSeeing->n; i++) {
-        psVectorAppend(fwhmValues, options->targetSeeing->data.F32[i]);
+
+    bool overLimit = psMetadataLookupBool(NULL, readoutOut->analysis, "DECONV.OVERLIMIT");
+    if (!overLimit) {
+        for (int i = 0; i < options->targetSeeing->n; i++) {
+            psVectorAppend(fwhmValues, options->targetSeeing->data.F32[i]);
+        }
     }
     psMetadataAddVector(readoutSrc->analysis, PS_LIST_TAIL, "STACK.PSF.FWHM.VALUES", PS_META_REPLACE, "PSF sizes", fwhmValues);