论文题目: 基于脑磁图的电流偶极子反演研究
论文类型: 博士论文
论文专业: 光学工程
作者: 朱红毅
导师: 何赛灵
关键词: 脑磁图,多电流偶极子,球模型,真实头模型,遗传算法,算法,边界元方法
文献来源: 浙江大学
发表年度: 2005
论文摘要: 脑磁图（Magnetoencephalography, MEG）是一种通过测量神经电流在头外产生的磁场来推断头内脑活动的无损脑功能,脑成像和脑疾病诊断的技术。脑磁图的时间分辨率小于1毫秒,定位皮层神经活动的空间分辨率在2-4毫米左右。利用测量的脑磁场(约50-500飞特斯拉;1fT=10-15T)分布,在合适的源模型下,可以唯一地确定皮层上因受到刺激而产生活动的神经元位置。脑磁图可用以研究包括脑的自发和诱发活动在内的多种脑神经活动过程,在临床上用以对癫痫病灶的定位。脑磁图的研究对于人类认识脑,开发脑,治愈脑疾病都有重要意义。 脑磁图研究中主要有两种源模型,一种是电流分布模型,一种是多偶极子模型,本文主要研究了多偶极子模型。 本论文由九章组成,前四章分别对脑磁图的研究概况,脑磁信号的神经生理学基础,脑磁正问题和脑磁逆问题做了较详细的介绍。在其中也包含了在本文中使用的边界元方法,遗传算法,Marquardt算法和各种头模型的构建方法的介绍。 在第五章中,我们首先介绍了一种基于对易理论快速求解脑磁引导场的方法。相比标准边界元方法,这种方法在计算速度和所需的存储空间上都有较大的优势。此外,采用这种方法,对于较浅的源,仍然能够保持很高的计算精度,而采用传统的边界元方法,对于这种源如果不对它附近的表面采用局部细化的网格划分的话（refine）,计算将产生很大的误差。然后我们利用上述计算引导场的方法建立了一个引导场网格,并通过三次插值获得不处于网格节点处的引导场值。基于这种引导场网格,可将脑磁正逆问题的求解速度相对标准边界元方法提高104量级。（不考虑只需要计算一次的引导场网格的计算时间） 在第六章中我们介绍了一种联合使用全局优化方法（遗传算法）与局部优化算法（Marquardt）的混合优化方法。首先通过遗传算法进行粗略搜索,并将遗传算法求解的值作为Marquardt算法的初值进行进一步的详细搜索,获得最终的解。采用这种方法,可以成功的解决局域优化方法会陷入局部极小值的问题,同时又可以解决遗传算法在后期收敛很慢的问题,提高搜索效率。为了进一步提高求解的速度,我们尝试了在需要较多计算时间的全局优化阶段采用计算速度很快的球模型,在局域优化阶段采用更精确的真实头模型。通过测试,我们发现对于源较少的情况这种方法可以获得很好的效果,但对于源较多的情况,由于情况太复杂,这种混合优化混合头模型的方法无法在获得令人满意的定位结果,但混合优化,纯真实头模型仍然可以较好的工作。此外,在算法中我们通过引入”局部最优偶极子”,将偶极子的强度用引导场表示,从而在目标函数中只含有位置参数,减少了需要求解的参数。我们也提出了一种迭代方法去求解“局部最优偶极子”。我们比较了这种方法和常用的奇异值分解方法的抗噪声能力,结果显示我们的迭代方法在强噪声下具有更好的稳定性。 在第七章中我们继续研究了混合头模型的使用。通过逐点比较源处于脑内不同位置
论文目录:
摘要
ABSTRACT
目录
第1章 绪论
§1.1 引言
§1.2 脑磁图研究概况
§1.3 脑磁图的医学应用
§1.4 本论文的思路和主要工作
第2章 脑磁产生的神经生理学基础
§2.1 大脑的组成
§2.2 脑神经
§2.3 离子动力学过程
§2.4 突触后电位(POSTSYNAPTIC POTENTIALS, PSP)
§2.5 动作电位(ACTION POTENTIALS)
§2.6 人脑结构及功能区
第3章 脑磁正问题基本电磁理论
§3.1 脑磁正问题的解析解
§3.2 球模型下脑磁正问题的解
§3.3 真实头模型下的正向计算
§3.4 基于时空数据的矩阵表示
第4章 脑磁逆问题
§4.1 重建方法
§4.1.1 最小范数解(Minimum Norm Solution, MNS)
§4.1.2 加权最小范数解(Weighted Minimum Norm Solution, WMNS)
§4.1.3 最大熵方法(Maximum Entropy Method, MEM)
§4.2 偶极子定位法
§4.2.1 构造目标函数
§4.2.2 优化方法
§4.2.3 MEG-MUSIC算法
第5章 根据对易理论求解脑磁图引导场
§5.1 引言
§5.2 理论导出
§5.3 引导场方程的数值解
§5.4 模拟结果与讨论
§5.4.1 所需存储空间
§5.4.2 计算精度比较
5.4.3计算精度比较
§5.5 网格法
§5.5.1 方法介绍
§5.5.2 计算结果
§5.6 本章小结
第6章 利用结合局域优化的遗传算法求解脑磁逆问题
§6.1.遗传算法
§6.1.1 遗传算法的基本原理
§6.1.2 遗传算法的特点
§6.2 结合局域优化的遗传算法
§6.2.1 对于瞬时脑磁的仿真
§6.2.2 基于时空数据的仿真计算
§6.3 一种决定偶极子最佳指向与瞬时行为的迭代法方法
§6.3.1 方法
§6.3.2 结果和讨论
§6.4 本章小结
第7章 利用混合头模型求解脑磁逆问题
§7.1 引言
§7.2 混合模型
§7.3 LDA中的快速算法
§7.4 偶极子源定位
§7.5 本章小结
第8章 多重信号分类方法在脑磁图中的应用
§8.1 引言
§8.2 算法
§8.3 仿真计算与讨论
§8.4 本章小结
第9章 总结与展望
参考文献
博士在读期间发表和撰写的论文
致谢
发布时间: 2005-07-27
参考文献
- [1].脑磁逆问题中的磁源重建算法研究[D]. 胡净.浙江大学2002
相关论文
- [1].生物医学电磁正反问题数值计算方法研究[D]. 寿国法.浙江大学2009
- [2].EEG/MEG反问题数值方法研究[D]. 王红艳.吉林大学2010
- [3].脑磁逆问题中的磁源重建算法研究[D]. 胡净.浙江大学2002
- [4].磁源性影像关于癫痫灶定位及脑功能区定位的研究[D]. 孙吉林.河北医科大学2002
- [5].棋盘格翻转图形刺激的脑磁图研究[D]. 张淑倩.河北医科大学2003
- [6].脑电逆问题求解的数值计算方法研究[D]. 李颖.河北工业大学2003
- [7].基于非线性优化及独立分量分析的脑电源定位研究[D]. 邹凌.浙江大学2005
- [8].有限元素法在三维各向异性脑电正问题及皮层成像技术中的应用[D]. 张迎春.浙江大学2005
- [9].脑磁图及磁源性影像在体感、听觉及癫痫灶定位中的应用[D]. 孙吉林.中国人民解放军军医进修学院2006
- [10].脑磁源定位技术研究[D]. 霍小林.浙江大学2001