Changeset 6056 for trunk/psphot/src/psphotSourceFits.c
- Timestamp:
- Jan 19, 2006, 2:55:33 PM (21 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotSourceFits.c (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotSourceFits.c
r5993 r6056 4 4 // given a source with an existing modelPSF, attempt a full PSF fit, subtract if successful 5 5 6 bool psphotFitPSF ( eamReadout *imdata, pmSource *source) {6 bool psphotFitPSF (pmReadout *readout, pmSource *source) { 7 7 8 8 float x, y; … … 12 12 13 13 // extend source radius as needed 14 psphotCheckRadiusPSF ( imdata, source, PSF);14 psphotCheckRadiusPSF (readout, source, PSF); 15 15 16 16 x = PSF->params->data.F32[2]; … … 61 61 } 62 62 63 bool psphotFitBlob ( eamReadout *imdata, pmSource *source, psArray *sources) {63 bool psphotFitBlob (pmReadout *readout, pmSource *source, psArray *sources) { 64 64 65 65 bool okEXT, okDBL; … … 79 79 pmSource *tmpSrc = pmSourceAlloc (); 80 80 81 pmModel *EXT = psphotFitEXT ( imdata, source);81 pmModel *EXT = psphotFitEXT (readout, source); 82 82 okEXT = psphotEvalEXT (tmpSrc, EXT); 83 83 chiEXT = EXT->chisq / EXT->nDOF; 84 84 85 psArray *DBL = psphotFitDBL ( imdata, source);85 psArray *DBL = psphotFitDBL (readout, source); 86 86 okDBL = psphotEvalDBL (tmpSrc, DBL->data[0]); 87 87 okDBL &= psphotEvalDBL (tmpSrc, DBL->data[1]); … … 129 129 130 130 // fit a double PSF source to an extended blob 131 psArray *psphotFitDBL ( eamReadout *imdata, pmSource *source) {131 psArray *psphotFitDBL (pmReadout *readout, pmSource *source) { 132 132 133 133 float x, y, dx, dy; … … 149 149 // save the PSF model from the Ensemble fit 150 150 PSF = source->modelPSF; 151 psphotCheckRadiusPSF ( imdata, source, PSF);151 psphotCheckRadiusPSF (readout, source, PSF); 152 152 153 153 modelSet = psArrayAlloc (2); … … 176 176 } 177 177 178 pmModel *psphotFitEXT ( eamReadout *imdata, pmSource *source) {178 pmModel *psphotFitEXT (pmReadout *readout, pmSource *source) { 179 179 180 180 float x, y; … … 183 183 pmModel *EXT = pmSourceModelGuess (source, modelTypeEXT); 184 184 185 psphotCheckRadiusEXT ( imdata, source, EXT);185 psphotCheckRadiusEXT (readout, source, EXT); 186 186 187 187 x = EXT->params->data.F32[2]; … … 196 196 } 197 197 198 bool psphotFitBlend ( eamReadout *imdata, pmSource *source) {198 bool psphotFitBlend (pmReadout *readout, pmSource *source) { 199 199 200 200 float x, y, dR; … … 202 202 // if this source is not a possible blend, just fit as PSF 203 203 if ((source->blends == NULL) || (source->mode & PM_SOURCE_SATSTAR)) { 204 bool status = psphotFitPSF ( imdata, source);204 bool status = psphotFitPSF (readout, source); 205 205 return status; 206 206 } … … 212 212 213 213 // extend source radius as needed 214 psphotCheckRadiusPSF ( imdata, source, PSF);214 psphotCheckRadiusPSF (readout, source, PSF); 215 215 216 216 x = PSF->params->data.F32[2];
Note:
See TracChangeset
for help on using the changeset viewer.
