- Timestamp:
- Jan 19, 2009, 4:57:24 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branch_20081230/psphot/src/psphotMakeResiduals.c
r21090 r21144 194 194 resid->Rx->data.F32[oy][ox] = 0.0; 195 195 resid->Ry->data.F32[oy][ox] = 0.0; 196 resid->mask->data.P S_TYPE_IMAGE_MASK_DATA[oy][ox] = badMask;196 resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = badMask; 197 197 continue; 198 198 } … … 228 228 229 229 if (fabs(resid->Ro->data.F32[oy][ox]) < pixelSN*fluxStats->sampleStdev/sqrt(nKeep)) { 230 resid->mask->data.P S_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.P S_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]); 234 234 235 235 } else { … … 268 268 269 269 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]); 271 272 272 273 if (fabs(resid->Ro->data.F32[oy][ox]) < pixelSN*dRo/sqrt(nKeep)) { 273 resid->mask->data.P S_TYPE_IMAGE_MASK_DATA[oy][ox] = 1;274 resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = 1; 274 275 } 275 276 //resid->weight->data.F32[oy][ox] = XXX;
Note:
See TracChangeset
for help on using the changeset viewer.
