Median Filter Implementación super eficiente

Question

Estoy buscando una implementación de Ansi C de un filtro mediano rápido/eficiente. ¿Alguna sugerencia?

Hasta ahora, he encontrado la implementación following, que es buena pero tengo curiosidad por las más rápidas. Solo necesito 1 dimensión.

Answer 1

Necesitaba extraer una señal de datos de consumo de CPU muy ruidosos. Aquí está el filtro mediano de Jeff McClintock ...

Inicialice el promedio y la mediana a cero, luego, para cada muestra 'en pulgadas', la mediana hacia la muestra de entrada en un pequeño incremento. Eventualmente se establecerá en un punto donde aproximadamente el 50% de las muestras de entrada son mayores y el 50% son menores que la mediana. El tamaño del incremento debe ser proporcional a la mediana real. Como no conocemos la mediana real, utilizo el promedio como una estimación aproximada. El tamaño del paso se calcula como 0.01 veces el estimado. Los tamaños de pasos más pequeños son más precisos, pero tardan más en resolverse.

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); 
}