high_density_parcellation.py

import mne
import os.path as op
import glob
import numpy as np

subjects = ['CC110033', 'CC110037', 'CC110045']
subject = subjects[0]

# For Desktop
# https://www.dropbox.com/sh/zkq03rj76usd6b5/AAC4pvtWoNMUkxNwrCrzUVNRa?dl=0
data_path = '/cluster/transcend/sheraz/Dropbox/mne-camcan-data/'

subjects_dir = op.join(data_path,'recons')
subject_dir = op.join(subjects_dir,subject)
bem_dir = op.join(subject_dir,'bem')
trans_file = op.join(data_path, 'trans',subject + '-trans.fif')
labels_fname  = glob.glob(op.join(data_path, 'labels', '*.label'))
labels = [mne.read_label(label, subject='fsaverageSK', color='r')
          for label in labels_fname]
for index, label in enumerate(labels):
    label.values.fill(1.0)
    labels[index] = label

labels = [label.morph('fsaverageSK', subject, subjects_dir=subjects_dir) for label in labels]