nelpy/ghost

Spectral analysis tools for neuroscience data

ghost's Issues

Adaptive sampling of wavelets

At the lower frequencies, we need more samples to get an accurate representation of the wavelet, but at higher frequencies, we need fewer. I need to figure out a better way of sampling the wavelets rather than the default of 16000 samples or so.

scipy: if _is_safe_size(n):

Hi,
This is my code with Accz and Time are array of float:

cwt = ContinuousWaveletTransform(wavelet= morse.Morse())
CPEF=cwt.transform(data =Accz, fs =1000, timestamps =Time, freq_limits =[0,300])

And the answer:

File ".......\lib\site-packages\scipy\fftpack\basic.py", line 87, in fft
if _is_safe_size(n):

TypeError: c2c() missing 1 required keyword-only argument: 'plan'

scipy fftpack missing

Once I fix the logsumexp issue, I get a new error:

from ghost.wave import ContinuousWaveletTransform

Now I get:

---------------------------------------------------------------------------

ModuleNotFoundError                       Traceback (most recent call last)

<ipython-input-5-e9048c87e447> in <module>()
----> 1 from ghost.wave import ContinuousWaveletTransform

5 frames

/content/ghost/ghost/sigtools/fourier.py in <module>()
      3 import logging
      4 import numpy as np
----> 5 import scipy.fftpack._fftpack as sff
      6 from . import convolution
      7 

ModuleNotFoundError: No module named 'scipy.fftpack._fftpack'

Too many points in wavelet spectrogram plot

If there are too many points in the wavelet spectrogram plot, some parts of the image will not render properly. I need to figure out a way to downsample correctly without messing up the temporal resolution.

Can't import ContinuousWaveletTransform with newest scipy

On google collaboratory, if I install ghost from git, and then run

from ghost.wave import ContinuousWaveletTransform

I get:

---------------------------------------------------------------------------

ImportError                               Traceback (most recent call last)

<ipython-input-8-e9048c87e447> in <module>()
----> 1 from ghost.wave import ContinuousWaveletTransform

2 frames

/content/ghost/ghost/wave/morlet.py in <module>()
      3 import numpy as np
      4 from .wavelet import Wavelet
----> 5 from scipy.misc import logsumexp
      6 import copy
      7 

ImportError: cannot import name 'logsumexp'

Handle timestamp info for wavelet spectrogram plot

Currently the timestamps are generated from the data but there should be a way to incorporate optional passed-in timestamps for plotting a wavelet spectrogram.