import scipy.io import mne from mne.bem import make_watershed_bem # Load .mat files inner_skull = scipy.io.loadmat('E:\MATLABproject\data\MRI\Visit1_040318\\tess_mri_COR_MPRAGE_RECON-mocoMEMPRAGE_FOV_220-298665.inner_skull.mat') outer_skull = scipy.io.loadmat('E:\MATLABproject\data\MRI\Visit1_040318\\tess_mri_COR_MPRAGE_RECON-mocoMEMPRAGE_FOV_220-298665.outer_skull.mat') scalp = scipy.io.loadmat('E:\MATLABproject\data\MRI\Visit1_040318\\tess_mri_COR_MPRAGE_RECON-mocoMEMPRAGE_FOV_220-298665.scalp.mat') print(inner_skull.keys()) # Assuming these .mat files contain triangulated surfaces, we will extract vertices and triangles # This might need adjustment based on the actual structure of your .mat files inner_skull_vertices = inner_skull['Vertices'] inner_skull_triangles = inner_skull['Faces'] outer_skull_vertices = outer_skull['Vertices'] outer_skull_triangles = outer_skull['Faces'] scalp_vertices = scalp['Vertices'] scalp_triangles = scalp['Faces'] # Prepare surfaces for MNE surfs = [ mne.bem.BEMSurface(inner_skull_vertices, inner_skull_triangles, sigma=0.01, id=4), # brain mne.bem.BEMSurface(outer_skull_vertices, outer_skull_triangles, sigma=0.016, id=3), # skull mne.bem.BEMSurface(scalp_vertices, scalp_triangles, sigma=0.33, id=5), # skin ] # Create BEM model model = mne.bem.BEM(surfs, conductivity=[0.3, 0.006, 0.3], is_sphere=False) model.plot(show=False) # Create BEM solution solution = mne.make_bem_solution(model) 运行代码时报错; Traceback (most recent call last): File "E:\pythonProject\MEG\头模型.py", line 24, in <module> mne.bem.BEMSurface(inner_skull_vertices, inner_skull_triangles, sigma=0.01, id=4), # brain AttributeError: module 'mne.bem' has no attribute 'BEMSurface'

根据你的代码和报错信息来看，mne 模块确实没有 `BEMSurface` 这个属性。根据官方文档，在 mne 的最新版本中，使用 `BEMSurface` 已经被废弃，而改用了 `make_watershed_bem` 函数来创建 BEM 模型。

你可以尝试修改你的代码如下：

import scipy.io
import mne
from mne.bem import make_watershed_bem

# Load .mat files
inner_skull = scipy.io.loadmat('E:\MATLABproject\data\MRI\Visit1_040318\\tess_mri_COR_MPRAGE_RECON-mocoMEMPRAGE_FOV_220-298665.inner_skull.mat')
outer_skull = scipy.io.loadmat('E:\MATLABproject\data\MRI\Visit1_040318\\tess_mri_COR_MPRAGE_RECON-mocoMEMPRAGE_FOV_220-298665.outer_skull.mat')
scalp = scipy.io.loadmat('E:\MATLABproject\data\MRI\Visit1_040318\\tess_mri_COR_MPRAGE_RECON-mocoMEMPRAGE_FOV_220-298665.scalp.mat')
print(inner_skull.keys())

# Assuming these .mat files contain triangulated surfaces, we will extract vertices and triangles
# This might need adjustment based on the actual structure of your .mat files

inner_skull_vertices = inner_skull['Vertices']
inner_skull_triangles = inner_skull['Faces']

outer_skull_vertices = outer_skull['Vertices']
outer_skull_triangles = outer_skull['Faces']

scalp_vertices = scalp['Vertices']
scalp_triangles = scalp['Faces']

# Prepare surfaces for MNE
surfs = [
    mne.make_bem_model(inner_skull_vertices, inner_skull_triangles, conductivity=0.01),  # brain
    mne.make_bem_model(outer_skull_vertices, outer_skull_triangles, conductivity=0.016),  # skull
    mne.make_bem_model(scalp_vertices, scalp_triangles, conductivity=0.33),  # skin
]

# Create BEM solution
model = make_watershed_bem(surfs)
solution = mne.make_bem_solution(model)

请注意，我还对 `make_bem_model` 的调用进行了修改，以适应新的 API。希望这能解决你的问题！如果还有其他疑问，请随时提问。