中值滤波器超高效实现

Question

我正在寻找快速/高效中值滤波器的Ansi C实现。任何指针？

到目前为止，我找到了following的实现，这很好，但我对更快的实现感到好奇。我只需要1维。

Answer 1

我需要从非常嘈杂的CPU消耗数据中提取信号。这里是Jeff McClintock中值滤波器...

将平均值和中值初始化为零，然后对于每个样本的“英寸”，向输入样本的中值增加一个小增量。最终，它将在大约50％的投入样本更大，50％比中位数少的点上解决。 增量的大小应该与实际的中位数成正比。由于我们不知道实际的中位数，我用平均数作为粗略估计。步长计算为估计值的0.01倍。越小的步长越准确，但需要更长时间才能解决。

float median = 0.0f; 
float average = 0.0f; 

// for each sample 
{ 
    average += (abs(sample) - average) * 0.1f; // rough running average magnitude. 
    median += _copysign(average * 0.01, sample - median); 
}