- Timestamp:
- Apr 20, 2011, 11:03:36 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20110404/psphot/src/psphotMagnitudes.c
r31313 r31328 167 167 pmSource *source = (pmSource *) sources->data[i]; 168 168 169 bool saveTest = false; 170 psImage *testImage = NULL; 171 # if (0) 172 if ((fabs(source->peak->xf-3518) < 5) && (fabs(source->peak->yf-3178) < 5)) { 173 saveTest = true; 174 psRegion subregion = psRegionSet (source->peak->xf - 200, source->peak->xf + 200, source->peak->yf - 200, source->peak->yf + 200); 175 testImage = psImageSubset ((psImage *) source->pixels->parent, subregion); 176 } 177 # endif 178 179 if (saveTest) { 180 psphotSaveImage(NULL, testImage, "test.image.1.fits"); 181 } 182 169 183 // replace object in image 170 184 if (source->tmpFlags & PM_SOURCE_TMPF_SUBTRACTED) { 171 185 pmSourceAdd (source, PM_MODEL_OP_FULL, maskVal); 172 186 } 187 188 if (saveTest) { 189 psphotSaveImage(NULL, testImage, "test.image.2.fits"); 190 } 173 191 174 192 // clear the mask bit and set the circular mask pixels … … 188 206 // re-subtract the object, leave local sky 189 207 pmSourceSub (source, PM_MODEL_OP_FULL, maskVal); 208 209 if (saveTest) { 210 psphotSaveImage(NULL, testImage, "test.image.3.fits"); 211 } 190 212 191 213 if (backModel) {
Note:
See TracChangeset
for help on using the changeset viewer.
