心脏磁场分布电流源重构及其精度分析.pdfVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 14 (2013) 148703 心脏磁场分布电流源重构及其精度分析* 1 1 1 2 2 1† 王伟远 赵晨 林玉章 张树林 谢晓明 蒋式勤 1) ( 同济大学电子与信息工程学院控制科学与控制工程系, 上海 201804 ) 2) ( 中国科学院上海微系统与信息技术研究所, 上海 200050 ) ( 2013年2月25 日收到; 2013年3月27 日收到修改稿) 用测量到的心脏磁场信号重构其电流源是一种无创揭示心脏电活动的方法. 心脏电活动的时空信息可视化, 将 有助于心脏功能的研究和心脏疾病诊断. 本文通过仿真实验研究了一种减时窗波束形成器对分布时变电流源的重 构能力, 以及源重构结果与心室兴奋传播的关系. 采用元胞自动机模拟心室的兴奋传播, 产生分布随时间变化的电 流源, 并用边界元法构建了一个心脏- 躯干模型, 模拟体电导的作用. 仿真结果表明, 这种减时窗波束形成器能够重 构分布时变电流源, 并达到较好的精度. 将该方法和无穷大均匀介质导联矩阵用于一例正常人的心脏磁场信号分析, 其结果可以反映心室兴奋传播的基本特征. 关键词: 心磁图, 源重构, 边界元, 波束形成器 PACS: 87.85.−d, 87.85.Pq, 87.85.Tu, 87.85.Ng DOI: 10.7498/aps.62.148703 法, 曾被用于脑电图(electroencephalography, EEG) 1 引言 和脑磁图 (magnetoencephalography, MEG) 的源估 11 计. 1997 年, van Veen 等 将线性约束最小方差 用心脏磁场测量信号重构心脏的电流源, 也称 (linear constraint minimum variance, LCMV) 波束形 之为心脏电流源重构(current source reconstruction, 成器引入脑电研究. 1999 年, Robinson 和Vrba12 CSR), 是一种无创揭示心脏电活动的方法1−3 . 其 在LCMV 波束形成器的基础上, 利用测量噪声的统 目的是通过求解心脏磁场逆问题, 使心脏电活动的 计特性对测量磁场的协方差矩阵进行处理, 提出了 时空信息可视化, 从而无创研究心脏功能和诊断心 一种合成孔径磁场测定法(synthetic aperture mag- 脏疾病. netometry, SAM), 并用于脑的功能源定位. 2001 年, 目前, 我国和许多发达国家已经能够通过多通 13 Sekihara 等 在研究LCMV 和 Borgiotti-Kaplan 道超导量子干涉器(superconducting quantum inter- 波束形成器的基础上, 提出了一种特征空间投影 ference device, SQUID) 完全无创地、非接触地测 (eigenspace projection, EP) 的方法, 能够有效地抑 量微弱的心脏磁场(约为10−12 T)3−6 , 得到记录体 14