CN103142211B - 一种基于极值圆的心脏磁信号处理方法 - Google Patents
一种基于极值圆的心脏磁信号处理方法 Download PDFInfo
- Publication number
- CN103142211B CN103142211B CN201110404119.1A CN201110404119A CN103142211B CN 103142211 B CN103142211 B CN 103142211B CN 201110404119 A CN201110404119 A CN 201110404119A CN 103142211 B CN103142211 B CN 103142211B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- extreme value
- heart
- magnetic
- section
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种基于极值圆的心脏磁信号处理方法,包括以下步骤:1)检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度,并同步测量心电图;2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据;3)对ST-T段内每一时刻6×6阵列检测点的磁感应强度进行三次样条插值处理,获取高分辨率的等磁场线图;4)根据等磁场线图获取极值圆和零磁场线;5)计算电生理参数L和θ。与现有技术相比,本发明只需要利用SQUID传感器检测到的ST-T段内的心脏磁场数据,在极值圆的范围内,就可以快速计算出辅助心脏疾病临床诊断的两个电生理参数及其呈阳性的参数大小和范围。
Description
技术领域
本发明涉及一种磁信号处理方法,尤其是涉及一种基于极值圆的心脏磁信号处理方法。
背景技术
七十年代,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%以上。人们希望通过相关理论的研究,不断探索新的适合临床应用的,具有高灵敏度和特异性的,计算速度快的方法。在一种基于极值圆的心脏磁信号处理新方法中,将计算两个电生理参数作为中间结果,辅助判断冠心病。
发明内容
本发明的目的是提供一种新的基于极值圆的可快速获取电生理参数的心脏磁信号处理方法。
本发明的目的可以通过以下技术方案来实现:
一种基于极值圆的心脏磁信号处理方法,包括以下步骤:
1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度,并同步测量心电图;
2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据;
3)对ST-T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理,获取高分辨率的等磁场线图;
4)处理等磁场线图确定极值圆和零磁场线。该极值圆的直径为等磁场线图中磁感应强度的极大值点与极小值点的连线,依次连接等磁场线图内的磁感应强度为零的点,得到零磁场线;
5)根据公式计算电生理参数L,其中L+为极值圆内磁感应强度大于零的面积,L-为极值圆内磁感应强度小于零的面积;
6)根据公式计算电生理参数θ,其中θn为零磁场线上相邻线段延长线的夹角。
所述的6×6阵列检测点中相邻检测点的间距为4cm,整个磁场检测平面大小为20×20cm。
所述的ST-T段的截止时候为心电图的T波峰处,ST-T段的起始时刻为所述T波峰向前1/3幅值处。
所述的步骤3)中三次样条插值后得到81×81点的等磁场线图。
与现有技术相比,本发明只需要利用SQUID传感器检测到的ST-T段内的心脏磁场数据,确定磁信号在ST-T段内数据的极值圆,并在极值圆的范围内快速计算出辅助心脏疾病诊断的两个电生理参数及其呈阳性的参数大小和范围。
附图说明
图1为本发明的流程图;
图2为本发明6×6阵列检测点的示意图;
图3为心电图中心电信号的ST-T段的曲线图;
图4为与心电图同步测量的心磁信号的ST-T段曲线图;
图5为三次样条插值后的等磁场线图,以及该等磁场线图中的极值圆;
图6为极值圆中零磁场线的示意图;
图7为正常人的ST-T段的电生理参数L的曲线图;
图8为正常人的ST-T段的电生理参数θ的曲线图;
图9为冠心病人的ST-T段的电生理参数L的曲线图;
图10为冠心病人的ST-T段的电生理参数θ的曲线图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
一种基于极值圆的心脏磁信号处理方法,该方法的流程如图1所示,包括以下步骤:
步骤S1:同步测量心磁信号和心电图,测量心磁信号时,受检者仰卧,多通道SQUID传感器阵列检测受检者胸腔表面如图2所示的6×6阵列检测点上的磁感应强度,由于每个相邻检测点的间距为4cm,整个检测平面大小为20×20cm。通过多通道SQUID传感器阵列可以实时记录垂直于人体胸腔平面Z方向上的磁感应强度Bz。
步骤S2:根据心电图截取心磁信号在ST-T段的数据。以心电图的T波峰值时刻作为对应心磁数据ST-T段的截止时刻tmax,以T波峰值向前的1/3幅值处作为心磁数据ST-T段的起始时刻tmin,如图3和图4所示。
步骤S3:对ST-T段内每一时刻6×6阵列检测点的磁感应强度进行三次样条插值处理,以获取分辨率较高的等磁场线图,方便后续处理。在本实施例中,插值点为81×81。
步骤S4:对等磁场线图进行处理。取磁感应强度的极大值点与极小值点的连线为直径,生成极值圆如图5所示。然后依次连接磁感应强度为零的点,获得零磁场线。
步骤S5:根据生成的极值圆和零磁场线,确定极值圆中磁感应强度大于零的面积L+,以及磁感应强度小于零的面积L-,通过公式计算出电生理参数L;极值圆中设有如图6所示的5个磁感应强度为零的点Z1、Z2、Z3、Z4和Z5,依次连接Z1、Z2、Z3、Z4和Z5,可以得到各线段延长线的夹角θ1、θ2、θ3,通过公式计算出电生理参数θ。
通过本发明处理通过计算机处理测得的心磁信号和心电图,步骤S2至步骤S4的成图过程均可通过制图软件自动完成,并自动获取图像的信息,处理得到电生理参数L和θ,可用于辅助判断心脏机能。正常人的ST-T段的电生理参数L和θ的曲线图分别如图7和图8所示。冠心病人的ST-T段的电生理参数L和θ的曲线图分别如图9和图10所示。正常人的L参数和θ参数小于冠心病患者。
Claims (4)
1.一种基于极值圆的心脏磁信号处理方法,其特征在于,包括以下步骤:
1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度,并同步测量心电图;
2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据;
3)对ST-T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理,获取高分辨率的等磁场线图;
4)根据等磁场线图获取极值圆和零磁场线,该极值圆的直径为等磁场线图中磁感应强度的极大值点与极小值点的连线,依次连接等磁场线图内磁感应强度为零的点,得到零磁场线;
5)根据公式计算电生理参数L,其中L+为极值圆内磁感应强度大于零的面积,L-为极值圆内磁感应强度小于零的面积;
6)根据公式计算电生理参数θ,其中θn为零磁场线上相邻线段延长线的夹角。
2.根据权利要求1所述的一种基于极值圆的心脏磁信号处理方法,其特征在于,所述的6×6阵列检测点中相邻检测点的间距为4cm。
3.根据权利要求1所述的一种基于极值圆的心脏磁信号处理方法,其特征在于,所述的ST-T段的截止时刻为心电图的T波峰处,ST-T段的起始时刻为所述T波峰向前1/3幅值处。
4.根据权利要求1所述的一种基于极值圆的心脏磁信号处理方法,其特征在于,所述的步骤3)中三次样条插值后得到81×81点的等磁场线图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110404119.1A CN103142211B (zh) | 2011-12-07 | 2011-12-07 | 一种基于极值圆的心脏磁信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110404119.1A CN103142211B (zh) | 2011-12-07 | 2011-12-07 | 一种基于极值圆的心脏磁信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103142211A CN103142211A (zh) | 2013-06-12 |
CN103142211B true CN103142211B (zh) | 2015-02-18 |
Family
ID=48540763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110404119.1A Expired - Fee Related CN103142211B (zh) | 2011-12-07 | 2011-12-07 | 一种基于极值圆的心脏磁信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103142211B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103810381A (zh) * | 2014-01-26 | 2014-05-21 | 同济大学 | 一种基于多通道心磁图的冠心病特征参数提取方法 |
CN104794474B (zh) * | 2015-03-26 | 2017-12-05 | 同济大学 | 一种基于心磁图的综合最大电流密度向量特征提取方法 |
CN105105710A (zh) * | 2015-07-24 | 2015-12-02 | 中国科学院上海微系统与信息技术研究所 | Squid 全张量测量模块、心磁信号探测装置及方法 |
CN105212897B (zh) * | 2015-08-25 | 2017-10-03 | 中国人民解放军第三〇九医院 | 一种基于心磁特征确定心脏磁场极大值点的方法 |
CN105606146A (zh) * | 2016-03-10 | 2016-05-25 | 江苏泰斯特电子设备制造有限公司 | 一种非线性修正技术 |
CN108937907A (zh) * | 2017-05-26 | 2018-12-07 | 北京小米移动软件有限公司 | 心率的采集方法及装置 |
CN109864733A (zh) * | 2019-01-16 | 2019-06-11 | 漫迪医疗仪器(上海)有限公司 | 心脑异常的检测方法、系统、介质及设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102077108B (zh) * | 2008-04-28 | 2015-02-25 | 康奈尔大学 | 分子mri中的磁敏度精确量化 |
-
2011
- 2011-12-07 CN CN201110404119.1A patent/CN103142211B/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
张磊 等.基于单电流偶极子的心脏磁场反演计算.《现代科学仪器》.2005,第2005年卷(第4期), * |
王春宇.心磁图在缺血性心脏病的诊断价值探讨.《中国优秀博硕士学位论文全文数据库(博士) 医药卫生科技辑》.2006,第2006年卷(第12期),E062-10. * |
Also Published As
Publication number | Publication date |
---|---|
CN103142211A (zh) | 2013-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103142211B (zh) | 一种基于极值圆的心脏磁信号处理方法 | |
WO2019015017A1 (zh) | 一种心电动力学数据量化分析方法 | |
CN107361760B (zh) | 心磁图仪诊断系统及应用该系统的心磁图仪 | |
CN108577825B (zh) | 一种心肌缺血病变位置定位方法及系统、存储介质及终端 | |
CN112842343B (zh) | 一种基于serf原子磁强计的心磁三维源位置估计方法 | |
CN109498013A (zh) | 基于单层eit电极边缘效应的膀胱体积测量方法 | |
CN103142223B (zh) | 一种基于极值差的心脏磁信号处理方法 | |
CN104794474B (zh) | 一种基于心磁图的综合最大电流密度向量特征提取方法 | |
JP3712348B2 (ja) | 心房粗動および心房細動の心臓磁界診断装置およびその作動方法 | |
CN103810381A (zh) | 一种基于多通道心磁图的冠心病特征参数提取方法 | |
JP2007268034A (ja) | 生体信号計測方法及び装置 | |
CN114403812B (zh) | 脑损伤状况的辅助分析方法、装置、系统及存储介质 | |
US7805179B2 (en) | Method of examining dynamic cardiac electromagnetic activity and detection of cardiac functions using results thereof | |
CN107374610B (zh) | 心磁图生成方法及生成系统 | |
CN113288102B (zh) | 一种无创监测脑血流的系统 | |
Udovychenko et al. | Current density distribution maps threshold processing | |
CN106236093B (zh) | 高精度宽频带分布式医用电阻抗层析成像装置及其应用 | |
CN114287943A (zh) | 一种基于可穿戴式心磁三维测量装置的三维成像方法 | |
JP2002028145A (ja) | 心室遅延電位の心臓磁界診断装置および心筋内興奮不均一伝播部位の同定方法 | |
JP4791797B2 (ja) | 生体磁場計測装置 | |
CN110916724A (zh) | 一种基于闭环最短路径的b超图像胎儿头围检测方法 | |
CN105559770B (zh) | 内探式生物组织磁分布图测量装置 | |
Zhao et al. | An integrated maximum current density approach for noninvasive detection of myocardial infarction | |
JP4342644B2 (ja) | 心内電気現象診断装置 | |
JP2795211B2 (ja) | 生体磁気計測装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150218 Termination date: 20171207 |
|
CF01 | Termination of patent right due to non-payment of annual fee |