Changeset 9807 for trunk/psphot/src/psphotGrowthCurve.c
- Timestamp:
- Oct 31, 2006, 9:37:39 AM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotGrowthCurve.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotGrowthCurve.c
r9595 r9807 5 5 // the 'center' option 6 6 7 bool psphotGrowthCurve (pmReadout *readout, pmPSF *psf) { 7 // XXX add a option to turn off the curve-of-growth (ie, make the apMag = fitMag everywhere); 8 9 bool psphotGrowthCurve (pmReadout *readout, pmPSF *psf, bool ignore) { 8 10 9 11 // bool status; … … 31 33 pmModel *model = pmModelFromPSF (modelRef, psf); 32 34 35 // measure the fitMag for this model 36 pmSourcePhotometryModel (&fitMag, model); 37 psf->growth->fitMag = fitMag; 38 33 39 // generate working image for this source 34 40 psRegion region = {xc - dx, xc + dx, yc - dy, yc + dy}; … … 46 52 47 53 psImageKeepCircle (mask, xc, yc, radius, "OR", PM_MASK_MARK); 54 pmModelAdd (image, mask, model, false, false); 55 pmSourcePhotometryAper (&apMag, model, image, mask); 56 psImageKeepCircle (mask, xc, yc, radius, "AND", PS_NOT_U8(PM_MASK_MARK)); 48 57 49 pmModelAdd (image, mask, model, false, false); 50 51 pmSourcePhotometryAper (&apMag, model, image, mask);52 53 psImageKeepCircle (mask, xc, yc, radius, "AND", PS_NOT_U8(PM_MASK_MARK));54 psf->growth->apMag->data.F32[i] = apMag; 58 // the 'ignore' mode is for testing 59 if (ignore) { 60 psf->growth->apMag->data.F32[i] = fitMag; 61 } else { 62 psf->growth->apMag->data.F32[i] = apMag; 63 } 55 64 } 56 57 pmSourcePhotometryModel (&fitMag, model);58 psf->growth->fitMag = fitMag;59 65 60 66 psf->growth->apRef = psVectorInterpolate (psf->growth->radius, psf->growth->apMag, psf->growth->refRadius);
Note:
See TracChangeset
for help on using the changeset viewer.
