Changeset 6862 for trunk/psphot/src/psphotFindPeaks.c
- Timestamp:
- Apr 14, 2006, 11:58:44 AM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotFindPeaks.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotFindPeaks.c
r6495 r6862 2 2 3 3 // 2006.02.02 : no leaks 4 psArray *psphotFindPeaks (pmReadout *readout, psMetadata * config) {4 psArray *psphotFindPeaks (pmReadout *readout, psMetadata *recipe) { 5 5 6 6 bool status = false; … … 12 12 psTimerStart ("psphot"); 13 13 14 SIGMA = psMetadataLookupF32 (&status, config, "PEAKS_SMOOTH_SIGMA");15 NSIGMA = psMetadataLookupF32 (&status, config, "PEAKS_SMOOTH_NSIGMA");14 SIGMA = psMetadataLookupF32 (&status, recipe, "PEAKS_SMOOTH_SIGMA"); 15 NSIGMA = psMetadataLookupF32 (&status, recipe, "PEAKS_SMOOTH_NSIGMA"); 16 16 17 17 psImage *smooth_im = psImageCopy (NULL, readout->image, PS_TYPE_F32); … … 37 37 38 38 // set peak threshold 39 NSIGMA = psMetadataLookupF32 (&status, config, "PEAKS_NSIGMA_LIMIT");39 NSIGMA = psMetadataLookupF32 (&status, recipe, "PEAKS_NSIGMA_LIMIT"); 40 40 // threshold = NSIGMA*sky->sampleStdev + sky->sampleMean; 41 41 threshold = PS_SQR(NSIGMA); … … 49 49 50 50 // optional dump of all peak data 51 char *output = psMetadataLookupStr (&status, config, "PEAKS_OUTPUT_FILE");51 char *output = psMetadataLookupStr (&status, recipe, "PEAKS_OUTPUT_FILE"); 52 52 if (status && (output != NULL) && (output[0])) { 53 53 pmPeaksWriteText (peaks, output); … … 55 55 psLogMsg ("psphot", 3, "%d peaks: %f sec\n", peaks->n, psTimerMark ("psphot")); 56 56 57 char *breakPt = psMetadataLookupStr (&status, config, "BREAK_POINT");57 char *breakPt = psMetadataLookupStr (&status, recipe, "BREAK_POINT"); 58 58 if (!strcasecmp (breakPt, "PEAKS")) exit (0); 59 59
Note:
See TracChangeset
for help on using the changeset viewer.
