Read and plot coregistrationΒΆ

We’ll take a look at the coregistration here.

# Author: Denis A. Enegemann
# License: BSD 3 clause

import os.path as op

from hcp.viz import plot_coregistration

we assume our data is inside a designated folder under $HOME

storage_dir = op.expanduser('~/mne-hcp-data')

and we assume to have the downloaded data, the MNE/freesurfer style anatomy directory, and the MNE style MEG directory. these can be obtained from make_mne_anatomy(). See also Extract essential anatomy information.

hcp_params = dict(
    subject='105923',
    hcp_path=op.join(storage_dir, 'HCP'),
    subjects_dir=op.join(storage_dir, 'subjects'),
    recordings_path=op.join(storage_dir, 'hcp-meg'))

let’s plot two views

for azim in (0, 90):
    plot_coregistration(
        view_init=(('azim', azim), ('elev', 0)), **hcp_params)
  • ../_images/sphx_glr_plot_coregistration_001.png
  • ../_images/sphx_glr_plot_coregistration_002.png

Out:

Reading 4D PDF file /Users/dengemann/mne-hcp-data/HCP/105923/unprocessed/MEG/3-Restin/4D/c,rfDC...
Creating Neuromag info structure ...
... Setting channel info structure.
... no headshape file supplied, doing nothing.
Currently direct inclusion of 4D weight tables is not supported. For critical use cases please take into account the MNE command "mne_create_comp_data" to include weights as printed out by the 4D "print_table" routine.
Triangle file:  nvert = 5000 ntri = 9996
Reading 4D PDF file /Users/dengemann/mne-hcp-data/HCP/105923/unprocessed/MEG/3-Restin/4D/c,rfDC...
Creating Neuromag info structure ...
... Setting channel info structure.
... no headshape file supplied, doing nothing.
Currently direct inclusion of 4D weight tables is not supported. For critical use cases please take into account the MNE command "mne_create_comp_data" to include weights as printed out by the 4D "print_table" routine.
Triangle file:  nvert = 5000 ntri = 9996

Total running time of the script: ( 0 minutes 0.451 seconds)

Generated by Sphinx-Gallery