专利名称：心磁信号噪声抑制中的信号选择性平均方法
技术领域：
本发明涉及一种心磁信号噪声抑制中的信号选择性平均方法，更确切地说，本发明涉及一种应用于心磁信号噪声抑制的信号选择性平均方法。
背景技术：
超导量子干涉器件(SuperconductingQuantum Interference Device, SQUID)是目前已知最灵敏的磁通传感器，典型低温超导SQUID的磁场灵敏度为3_5fT/sqrt (Hz)。作为SQUID的一个重要应用领域，经临床研究证明，心磁系统在心脏疾病诊断、功能研究等方面有独特的应用潜力[V. Pizzela etal, Supercond. Sc1. Technol. 14(2001)R79-R114。心磁系统采集获得的信号包括有用信号及噪声干扰，其中噪声包括白噪声、非白噪声、外界环境干扰，以及一些来自信号本身的变化、局部信号异常等。经过预处理后，心磁信号所得的结果中仍然存在一些‘顽固’噪声，如脉冲、随机出现的扰动及局部显著的信号毛刺等。为了消除这些干扰噪声，多周期平均处理得到了广泛地使用。平均是生物医学信号处理过程中常用的提高周期或准周期信号质量的最常用方法之一，在心磁信号处理中被广泛使用。在心磁信号平均处理中，通过采集同步心电信号，以心电信号的R峰为基准进行截断平均。目前，最常用的方法主要是对全段信号进行平均处理，该方法使得这些顽固噪声引入到最后的平均结果中，造成信号失真。除此之外还有聚类分析法DDiPietroPaolo, H-P Muller and S N Erne, Phys. Med.Biol. 50(2005) 2415-2426，它是将每一个信号段按照信号特征划分为树形结构，但是由于过于精细造成复杂度较大，效率较低。为此，探寻高效率、强实用性、易操作的心磁信号平均方法是本领域技术人员努力追求的目标之一，同时在实际应用和分析方面具有十分重要的意义，构筑成本发明的构思。
发明内容
本发明的目的在于提供一种心磁信号噪声抑制中的信号选择性平均方法，所提供的方法是一种在信号噪声抑制中的信号选择性平均方法，所述的方法基于模版匹配的方式，通过计算模版信号与被检信号之间的相似系数，经阈值判断剔除部分干扰周波进行平均，以实现可选择信号平均。本发明所述的信号选择性平均方法是(I)查找与待处理心磁信号对应的参考心电信号R峰，以对应心电信号R峰为基准，将心磁信号分割为各个独立周波；进行心磁的单周波分段；(2)将待处理信号进行全平均处理，获取的平均信号结果作为初级模版；(3)计算每一原始分段周波信号与初级模版之间的相似系数，与设定的初阈值进行比较，识别并剔除波形变化较大的周波段，将保留下来的周波段进行平均，获得最终模版信号;相似系数Y代表两者之间的相似程度，定义为
y =1- δδ为差值系数，反映了当前信号与模版信号在各时间点上的形状差异
权利要求
1.一种心磁信号噪声抑制中的信号选择性平均方法，其特征在于所述的方法基于模版匹配的方式，通过计算模版信号与被检信号之间的相似系数，经阀值判断剔除部分干扰周波进行平均，以实现可选择信号平均。
2.按权利要求1所述的方法，其特征在于具体步骤是(1)查找与待处理心磁信号对应的参考心电信号R峰，以对应心电信号R峰为基准，将心磁信号分割为各个独立周波；进行心磁的单周波分段；(2)将待处理信号进行全平均处理，获取的平均信号结果作为初级模版；(3)计算每一原始分段周波信号与初级模版之间的相似系数，与设定的初阈值进行比较，识别并剔除波形变化较大的周波段，将保留下来的周波段进行平均，获得最终模版信号;(4)重复步骤(3)，模版匹配以最终模版代替初级模版获取最终的平均心磁信号。
3.按权利要求1或2所述的方法，其特征在于所述的相似系数Y代表两者之间的相似程度，定义为Y =1- δδ为差值系数，反映了当前信号与模版信号在各时间点上的形状差异
4.按权利要求3所述的方法，其特征在于在波形识别过程中，相似系数对信号的幅度变化及波宽更为敏感，对于波形变化的识别能力更强。
5.按权利要求2所述的方法，其特征在于设定的初阀值大小决定信号选择剔除的程度，直接影响模版信号的优劣，采用相对阈值，定义为Gate = μ X Y mean, O ^ μ ^ I式中，Y m_为相似系数序列的均值，μ为强度因子，根据实际数据的情况进行选择。
6.按权利要求2所述的方法，其特征在于步骤(2)对待处理信号进行全平均处理即对应点叠加平均。
全文摘要
本发明涉及一种心磁信号噪声抑制中的信号选择性平均方法，其特征在于所述的方法基于模版匹配的方式，通过计算模版信号与被检信号之间的相似系数，经阀值判断剔除部分干扰周波进行平均，以实现可选择信号平均。本方法的特点是利用模版匹配剔除信号中存在噪声干扰或异常信号的周波，避免对平均信号结果引入不必要的噪声引发信号失真，灵活度大，适用于实际采集心磁信号过程中遇到的各种干扰问题。
文档编号A61B5/0402GK102988041SQ20121046549
公开日2013年3月27日 申请日期2012年11月16日 优先权日2012年11月16日
发明者刘明, 张树林, 曾佳, 邱阳, 李华, 王永良, 孔祥燕, 谢晓明 申请人:中国科学院上海微系统与信息技术研究所