efficient running median

[Peter Otten]

Janto Dreijer wrote:

> I'm looking for code that will calculate the running median of a
> sequence, efficiently. (I'm trying to subtract the running median from
> a signal to correct for gradual drift).
> 
> My naive attempt (taking the median of a sliding window) is
> unfortunately too slow as my sliding windows are quite large (~1k) and
> so are my sequences (~50k). On my PC it takes about 18 seconds per
> sequence. 17 of those seconds is spent in sorting the sliding windows.
> 
> I've googled around and it looks like there are some recent journal
> articles on it, but no code. Any suggestions?

If you aren't using numpy, try that. Showing your code might also be a good 
idea...

Peter