more general filtering (e.g. band pass filters)

serazing / xscale

Xscale a library of multi-dimensional signal processing tools using parallel computing

Apache License 2.0

25 stars 12 forks source link

more general filtering (e.g. band pass filters) #14

Open apatlpo opened 6 years ago

apatlpo commented 6 years ago

I've got the feeling this is not implemented yet, while scipy.signal.firwin allows in principle. Does anybody confirm?

serazing commented 6 years ago

Yes, you can only do a low-pass filtering for the moment. I usually estimate the high-pass component by subtracting the low-pass component to the original field. Implementing high-pass filtering would just require to set 1 - low-pass kernel. The band-pass would require to make a composition of a high-pass and a low pass-kernel. Feel free to modify the code in xscale.window.linearfilters.Window.set and submit a pull request if you want this functionality added to xscale.

stephanieleroux commented 5 years ago

hi Guillaume, how about a lanczos "window" like you had added in oocgcm?