Accuracy and run-time comparison for different potential approaches and iterative solvers in finite element method based EEG source analysis.
Accuracy and run-time comparison for different potential approaches and iterative solvers in finite element method based EEG source analysis.
复制标题
DOI：
10.1016/j.apnum.2009.02.006
复制
发表时间：
2009-08
影响因子：
2.8
通讯作者：
MacLeod, R. S.
中科院分区：
文献类型：
--
作者：
Lew, S.;Wolters, C. H.;Dierkes, T.;Roeer, C.;MacLeod, R. S.
关键词：
ElectroencephalographySource reconstructionFinite element methodDipole singularityFull subtraction potential approachVenant potential approachPartial integration potential approachPreconditioned conjugate gradient methodAlgebraic multi-grid preconditionerIncomplete Cholesky preconditionerJacobi preconditionerConstrained Delaunay tetrahedralizationAnisotropic four-layer sphere model
准确性和运行时间在医学诊断和研究以及神经科学领域都发挥着重要作用。在脑电(EEG)源重建中，根据测量到的头部表面的电位非侵入性地重建人脑中的电流分布(EEG逆问题)。数值模拟技术被用来模拟人体大脑皮层偶极电流源的头部表面电位，即所谓的脑电正问题。本文比较了代数多重网格(AMG)、不完全Cholesky(IC)和Jacobi预条件子在迭代求解基于有限元(FE)方法的脑电正问题中的效率。研究了处理偶极子奇异性的三种求解器--全减法和两种直接位势法、维纳法和部分积分法的相互作用。检验是在具有各向异性头盖层的四室球体模型中进行的，其中准解析解允许精确地量化计算速度与数值误差之间的关系。特调约束Delaunay四面化(CDT)有限元网格在全减法和直接势方法中都有很高的精度。如果满足齐性条件，则采用全减法可获得最佳精度。结果表明，在减小网格尺寸的情况下，AMG-CG的计算速度比使用标准预条件的CG提高了一个数量级，并且增加了增益因子。我们的研究结果将为高分辨率高精度有限元体导体模型在震源分析中的应用提供更广阔的空间。
Accuracy and run-time play an important role in medical diagnostics and research as well as in the field of neuroscience. In Electroencephalography (EEG) source reconstruction, a current distribution in the human brain is reconstructed noninvasively from measured potentials at the head surface (the EEG inverse problem). Numerical modeling techniques are used to simulate head surface potentials for dipolar current sources in the human cortex, the so-called EEG forward problem. In this paper, the efficiency of algebraic multigrid (AMG), incomplete Cholesky (IC) and Jacobi preconditioners for the conjugate gradient (CG) method are compared for iteratively solving the finite element (FE) method based EEG forward problem. The interplay of the three solvers with a full subtraction approach and two direct potential approaches, the Venant and the partial integration method for the treatment of the dipole singularity is examined. The examination is performed in a four-compartment sphere model with anisotropic skull layer, where quasi-analytical solutions allow for an exact quantification of computational speed versus numerical error. Specifically-tuned constrained Delaunay tetrahedralization (CDT) FE meshes lead to high accuracies for both the full subtraction and the direct potential approaches. Best accuracies are achieved by the full subtraction approach if the homogeneity condition is fulfilled. It is shown that the AMG-CG achieves an order of magnitude higher computational speed than the CG with the standard preconditioners with an increasing gain factor when decreasing mesh size. Our results should broaden the application of accurate and fast high-resolution FE volume conductor modeling in source analysis routine.
登录
查看更多内容
影响因子：
5.7
作者：
Rullmann M;Anwander A;Dannhauer M;Warfield SK;Duffy FH;Wolters CH
通讯作者：
Wolters CH
DOI：
10.1002/nme.2318
发表时间：
2008-08-13
影响因子：
2.9
作者：
Si, H.
通讯作者：
Si, H.
影响因子：
3.1
作者：
Haase, G;Reitzinger, S
通讯作者：
Reitzinger, S
DOI：
10.1016/s0013-4694(96)95698-9
发表时间：
1997-04-01
期刊：
ELECTROENCEPHALOGRAPHY AND CLINICAL NEUROPHYSIOLOGY
影响因子：
--
作者：
Buchner, H;Knoll, G;Pesch, J
通讯作者：
Pesch, J
影响因子：
4.6
作者：
Awada, KA;Jackson, DR;Papanicolaou, AC
通讯作者：
Papanicolaou, AC