CN103142223A - 一种基于极值差的心脏磁信号处理方法 - Google Patents

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

Info

Publication number
CN103142223A
CN103142223A CN2011104041613A CN201110404161A CN103142223A CN 103142223 A CN103142223 A CN 103142223A CN 2011104041613 A CN2011104041613 A CN 2011104041613A CN 201110404161 A CN201110404161 A CN 201110404161A CN 103142223 A CN103142223 A CN 103142223A
Authority
CN
China
Prior art keywords
magnetic field
heart
section
magnetic
signal processing
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.)
Granted
Application number
CN2011104041613A
Other languages
English (en)
Other versions
CN103142223B (zh
Inventor
蒋式勤
赵晨
石明伟
周大方
顾嘉期
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN201110404161.3A priority Critical patent/CN103142223B/zh
Publication of CN103142223A publication Critical patent/CN103142223A/zh
Application granted granted Critical
Publication of CN103142223B publication Critical patent/CN103142223B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明涉及一种基于极值差的心脏磁信号处理方法，其特征在于，包括以下步骤：1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度，并同步测量心电图；2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据；3)对ST-T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图；4)根据等磁场线图获取极大值Bzmax和极小值Bzmin；5)根据公式计算电生理参数JS，并生成JS曲线。与现有技术相比，本发明只需要利用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％以上。人们希望通过相关理论的研究，不断探索新的适合临床应用的，具有高灵敏度和特异性的，计算速度快的信号处理方法。在一种基于极值差的心脏磁信号处理新方法中，需要通过计算电生理参数JS作为中间结果，辅助判断冠心病。
发明内容
本发明的目的是提供一种新的基于极值差的可快速获取电生理参数的心脏磁信号处理方法。
本发明的目的可以通过以下技术方案来实现：
一种基于极值差的心脏磁信号处理方法，包括以下步骤：
1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度，并同步测量心电图；
2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据；
3)对ST-T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图；
4)根据等磁场线图获取极大值和极小值，该等磁场线图中磁感应强度的极大值点与极小值点的连线为D；
5)根据公式
Figure BDA0000117101090000021
计算电生理参数JS，并生成JS曲线。
所述的6×6阵列检测点中相邻检测点的间距为4cm，整个磁场检测平面大小为20×20cm。
所述的ST-T段的截止时候为心电图的T波峰处，ST-T段的起始时刻为所述T波峰向前1/3幅值处。
所述的步骤3)中三次样条插值后得到81×81点的等磁场线图。
与现有技术相比，本发明只需要利用SQUID传感器检测到的ST-T段内的心脏磁场数据，根据等磁场线图中的极大值和极小值，就可以快速计算出辅助心脏疾病诊断的电生理参数JS，并分析心脏电活动的功能。
附图说明
图1为本发明的流程图；
图2为本发明6×6阵列检测点的示意图；
图3为心电图中心电信号的ST-T段的曲线图；
图4为与心电图同步测量的心磁信号的ST-T段曲线图；
图5为三次样条插值后的等磁场线图；
图6为磁场轮廓连线示意图；
图7为正常人的ST-T段的电生理参数JS的曲线图；
图8为冠心病人的ST-T段的电生理参数JS的曲线图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
一种基于极值差的心脏磁信号处理方法，该方法的流程如图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阵列检测点的磁感应强度进行三次样条插值处理，以获取分辨率较高的等磁场线图，如图5所示。在本实施例中，插值点为81×81。
步骤S4：根据等磁场线图获取极大值和极小值。通过磁感应强度的极大值点与极小值点的连线D，可以得到与磁场极大值和极小值有关的磁场轮廓线示意图。如图6所示。即通过图5中连线D的平面切割等磁场线图，可以得到如图6所示的BzσE曲线。当假设心脏磁场是由单个等效偶极子源产生时，根据Biot-Savart定律，该磁场可以分为两部分。其中一部分B∞(r)是由等效偶极子源在无穷大均匀介质中产生的；另一部分BσE是由体积电流σE产生的。在极值点处，由Z方向的测量值Bzmax，Bzmin可以计算得到BzσE。BzσE的大小表征了对应时刻心脏中体积电流产生的磁感应强度。BzσE曲线表示了该体积电流产生的磁感应强度的时变特性。
步骤S5：根据等磁场线图的极大值和极小值，和公式计算出电生理参数JS，并生成JS曲线。
通过本发明可快速处理得到电生理参数JS，用于辅助判断心脏机能。正常人的ST-T段的电生理参数JS的曲线图如图7所示。冠心病人的ST-T段的电生理参数JS的曲线图如图8所示。正常人与冠心病患者的JS参数有明显的差别。

Claims (4)

1.一种基于极值差的心脏磁信号处理方法，其特征在于，包括以下步骤：
1)SQUID传感器检测心脏磁场垂直于胸腔平面的6×6阵列检测点的磁感应强度，并同步测量心电图；
2)根据心电图截取心脏磁场的磁感应强度在ST-T段的数据；
3)对ST-T段内每一时刻6×6阵列检测点上的磁感应强度进行三次样条插值处理，获取高分辨率的等磁场线图；
4)根据等磁场线图获取极大值Bzmax和极小值Bzmin；
5)根据公式
Figure FDA0000117101080000011
计算电生理参数JS，并生成JS曲线。
2.根据权利要求1所述的一种基于极值差的心脏磁信号处理方法，其特征在于，所述的6×6阵列检测点中相邻检测点的间距为4cm。
3.根据权利要求1所述的一种基于极值差的心脏磁信号处理方法，其特征在于，所述的ST-T段的截止时刻为心电图的T波峰处，ST-T段的起始时刻为所述T波峰向前1/3幅值处。
4.根据权利要求1所述的一种基于极值差的心脏磁信号处理方法，其特征在于，所述的步骤3)中三次样条插值后得到81×81点的等磁场线图。
CN201110404161.3A 2011-12-07 2011-12-07 一种基于极值差的心脏磁信号处理方法 Expired - Fee Related CN103142223B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110404161.3A CN103142223B (zh) 2011-12-07 2011-12-07 一种基于极值差的心脏磁信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110404161.3A CN103142223B (zh) 2011-12-07 2011-12-07 一种基于极值差的心脏磁信号处理方法

Publications (2)

Publication Number Publication Date
CN103142223A true CN103142223A (zh) 2013-06-12
CN103142223B CN103142223B (zh) 2015-02-18

Family

ID=48540775

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110404161.3A Expired - Fee Related CN103142223B (zh) 2011-12-07 2011-12-07 一种基于极值差的心脏磁信号处理方法

Country Status (1)

Country Link
CN (1) CN103142223B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103810381A (zh) * 2014-01-26 2014-05-21 同济大学 一种基于多通道心磁图的冠心病特征参数提取方法
CN109864733A (zh) * 2019-01-16 2019-06-11 漫迪医疗仪器（上海）有限公司 心脑异常的检测方法、***、介质及设备
CN110074774A (zh) * 2019-04-28 2019-08-02 漫迪医疗仪器（上海）有限公司 基于心磁图的心脏室间隔异常的分析方法、***、介质及终端
CN112932492A (zh) * 2021-01-29 2021-06-11 上海跃磁生物科技有限公司 心磁图仪多个扫描数据同步方法、***、介质及装置

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11253412A (ja) * 1998-03-17 1999-09-21 Hitachi Ltd 生体磁場計測装置
JP2001095771A (ja) * 2000-08-09 2001-04-10 Hitachi Ltd 磁場計測装置
US6275719B1 (en) * 1998-09-09 2001-08-14 Hitachi, Ltd. Biomagnetic field measurement apparatus
US6522908B1 (en) * 1999-10-06 2003-02-18 Hitachi, Ltd. Biomagnetic field measuring apparatus
US20030149354A1 (en) * 2001-08-23 2003-08-07 Bakharev Alexander A. Ischemia identification, quantification and partial localization MCG
CN1471374A (zh) * 2000-08-29 2004-01-28 心磁图检查中的局部缺血识别、量化和局部定位
CN1552285A (zh) * 2003-05-26 2004-12-08 中国科学院物理研究所 一种心磁波组图及其分析方法
CN1556687A (zh) * 2001-09-21 2004-12-22 使用小波变换的心磁图非线性噪声降低
US20050148844A1 (en) * 2003-09-10 2005-07-07 Kuniomi Ogata Biomagnetic measurement apparatus
JP2008200331A (ja) * 2007-02-21 2008-09-04 Hitachi High-Technologies Corp 生体磁場計測装置
CN101361655A (zh) * 2008-09-12 2009-02-11 北京大学 一种心磁测量区域极值点的测定方法
CN101836860A (zh) * 2010-06-10 2010-09-22 北京大学 一种心脏磁场最大值和最小值的点位确定方法
CN101862183A (zh) * 2010-06-13 2010-10-20 北京大学 一种基于心磁特征快速确定心磁测量范围的方法
CN101897576A (zh) * 2010-07-18 2010-12-01 北京大学 基于心脏单电流偶极子模型确定心磁测量范围的方法

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11253412A (ja) * 1998-03-17 1999-09-21 Hitachi Ltd 生体磁場計測装置
US6275719B1 (en) * 1998-09-09 2001-08-14 Hitachi, Ltd. Biomagnetic field measurement apparatus
US6522908B1 (en) * 1999-10-06 2003-02-18 Hitachi, Ltd. Biomagnetic field measuring apparatus
JP2001095771A (ja) * 2000-08-09 2001-04-10 Hitachi Ltd 磁場計測装置
CN100337584C (zh) * 2000-08-29 2007-09-19 心脏磁力成像公司 心磁图检查中的局部缺血识别、量化和局部定位的***
CN1471374A (zh) * 2000-08-29 2004-01-28 心磁图检查中的局部缺血识别、量化和局部定位
US20030149354A1 (en) * 2001-08-23 2003-08-07 Bakharev Alexander A. Ischemia identification, quantification and partial localization MCG
CN1556687A (zh) * 2001-09-21 2004-12-22 使用小波变换的心磁图非线性噪声降低
CN1552285A (zh) * 2003-05-26 2004-12-08 中国科学院物理研究所 一种心磁波组图及其分析方法
US20050148844A1 (en) * 2003-09-10 2005-07-07 Kuniomi Ogata Biomagnetic measurement apparatus
JP2008200331A (ja) * 2007-02-21 2008-09-04 Hitachi High-Technologies Corp 生体磁場計測装置
CN101361655A (zh) * 2008-09-12 2009-02-11 北京大学 一种心磁测量区域极值点的测定方法
CN101836860A (zh) * 2010-06-10 2010-09-22 北京大学 一种心脏磁场最大值和最小值的点位确定方法
CN101862183A (zh) * 2010-06-13 2010-10-20 北京大学 一种基于心磁特征快速确定心磁测量范围的方法
CN101897576A (zh) * 2010-07-18 2010-12-01 北京大学 基于心脏单电流偶极子模型确定心磁测量范围的方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
SHIQIN JIANG, ET AL.: "A new ECG obtained from MCG-recordings", 《ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY, 2004. IEMBS "04. 26TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE 》 *
刘芳 等: "心磁图的临床应用及相关比较研究", 《中国心血管杂志》 *
张磊: "基于单电流偶极子的心脏磁场反演计算", 《现代科学仪器》 *
易攀科 等: "心磁和心电信号的混沌时间序列分析", 《同济大学学报(自然科学版)》 *
李英梅 等: "冠心病患者的心磁图分析及应用", 《临床心血管病杂志》 *
李英梅 等: "心磁图在冠心病诊断及治疗后监测中的临床应用", 《心血管病学进展》 *
杜萱 等: "心磁图对老年冠心病的诊断价值", 《中华老年心脑血管病杂志》 *
王会武 等: "用单通道直流高温超导量子干涉器磁强计测量人体心磁", 《科学通报》 *
石明伟 等: "浅谈心脏磁场的检测及心脏源模型建模", 《硅谷》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103810381A (zh) * 2014-01-26 2014-05-21 同济大学 一种基于多通道心磁图的冠心病特征参数提取方法
CN109864733A (zh) * 2019-01-16 2019-06-11 漫迪医疗仪器（上海）有限公司 心脑异常的检测方法、***、介质及设备
CN110074774A (zh) * 2019-04-28 2019-08-02 漫迪医疗仪器（上海）有限公司 基于心磁图的心脏室间隔异常的分析方法、***、介质及终端
CN112932492A (zh) * 2021-01-29 2021-06-11 上海跃磁生物科技有限公司 心磁图仪多个扫描数据同步方法、***、介质及装置
CN112932492B (zh) * 2021-01-29 2022-08-05 上海跃磁生物科技有限公司 心磁图仪多个扫描数据同步方法、***、介质及装置

Also Published As

Publication number Publication date
CN103142223B (zh) 2015-02-18

Similar Documents

Publication Publication Date Title
CN103142211B (zh) 一种基于极值圆的心脏磁信号处理方法
US9125581B2 (en) Continuous modeling for dipole localization from 2D MCG images with unknown depth
CN103584847B (zh) 一种非接触磁感应心率和呼吸率同步检测方法及***
US20020173714A1 (en) Biomagnetic field measuring method and apparatus
WO2019015017A1 (zh) 一种心电动力学数据量化分析方法
JPS63279823A (ja) 生物電気信号の検出処理装置と方法
CN103142223B (zh) 一种基于极值差的心脏磁信号处理方法
CN104794474B (zh) 一种基于心磁图的综合最大电流密度向量特征提取方法
CN112515679A (zh) 无屏蔽式心磁图设备
JP2007268034A (ja) 生体信号計測方法及び装置
CN109157220B (zh) 一种基于多通道动态监测的呼吸指数提取***及工作方法
JP4276107B2 (ja) 心臓磁気計測装置
CN103810381A (zh) 一种基于多通道心磁图的冠心病特征参数提取方法
JP2002028143A (ja) 心房粗動および心房細動の心臓磁界診断装置および心房粗動および心房細動の電気的旋回路の同定方法
CN103654775A (zh) 可视化crbbb/clbbb心脏电活动的等效源重构方法
CN113081261A (zh) 血肿穿刺导航探头装置及脑部磁探测电阻抗成像***
JP2002028144A5 (zh)
CN113288102B (zh) 一种无创监测脑血流的***
CN106236093B (zh) 高精度宽频带分布式医用电阻抗层析成像装置及其应用
Udovychenko et al. Current density distribution maps threshold processing
Wu et al. Noninvasive diagnosis of coronary artery disease using two parameters extracted in an extrema circle of magnetocardiogram
JP4791797B2 (ja) 生体磁場計測装置
JP2002028145A (ja) 心室遅延電位の心臓磁界診断装置および心筋内興奮不均一伝播部位の同定方法
CN105559770B (zh) 内探式生物组织磁分布图测量装置
JPH07369A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150218

Termination date: 20171207