Changeset 28909
- Timestamp:
- Aug 13, 2010, 12:18:27 PM (16 years ago)
- Location:
- branches/eam_branches/ipp-20100621/psModules/src/objects
- Files:
-
- 14 edited
-
pmSource.c (modified) (1 diff)
-
pmSource.h (modified) (1 diff)
-
pmSourceIO_CMF_PS1_DV1.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_DV2.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_SV1.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_V1.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_V2.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_V3.c (modified) (2 diffs)
-
pmSourceIO_PS1_CAL_0.c (modified) (2 diffs)
-
pmSourceIO_PS1_DEV_0.c (modified) (2 diffs)
-
pmSourceIO_PS1_DEV_1.c (modified) (2 diffs)
-
pmSourceIO_RAW.c (modified) (2 diffs)
-
pmSourceIO_SMPDATA.c (modified) (2 diffs)
-
pmSourcePhotometry.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSource.c
r28860 r28909 140 140 source->sky = NAN; 141 141 source->skyErr = NAN; 142 source->pixWeight = NAN; 142 source->pixWeightNotBad = NAN; 143 source->pixWeightNotPoor = NAN; 143 144 144 145 source->psfChisq = NAN; -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSource.h
r28904 r28909 77 77 float apRadius; ///< radius for aperture magnitude 78 78 79 float p sfWeightNotBad; ///< PSF-weighted coverage of unmasked (not BAD) pixels80 float p sfWeightNotPoor; ///< PSF-weighted coverage of unmasked (not POOR) pixels79 float pixWeightNotBad; ///< PSF-weighted coverage of unmasked (not BAD) pixels 80 float pixWeightNotPoor; ///< PSF-weighted coverage of unmasked (not POOR) pixels 81 81 82 82 float psfChisq; ///< probability of PSF -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_CMF_PS1_DV1.c
r28860 r28909 202 202 psMetadataAdd (row, PS_LIST_TAIL, "PSF_MINOR", PS_DATA_F32, "PSF width (minor axis)", axes.minor); 203 203 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 204 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );204 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 205 205 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NDOF", PS_DATA_S32, "degrees of freedom", nDOF); 206 206 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NPIX", PS_DATA_S32, "number of pixels in fit", nPix); … … 342 342 } 343 343 344 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");344 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 345 345 source->crNsigma = psMetadataLookupF32 (&status, row, "CR_NSIGMA"); 346 346 source->extNsigma = psMetadataLookupF32 (&status, row, "EXT_NSIGMA"); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_CMF_PS1_DV2.c
r28860 r28909 202 202 psMetadataAdd (row, PS_LIST_TAIL, "PSF_MINOR", PS_DATA_F32, "PSF width (minor axis)", axes.minor); 203 203 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 204 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight); 204 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor (bad)", source->pixWeightNotBad); 205 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF_PERFECT", PS_DATA_F32, "PSF coverage/quality factor (poor)", source->pixWeightNotPoor); 205 206 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NDOF", PS_DATA_S32, "degrees of freedom", nDOF); 206 207 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NPIX", PS_DATA_S32, "number of pixels in fit", nPix); … … 348 349 } 349 350 350 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF"); 351 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 352 source->pixWeightNotPoor = psMetadataLookupF32 (&status, row, "PSF_QF_PERFECT"); 351 353 source->crNsigma = psMetadataLookupF32 (&status, row, "CR_NSIGMA"); 352 354 source->extNsigma = psMetadataLookupF32 (&status, row, "EXT_NSIGMA"); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_CMF_PS1_SV1.c
r28860 r28909 202 202 psMetadataAdd (row, PS_LIST_TAIL, "PSF_MINOR", PS_DATA_F32, "PSF width (minor axis)", axes.minor); 203 203 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 204 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );204 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 205 205 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NDOF", PS_DATA_S32, "degrees of freedom", nDOF); 206 206 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NPIX", PS_DATA_S32, "number of pixels in fit", nPix); … … 371 371 } 372 372 373 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");373 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 374 374 source->crNsigma = psMetadataLookupF32 (&status, row, "CR_NSIGMA"); 375 375 source->extNsigma = psMetadataLookupF32 (&status, row, "EXT_NSIGMA"); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_CMF_PS1_V1.c
r28860 r28909 192 192 psMetadataAdd (row, PS_LIST_TAIL, "PSF_MINOR", PS_DATA_F32, "PSF width (minor axis)", axes.minor); 193 193 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 194 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );194 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 195 195 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NDOF", PS_DATA_S32, "degrees of freedom", nDOF); 196 196 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NPIX", PS_DATA_S32, "number of pixels in fit", nPix); … … 328 328 } 329 329 330 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");330 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 331 331 source->crNsigma = psMetadataLookupF32 (&status, row, "CR_NSIGMA"); 332 332 source->extNsigma = psMetadataLookupF32 (&status, row, "EXT_NSIGMA"); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_CMF_PS1_V2.c
r28860 r28909 191 191 psMetadataAdd (row, PS_LIST_TAIL, "PSF_MINOR", PS_DATA_F32, "PSF width (minor axis)", axes.minor); 192 192 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 193 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );193 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 194 194 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NDOF", PS_DATA_S32, "degrees of freedom", nDOF); 195 195 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NPIX", PS_DATA_S32, "number of pixels in fit", nPix); … … 333 333 } 334 334 335 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");335 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 336 336 source->crNsigma = psMetadataLookupF32 (&status, row, "CR_NSIGMA"); 337 337 source->extNsigma = psMetadataLookupF32 (&status, row, "EXT_NSIGMA"); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_CMF_PS1_V3.c
r28860 r28909 193 193 psMetadataAdd (row, PS_LIST_TAIL, "PSF_MINOR", PS_DATA_F32, "PSF width (minor axis)", axes.minor); 194 194 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 195 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight); 195 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor (bad)", source->pixWeightNotBad); 196 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF_PERFECT", PS_DATA_F32, "PSF coverage/quality factor (poor)", source->pixWeightNotPoor); 196 197 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NDOF", PS_DATA_S32, "degrees of freedom", nDOF); 197 198 psMetadataAdd (row, PS_LIST_TAIL, "PSF_NPIX", PS_DATA_S32, "number of pixels in fit", nPix); … … 362 363 } 363 364 364 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF"); 365 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 366 source->pixWeightNotPoor = psMetadataLookupF32 (&status, row, "PSF_QF_PERFECT"); 365 367 source->crNsigma = psMetadataLookupF32 (&status, row, "CR_NSIGMA"); 366 368 source->extNsigma = psMetadataLookupF32 (&status, row, "EXT_NSIGMA"); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_PS1_CAL_0.c
r28860 r28909 188 188 psMetadataAdd (row, PS_LIST_TAIL, "PSF_WIDTH_Y", PS_DATA_F32, "PSF width in y coordinate", axes.minor); 189 189 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 190 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );190 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 191 191 192 192 // XXX not sure how to get this : need to load Nimages with weight? … … 300 300 source->peak->yf = PAR[PM_PAR_YPOS]; // more accurate position 301 301 302 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");302 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 303 303 304 304 // note that some older versions used PSF_PROBABILITY: this was not well defined. -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_PS1_DEV_0.c
r28860 r28909 125 125 psMetadataAdd (row, PS_LIST_TAIL, "PSF_WIDTH_Y", PS_DATA_F32, "PSF width in y coordinate", axes.minor); 126 126 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 128 128 // XXX not sure how to get this : need to load Nimages with weight 129 129 psMetadataAdd (row, PS_LIST_TAIL, "N_FRAMES", PS_DATA_U16, "Number of frames overlapping source center", nImageOverlap); … … 226 226 source->peak->yf = PAR[PM_PAR_YPOS]; // more accurate position 227 227 228 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");228 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 229 229 230 230 // XXX other values saved but not loaded? -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_PS1_DEV_1.c
r28860 r28909 149 149 psMetadataAdd (row, PS_LIST_TAIL, "PSF_WIDTH_Y", PS_DATA_F32, "PSF width in y coordinate", axes.minor); 150 150 psMetadataAdd (row, PS_LIST_TAIL, "PSF_THETA", PS_DATA_F32, "PSF orientation angle", axes.theta); 151 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeight );151 psMetadataAdd (row, PS_LIST_TAIL, "PSF_QF", PS_DATA_F32, "PSF coverage/quality factor", source->pixWeightNotBad); 152 152 153 153 // XXX not sure how to get this : need to load Nimages with weight? … … 269 269 source->peak->yf = PAR[PM_PAR_YPOS]; // more accurate position 270 270 271 source->pixWeight = psMetadataLookupF32 (&status, row, "PSF_QF");271 source->pixWeightNotBad = psMetadataLookupF32 (&status, row, "PSF_QF"); 272 272 273 273 // note that some older versions used PSF_PROBABILITY: this was not well defined. -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_RAW.c
r28643 r28909 128 128 source[0].peak->SN, 129 129 source[0].apRadius, 130 source[0].pixWeight ,130 source[0].pixWeightNotBad, 131 131 model[0].nDOF, 132 132 model[0].nIter); … … 190 190 source[0].peak->SN, 191 191 source[0].apRadius, 192 source[0].pixWeight ,192 source[0].pixWeightNotBad, 193 193 model[0].nDOF, 194 194 model[0].nIter); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceIO_SMPDATA.c
r28643 r28909 114 114 psMetadataAdd (row, PS_LIST_TAIL, "THETA", PS_DATA_F32, "", axes.theta); 115 115 psMetadataAdd (row, PS_LIST_TAIL, "DOPHOT", PS_DATA_U8, "", type); 116 psMetadataAdd (row, PS_LIST_TAIL, "WEIGHT", PS_DATA_U8, "", PS_MIN (255, PS_MAX(0, 255*source->pixWeight )));116 psMetadataAdd (row, PS_LIST_TAIL, "WEIGHT", PS_DATA_U8, "", PS_MIN (255, PS_MAX(0, 255*source->pixWeightNotBad))); 117 117 psMetadataAdd (row, PS_LIST_TAIL, "DUMMY", PS_DATA_U16, "", 0); 118 118 … … 197 197 source->apMag = psMetadataLookupF32 (&status, row, "MAG_AP") - ZERO_POINT; 198 198 199 source->pixWeight = psMetadataLookupU8 (&status, row, "WEIGHT")/255.0;199 source->pixWeightNotBad = psMetadataLookupU8 (&status, row, "WEIGHT")/255.0; 200 200 int dophot = psMetadataLookupU8 (&status, row, "DOPHOT"); 201 201 pmSourceSetDophotType (source, dophot); -
branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourcePhotometry.c
r28904 r28909 209 209 // if the aper mag is NAN, the flux < 0. this can happen for sources near the 210 210 // detection limits (esp near bright neighbors) 211 source->apMag = source->apMagRaw; 211 212 if (isfinite (source->apMag) && isPSF && psf) { 212 213 if (psf->growth && (mode & PM_SOURCE_PHOT_GROWTH)) { … … 311 312 bool pmSourcePixelWeight (float *pixWeightNotBad, float *pixWeightNotPoor, pmModel *model, psImage *mask, psImageMaskType maskVal, psImageMaskType markVal) 312 313 { 313 PS_ASSERT_PTR_NON_NULL(pixWeight, false); 314 PS_ASSERT_PTR_NON_NULL(pixWeightNotBad, false); 315 PS_ASSERT_PTR_NON_NULL(pixWeightNotPoor, false); 314 316 PS_ASSERT_PTR_NON_NULL(mask, false); 315 317 PS_ASSERT_PTR_NON_NULL(model, false);
Note:
See TracChangeset
for help on using the changeset viewer.
