IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jan 19, 2009, 4:57:24 PM (18 years ago)
Author:
eugene
Message:

bug fixes for mask types

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/eam_branch_20081230/psphot/src/psphotMakeResiduals.c

    r21090 r21144  
    194194                resid->Rx->data.F32[oy][ox] = 0.0;
    195195                resid->Ry->data.F32[oy][ox] = 0.0;
    196                 resid->mask->data.PS_TYPE_IMAGE_MASK_DATA[oy][ox] = badMask;
     196                resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = badMask;
    197197                continue;
    198198            }
     
    228228
    229229                if (fabs(resid->Ro->data.F32[oy][ox]) < pixelSN*fluxStats->sampleStdev/sqrt(nKeep)) {
    230                   resid->mask->data.PS_TYPE_IMAGE_MASK_DATA[oy][ox] = 1;
    231                 }
    232 
    233                 // fprintf (stderr, "res: %2d %2d : %6.4f  %6.4f  %6.4f   %3d  %1d\n", ox, oy, resid->Ro->data.F32[oy][ox], fluxStats->sampleStdev, fluxStats->sampleStdev/sqrt(nKeep), nKeep, resid->mask->data.PS_TYPE_IMAGE_MASK_DATA[oy][ox]);
     230                  resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = 1;
     231                }
     232
     233                // fprintf (stderr, "res: %2d %2d : %6.4f  %6.4f  %6.4f   %3d  %1d\n", ox, oy, resid->Ro->data.F32[oy][ox], fluxStats->sampleStdev, fluxStats->sampleStdev/sqrt(nKeep), nKeep, resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox]);
    234234
    235235            } else {
     
    268268
    269269                float dRo = sqrt(A->data.F32[0][0]);
    270                 // fprintf (stderr, "res: %2d %2d : %6.4f  %6.4f  %6.4f   %3d  %1d\n", ox, oy, resid->Ro->data.F32[oy][ox], dRo, dRo/sqrt(nKeep), nKeep, resid->mask->data.PS_TYPE_IMAGE_MASK_DATA[oy][ox]);
     270                // fprintf (stderr, "res: %2d %2d : %6.4f  %6.4f  %6.4f   %3d  %1d\n",
     271                // ox, oy, resid->Ro->data.F32[oy][ox], dRo, dRo/sqrt(nKeep), nKeep, resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox]);
    271272
    272273                if (fabs(resid->Ro->data.F32[oy][ox]) < pixelSN*dRo/sqrt(nKeep)) {
    273                   resid->mask->data.PS_TYPE_IMAGE_MASK_DATA[oy][ox] = 1;
     274                  resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = 1;
    274275                }
    275276                //resid->weight->data.F32[oy][ox] = XXX;
Note: See TracChangeset for help on using the changeset viewer.