Changeset 4949 for trunk/psphot/src/psphotFitGalaxies.c
- Timestamp:
- Sep 5, 2005, 8:44:00 PM (21 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotFitGalaxies.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotFitGalaxies.c
r4946 r4949 21 21 22 22 char *modelName = psMetadataLookupPtr (&status, config, "GAL_MODEL"); 23 p sModelType modelType = psModelSetType (modelName);24 p sModelRadius modelRadius = psModelRadius_GetFunction (modelType);23 pmModelType modelType = pmModelSetType (modelName); 24 pmModelRadius modelRadius = pmModelRadius_GetFunction (modelType); 25 25 26 26 psTimerStart ("psphot"); 27 27 28 28 for (int i = 0; i < sources->n; i++) { 29 p sSource *source = sources->data[i];29 pmSource *source = sources->data[i]; 30 30 31 31 // sources which should not be fitted … … 50 50 51 51 // recalculate the source moments using the larger galaxy moments radius 52 status = pmSourceMoments _EAM(source, GAL_MOMENTS_RAD);52 status = pmSourceMoments (source, GAL_MOMENTS_RAD); 53 53 if (!status) { 54 54 source->type = PS_SOURCE_DROP_GALAXY; // better choice? … … 57 57 58 58 // use the source moments, etc to guess basic model parameters 59 p sModel *model = pmSourceModelGuess (source, modelType);59 pmModel *model = pmSourceModelGuess (source, modelType); 60 60 x = model->params->data.F32[2]; 61 61 y = model->params->data.F32[3]; … … 75 75 // fit FLT (not PSF) model (set/unset the pixel mask) 76 76 psImageKeepCircle (source->mask, x, y, model->radius, OR, PSPHOT_MASK_KEEP); 77 status = pmSourceFitModel _EAM(source, model, false);77 status = pmSourceFitModel (source, model, false); 78 78 psImageKeepCircle (source->mask, x, y, model->radius, AND, ~PSPHOT_MASK_KEEP); 79 79 if (!status) {
Note:
See TracChangeset
for help on using the changeset viewer.
