[Mne_analysis] Strange behavior of mne.preprocessing.current source density

[Alexandre Gramfort]

        External Email - Use Caution        

hi Mat,

please open an issue here https://github.com/mne-tools/mne-python/issues
with enough information so we can reproduce.

I would start here by testing that the matlab code gives different results
on the data from this example:

https://mne.tools/dev/auto_examples/preprocessing/plot_eeg_csd.html
<https://mne.tools/dev/auto_examples/preprocessing/plot_eeg_csd.html#sphx-glr-auto-examples-preprocessing-plot-eeg-csd-py>

Alex


On Tue, Jun 2, 2020 at 11:09 PM <mathieu.servant at univ-fcomte.fr> wrote:

>         External Email - Use Caution
>
> Hi,
> The current source density implementation in mne.python seems to compute
> incorrect CSD estimates. The attached file illustrates the issue. The first
> row shows topoplots of monopolar data (average reference). The second row
> shows topoplots of CSD estimates of these data, computed using the
> mne.preprocessing.courrent_cource_density function using default parameters
> (lambda2=1e-05, stiffness=4, n_legendre_terms=50). the third row shows
> topoplots of CSD estimates of these data, computed using Mike Cohen's
> code and similar parameters ( lambda2=1e-05, stiffness=4,
> n_legendre_terms=50).
> Do you have any idea of what's going on? Note that for rows 2 and 3, I
> used the 'biosemi64' electrode coordinates provided by mne.
>
> Best,
> Mat
>
>
>
> _______________________________________________
> Mne_analysis mailing list
> Mne_analysis at nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20200603/bbf1b977/attachment.html