I have used the following method to visualize brain parcels, but I want to add parcel names in addition to parcels.
What I want is for example for “IFSa” parcel I want to see the label name of that parcel on IFSa parcel.
Which method I shall follow for the desired modification?
Code:
import mne
Brain = mne.viz.get_brain_class()
subjects_dir = mne.datasets.sample.data_path() / “subjects”
mne.datasets.fetch_hcp_mmp_parcellation(subjects_dir=subjects_dir, verbose=True)
mne.datasets.fetch_aparc_sub_parcellation(subjects_dir=subjects_dir, verbose=True)
labels = mne.read_labels_from_annot(
“fsaverage”, “HCPMMP1”, “lh”, subjects_dir=subjects_dir
)
brain = Brain(
“fsaverage”,
“lh”,
“inflated”,
subjects_dir=subjects_dir,
cortex=“low_contrast”,
background=“white”,
size=(800, 600),
)
brain.add_annotation(“HCPMMP1”)
parcels_to_visualize = [
“L_SFL_ROI-lh”,
“L_33pr_ROI-lh”,
“L_p32pr_ROI-lh”,
“L_a24_ROI-lh”,
“L_8BM_ROI-lh”,
“L_8Av_ROI-lh”,
“L_8Ad_ROI-lh”,
“L_8BL_ROI-lh”,
“L_8C_ROI-lh”,
“L_44_ROI-lh”,
“L_IFJa_ROI-lh”,
“L_IFSp_ROI-lh”,
“L_46_ROI-lh”,
“L_p9-46v_ROI-lh”,
“L_47s_ROI-lh”,
“L_s6-8_ROI-lh”,
“L_25_ROI-lh”,
“L_a24pr_ROI-lh”,
“L_FOP4_ROI-lh”,
“L_AVI_ROI-lh”,
“L_FOP5_ROI-lh”,
“L_a32pr_ROI-lh”,
“L_p24_ROI-lh”,
]
for parcel_name in parcels_to_visualize:
parcel_label = [label for label in labels if label.name == parcel_name][0]
brain.add_label(parcel_label, borders=False)
brain.show()
