一种基于极值差的心脏磁信号处理方法技术

本发明专利技术涉及一种基于极值差的心脏磁信号处理方法，其特征在于，包括以下步骤：1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度，并同步测量心电图；2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据；3)对ST-T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图；4)根据等磁场线图获取极大值Bzmax和极小值Bzmin；5)根据公式计算电生理参数JS，并生成JS曲线。与现有技术相比，本发明专利技术只需要利用SQUID传感器检测到的ST-T段内的心脏磁场数据的极大值和极小值，就可以快速计算出辅助心脏疾病临床诊断的电生理参数，并分析心脏电活动的功能。

【技术实现步骤摘要】

本专利技术涉及一种磁信号处理方法，尤其是涉及。
技术介绍
七十年代， D.Cohen等人在实验室首次用超导量子干涉器(SQUID)测量到人体心脏磁场。1976年，心磁仪的开发者相继提出了用心磁图和伪电流密度图(或箭头图，也称之为Hosaka-Cohen变换)诊断心肌缺血的理论与方法。九十年代，美国CMI公司推出了用单磁偶极子算法诊断心肌缺血等心脏疾病的9通道心磁仪.2005年德国J.W.Park等人提出了用心磁仪的测量数据预测冠心病(CAD)的方法。2006年，德国W.Haberkorn等在伪电流密度图的基础上，提出了具有电生理意义的伪电流密度成像方法。同年，K.Tolstrup等也提出了一种快速磁成像检测心肌缺血的方法。2007年，台湾与韩国合作研究机构提出了一种用心磁T波信号的二维传播成像及面积比方法，并用来诊断心肌缺血等疾病。同年，日本提出了一种用心磁图JT段积分值筛选冠心病的方法.A.Gapelyuk等也提出了一种用心脏磁场图探测CAD的方法。这些诊断方法的灵敏度和特异性在60%-80%左右。近年来，随着心脏磁场检测技术的发展，心磁测量设备性能的不断提高，无创诊断心脏疾病的方法也有了一定的进步。2008年,P.V.Leeuwen等定量分析比较了三种无创诊断冠心病的方法。2010年，Kwon等和A.Gapelyuk等分别用心磁图分类识别,用KL熵和剩余参数两种方法的组合，将这些诊断方法的灵敏度和特异性提高到80%以上。人们希望通过相关理论的研究，不断探索新的适合临床应用的，具有高灵敏度和特异性的，计算速度快的信号处理方法。在一种基于极值差的心脏磁信号处理新方法中，需要通过计算电生理参数JS作为中间结果，辅助判断冠心病。
技术实现思路
本专利技术的目的是提供一种新的基于极值差的可快速获取电生理参数的心脏磁信号处理方法。本专利技术的目的可以通过以下技术方案来实现:，包括以下步骤:I) SQUID传感器检测心脏磁场垂直于胸腔平面的6X6阵列检测点的磁感应强度，并同步测量心电图；2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据；3)对ST-T段内每一时刻6X6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图；4)根据等磁场线图获取极大值和极小值，该等磁场线图中磁感应强度的极大值点与极小值点的连线为D ；β _ β5)根据公式Μ =|,|Γ I)计算电生理参数Js，并生成Js曲线。所述的6X6阵列检测点中相邻检测点的间距为4cm，整个磁场检测平面大小为20X20cm。所述的ST-T段的截止时候为心电图的T波峰处，ST-T段的起始时刻为所述T波峰向前1/3幅值处。所述的步骤3)中三次样条插值后得到81X81点的等磁场线图。与现有技术相比，本专利技术只需要利用SQUID传感器检测到的ST-T段内的心脏磁场数据，根据等磁场线图中的极大值和极小值，就可以快速计算出辅助心脏疾病诊断的电生理参数JS,并分析心脏电活动的功能。附图说明图1为本专利技术的流程图；图2为本专利技术6X6阵列检测点的示意图；图3为心电图中心电信号的ST-T段的曲线图；图4为与心电图同步测量的心磁信号的ST-T段曲线图；图5为三次样条插值后的等磁场线图；图6为磁场轮廓连线示意图；图7为正常人的ST-T段的电生理参数JS的曲线图；图8为冠心病人的ST-T段的电生理参数JS的曲线图。具体实施例方式下面结合附图和具体实施例对本专利技术进行详细说明。实施例，该方法的流程如图1所示，包括以下步骤:步骤S1:同步测量心磁图和心电图，测量心磁信号时，受检者仰卧，多通道SQUID传感器阵列检测受检者胸腔表面如图2所示的6X6阵列检测点上的磁感应强度，由于每个相邻检测点的间距为4cm，整个检测平面大小为20X20cm。通过多通道SQUID传感器阵列可以实时记录垂直于人体胸腔平面Z方向上的磁感应强度Bz。步骤S2:根据心电图截取心磁信号在ST-T段的数据。以心电图的T波峰值时刻作为对应心磁数据ST-T段的截止时刻tmax，以T波峰值向前的1/3幅值处作为心磁数据ST-T段的起始时刻tmin，如图3和图4所示。步骤S3:对ST-T段内每一时刻6X6阵列检测点的磁感应强度进行三次样条插值处理，以获取分辨率较高的等磁场线图，如图5所示。在本实施例中，插值点为81X81。步骤S4:根据等磁场线图获取极大值和极小值。通过磁感应强度的极大值点与极小值点的连线D，可以得到与磁场极大值和极小值有关的磁场轮廓线示意图。如图6所示。即通过图5中连线D的平面切割等磁场线图，可以得到如图6所示的Bzm曲线。当假设心脏磁场是由单个等效偶极子源产生时， 根据Biot-Savart定律，该磁场可以分为两部分。其中一部分B00OO是由等效偶极子源在无穷大均匀介质中产生的；另一部分Ie是由体积电流σ E产生的。在极值点处，由Z方向的测量值Bzmax，Bzmin可以计算得到B_。Bz^的大小表征了对应时刻心脏中体积电流产生的磁感应强度。Βζμ曲线表示了该体积电流产生的磁感应强度的时变特性。步骤S5:根据等磁场线图的极大值和极小值，和公式本文档来自技高网...

【技术保护点】
一种基于极值差的心脏磁信号处理方法，其特征在于，包括以下步骤：1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度，并同步测量心电图；2)根据心电图截取心脏磁场的磁感应强度在ST?T段的数据；3)对ST?T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图；4)根据等磁场线图获取极大值Bzmax和极小值Bzmin；5)根据公式计算电生理参数JS，并生成JS曲线。FDA0000117101080000011.tif

【技术特征摘要】
1.一种基于极值差的心脏磁信号处理方法，其特征在于，包括以下步骤: DSQUID传感器检测心脏磁场垂直于胸腔平面的6X6阵列检测点的磁感应强度，并同步测量心电图； 2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据； 3)对ST-T段内每一时刻6X6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图； 4)根据等磁场线图获取极大值Bzmax和极小值Bzmin；5)根据公式2.根...

【专利技术属性】
技术研发人员：蒋式勤，赵晨，石明伟，周大方，顾嘉期，
申请(专利权)人：同济大学，
类型：发明
国别省市：