Changeset 15932
- Timestamp:
- Dec 26, 2007, 4:09:10 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotImageMedian.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotImageMedian.c
r15127 r15932 1 1 # include "psphotInternal.h" 2 2 static int npass = 0; 3 4 // we have 4 possibilities: (INTERNAL or I/O file) and (exists or not)5 // select model pixels (from output background model file, or create internal file)6 static pmReadout *get_model_readout(const char *name, // name of internal/external file7 const pmConfig *config, // configuration information8 const pmFPAview *view,9 pmFPA *inFPA,10 const psImageBinning *binning) {11 pmReadout *model = NULL;12 13 bool status = true;14 pmFPAfile *file = psMetadataLookupPtr(&status, config->files, name);15 if (file == NULL) {16 // we are not using PSPHOT.BACKMDL as an I/O file: define an internal version17 model = pmFPAfileDefineInternal (config->files, name, binning->nXruff, binning->nYruff, PS_TYPE_F32);18 } else {19 if (file->mode == PM_FPA_MODE_INTERNAL) {20 // we are not using PSPHOT.BACKMDL as an I/O file: already defined above21 model = file->readout;22 } else {23 // we are using PSPHOT.BACKMDL as an I/O file: select readout or create24 model = pmFPAviewThisReadout (view, file->fpa);25 if (model == NULL) {26 // readout does not yet exist: create from input27 // XXX we have an inconsistency in this calculation here and in pmFPACopy28 // XXX use the psImageBinning functions to set the output image size29 pmFPAfileCopyStructureView (file->fpa, inFPA, binning->nXbin, binning->nYbin, view);30 model = pmFPAviewThisReadout (view, file->fpa);31 PS_ASSERT (binning->nXruff == model->image->numCols, false);32 PS_ASSERT (binning->nYruff == model->image->numRows, false);33 }34 }35 }36 37 return model;38 }39 3 40 4 // generate the median in NxN boxes, clipping heavily … … 165 129 166 130 // we save the binning structure for use in psphotMagnitudes 167 pmReadout *model = get_model_readout("PSPHOT.BACKMDL", config, view,inFPA, binning);168 pmReadout *modelStdev = get_model_readout("PSPHOT.BACKMDL.STDEV", config, view, inFPA, binning);131 pmReadout *model = pmFPAGenerateReadout (config, view, "PSPHOT.BACKMDL", inFPA, binning); 132 pmReadout *modelStdev = pmFPAGenerateReadout (config, view, "PSPHOT.BACKMDL.STDEV", inFPA, binning); 169 133 170 134 psF32 **modelData = model->image->data.F32;
Note:
See TracChangeset
for help on using the changeset viewer.
