[Numpy-discussion] Functions for finding the relative extrema of numeric data

[Jacob Silterra]

Hello all,

I just opened a pull request to add scipy.signal._peak_finding.find_peaks,
which finds the relative maxima in a 1d ndarray. The algorithm behind this
function has been discussed on the scipy list, but here's a quick recap: The
idea is basically convolve the function with a wavelet of varying widths
(default is ricker, or "mexican hat" wavelet) and identify as peaks those
locations where there are relative maxima (meaning they are larger than
their nearest neighbor on each side) at each of several widths, and have a
sufficiently high SNR. More details are provided in the docstrings of the
functions. Pull request at https://github.com/scipy/scipy/pull/85.

I believe the code quality/docstrings/test coverage are pretty good, however
it would be of great benefit for someone with experience with similar
algorithms to review the code.

-Jacob

PS  I've sent this to both the scipy and numpy lists, apologies to those who
get it twice.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20110930/2914ec4e/attachment.html>