xtyangpsp / StackMaster

A collection of methods for data stacking
MIT License
20 stars 5 forks source link

Add DOST-based tfpws algorithm #1

Closed jaredbryan881 closed 1 year ago

jaredbryan881 commented 1 year ago

Time-frequency phase weighted stacking based on the full stockwell transform can be too costly to compute for large datasets. This pull request adds an additional stacking method (tfpws_fast) based on the discrete orthonormal stockwell transform (DOST), which is essentially a short-time fourier transform with frequency-dependent bandwidth. The DOST implementation is described in [1] and is essentially a python translation of the Matlab implementation in [2].

[1] U. Battisti, L. Riba, "Window-dependent bases for efficient representations of the Stockwell transform", Applied and Computational Harmonic Analysis, 23 February 2015, http://dx.doi.org/10.1016/j.acha.2015.02.002. [2] Luigi Riba (2020). FFT-fast S-transforms: DOST, DCST, DOST2 and DCST2 (https://www.mathworks.com/matlabcentral/fileexchange/53910-fft-fast-s-transforms-dost-dcst-dost2-and-dcst2), MATLAB Central File Exchange. Retrieved November 28, 2020.