基于贪婪稀疏方法的心脏磁场源重构 - 物理学报.PDF

物理学报 Acta Phys. Sin. Vol. 62, No. 11 (2013) 118703 基于贪婪稀疏方法的心脏磁场源重构* 1 1 2 1† 邴璐 王伟远 王永良 蒋式勤 1) ( 同济大学电子与信息工程学院, 控制科学与控制工程系, 上海 201804 ) 2) ( 中国科学院上海微系统与信息技术研究所, 上海 200050 ) ( 2012 年12 月29 日收到) 心脏磁场源重构, 即通过人体胸腔表面的磁场阵列测量信号反演产生该磁场的电流源分布, 是一种无创地研究 心脏电活动的方法. 本文用线性化方程描述测量磁场与心脏内部电流源的关系, 并通过一种贪婪优化的方法得到了 源重构的稀疏解. 通过近似正交化过程和改变迭代算法中原子的选择方式, 降低了优化算法的复杂度, 可在保证源重 构精度的情况下, 快速搜索源空间中强度比较大的位置. 文中通过一组正常人的心磁测量数据, 说明了源重构方法的 有效性. 这组数据的研究结果表明, 强度大于65% 的重构电流源的空间移动轨迹与心脏除极和复极电活动的传导过 程基本符合. 其QRS 和ST-T 段的实测心磁图与重构电流源产生的磁场图的平均拟合优度分别为99.36%和99.78%. 关键词: 心磁图, 逆问题, 稀疏解, 正交匹配追踪 PACS: 87.85.−d, 87.85.Pq, 87.85.Tu, 87.85.Ng DOI: 10.7498/aps.62.118703 磁场的激励源. 本文采用固定位置的分布源模型, 1 引言 将计算磁场源参数转化为求解偶极矩的逆问题. 但 是, 这种基于分布源模型的线性方程中, 未知源的 研究心脏电活动的方法之一是用人体胸腔表 数量往往超过磁场测点的个数, 理论上没有唯一 面测量到的心脏磁场数据重构产生该磁场的电 解. 因此, 需要施加各种约束条件, 并采用优化方法 流源. 这种方法也被称作心脏电流源重构(cardiac 67 求解 . 1 source reconstruction) 或心电活动的三维磁成像 . 稀疏信号重构是一种通过有限数据求局部能 目的是通过心脏电活动的可视图像, 无创地研究其 8 量的优化方法 . 被广泛应用于各种信号估计, 信 心脏功能与诊断心脏疾病. 号分类, 信号处理, 图像压缩和生物医学成像. 简单 心脏磁场与脑磁场的源重构具有相似的理 说, 稀疏信号(sparse signal) 是指一个信号向量中只 论基础. 一些重要的源重构方法, 比如, 空间滤 有有限个非零元素. 稀疏分解方法(sparse decompo- 23 4 波法 和稀疏分解方法 , 最早均用于脑电/磁 sition) 可以将一组信号在给定变换域(过完备字典) (EEG/MEG) 研究. 脑的外部磁场信号是由头颅内特 中, 通过其中向量的线性组合, 得到信号最稀疏的 殊部位的细胞电流产生的. 这些特殊部位数量不多, 表达. 体积很小. 然而, 心脏内部的细胞电流较为复杂, 电 1995 年