1 / 19

A “Peak” at the Algorithm Behind “Peaklet Analysis” Software Bruce Kessler Western Kentucky University for KY MAA Annual

A “Peak” at the Algorithm Behind “Peaklet Analysis” Software Bruce Kessler Western Kentucky University for KY MAA Annual Meeting March 26, 2011. * This work was partially supported by KSEF grant KSEF-1653-RDE-011, and is currently supported by KCF grant COMMFUND-1280-RFP-011.

reya
Download Presentation

A “Peak” at the Algorithm Behind “Peaklet Analysis” Software Bruce Kessler Western Kentucky University for KY MAA Annual

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A “Peak” at the Algorithm Behind“Peaklet Analysis” SoftwareBruce KesslerWestern Kentucky UniversityforKY MAA Annual MeetingMarch 26, 2011 * This work was partially supported by KSEF grant KSEF-1653-RDE-011, and is currently supported by KCF grant COMMFUND-1280-RFP-011. * Patent pending on the algorithm, application #12/756,991.

  2. Wavelet Basics Wavelets are signal processing and analysis tools, analogous to trigonometric polynomials in the Fourier transform. • Wavelet bases are local, with coefficients representing different translates of the basis. Thus, coefficients reflect not just if a feature occurs, but give some idea of where a feature occurs. • Wavelet transforms ignore polynomial components of the signal up to the approximation order of the basis. • There are a wide selection of wavelet bases that can be used, depending on the needs or application of the user. • Multiwavelets can be simultaneously orthogonal (quick decompositions) and have symmetry properties. There is only one single wavelet construction that can have both.

  3. Wavelet Basics

  4. Wavelet Basics Then we may decompose a parent space into a smoother space and all of the included wavelet spaces. This type of decomposition has numerous applications in denoising, edge detection, data/image compression, etc., and now spectral analysis!

  5. Example with r = 1 (Haar) What is approximation order? If a scaling vector  has approximation orderk, then all of the spaces Vj generated by  will contain polynomials of degree k – 1. Approximation order 1 basis – just showed it to you.

  6. Approximation order 2 bases(Hardin, Geronimo, Massopust) orthogonalto integertranslates Another approximation order 2 basis (Daubechies) Also orthogonal to itsinteger translates, butnot when restricted toa shorter interval with integer endpoints.

  7. Approximation order 3 bases (Donovan, Hardin, Massopust)

  8. Approximation order 4 bases (DHM)

  9. Approximation order 4 bases (with one derivative, K.)

  10. When to use (multi)wavelets – Quantitative data. • Background noise. Wavelet decompositions ignore signal components up to the polynomial approximation order of the basis, so a wavelet analysis to look for patterns over the top of background noise. • Random noise. Gaussian noise appears as wavelet coefficients that are very close to 0, so they have little effect (and can be threshholded out if we desire). • Non-sparse data. Multiwavelets allow us to build boundary basis functions, allowing us to ignore data outside our analysis window.

  11. Background on the problem: There are several different ways of interrogating material that generate counts at different energy- levels or frequencies, with peak locations indicative of different elements.

  12. The problem is . . . the actual readings never look that clean. Background subtraction is typical. Some other numerical techniques may be used, but results are sometimes subjective.

  13. Wavelet-Based Algorithm  The user inputs peak information (center, standard deviation) for each element for which they are searching, along with other peaks in the same horizontal range. The user has the option of combining area under multiple peaks into one combined measure, and of setting relative areas.  The spectral data is input by the user. The peak centers may be adjusted slightly, within a user-set tolerance, to match the apparent calibration of the spectral data.  Wavelet decompositions are generated, to a level specified by the user, for each of the user-defined Gaussian peaks, using the adjusted centers.

  14. Wavelet-Based Algorithm  The wavelet decomposition of the spectral data is calculated to the user-defined level, and the wavelets are replaced with a variable linear combination of the wavelets from the elemental decompositions. The spectral data is then reconstructed.  The variables are set to amounts that give a least-squares fit to the original data.  (Optional) The user may input a library of elemental ratios for substances for which they are searching, and an angle will be calculated from the fit vector and each library vector, or combination of vectors.

  15. Example 1:

  16. Example 1:

  17. Example 1:  If relative heights are given, then the n-vector of amounts is compared to the n-vectors of substances in a user-input library.

  18. Example 1:

  19. Next steps: • The algorithm is partially coded into JAVA, but we may be starting over in a partnership with a software company, to provide a “try-before-you-buy” webpage. • We have our first potential customers waiting to test the algorithm on their data. • Licensing? Nobel? Mission Critical: I would like to thank you for your kind attention today. “And that’s all I have to say about that.” – Forrest Gump Thanks!works.bepress.com/bruce_kesslerbruce.kessler@wku.edu

More Related