Changeset 4949 for trunk/psphot/src/psphotApplyPSF.c
- Timestamp:
- Sep 5, 2005, 8:44:00 PM (21 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotApplyPSF.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotApplyPSF.c
r4946 r4949 31 31 32 32 // this function specifies the radius at this the model hits the given flux 33 p sModelRadius modelRadius = psModelRadius_GetFunction (psf->type);33 pmModelRadius modelRadius = pmModelRadius_GetFunction (psf->type); 34 34 35 35 for (int i = 0; i < sources->n; i++) { 36 36 37 p sSource *source = sources->data[i];37 pmSource *source = sources->data[i]; 38 38 39 39 // skip non-astronomical objects (very likely defects) … … 43 43 44 44 // use the source moments, etc to guess basic model parameters 45 p sModel *modelFLT = pmSourceModelGuess (source, psf->type);45 pmModel *modelFLT = pmSourceModelGuess (source, psf->type); 46 46 47 47 // set PSF parameters for this model 48 p sModel *model = psModelFromPSF (modelFLT, psf);48 pmModel *model = pmModelFromPSF (modelFLT, psf); 49 49 x = model->params->data.F32[2]; 50 50 y = model->params->data.F32[3]; … … 66 66 67 67 // fit PSF model (set/unset the pixel mask) 68 psImageKeepCircle (source->mask, x, y, model->radius, OR, PSPHOT_MASK_KEEP);68 psImageKeepCircle (source->mask, x, y, model->radius, "or", PSPHOT_MASK_KEEP); 69 69 status = pmSourceFitModel (source, model, true); 70 psImageKeepCircle (source->mask, x, y, model->radius, AND, ~PSPHOT_MASK_KEEP);70 psImageKeepCircle (source->mask, x, y, model->radius, "and", ~PSPHOT_MASK_KEEP); 71 71 if (!status || (model->params->data.F32[1] < 0)) { 72 72 psLogMsg ("psphot", 3, "PSF fit failed for %f, %f (%d iterations)\n", x, y, model->nIter);
Note:
See TracChangeset
for help on using the changeset viewer.
