Changeset 25645
- Timestamp:
- Sep 28, 2009, 5:53:45 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/20090715/psphot/src/psphotMakeResiduals.c
r25531 r25645 174 174 bool offImage = false; 175 175 if (psImageInterpolate (&flux, &dflux, &mflux, ix, iy, interp) == PS_INTERPOLATE_STATUS_OFF) { 176 // fprintf (stderr, "off image: %f %f : %f %f\n", ix, iy, flux, dflux);177 176 // This pixel is off the image 178 177 offImage = true; … … 256 255 resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = 1; 257 256 } 258 259 // 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]);260 261 257 } else { 262 258 assert (SPATIAL_ORDER == 1); … … 294 290 295 291 float dRo = sqrt(A->data.F32[0][0]); 296 // fprintf (stderr, "res: %2d %2d : %6.4f %6.4f %6.4f %3d %1d\n",297 // ox, oy, resid->Ro->data.F32[oy][ox], dRo, dRo/sqrt(nKeep), nKeep, resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox]);298 292 299 293 if (fabs(resid->Ro->data.F32[oy][ox]) < pixelSN*dRo/sqrt(nKeep)) { 300 294 resid->mask->data.PM_TYPE_RESID_MASK_DATA[oy][ox] = 1; 301 295 } 302 //resid->variance->data.F32[oy][ox] = XXX;303 296 } 304 297 }
Note:
See TracChangeset
for help on using the changeset viewer.
