Changeset 10801 for trunk/psphot/src/psphotFindPeaks.c
- Timestamp:
- Dec 17, 2006, 12:49:47 AM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotFindPeaks.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotFindPeaks.c
r10295 r10801 20 20 psImage *smooth_im = psImageCopy (NULL, readout->image, PS_TYPE_F32); 21 21 psImageSmooth (smooth_im, SIGMA, NSIGMA); 22 psLogMsg ("psphot", 4, "smooth image: %f sec\n", psTimerMark ("psphot"));22 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth image: %f sec\n", psTimerMark ("psphot")); 23 23 24 24 psImage *smooth_wt = psImageCopy (NULL, readout->weight, PS_TYPE_F32); … … 26 26 smooth_wt = (psImage*)psBinaryOp(smooth_wt, smooth_wt, "*", 27 27 psScalarAlloc(1/(4*M_PI*PS_SQR(SIGMA)), PS_TYPE_F32)); 28 psLogMsg ("psphot", 4, "smooth weight: %f sec\n", psTimerMark ("psphot"));28 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth weight: %f sec\n", psTimerMark ("psphot")); 29 29 30 30 psImage *mask = readout->mask; … … 36 36 } 37 37 38 psTimerStart ("psphot");39 38 for (int j = 0; j < smooth_im->numRows; j++) { 40 39 for (int i = 0; i < smooth_im->numCols; i++) { … … 47 46 } 48 47 } 49 psLogMsg ("psphot", 4, "built S/N map: %f sec\n", psTimerMark ("psphot")); 50 51 psTimerStart ("psphot"); 48 psLogMsg ("psphot", PS_LOG_MINUTIA, "built S/N map: %f sec\n", psTimerMark ("psphot")); 52 49 53 50 // set peak threshold … … 55 52 // threshold = NSIGMA*sky->sampleStdev + sky->sampleMean; 56 53 threshold = PS_SQR(NSIGMA); 57 psLogMsg ("psphot", 4, "threshold: %f DN\n", threshold);54 psLogMsg ("psphot", PS_LOG_DETAIL, "threshold: %f DN\n", threshold); 58 55 59 56 // find the peaks in the smoothed image … … 68 65 pmPeaksWriteText (peaks, output); 69 66 } 70 psLogMsg ("psphot", 3, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("psphot"));67 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("psphot")); 71 68 72 69 return (peaks);
Note:
See TracChangeset
for help on using the changeset viewer.
