CN103142223A - Heart magnetic signal processing method based on extreme value difference - Google Patents
Heart magnetic signal processing method based on extreme value difference Download PDFInfo
- 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
Links
Images
Abstract
The invention relates to a heart magnetic signal processing method based on extreme value difference, and the method is characterized by comprising the following steps of 1) detecting magnetic induction intensity of 6*6 array detection points of a heart magnetic field which is vertical to a chest plane through a superconducting quantum interference device (SQUID) sensor, and synchronously measuring an electrocardiogram; 2) intercepting data of the magnetic induction intensity of the heart magnetic field on an ST-T section according to the electrocardiogram; 3) conducting three times of spline interpolation for the magnetic induction intensity on the 6*6 array detection points at every moment win the ST-T section, and acquiring a high-resolution isogamme chart; 4) acquiring a maximal value Bzmax and a minimal value Bzmin according to the isogamme chart; and 5) calculating electrophysiological parameters JS according to a formula, and generating a JS curve. Compared with the prior art, the electrophysiological parameters capable of assisting the clinical diagnosis of heart diseases can be rapidly calculated only by utilizing the SQUID sensor to detect the maximal value and the minimal value of the heart magnetic field data in the ST-T section, and the electrical activity function of the heart can be analyzed.
Description
Technical field
The present invention relates to a kind of magnetic signal processing method, especially relate to a kind of based on the poor heart magnetic signal processing method of extreme value.
Background technology
The seventies, the people such as D.Cohen use superconducting quantum interference device (SQUID) to measure human heart magnetic field at laboratory first.1976, the developer of heart magnetic instrument proposed theory and the method with magnetocardiogram and pseudo-electric current density figure (or arrow plot, also be referred to as the Hosaka-Cohen conversion) diagnosis of myocardial ischemia in succession.The nineties, the people such as 9 passage heart magnetic instrument .2005 Germany J.W.Park that U.S. CMI company has released with heart diseases such as single magnetic dipole algorithm diagnosis of myocardial ischemia have proposed the method for the measurement data prediction coronary heart disease (CAD) of magnetic instrument diligently.2006, German W.Haberkorn etc. proposed to have the pseudo-Current density imaging method of electric physiological significance on the basis of pseudo-electric current density figure.In the same year, K.Tolstrup etc. have also proposed a kind of method that quick magnetic imaging detects myocardial ischemia.2007, Taiwan and Cooperation in Korea research institution proposed a kind of two dimension of motive magnetic T ripple signal and have propagated imaging and area ratio method, and were used for the disease such as diagnosis of myocardial ischemia.In the same year, Japan has proposed a kind of method .A.Gapelyuk with magnetocardiogram JT section integrated value screening coronary heart disease etc. and has also proposed a kind of method of surveying CAD with the heart magnetic chart.The sensitivity of these diagnostic methods and specificity are in the 60%-80% left and right.In recent years, along with the development of cardiac magnetic field detection technique, the improving constantly of magnetocardiogram measurement equipment performance, the method for non-invasive diagnosis heart disease has also had certain progress.2008, the method for three kinds of non-invasive diagnosis coronary heart disease was compared in the quantitative analyses such as P.V.Leeuwen.2010, Kwon etc. and A.Gapelyuk philosophy magnetocardiogram Classification and Identification with the combination of KL entropy and two kinds of methods of rest parameter, brought up to sensitivity and the specificity of these diagnostic methods more than 80%.People wish the research by correlation theory, constantly explore new suitable clinical practice, have high sensitivity and specific, a signal processing method that computational speed is fast.In a kind of heart magnetic signal New Method for Processing poor based on extreme value, need by calculating electric-physiology parameter JS as intermediate object program, auxiliary judgment coronary heart disease.
Summary of the invention
But the purpose of this invention is to provide a kind of new heart magnetic signal processing method based on the poor quick obtaining electric-physiology parameter of extreme value.
Purpose of the present invention can be achieved through the following technical solutions:
A kind of based on the poor heart magnetic signal processing method of extreme value, comprise the following steps:
1) the SQUID sensor detects cardiac magnetic field perpendicular to the magnetic induction of 6 * 6 array detection points on plane, thoracic cavity, and the synchro measure electrocardiogram;
2) intercept the magnetic induction of cardiac magnetic field in the data of ST-T section according to electrocardiogram;
3) magnetic induction on each moment 6 * 6 array detection point in the ST-T section is carried out cubic spline interpolation and process, obtain the high-resolution magnetic field line figure that waits;
4) according to waiting magnetic field line figure to obtain maximum and minimum, in these magnetic field line charts, the line of the maximum point of magnetic induction and minimum point is D;
5) according to formula
Calculate electric-physiology parameter JS, and generate the JS curve.
In described 6 * 6 array detection points, the spacing of adjacent test point is 4cm, and whole magnetic field detection plane sizes is 20 * 20cm.
Be Electrocardiographic T crest place in the time of the cut-off of described ST-T section, the initial time of ST-T section is described T crest 1/3 amplitude place forward.
Described step 3) obtain the figure such as magnetic field line such as grade of 81 * 81 in after cubic spline interpolation.
Compared with prior art, the interior cardiac magnetic field data of ST-T section that the present invention only need to utilize the SQUID sensor to detect, according to the maximum and the minimum that wait in magnetic field line figure, just can calculate fast the electric-physiology parameter JS of accessory heart medical diagnosis on disease, and the function of analysis of cardiac electrical activity.
Description of drawings
Fig. 1 is flow chart of the present invention;
Fig. 2 is the schematic diagram of the present invention's 6 * 6 array detection points;
Fig. 3 is the curve chart of the ST-T section of electrocardiosignal in electrocardiogram;
Fig. 4 is the ST-T section curve chart with the mcg-signals of electrocardiogram synchro measure;
Fig. 5 is the figure such as magnetic field line such as grade after cubic spline interpolation;
Fig. 6 is magnetic field profile line schematic diagram;
Fig. 7 is the curve chart of electric-physiology parameter JS of normal person's ST-T section;
Fig. 8 is the curve chart of electric-physiology parameter JS of the ST-T section of Coronary Heart Disease Patients.
The specific embodiment
The present invention is described in detail below in conjunction with the drawings and specific embodiments.
Embodiment
A kind of based on the poor heart magnetic signal processing method of extreme value, the flow process of the method comprises the following steps as shown in Figure 1:
Step S1: synchro measure magnetocardiogram and electrocardiogram, when measuring mcg-signals, the person under inspection lies on the back, multichannel SQUID sensor array detects the magnetic induction on 6 * 6 array detection points as shown in Figure 2 of surface, person under inspection thoracic cavity, because the spacing of each adjacent test point is 4cm, whole detection plane size is 20 * 20cm.By multichannel SQUID sensor array can real time record perpendicular to the magnetic induction density B z on the Z direction of torso model plane.
Step S2: according to the data of electrocardiogram intercepting mcg-signals in the ST-T section.With the cut-off time t of the Electrocardiographic T crest value moment as corresponding heart magnetic data ST-T section
max, with the initial time t as heart magnetic data ST-T section of T crest value 1/3 amplitude place forward
min, as shown in Figure 3 and Figure 4.
Step S3: the magnetic induction to each moment 6 * 6 array detection point in the ST-T section carries out the cubic spline interpolation processing, to obtain the higher figure such as magnetic field line such as grade of resolution, as shown in Figure 5.In the present embodiment, interpolation point is 81 * 81.
Step S4: according to waiting magnetic field line figure to obtain maximum and minimum.By the maximum point of magnetic induction and the line D of minimum point, can obtain the magnetic field contour line schematic diagram relevant with minimum with magnetic field maximum.As shown in Figure 6.Namely by the magnetic field line figure such as plane cutting of line D in Fig. 5, can obtain B as shown in Figure 6
Z σ ECurve.When the hypothesis cardiac magnetic field is when being produced by single equivalent dipole source, according to the Biot-Savart law, this magnetic field can be divided into two parts.A part of B wherein
∞(r) produced in infinitely great uniform dielectric by the equivalent dipole source; Another part B
σ EProduced by volume electric current σ E.At the extreme point place, by the measured value B of Z direction
Zmax, B
ZminCan calculate B
Z σ EB
Z σ ESize characterized the magnetic induction that volume electric current in corresponding moment heart produces.B
Z σ ECurve representation the time-varying characteristics of the magnetic induction that produces of this volume electric current.
Step S5: according to the maximum and the minimum that wait magnetic field line figure, and formula
Calculate electric-physiology parameter JS, and generate the JS curve.
But obtain electric-physiology parameter JS by the present invention's fast processing, be used for the auxiliary judgment heart body.The curve chart of the electric-physiology parameter JS of normal person's ST-T section as shown in Figure 7.The curve chart of the electric-physiology parameter JS of the ST-T section of Coronary Heart Disease Patients as shown in Figure 8.The JS parameter of normal person and patients with coronary heart disease has obvious difference.
Claims (4)
1. one kind based on the poor heart magnetic signal processing method of extreme value, it is characterized in that, comprises the following steps:
1) the SQUID sensor detects cardiac magnetic field perpendicular to the magnetic induction of 6 * 6 array detection points on plane, thoracic cavity, and the synchro measure electrocardiogram;
2) intercept the magnetic induction of cardiac magnetic field in the data of ST-T section according to electrocardiogram;
3) magnetic induction on each moment 6 * 6 array detection point in the ST-T section is carried out cubic spline interpolation and process, obtain the high-resolution magnetic field line figure that waits;
4) according to waiting magnetic field line figure to obtain maximum B
ZmaxWith minimum B
Zmin
5) according to formula
Calculate electric-physiology parameter JS, and generate the JS curve.
2. according to claim 1 a kind ofly it is characterized in that based on the poor heart magnetic signal processing method of extreme value, in described 6 * 6 array detection points, the spacing of adjacent test point is 4cm.
3. according to claim 1 a kind ofly it is characterized in that based on the poor heart magnetic signal processing method of extreme value, the cut-off time of described ST-T section is Electrocardiographic T crest place, and the initial time of ST-T section is described T crest 1/3 amplitude place forward.
4. according to claim 1 a kind of based on the poor heart magnetic signal processing method of extreme value, it is characterized in that described step 3) in obtain after cubic spline interpolation 81 * 81 wait magnetic field line figure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110404161.3A CN103142223B (en) | 2011-12-07 | 2011-12-07 | Heart magnetic signal processing method based on extreme value difference |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110404161.3A CN103142223B (en) | 2011-12-07 | 2011-12-07 | Heart magnetic signal processing method based on extreme value difference |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103142223A true CN103142223A (en) | 2013-06-12 |
| CN103142223B CN103142223B (en) | 2015-02-18 |
Family
ID=48540775
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110404161.3A Expired - Fee Related CN103142223B (en) | 2011-12-07 | 2011-12-07 | Heart magnetic signal processing method based on extreme value difference |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103142223B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103810381A (en) * | 2014-01-26 | 2014-05-21 | 同济大学 | Coronary heart disease characteristic parameter extracting method based on multichannel magnetocardiogram |
| CN109864733A (en) * | 2019-01-16 | 2019-06-11 | 漫迪医疗仪器（上海）有限公司 | Detection method, system, medium and the equipment of heart and brain exception |
| CN110074774A (en) * | 2019-04-28 | 2019-08-02 | 漫迪医疗仪器（上海）有限公司 | Analysis method, system, medium and the terminal of cardiac septum exception based on magnetocardiogram |
| CN112932492A (en) * | 2021-01-29 | 2021-06-11 | 上海跃磁生物科技有限公司 | Method, system, medium, and apparatus for synchronizing a plurality of scan data of a magnetocardiograph |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11253412A (en) * | 1998-03-17 | 1999-09-21 | Hitachi Ltd | Boimagnetic field measuring device |
| JP2001095771A (en) * | 2000-08-09 | 2001-04-10 | Hitachi Ltd | Magnetic field measuring apparatus |
| 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 (en) * | 2000-08-29 | 2004-01-28 | Ischemia identification, quantification and partial localization in MCG | |
| CN1552285A (en) * | 2003-05-26 | 2004-12-08 | 中国科学院物理研究所 | Magnetocardiogram and its analyzing method |
| CN1556687A (en) * | 2001-09-21 | 2004-12-22 | Nonlinear noise reduction for magnetocardiograms using wavelet transforms | |
| US20050148844A1 (en) * | 2003-09-10 | 2005-07-07 | Kuniomi Ogata | Biomagnetic measurement apparatus |
| JP2008200331A (en) * | 2007-02-21 | 2008-09-04 | Hitachi High-Technologies Corp | Biological magnetic field measuring apparatus |
| CN101361655A (en) * | 2008-09-12 | 2009-02-11 | 北京大学 | Measuring method of extreme point in heart magnetic measurement zone |
| CN101836860A (en) * | 2010-06-10 | 2010-09-22 | 北京大学 | Method for determining point positions of maximum value and minimum value of cardiac magnetic field |
| CN101862183A (en) * | 2010-06-13 | 2010-10-20 | 北京大学 | Method for quickly determining magnetocardiogram measurement range based on magnetocardiogram characteristics |
| CN101897576A (en) * | 2010-07-18 | 2010-12-01 | 北京大学 | Method for determining magnetocardiogram measurement range based on cardiac single current dipole model |
-
2011
- 2011-12-07 CN CN201110404161.3A patent/CN103142223B/en not_active Expired - Fee Related
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11253412A (en) * | 1998-03-17 | 1999-09-21 | Hitachi Ltd | Boimagnetic field measuring device |
| 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 (en) * | 2000-08-09 | 2001-04-10 | Hitachi Ltd | Magnetic field measuring apparatus |
| CN100337584C (en) * | 2000-08-29 | 2007-09-19 | 心脏磁力成像公司 | Ischemia identification, quantification and partial localization in MCG |
| CN1471374A (en) * | 2000-08-29 | 2004-01-28 | Ischemia identification, quantification and partial localization in MCG | |
| US20030149354A1 (en) * | 2001-08-23 | 2003-08-07 | Bakharev Alexander A. | Ischemia identification, quantification and partial localization MCG |
| CN1556687A (en) * | 2001-09-21 | 2004-12-22 | Nonlinear noise reduction for magnetocardiograms using wavelet transforms | |
| CN1552285A (en) * | 2003-05-26 | 2004-12-08 | 中国科学院物理研究所 | Magnetocardiogram and its analyzing method |
| US20050148844A1 (en) * | 2003-09-10 | 2005-07-07 | Kuniomi Ogata | Biomagnetic measurement apparatus |
| JP2008200331A (en) * | 2007-02-21 | 2008-09-04 | Hitachi High-Technologies Corp | Biological magnetic field measuring apparatus |
| CN101361655A (en) * | 2008-09-12 | 2009-02-11 | 北京大学 | Measuring method of extreme point in heart magnetic measurement zone |
| CN101836860A (en) * | 2010-06-10 | 2010-09-22 | 北京大学 | Method for determining point positions of maximum value and minimum value of cardiac magnetic field |
| CN101862183A (en) * | 2010-06-13 | 2010-10-20 | 北京大学 | Method for quickly determining magnetocardiogram measurement range based on magnetocardiogram characteristics |
| CN101897576A (en) * | 2010-07-18 | 2010-12-01 | 北京大学 | Method for determining magnetocardiogram measurement range based on cardiac single current dipole model |
Non-Patent Citations (9)
| 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)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103810381A (en) * | 2014-01-26 | 2014-05-21 | 同济大学 | Coronary heart disease characteristic parameter extracting method based on multichannel magnetocardiogram |
| CN109864733A (en) * | 2019-01-16 | 2019-06-11 | 漫迪医疗仪器（上海）有限公司 | Detection method, system, medium and the equipment of heart and brain exception |
| CN110074774A (en) * | 2019-04-28 | 2019-08-02 | 漫迪医疗仪器（上海）有限公司 | Analysis method, system, medium and the terminal of cardiac septum exception based on magnetocardiogram |
| CN112932492A (en) * | 2021-01-29 | 2021-06-11 | 上海跃磁生物科技有限公司 | Method, system, medium, and apparatus for synchronizing a plurality of scan data of a magnetocardiograph |
| CN112932492B (en) * | 2021-01-29 | 2022-08-05 | 上海跃磁生物科技有限公司 | Method, system, medium, and apparatus for synchronizing a plurality of scan data of a magnetocardiograph |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103142223B (en) | 2015-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103142211B (en) | Heart magnetic signal processing method based on extreme value circle | |
| US9125581B2 (en) | Continuous modeling for dipole localization from 2D MCG images with unknown depth | |
| US20010009975A1 (en) | Biomagnetic field measuring method and apparatus | |
| JPS63279823A (en) | Apparatus and method for detecting processing bio-electric signal | |
| WO2019015017A1 (en) | Quantitative analysis method for electrocardio dynamics data | |
| CN103142223B (en) | Heart magnetic signal processing method based on extreme value difference | |
| CN112515679A (en) | Unshielded magnetocardiogram device | |
| CN109157220B (en) | Respiratory index extraction system based on multi-channel dynamic monitoring and working method | |
| JP4276107B2 (en) | Cardiac magnetic measuring device | |
| CN103810381A (en) | Coronary heart disease characteristic parameter extracting method based on multichannel magnetocardiogram | |
| JP2002028143A (en) | Cardiac magnetic field diagnostic apparatus of atrial flutter and atrial fibrillation, and method for identifying electrical re-entry circuit of atrial flutter and atrial fibrillation | |
| JP2007268034A (en) | Method and device for measuring biological signal | |
| CN103654775A (en) | Equivalent source reconstruction method of visual CRBBB/CLBBB heart electrical activity | |
| JP2002028144A5 (en) | ||
| CN113288102B (en) | System for monitoring cerebral blood flow without wound | |
| CN107374610B (en) | Magnetocardiogram generation method and system | |
| CN106236093B (en) | High-precision broadband distributed medical electrical impedance tomography device and application thereof | |
| Udovychenko et al. | Current density distribution maps threshold processing | |
| CN114287943A (en) | Three-dimensional imaging method based on wearable magnetocardiogram three-dimensional measuring device | |
| CN113081261A (en) | Hematoma puncture navigation probe device and brain magnetic detection electrical impedance imaging system | |
| CN114468996A (en) | Method for analyzing breast signs based on orderliness, multimodality and symmetry deficiency | |
| Wu et al. | Noninvasive diagnosis of coronary artery disease using two parameters extracted in an extrema circle of magnetocardiogram | |
| JP4791797B2 (en) | Biomagnetic field measurement device | |
| CN105559770B (en) | Nei Tan formulas biological tissue magnetic distribution map measuring device | |
| JPH07369A (en) | High-speed imaging method of internal impedance distribution |
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 |