Changeset 17350
- Timestamp:
- Apr 6, 2008, 2:42:06 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branch_20080324/psphot/src/psphotFindPeaks.c
r16820 r17350 11 11 12 12 // smooth the image and weight map 13 psTimerStart ("p sphot");13 psTimerStart ("peaks"); 14 14 15 15 // XXX if we have been supplied a PSF, we can use that to set the FWHM_X,FWHM_Y values … … 32 32 psImage *smooth_im = psImageCopy (NULL, readout->image, PS_TYPE_F32); 33 33 psImageSmoothMaskF32 (smooth_im, readout->mask, maskVal, SIGMA_SMTH, NSIGMA_SMTH); 34 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth image: %f sec\n", psTimerMark ("p sphot"));34 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth image: %f sec\n", psTimerMark ("peaks")); 35 35 36 36 // smooth the weight, applying the mask as we go 37 37 psImage *smooth_wt = psImageCopy (NULL, readout->weight, PS_TYPE_F32); 38 38 psImageSmoothMaskF32 (smooth_wt, readout->mask, maskVal, SIGMA_SMTH/sqrt(2), NSIGMA_SMTH); 39 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth weight: %f sec\n", psTimerMark ("p sphot"));39 psLogMsg ("psphot", PS_LOG_MINUTIA, "smooth weight: %f sec\n", psTimerMark ("peaks")); 40 40 41 41 psImage *mask = readout->mask; … … 61 61 } 62 62 } 63 psLogMsg ("psphot", PS_LOG_INFO, "built smoothed signficance image: %f sec\n", psTimerMark ("p sphot"));63 psLogMsg ("psphot", PS_LOG_INFO, "built smoothed signficance image: %f sec\n", psTimerMark ("peaks")); 64 64 65 65 // optionally save example images under trace … … 70 70 } 71 71 72 psTimerStart ("p sphot");72 psTimerStart ("peaks"); 73 73 // set peak threshold 74 74 … … 140 140 pmPeaksWriteText (peaks, output); 141 141 } 142 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("p sphot"));142 psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("peaks")); 143 143 144 144 // If they asked us to return a set of pmFootprints, not a raw pmPeak list, … … 167 167 peaks = footprints; // well, you know what I mean 168 168 169 psLogMsg ("psphot", PS_LOG_INFO, "%ld footprints: %f sec\n", peaks->n, psTimerMark ("p sphot"));169 psLogMsg ("psphot", PS_LOG_INFO, "%ld footprints: %f sec\n", peaks->n, psTimerMark ("peaks")); 170 170 } 171 171
Note:
See TracChangeset
for help on using the changeset viewer.
