(resubmitted) Problem with plotting VolSourceEstimate

✉️ Mailing List Archive (read-only)
1

External Email - Use Caution

  Hello, Alexandre,

  Here is
<https://drive.google.com/file/d/1e6-p5S3XXZRrIFvpLl6WNrLJsaI4RHPz/view?usp=sharing&gt;
the code you asked for, I hope this helps.

  Sincerely,
GM.

        External Email - Use Caution

Dear Gregory,

I mean a script I can run on my machine so I can replicate this issue.

for this I need a script that operates on a public dataset so I suggest
you use
the sample dataset we use in all our examples online.

thanks

Alex

>
> External Email - Use Caution
>
> Did you mean code with all preprocessing, events creation and other
steps? If that's the case, I have uploaded it here. LINK
>
> GM.
> _______________________________________________
> Mne_analysis mailing list
> Mne_analysis at nmr.mgh.harvard.edu
> Mne_analysis Info Page

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

_______________________________________________
Mne_analysis mailing list
Mne_analysis at nmr.mgh.harvard.edu
Mne_analysis Info Page

End of Mne_analysis Digest, Vol 143, Issue 23
*********************************************

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20191224/26e734fe/attachment.html

2

External Email - Use Caution

hi,

I can replicate the pb. There is a bug in setup_volume_source_space

this:
src = mne.setup_volume_source_space(bem=bem)

returns a discrete source space which is not expected. It should be a volume
source space.

can you open an issue on github?

thanks

Alex

3

External Email - Use Caution

hi,

I looked into the issue. The problem is following:

src = mne.setup_volume_source_space(bem=bem)

returns a discrete source space as it does not know what subject to
which the bem corresponds to.
Basically it does not know where to get the MRI from.

to fix your problem you should use:

src = mne.setup_volume_source_space(bem=bem, subject='sample',
subjects_dir=subjects_dir)

I agree that this is confusing as you needed to pass subject and
subjects_dir in mne.make_bem_model
to technically you provided this information.

Eric what do you think?

Alex

4

External Email - Use Caution

This sounds right.

We could try to infer the subject (it might be in the BEM somewhere), but
then we still wouldn't know subjects_dir. Maybe we should add a warning
saying that no MRI was provided and this the source space will be discrete,
since it seems rare that you'd want to use a BEM, spacing as a float (i.e.,
grid), but have no MRI associated with it.

Eric