C#实现局部峰值查找,功能对应Matlab中的findpeaks.m

相关算法的原理参考Ronny，地址：图像分析：投影曲线的波峰查找，这里感谢下原作者。

参照C++的代码实现，我用C#翻译了下，其实原理也很简单的，下面放相关实现代码：

        private double[] oneDiff(double[] data)
        {
            double[] result = new double[data.Length - ];
            for (int i = ; i < result.Length; i++)
            {
                result[i] = data[i + ] - data[i];
            }
            return result;
        }
        private int[] trendSign(double[] data)
        {
            int[] sign = new int[data.Length];
            for (int i = ; i < sign.Length; i++)
            {
                if (data[i] > ) sign[i] =;
                else if (data[i] == ) sign[i] = ;
                else sign[i] = -;
            }

            for (int i = sign.Length - ; i >= ; i--)
            {
                if (sign[i] ==  && i ==sign.Length - )
                {
                    sign[i] = ;
                }
                else if (sign[i] == )
                {
                    if (sign[i+] >= )
                    {
                        sign[i] = ;
                    }
                    else
                    {
                        sign[i] = -;
                    }
                }
            }
            return sign;
        }
        private int[] getPeaksIndex(int[] diff)
        {
            List<int> data = new List<int>();
            for (int i = ; i != diff.Length - ; i++)
            {
                if (diff[i+] - diff[i] == -)
                {
                    data.Add(i + );
                }
            }
            int[] result = new int[data.Count];
            for (int i = ; i < result.Length; i++)
            {
                result[i] = data[i];
            }
            return result;//相当于原数组的下标
        }

调用方法：

double[] data = {, , , , , , , , , , };
int[] index = getPeaksIndex(trendSign(oneDiff(Constant.data)));

返回的int[]数组内容是对应原数组中峰值对应的索引(从0开始)