Changeset 21519 for trunk/psphot/src/psphotSourceSize.c
- Timestamp:
- Feb 16, 2009, 12:34:24 PM (17 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotSourceSize.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotSourceSize.c
r21366 r21519 56 56 57 57 // source must have been subtracted 58 if (!(source->mode & PM_SOURCE_MODE_SUBTRACTED)) { 58 if (!(source->tmpFlags & PM_SOURCE_TMPF_SUBTRACTED)) { 59 source->mode |= PM_SOURCE_MODE_SIZE_SKIPPED; 59 60 psTrace("psphot", 7, "Not calculating extNsigma,crNsigma since source is not subtracted\n"); 60 61 continue; … … 63 64 psF32 **resid = source->pixels->data.F32; 64 65 psF32 **variance = source->variance->data.F32; 65 psImageMaskType **mask = source->maskObj->data.PS_TYPE_IMAGE_MASK_DATA;66 psImageMaskType **mask = source->maskObj->data.PS_TYPE_IMAGE_MASK_DATA; 66 67 67 68 // check for extendedness: measure the delta flux significance at the 1 sigma contour … … 82 83 if (xPeak < 1 || xPeak > source->pixels->numCols - 2 || 83 84 yPeak < 1 || yPeak > source->pixels->numRows - 2) { 85 source->mode |= PM_SOURCE_MODE_SIZE_SKIPPED; 84 86 psTrace("psphot", 7, "Not calculating crNsigma due to edge\n"); 85 87 continue; … … 98 100 if (!keep) { 99 101 psTrace("psphot", 7, "Not calculating crNsigma due to masked pixels\n"); 102 source->mode |= PM_SOURCE_MODE_SIZE_SKIPPED; 100 103 continue; 101 104 } … … 156 159 } 157 160 source->crNsigma = (nCR > 0) ? fCR / nCR : 0.0; 158 if (!isfinite(source->crNsigma)) continue; 161 if (!isfinite(source->crNsigma)) { 162 continue; 163 } 159 164 160 165 // this source is thought to be a cosmic ray. flag the detection and mask the pixels … … 264 269 // replace the source flux 265 270 pmSourceAdd (source, PM_MODEL_OP_FULL, maskVal); 266 source-> mode &= ~PM_SOURCE_MODE_SUBTRACTED;271 source->tmpFlags &= ~PM_SOURCE_TMPF_SUBTRACTED; 267 272 268 273 // flag this as a CR
Note:
See TracChangeset
for help on using the changeset viewer.
