Changeset 25613
- Timestamp:
- Sep 27, 2009, 11:10:45 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/20090715/psphot/src/psphotReadout.c
r25533 r25613 124 124 return psphotReadoutCleanup (config, readout, recipe, detections, psf, sources); 125 125 } 126 127 126 psphotVisualShowPSFModel (readout, psf); 128 127 … … 130 129 psphotLoadExtSources (config, view, sources); 131 130 132 // construct an initial model for each object 131 // construct an initial model for each object, set the radius to fitRadius, set circular fit mask 133 132 psphotGuessModels (config, readout, sources, psf); 134 133 135 // XXX test output of models 136 // psphotTestSourceOutput (readout, sources, recipe, psf); 137 138 // linear PSF fit to source peaks 134 // linear PSF fit to source peaks, subtract the models from the image (in PSF mask) 139 135 psphotFitSourcesLinear (readout, sources, recipe, psf, FALSE); 140 136 … … 151 147 152 148 // non-linear PSF and EXT fit to brighter sources 149 // replace model flux, adjust mask as needed, fit, subtract the models (full stamp) 153 150 psphotBlendFit (config, readout, sources, psf); 154 151 … … 156 153 psphotReplaceAllSources (sources, recipe); 157 154 158 // linear fit to include all sources 155 // linear fit to include all sources (subtract again) 159 156 psphotFitSourcesLinear (readout, sources, recipe, psf, TRUE); 160 157 … … 163 160 goto pass1finish; 164 161 } 165 166 // XXX for the moment, drop the re-calc of the background (prove this works) 167 // replace background in residual image 168 // psphotSkyReplace (config, view); 169 // re-measure background model (median, smoothed image) 170 // psphotImageMedian (config, view); 162 // NOTE: possibly re-measure background model here with objects subtracted 171 163 172 164 // add noise for subtracted objects … … 188 180 } 189 181 190 // create full input models 182 // create full input models, set the radius to fitRadius, set circular fit mask 191 183 psphotGuessModels (config, readout, newSources, psf); 192 184 … … 212 204 finish: 213 205 214 // XXX drop this : test dev version215 // psphotPetrosianAnalysis (readout, sources, recipe);216 217 206 // plot positive sources 218 207 // psphotSourcePlots (readout, sources, recipe);
Note:
See TracChangeset
for help on using the changeset viewer.
