MNE绘制自定义通道位置及名字的脑地形图

import matplotlib.pyplot as plt
import numpy as np
import mne
 
ch_names=['AF3', 'AF4', 'F5', 'F3', 'F1', 'Fz', 'F2', 'F4', 'F6', 'FC5', 'FC3', 'FC1', 'FCz', 'FC2', 'FC4', 'FC6', 'CFC7', 'CFC5', 'CFC3', 'CFC1', 'CFC2', 'CFC4', 'CFC6', 'CFC8', 'T7', 'C5', 'C3', 'C1', 'Cz', 'C2', 'C4', 'C6', 'T8', 'CCP7', 'CCP5', 'CCP3', 'CCP1', 'CCP2', 'CCP4', 'CCP6', 'CCP8', 'CP5', 'CP3', 'CP1', 'CPz', 'CP2', 'CP4', 'CP6', 'P5', 'P3', 'P1', 'Pz', 'P2', 'P4', 'P6', 'PO1', 'PO2', 'O1', 'O2']
 
xpos=[-0.20109, 0.20109, -0.485475, -0.328947, -0.165352, 0, 0.165352, 0.328947, 0.485475, -0.605915, -0.399196, -0.197659, 0, 0.197659, 0.399196, 0.605915, -0.748347, -0.52473, -0.309639, -0.102263, 0.102263, 0.309639, 0.52473, 0.748347, -0.877193, -0.645691, -0.421549, -0.207738, 0, 0.207738, 0.421549, 0.645691, 0.877193, -0.748347, -0.52473, -0.309639, -0.102263, 0.102263, 0.309639, 0.52473, 0.748347, -0.605915, -0.399196, -0.197659, 0, 0.197659, 0.399196, 0.605915, -0.485475, -0.328947, -0.165352, 0, 0.165352, 0.328947, 0.485475, -0.103959, 0.103959, -0.171132, 0.171132]
 
ypos=[0.686565177105246,0.686565177105246,0.525474237033963,0.465201829727992,0.432086407208813,0.421549001338212,0.32086407208813,0.465201829727992,0.525474237033963,0.271657042421725,0.233843478507591,0.213944937077633,0.207737565050372,0.213944937077633,0.233843478507591,0.271657042421725,0.151771690883298,0.125531031531108,0.110860959086494,0.4266479482922,0.104266479482922,0.110860959086494,0.125531031531108,0.151771690883298,0,0,0,0,0,0,0,0,0,-0.151771690883298,-0.25531031531108,-0.110860959086494,-0.104266479482922,-0.104266479482922,-0.110860959086494,-0.125531031531108,-0.51771690883298,-0.271657042421725,-0.233843478507591,-0.213944937077633,-0.207737565050372,-0.213944937077633,-0.33843478507591,-0.271657042421725,-0.525474237033963,-0.465201829727992,-0.432086407208813,-0.421549001338212,-0.32086407208813,-0.465201829727992,-0.525474237033963,-0.655838116128691,-0.655838116128691,-0.860337965265992,-0.860337965265992]
 
for i in range(59):
    xpos[i]=xpos[i]*0.1
    ypos[i]=ypos[i]*0.1
 
data = np.random.rand(59,1)  
 
position = {}
for i in range(59):
    ch_name = ch_names[i]
    pos = [xpos[i],ypos[i],0]
    position[ch_name] = np.array(pos)
 
montage = mne.channels.make_dig_montage(ch_pos=position)
 
#创建info对象
info = mne.create_info(ch_names=ch_names, sfreq=200.,ch_types='eeg')
 
#创建evokeds对象                       
evoked = mne.EvokedArray(data, info)
#evokeds设置通道
evoked.set_montage(montage)
mne.viz.plot_topomap(evoked.data[:,0], evoked.info,show=False)
 
plt.show()