Changeset 16911 for trunk/psModules/src/objects/pmPSFtry.c

Timestamp:

Mar 10, 2008, 2:49:14 PM (18 years ago)

Author:

Paul Price

Message:

The NANs in the pmPSFtry->metric vector were coming from
pmSourceMagnitudes with mode=PM_SOURCE_PHOT_INTERP, and some bad
pixels in the source postage stamp. The NANs were effectively grown
by the call to psImageShift, while the mask was not grown. This
results in NAN pixels being accumulated in pmSourcePhotometryAper,
resulting in a NAN aperMag, so that the metric is also NAN.

Another possible error I found while debugging is that
pmSourcePixelWeight and pmSourcePhotometryAper should be run with the
mask value set to "maskVal | mark", so that only the circular aperture
is used.

File:

• trunk/psModules/src/objects/pmPSFtry.c (modified) (3 diffs)

trunk/psModules/src/objects/pmPSFtry.c

@@ -5 +5 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.54 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-03-10 21:03:30 $
+ *  @version $Revision: 1.55 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-03-11 00:49:14 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -193 +193 @@
         psfTry->metric->data.F32[i] = source->apMag - source->psfMag;
         psfTry->metricErr->data.F32[i] = source->errMag;
+
         Npsf ++;
     }
@@ -281 +282 @@
 
     for (int i = 0; i < psfTry->sources->n; i++) {
-        if (!isfinite(psfTry->metric->data.F32[i]) || !isfinite(psfTry->metricErr->data.F32[i])) {
-            psfTry->mask->data.U8[i] = PSFTRY_MASK_ALL;
-            continue;
-        }
         if (psfTry->mask->data.U8[i] & PSFTRY_MASK_ALL)
             continue;