Plotting sEEG contacts

šŸ’¬ Support & Discussions
,
1

Hi, I may have a similar question. I did freesurferā€™s recon-all. And I want to plot the SEEG contacts(in MNI space) in lh.pial and lr.pial of this subject.
I donā€™t know if I understand right. Iā€™m following the tutorial here How MNE uses FreeSurferā€™s outputs ā€” MNE 0.24.0 documentation
This is what I did.

import mne
import nibabel as nib
mri_mni_trans = mne.read_talxfm(subject, subjects_dir)
ver, face = nib.freesurfer.read_geometry(lh_pial_path)
ver_mni = apply_trans(mri_mni_trans, ver)

with ver, so I could plot this pial in MNI together with SEEG contacts.
However, the result seems not right. Where did I do wrong?

image
image821Ɨ647 68.7 KB

2

@BarryLiu97 Iā€™ve moved your message to a new topic. Please always create a new topic for a new question to make conversation easier. Thank you!

cc @alexrockhill

3

That case is handled exactly in Working with sEEG data ā€” MNE 1.0.dev0 documentation, is there something there that is unclear?

4

Also your code doesnā€™t cover how you plotted the contacts and clearly there is a coordinate frame mismatch so that information is needed too.

5

My contacts are in MNI space, just like the form in the tutorial (x, y, z) in mm. when I use ā€˜fsaverageā€™, the figure was right. So I was thinking whether I transform the pial in a wrong way.

6

If you put your contact locations in a montage and assigned it to a raw object with fiducials, that would transform the contacts into ā€œheadā€ space because thatā€™s how MNE stores coordinates internally. So itā€™s a bit more complicated and without the full code to reproduce your figure, I really canā€™t guess what is going wrong.

1 Like
7

Hi, I would post my code and data here after organizing it.