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开始)