Filtr Median Bardzo wydajna implementacja

Question

Poszukuję implementacji Ansi C szybkiego/wydajnego filtru median. Jakieś wskazówki?

Do tej pory znalazłem implementację following, która jest dobra, ale jestem ciekawa na szybsze. Potrzebuję tylko jednego wymiaru.

Answer 1

Potrzebowałem wyodrębnić sygnał z bardzo hałaśliwych danych o zużyciu procesora. Oto filtr mediana Jeffa McClintocka ...

Inicjalizuj średnią i medianę na zero, następnie dla każdej próbki "cal" medianę w kierunku próbki wejściowej o mały przyrost. Ostatecznie ustabilizuje się w punkcie, w którym około 50% próbek wejściowych jest większych, a 50% jest mniejszych niż mediana. Wielkość inkrementu powinna być proporcjonalna do rzeczywistej mediany. Ponieważ nie znamy rzeczywistej mediany, używam średniej jako przybliżonej. Wielkość kroku jest obliczana jako 0,01 wartości szacunkowej. Mniejsze rozmiary kroku są bardziej dokładne, ale wymagają więcej czasu.

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