Skip to content

Commit

Permalink
Add more details to wavelet transform (#74)
Browse files Browse the repository at this point in the history
* Add more details to wavelet transform

* Fixed flake errors

* Added blank line for docs build

* Try 2

* If this doesn't work I quit
  • Loading branch information
jthermiz committed Jun 24, 2022
1 parent 8369e58 commit ffcab29
Showing 1 changed file with 11 additions and 1 deletion.
12 changes: 11 additions & 1 deletion src/process_nwb/wavelet_transform.py
Expand Up @@ -240,7 +240,17 @@ def wavelet_transform(X, rate, filters='rat', hg_only=True, X_fft_h=None, npad='
"""Apply a wavelet transform using a prespecified set of filters.
Calculates the center frequencies and bandwidths for the wavelets and applies them along with
a heavyside function to the fft of the signal before performing an inverse fft.
a heavyside function to the fft of the signal before performing an inverse fft. Here are additional details:
**1.** Computes the FFT of the signal and applies 2u(f), where u(f) = heaviside function necessary for
calculating analytic signal
**2.** Filters in the frequency domain by multiplying with a gaussian kernel
(or equivalently a complex morlet wavelet in the time domain). The gaussian kernel location = the center frequency
and the standard deviation = the center frequency / Q. For the 'rat' and 'human' filters, Q is constant with a
default value of 8
**3.** Computes the IFFT returning the analytic bandpassed (via a complex morlet wavelet) signal for each filter
Parameters
----------
Expand Down

0 comments on commit ffcab29

Please sign in to comment.