Changeset 21183 for trunk/psphot/src/psphotFitSourcesLinear.c
- Timestamp:
- Jan 26, 2009, 8:40:07 PM (17 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotFitSourcesLinear.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotFitSourcesLinear.c
r20938 r21183 10 10 // the analysis is performed wrt the simulated pixel values 11 11 12 static bool SetBorderMatrixElements (psSparseBorder *border, pmReadout *readout, psArray *sources, bool constant_weights, int SKY_FIT_ORDER, ps MaskType markVal);12 static bool SetBorderMatrixElements (psSparseBorder *border, pmReadout *readout, psArray *sources, bool constant_weights, int SKY_FIT_ORDER, psImageMaskType markVal); 13 13 14 14 bool psphotFitSourcesLinear (pmReadout *readout, psArray *sources, psMetadata *recipe, pmPSF *psf, bool final) { … … 23 23 24 24 // bit-masks to test for good/bad pixels 25 ps MaskType maskVal = psMetadataLookupU8(&status, recipe, "MASK.PSPHOT");25 psImageMaskType maskVal = psMetadataLookupImageMask(&status, recipe, "MASK.PSPHOT"); 26 26 assert (maskVal); 27 27 28 28 // bit-mask to mark pixels not used in analysis 29 ps MaskType markVal = psMetadataLookupU8(&status, recipe, "MARK.PSPHOT");29 psImageMaskType markVal = psMetadataLookupImageMask(&status, recipe, "MARK.PSPHOT"); 30 30 assert (markVal); 31 31 … … 233 233 // on the pixels which correspond to all of the sources of interest. These elements fill in 234 234 // the border matrix components in the sparse matrix equation. 235 static bool SetBorderMatrixElements (psSparseBorder *border, pmReadout *readout, psArray *sources, bool constant_weights, int SKY_FIT_ORDER, ps MaskType markVal) {235 static bool SetBorderMatrixElements (psSparseBorder *border, pmReadout *readout, psArray *sources, bool constant_weights, int SKY_FIT_ORDER, psImageMaskType markVal) { 236 236 237 237 // generate the image-wide weight terms … … 248 248 float x = model->params->data.F32[PM_PAR_XPOS]; 249 249 float y = model->params->data.F32[PM_PAR_YPOS]; 250 psImageMaskCircle (source->maskView, x, y, model->radiusFit, "AND", PS_NOT_ U8(markVal));250 psImageMaskCircle (source->maskView, x, y, model->radiusFit, "AND", PS_NOT_IMAGE_MASK(markVal)); 251 251 } 252 252 … … 254 254 psF32 **image = readout->image->data.F32; 255 255 psF32 **weight = readout->weight->data.F32; 256 ps U8 **mask = readout->mask->data.U8;256 psImageMaskType **mask = readout->mask->data.PS_TYPE_IMAGE_MASK_DATA; 257 257 258 258 double w, x, y, x2, xy, y2, xc, yc, wt, f, fo, fx, fy; … … 288 288 289 289 // turn off MARK for all image pixels 290 psImageMaskRegion (readout->mask, fullArray, "AND", PS_NOT_ U8(markVal));290 psImageMaskRegion (readout->mask, fullArray, "AND", PS_NOT_IMAGE_MASK(markVal)); 291 291 292 292 // set the Border T elements
Note:
See TracChangeset
for help on using the changeset viewer.
