CN103142223B - 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
- CN103142223B CN103142223B CN201110404161.3A CN201110404161A CN103142223B CN 103142223 B CN103142223 B CN 103142223B CN 201110404161 A CN201110404161 A CN 201110404161A CN 103142223 B CN103142223 B CN 103142223B
- Authority
- CN
- China
- Prior art keywords
- heart
- magnetic field
- 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.)
- Expired - Fee Related
Links
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 heart magnetic signal processing method based on extreme value difference.
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 the theory and means by magnetocardiogram and pseudo-electric current density figure (or arrow plot, be also referred to as Hosaka-Cohen conversion) diagnosis of myocardial ischemia in succession.The nineties, CMI company of the U.S. is proposed the method proposing measurement data prediction coronary heart disease (CAD) of magnetic instrument diligently with 9 passage heart magnetic instrument .2005 German such as J.W.Park of the heart diseases such as single magnetic dipole algorithm diagnosis of myocardial ischemia.2006, German W.Haberkorn etc., on the basis of pseudo-electric current density figure, proposed the pseudo-Current density imaging method with electro physiology meaning.In the same year, K.Tolstrup etc. it is also proposed a kind of method that quick magnetic imaging detects myocardial ischemia.2007, the two dimension that Taiwan and Cooperation in Korea research institution propose a kind of motive magnetic T ripple signal propagated imaging and area ratio method, and was used for the diseases such as diagnosis of myocardial ischemia.In the same year, the method .A.Gapelyuk etc. that Japan proposes a kind of magnetocardiogram JT section integrated value screening coronary heart disease it is also proposed a kind of method that heart magnetic chart detects CAD.The sensitivity of these diagnostic methods and specificity are at about 60%-80%.In recent years, along with the development of cardiac magnetic field detection technique, improving constantly of magnetocardiogram measurement equipment performance, the method for non-invasive diagnosis heart disease there has also been certain progress.2008, the quantitative analyses such as P.V.Leeuwen compared the method for three kinds of non-invasive diagnosis coronary heart disease.2010, the magnetic chart Classification and Identification diligently respectively such as Kwon etc. and A.Gapelyuk, with the combination of KL entropy and rest parameter two kinds of methods, brought up to more than 80% by the sensitivity of these diagnostic methods and specificity.People wish the research by correlation theory, constantly to explore new applicable clinical practice, have high sensitivity and specific, the signal processing method that computational speed is fast.A kind of based in the heart magnetic signal New Method for Processing of extreme value difference, need by calculating electric-physiology parameter JS as intermediate object program, auxiliary judgment coronary heart disease.
Summary of the invention
The object of this invention is to provide a kind of newly based on extreme value difference can the heart magnetic signal processing method of quick obtaining electric-physiology parameter.
Object of the present invention can be achieved through the following technical solutions:
Based on a heart magnetic signal processing method for extreme value difference, comprise the following steps:
1) SQUID sensor detects the magnetic induction of cardiac magnetic field perpendicular to 6 × 6 array detection points of thoracic cavity plane, and synchro measure electrocardiogram;
2) data of magnetic induction in ST-T section of cardiac magnetic field are intercepted according to electrocardiogram;
3) cubic spline interpolation process is carried out to the magnetic induction on each moment 6 × 6 array detection point in ST-T section, obtain high-resolution figure such as magnetic field line such as grade;
4) according to waiting magnetic field line figure to obtain maximum and minimum, in these magnetic field line charts, the maximum point of magnetic induction and the line of minimum point are D;
5) according to formula
calculate electric-physiology parameter JS, and generate JS curve.
In 6 × 6 described 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 when 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) in obtain after cubic spline interpolation 81 × 81 wait magnetic field line figure.
Compared with prior art, cardiac magnetic field data in the ST-T section that the present invention only needs to utilize SQUID sensor to detect, according to the maximum waited in magnetic field line figure and minimum, the electric-physiology parameter JS of accessory heart medical diagnosis on disease just can be calculated fast, and the function of analysis of cardiac electrical activity.
Accompanying drawing explanation
Fig. 1 is flow chart of the present invention;
Fig. 2 is the schematic diagram of the present invention 6 × 6 array detection point;
Fig. 3 is the curve chart of the ST-T section of electrocardiogram center telecommunications number;
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 field profile line schematic diagram;
Fig. 7 is the curve chart of the electric-physiology parameter JS of the ST-T section of normal person;
Fig. 8 is the curve chart of the electric-physiology parameter JS of the ST-T section of Coronary Heart Disease Patients.
Detailed description of the invention
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.
Embodiment
Based on a heart magnetic signal processing method for extreme value difference, the flow process of the method as shown in Figure 1, comprises the following steps:
Step S1: synchro measure magnetocardiogram and electrocardiogram, when measuring mcg-signals, supine, 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, spacing due to each adjacent test point is 4cm, and whole detection plane size is 20 × 20cm.By multichannel SQUID sensor array can real time record perpendicular to the magnetic induction density B z in torso model plane Z-direction.
Step S2: intercept the data of mcg-signals in ST-T section according to electrocardiogram.Cut-off time t using the Electrocardiographic T crest value moment as corresponding heart magnetic data ST-T section
max, using T crest value 1/3 amplitude place forward as the initial time t of heart magnetic data ST-T section
min, as shown in Figure 3 and Figure 4.
Step S3: cubic spline interpolation process is carried out to the magnetic induction of each moment 6 × 6 array detection point in ST-T section, 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: obtain maximum and minimum according to waiting magnetic field line figure.By the maximum point of magnetic induction and the line D of minimum point, the Magnetic field contours schematic diagram relevant with minimum with magnetic field maximum can be obtained.As shown in Figure 6.Namely pass through the magnetic field line figure such as the plane cutting of line D in Fig. 5, B as shown in Figure 6 can be obtained
z σ Ecurve.When supposing that cardiac magnetic field is by single equivalent dipole source generation, according to Biot-Savart law, this magnetic field can be divided into two parts.Wherein a part of B
∞r () is produced in infinitely great uniform dielectric by equivalent dipole source; Another part B
σ Eproduced by volumetric energy σ E.At extreme point place, by the measured value B of Z-direction
zmax, B
zminb can be calculated
z σ E.B
z σ Esize characterize the magnetic induction that in corresponding moment heart, volumetric energy produces.B
z σ Ecurve illustrates the time-varying characteristics of the magnetic induction that this volumetric energy produces.
Step S5: according to the maximum and the minimum that wait magnetic field line figure, and formula
calculate electric-physiology parameter JS, and generate JS curve.
Fast processing electric-physiology parameter JS can be obtained, for auxiliary judgment heart body by the present invention.The curve chart of the electric-physiology parameter JS of the ST-T section of normal person 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 (3)
1., based on a heart magnetic signal processing method for extreme value difference, it is characterized in that, comprise the following steps:
1) SQUID sensor detects the magnetic induction of cardiac magnetic field perpendicular to 6 × 6 array detection points of thoracic cavity plane, and synchro measure electrocardiogram;
2) data of magnetic induction in ST-T section of cardiac magnetic field are intercepted according to electrocardiogram;
3) cubic spline interpolation process is carried out to the magnetic induction on each moment 6 × 6 array detection point in ST-T section, obtain high-resolution figure such as magnetic field line such as grade;
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 JS curve;
Described step 3) in obtain after cubic spline interpolation 81 × 81 wait magnetic field line figure.
2. a kind of heart magnetic signal processing method based on extreme value difference according to claim 1, is characterized in that, in 6 × 6 described array detection points, the spacing of adjacent test point is 4cm.
3. a kind of heart magnetic signal processing method based on extreme value difference according to claim 1, is characterized in that, 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.
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 CN103142223A (en) | 2013-06-12 |
| CN103142223B true 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) |
Families Citing this family (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 |
| CN110074774B (en) * | 2019-04-28 | 2022-04-12 | 漫迪医疗仪器（上海）有限公司 | Analysis method, system, medium and terminal for heart ventricular septal abnormality based on magnetocardiogram |
| CN112932492B (en) * | 2021-01-29 | 2022-08-05 | 上海跃磁生物科技有限公司 | Method, system, medium, and apparatus for synchronizing a plurality of scan data of a magnetocardiograph |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 | |
| 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 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3546686B2 (en) * | 1998-03-17 | 2004-07-28 | 株式会社日立製作所 | Biomagnetic field measurement device |
| JP3424664B2 (en) * | 2000-08-09 | 2003-07-07 | 株式会社日立製作所 | Magnetic field measurement device |
| US20030149354A1 (en) * | 2001-08-23 | 2003-08-07 | Bakharev Alexander A. | Ischemia identification, quantification and partial localization MCG |
| JP4027867B2 (en) * | 2003-09-10 | 2007-12-26 | 株式会社日立ハイテクノロジーズ | Biomagnetic field measurement device |
| JP5085155B2 (en) * | 2007-02-21 | 2012-11-28 | 株式会社日立ハイテクノロジーズ | Biomagnetic field measurement device |
-
2011
- 2011-12-07 CN CN201110404161.3A patent/CN103142223B/en not_active Expired - Fee Related
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN1471374A (en) * | 2000-08-29 | 2004-01-28 | Ischemia identification, quantification and partial localization in MCG | |
| CN100337584C (en) * | 2000-08-29 | 2007-09-19 | 心脏磁力成像公司 | Ischemia identification, quantification and partial localization in 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 |
| 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 》.2004,1945-1948. * |
| 冠心病患者的心磁图分析及应用;李英梅 等;《临床心血管病杂志》;20040425;第20卷(第4期);229-231 * |
| 基于单电流偶极子的心脏磁场反演计算;张磊;《现代科学仪器》;20050825(第4期);31-34 * |
| 心磁和心电信号的混沌时间序列分析;易攀科 等;《同济大学学报(自然科学版)》;20080115;第36卷(第1期);108-112 * |
| 心磁图在冠心病诊断及治疗后监测中的临床应用;李英梅 等;《心血管病学进展》;20040525;第25卷(第3期);193-197 * |
| 心磁图对老年冠心病的诊断价值;杜萱 等;《中华老年心脑血管病杂志》;20060215;第8卷(第2期);90-92 * |
| 心磁图的临床应用及相关比较研究;刘芳 等;《中国心血管杂志》;20051225;第10卷(第6期);472-475 * |
| 浅谈心脏磁场的检测及心脏源模型建模;石明伟 等;《硅谷》;20110208(第3期);44-45 * |
| 用单通道直流高温超导量子干涉器磁强计测量人体心磁;王会武 等;《科学通报》;20020915;第47卷(第17期);1296-1298 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103142223A (en) | 2013-06-12 |
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 | |
| CN103584847B (en) | Non-contact magnetic induction heart rate and respiration rate synchronous detection method and system | |
| WO2019015017A1 (en) | Quantitative analysis method for electrocardio dynamics data | |
| JPS63279823A (en) | Apparatus and method for detecting processing bio-electric signal | |
| CN103142223B (en) | Heart magnetic signal processing method based on extreme value difference | |
| CN110148466A (en) | A kind of heart impact signal atrial fibrillation computer aided diagnosing method based on transfer learning | |
| JP2007268034A (en) | Method and device for measuring biological signal | |
| JP4276107B2 (en) | Cardiac magnetic measuring device | |
| CN103810381A (en) | Coronary heart disease characteristic parameter extracting method based on multichannel magnetocardiogram | |
| CN109157220B (en) | Respiratory index extraction system based on multi-channel dynamic monitoring and working method | |
| Calder et al. | A simulated anatomically accurate investigation into the effects of biodiversity on electrogastrography | |
| CN107374610B (en) | Magnetocardiogram generation method and system | |
| Udovychenko et al. | Current density distribution maps threshold processing | |
| CN106236093B (en) | High-precision broadband distributed medical electrical impedance tomography device and application thereof | |
| CN113288102B (en) | System for monitoring cerebral blood flow without wound | |
| Apama et al. | Interpretation of ResNet50 model for MI related cardiac events using Explainable Grad-CAM approach | |
| CN114287943A (en) | Three-dimensional imaging method based on wearable magnetocardiogram three-dimensional measuring device | |
| CN114468996A (en) | Method for analyzing breast signs based on orderliness, multimodality and symmetry deficiency | |
| JP4791797B2 (en) | Biomagnetic field measurement device | |
| Dossel et al. | Optimization of electrode positions for multichannel electrocardiography with respect to electrical imaging of the heart | |
| Wu et al. | Noninvasive diagnosis of coronary artery disease using two parameters extracted in an extrema circle of magnetocardiogram | |
| CN105559770B (en) | Nei Tan formulas biological tissue magnetic distribution map measuring device | |
| TWI802888B (en) | Cardiovascular Function Assessment System | |
| JP3518493B2 (en) | Calculation method of isometric diagram of biomagnetic field |
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 |