No implementation for csd_welch?

đź’¬ Support & Discussions
1

MNE version: 1.7.1
OS: Windows 10

In the J. Gross et al. 2001 DICS paper, the welch method is used to compute CSD

The complex cross spectral density C for signals x (t ) and y (t ) is computed by using Welch’s method of spectral density estimation. A Hanning window is applied to the segments of signals x (t ) and y (t ).

But I can’t find implementation for computing welch CSD on epochs object directly like csd_morlet or csd_fourier (well you can sort of make it into small windows and apply a hann window to get welch spectral estimate, but that’s extra work). Is there any particular reason not having this function?

I found scipy.signal.csd computes CSD using welch by default. Would you recommend to compute CSD for DICS using this function? For example:
scipy.signal.csd(epoch.get_data(), epoch.get_data())

I can contribute to write a csd_welch function if you think this is a good idea.

Thanks,
ST

2

the authors of that paper later switched to using Morlet wavelets in their DICS analysis (see van Vliet et al. 2018 for a description of their current pipeline). This is why I didn’t implement Welch.

1 Like
3

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.