Peak finding and measurement (version 2)

Function to locate and measure the positive peaks in noisy data sets.

Description

A fast customizable function for locating and measuring the peaks in noisy time-series signals. Adjustable parameters allow discrimination of "real" signal peaks from noise and background. Determines the position, height, and width of each peak by least-squares curve-fitting. It can find and measure 1000 peaks in a 1,000,000 point signal in 8 seconds. Includes two interactive versions, one with mouse-controlled sliders and one with keyboard control, for adjusting the peak finding criteria in real-time. Self-contained demos show how it works. See http://www.wam.umd.edu/~toh/spectrum/PeakFindingandMeasurement.htm for details.

• Acknowledgements: This submission has inspired the following:; iPeak: Interactive Peak Detector (Version 1.1)
• MATLAB release: MATLAB 6.5 (R13)

Other Files

InteractivePeakfinder2/AmpTSlider.m,
InteractivePeakfinder2/BGSlider.m,
InteractivePeakfinder2/condense.m,
InteractivePeakfinder2/DemoFindPeak.m,
InteractivePeakfinder2/DemoFindPeakSliders.m,
InteractivePeakfinder2/DemoFindPeakSliders2.m,
InteractivePeakfinder2/DemoFindPeakSlidersG.m,
InteractivePeakfinder2/deriv.m,
InteractivePeakfinder2/fastsmooth.m,
InteractivePeakfinder2/findpeaks.m,
InteractivePeakfinder2/FindPeakSliders.m,
InteractivePeakfinder2/FindPeakSlidersG.m,
InteractivePeakfinder2/findpeakslidersRedraw.m,
InteractivePeakfinder2/FitSlider.m,
InteractivePeakfinder2/gaussian.m,
InteractivePeakfinder2/lorentzian.m,
InteractivePeakfinder2/Peakfinder info.txt,
InteractivePeakfinder2/PeakSlider.m,
InteractivePeakfinder2/RedrawPeak.m,
InteractivePeakfinder2/rtslid.m,
InteractivePeakfinder2/secderiv.m,
InteractivePeakfinder2/SlopeTSlider.m,
InteractivePeakfinder2/SmoothSlider.m,
InteractivePeakfinder2/TestPeakfind.m,
InteractivePeakfinder2/Thumbs.db,
InteractivePeakfinder2/val2ind.m,
license.txt

[Frank Anold]

It is pretty good and easy to use. Yet some commercial software,such as origin, can do faster.

[Tom O'Haver]

Version 1.2, dated August 24, 2006, is even faster than the previous version It can find and measure 1000 peaks in a 1,000,000 point signal in 13 seconds.

[K. Tampa]

Get EXTREME.M in this site and you'll find the ALL peaks (minima and maxima!) of
rand(10^6,1) in less than a second, 10^7 in about 10 sec, and from matrix too, so...

[mcs swamy]

 

[Tom O'Haver]

K. Tampa is correct; extreme.m is the one to use if you want to find ALL the minima and maxima. But "Peak finding and measurement" is intended to distinguish real peaks from noise in experimentally measured signals and to measure their height, position, and width. Two completely different purposes.

[BHARATH KUMAR]

 

[Dmitry Kazachk]

Would be greate if it will do deconvolution and gave the areas for deconvolved peaks.
In general - great tool!

[Tom Driscoll]

Perfect for scientific data processing. Good performance without being overly complicated. You just saved me a week of my time.

[Matthew Butcher]

You saved my MEng project with this one. I credited you thogh, don't worry.
Many thanks Tom,
Matt.

[Li XinBo]

Your algorithm is excellent.But I find a problem in it.
Your algorithm is limited only when the data are positive number.If all the data are negative and there is peak points.Your algorithm can't detect the peak points.
But I can translate the coordinate axis.
It is to make all the y data into positive number. Then using your algorithm we can find the peak points.
I think you can modify your algorithm and make it adapt to the situation that the data is all negative number.
If you can modify your algorithm,please send the message to my email.([e-mail])
Thank your.
A student from BIT of China

[Karan Jain]

This is a good algorithm but can it be changed to give peeks over a certain datum instead of the trough values.

[Hugo Banziger]

Please gat back to me with your private phone number.

Regards,
Dr. Hugo.

[Julius Kusuma]

Very useful interactive software! Thanks for sharing this.

[Ralph Ernstorfer]

fast and reliable tool ... thanks for sharing!

[Christof Devriendt]

Just what I needed! Thank you very much!

[Tom O'Haver]

Latest version just uploaded corrects a bug in the FindPeakSlidersG script that caused a malfunction when individual peaks are viewed if the x-axis data series does not start with zero.

[Michael Coen]

Very nifty idea but it crashes quickly in Matlab 2008b, 2009a, and 2009b. There are clearly some GUI incompatibilities.

[Michael Coen]

Ah, commenting out line 337 of rtslid.m:

"figure(params.fig);"

eliminates crashes for later versions of Matlab...

[Michael Coen]

And comment out line 372 of rtslid.m as well... :)