In February 2008, Natalia Stefanova asked Ferret users if anyone had a median filter script to remove spikes from a time series. I, too, need one, so I wrote one. It works, but it is slow if your series is long. The Ferret group may add a built-in @MED:N smoothing function to Ferret someday, but in the meantime, here's mine. I attach 3 files: (1) median_test.jnl: A test file that gives 2 examples of how to use the script. (2) median_l.jnl: The script for computing the median in l, the time axis index. (3) median_l_aux.jnl: A sub-procedure called by median_l.jnl. ***** |
Attachment:
median_test.jnl
Description: Binary data
***** |
Attachment:
median_l.jnl
Description: Binary data
***** |
Attachment:
median_l_aux.jnl
Description: Binary data
***** ------------------------------------------------------------------------------------------ Edward D. (Ned) Cokelet, Ph.D. Oceanographer NOAA/PMEL off: (206) 526-6820 7600 Sand Point Way NE fax: (206) 526-6485 Seattle, WA 98115-6349 The contents of this message are mine personally and do not necessarily reflect any position of the Government or the National Oceanic and Atmospheric Administration. |