基于贪婪稀疏方法的心脏磁场源重构

心脏磁场源重构, 即通过人体胸腔表面的磁场阵列测量信号反演产生该磁场的电流源分布, 是一种无创地研究心脏电活动的方法. 本文用线性化方程描述测量磁场与心脏内部电流源的关系, 并通过一种贪婪优化的方法得到了源重构的稀疏解. 通过近似正交化过程和改变迭代算法中原子的选择方式, 降低了优化算法的复杂度, 可在保证源重构精度的情况下, 快速搜索源空间中强度比较大的位置. 文中通过一组正常人的心磁测量数据, 说明了源重构方法的有效性. 这组数据的研究结果表明, 强度大于65%的重构电流源的空间移动轨迹与心脏除极和复极电活动的传导过程基本符合. 其QRS 和ST-T 段的实测心磁图与重构电流源产生的磁场图的平均拟合优度分别为99.36%和99.78%. 关键词： 心磁图 逆问题 稀疏解 正交匹配追踪

MCG source reconstruction based on greedy sparse method

Current source reconstruction, i.e., reconstructing current dipole distribution through measured array signals of cardiac magnetic field on body surface, is a method for non-invasively study on the heart electrical activity. In this paper, the relationship between measured magnetic signals and current dipole distribution is described by a linear equation, and a sparse solution of current source reconstruction is achieved using a fast greedy method. This method can significantly decrease the computational complexity of or- thogonal matching pursuit (OMP) algorithm by means of approximating orthogonalisation and improving the selection vector strategy per iteration. Thereby, the sources with large dipole strength can be fast searched out with high accuracy. A set of magnetocardiogram (MCG) data of normal subject is used to demonstrate the effectiveness of this method that the trajectory of reconstructed dominant sources, whose strengths are more than 65%, is almost consistent with conduction process in depolarization and repolarization. The average goodness of fit (GOFs) of measured MCG and the magnetic field map generated by the reconstructed current sources during QRS complex and ST-T segment are 99.36% and 99.78%, respectively.
Keywords： magnetocardiography (MCG) inverse problem sparse solution orthogonal matching pursuit (OMP)