CN108245152B - Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal - Google Patents

Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal Download PDF

Info

Publication number
CN108245152B
CN108245152B CN201810027163.7A CN201810027163A CN108245152B CN 108245152 B CN108245152 B CN 108245152B CN 201810027163 A CN201810027163 A CN 201810027163A CN 108245152 B CN108245152 B CN 108245152B
Authority
CN
China
Prior art keywords
channel
interval
magnetocardiogram
magnetic pole
myocardial ischemia
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.)
Active
Application number
CN201810027163.7A
Other languages
Chinese (zh)
Other versions
CN108245152A (en
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.)
Shanghai Institute of Microsystem and Information Technology of CAS
Original Assignee
Shanghai Institute of Microsystem and Information Technology of CAS
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 Shanghai Institute of Microsystem and Information Technology of CAS filed Critical Shanghai Institute of Microsystem and Information Technology of CAS
Priority to CN201810027163.7A priority Critical patent/CN108245152B/en
Publication of CN108245152A publication Critical patent/CN108245152A/en
Application granted granted Critical
Publication of CN108245152B publication Critical patent/CN108245152B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/242Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
    • A61B5/243Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetocardiographic [MCG] signals

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Cardiology (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

The invention provides a myocardial ischemia characteristic parameter extraction method and system, a storage medium and a terminal, and the method comprises the following steps: acquiring an isomagnetic map and a current density map of a T wave band based on the magnetic field intensity and channel position acquired by each channel of the multi-channel magnetocardiogram instrument; extracting time domain characteristic parameters based on the isomagnetic map and the current density map; performing discrete wavelet transformation on T-wave band magnetocardiogram signals acquired by a multi-channel magnetocardiogram instrument to acquire low-frequency signal components, and performing inverse transformation on the low-frequency signal components to acquire reconstructed signals; extracting frequency domain feature parameters based on the low frequency signal component and the reconstructed signal; and extracting information theory characteristic parameters based on T wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument. The myocardial ischemia characteristic parameter extraction method and system, the storage medium and the terminal can extract time domain magnetocardiogram parameters, frequency domain parameters and information theory parameters, and have strong anti-interference performance and high accuracy.

Description

Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal
Technical Field
The invention relates to the technical field of data processing, in particular to a myocardial ischemia characteristic parameter extraction method and system, a storage medium and a terminal.
Background
The back of the human body life activity is hidden with abundant electromagnetic information. The magnetocardiogram is a new type of heart disease diagnosis method for imaging analysis by detecting the space magnetic field generated by the electrical activity of human heart. Similar to the conventional electrocardiogram, the magnetocardiogram reflects the electrophysiological activity of the heart, and is a functional imaging method. Because the electrocardiogram is completely non-invasive, non-radiative, non-contact, less influenced by body fluid, bones and the like, and can respond to the annular vortex current, the electrocardiogram contains a lot of electrophysiological information which can not be reflected by the traditional electrocardiogram, thereby presenting better sensitivity and early diagnosis capability. Clinical research shows that the magnetocardiogram has good application potential in the aspects of coronary heart disease, myocardial ischemia and the like, thereby having extremely high clinical research and application values.
The sensitivity and reliability of magnetocardiogram prediction for diseases such as myocardial ischemia directly depend on data interpretation. Currently, interpretation of data for magnetocardiogram is based on extracting time-domain characteristic parameters of magnetocardiogram, such as: the angle and position of the current density vector of the heart, the positions and the intensities of the positive and negative magnetic poles on the isomagnetic map, the morphology and the duration of a QRS wave or a T wave and the like to identify normal people and myocardial ischemia patients.
However, the conventional diagnosis scheme of the time domain characteristic parameter has the following problems:
(1) the signal to noise ratio of the signal is high, and the signal is easily influenced by environmental noise and measurement interference;
(2) the interpretation dimension of the magnetocardiogram signals is single, and missed diagnosis and misdiagnosis are easy to happen.
Disclosure of Invention
In view of the above drawbacks of the prior art, an object of the present invention is to provide a method and a system for extracting myocardial ischemia feature parameters, a storage medium, and a terminal, which are capable of extracting time-domain magnetocardiogram parameters, frequency-domain parameters, and information theory parameters, and have strong anti-interference performance and high accuracy.
In order to achieve the above objects and other related objects, the present invention provides a myocardial ischemia feature parameter extraction method based on a multi-channel magnetocardiogram apparatus, comprising the following steps: acquiring an isomagnetic map and a current density map of a T wave band based on the magnetic field intensity and channel position acquired by each channel of the multi-channel magnetocardiogram instrument; extracting time domain characteristic parameters based on the isomagnetic map and the current density map; performing discrete wavelet transformation on T-wave band magnetocardiogram signals acquired by a multi-channel magnetocardiogram instrument to acquire low-frequency signal components, and performing inverse transformation on the low-frequency signal components to acquire reconstructed signals; extracting frequency domain feature parameters based on the low frequency signal component and the reconstructed signal; and extracting information theory characteristic parameters based on T wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument.
In an embodiment of the present invention, the number of the time domain characteristic parameters is 18, the number of the frequency domain characteristic parameters is 3 × N, and the number of the information theory characteristic parameters is N +2, where N is the number of channels of the multi-channel magnetocardiograph.
In an embodiment of the invention, the time-domain characteristic parameters include a maximum current angle of a T-peak, a maximum current angle of a TT-interval, a minimum current angle of a TT-interval, a maximum current angle variation value of a TT-interval, a maximum current position variation standard deviation of a TT-interval, a maximum current position displacement sum of a TT-interval, magnetic field angles of positive and negative poles of a T-peak, a positive and negative pole intensity ratio of a T-peak, a maximum magnetic field angle of a TT-interval, a minimum magnetic field angle of a TT-interval, a magnetic field angle variation value of a TT-interval, a magnetic pole distance displacement variation value of a TT-interval, a negative magnetic pole area variation standard deviation of a TT-interval, a negative magnetic pole position variation sum of a TT-interval, a positive magnetic pole.
In an embodiment of the present invention, each channel of the multi-channel magnetocardiograph includes the following three frequency domain characteristic parameters: frequency domain signal energy, frequency domain signal power, and frequency domain signal to original signal power ratio.
In an embodiment of the present invention, the characteristic parameters of the information theory include a TT interval shannon entropy, a standard deviation of TT interval kini indexes, and a TT interval singular value entropy corresponding to each cardiac magnetic channel.
Correspondingly, the invention provides a myocardial ischemia characteristic parameter extraction system based on a multi-channel magnetocardiogram instrument, which comprises a first extraction module, a second extraction module and a third extraction module;
the first extraction module is used for acquiring an isomagnetic map and a current density map based on the magnetic field intensity and channel position acquired by each channel of the multi-channel magnetocardiogram instrument; extracting time domain characteristic parameters based on the isomagnetic map and the current density map;
the second extraction module is used for performing discrete wavelet transformation on T-wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument to acquire low-frequency signal components, and performing inverse transformation on the low-frequency signal components to acquire reconstructed signals; extracting frequency domain feature parameters based on the low frequency signal component and the reconstructed signal;
the third extraction module is used for extracting information theory characteristic parameters based on T wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument.
The invention provides a storage medium, which stores a computer program, and the program is executed by a processor to realize the myocardial ischemia characteristic parameter extraction method based on the multi-channel magnetocardiograph.
The present invention provides a terminal, including: a processor and a memory;
the memory is used for storing a computer program;
the processor is used for executing the computer program stored in the memory so as to enable the terminal to execute the myocardial ischemia characteristic parameter extraction method based on the multi-channel magnetocardiogram instrument.
Finally, the invention provides a myocardial ischemia characteristic parameter extraction system based on a multi-channel magnetocardiogram instrument, which comprises the terminal and the multi-channel magnetocardiogram instrument.
In an embodiment of the present invention, the multi-channel magnetocardiograph includes 36 channels.
As described above, the myocardial ischemia feature parameter extraction method and system, the storage medium, and the terminal according to the present invention have the following advantageous effects:
(1) time domain magnetocardiogram parameters, frequency domain parameters and information theory parameters can be extracted, the anti-interference performance is strong, and the accuracy is high;
(2) the sensitivity and the reliability of early diagnosis of myocardial ischemia are effectively improved, and the influence of signal noise and environmental interference on the diagnosis result is reduced.
Drawings
FIG. 1 is a flowchart illustrating a myocardial ischemia feature parameter extracting method based on a multi-channel magnetocardiogram apparatus according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of a myocardial ischemia feature parameter extraction system based on a multi-channel magnetocardiogram apparatus according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a terminal according to an embodiment of the invention;
fig. 4 is a schematic structural diagram of a myocardial ischemia feature parameter extraction system based on a multi-channel magnetocardiogram apparatus according to another embodiment of the present invention.
Description of the element reference numerals
21 first extraction module
22 second extraction module
23 third extraction Module
31 processor
32 memory
41 terminal
42 multi-channel magnetocardiogram instrument
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It is to be noted that the features in the following embodiments and examples may be combined with each other without conflict.
It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention, and the components related to the present invention are only shown in the drawings rather than drawn according to the number, shape and size of the components in actual implementation, and the type, quantity and proportion of the components in actual implementation may be changed freely, and the layout of the components may be more complicated.
The myocardial ischemia characteristic parameter extraction method and system based on the multi-channel magnetocardiogram instrument, the storage medium and the terminal can extract time domain magnetocardiogram parameters, frequency domain parameters and information theory parameters, improve the clinical application value of the magnetocardiogram instrument, and have strong anti-interference performance and high accuracy.
As shown in fig. 1, in an embodiment, the method for extracting myocardial ischemia feature parameters based on a multi-channel magnetocardiogram apparatus of the present invention includes the following steps:
s1, acquiring an isomagnetic map and a current density map of a T wave band based on the magnetic field intensity and the channel position acquired by each channel of the multi-channel magnetocardiogram instrument; and extracting time domain characteristic parameters based on the isomagnetic map and the current density map.
Specifically, a magnetic field isoparameter in the T wave band above the thorax is plotted according to the magnetic field intensity recorded in each cardiac magnetic channel and the positional information of the channel. Assuming that the coordinates on the iso-magnetic map of a point above the thorax are (x, y), the magnetic field strength corresponding to this point is F (x,y) the current density of the spot
Figure BDA0001545312640000041
Therefore, a current density map can be obtained by calculating the current densities of all points on the iso-magnetic map.
Then, time-domain feature parameters are extracted according to the isomagnetic map and the current density map. In an embodiment of the present invention, the number of the time domain feature parameters is 18, which is specifically as follows:
1. maximum current angle of T peak
And drawing an isomagnetic graph of the T wave crest moment, recording the maximum current density value on the isomagnetic graph, and calculating the angle by an arctan function.
2. TT interval maximum current angle
The whole T wave band is evenly divided into N equal parts, N equal magnetic graphs are respectively drawn, the maximum current angle corresponding to each equal magnetic graph is recorded, and then the maximum value of the N maximum current angles is selected.
3. Minimum current angle of TT interval
The whole T wave band is evenly divided into N equal parts, N equal magnetic graphs are respectively drawn, the maximum current angle corresponding to each equal magnetic graph is recorded, and then the minimum value of the N maximum current angles is taken.
4. Maximum current angle change value of TT interval
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the maximum current angle corresponding to each equal magnetic map is recorded, and the difference value between the maximum value and the minimum value of the N values is calculated.
5. Standard deviation of position variation value of maximum current in TT interval
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, so that the maximum current displacement on the equal magnetic maps at adjacent moments is recorded, N-1 values are obtained in total, and the standard deviation of the N-1 values is recorded.
6. TT interval maximum current position displacement sum
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the position of the maximum current corresponding to each equal magnetic map is recorded, the moving track of the maximum current at different moments is drawn, and the length of the track is recorded.
7. Magnetic field angle of T peak positive and negative magnetic poles
Wherein, an isomagnetic map of the T wave crest moment is drawn, and the included angle between the positive and negative magnetic pole connecting line and the X axis is calculated.
8. T peak positive and negative magnetic pole strength ratio
Wherein, an isomagnetic map of the T wave crest moment is drawn, and the intensity ratio of the positive and negative magnetic poles is calculated.
9. Maximum value of magnetic field angle in TT interval
The whole T-wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the magnetic field angle (namely the included angle between the positive and negative magnetic pole connecting line and the X axis) corresponding to each equal magnetic map is recorded, and then the maximum value of the N magnetic field angles is selected.
10. Minimum value of magnetic field angle in TT interval
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the magnetic field angle corresponding to each equal magnetic map is recorded, and then the minimum value of the N magnetic field angles is selected.
11. Angle change value of TT interval magnetic field
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the magnetic field angle corresponding to each equal magnetic map is recorded, and then the difference value between the maximum value and the minimum value of the N magnetic field angles is taken.
12. TT interval magnetic pole distance displacement variation value
The whole T-wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the position coordinates of the middle points of the connecting lines of the positive and negative magnetic poles are recorded, and the distance between the farthest two points in the N coordinates is calculated.
13. Standard deviation of variation of TT interval negative magnetic pole area value
The method comprises the steps of uniformly dividing the whole T-wave band by N, respectively drawing N equal magnetic maps, dividing an area with the magnetic field intensity being more than or equal to 80% multiplied by the magnetic field intensity of a negative magnetic pole on the equal magnetic maps into negative magnetic pole areas, calculating the area of the negative magnetic pole areas, and then taking the standard deviation of the area change values of the N negative magnetic poles.
14. TT interval negative magnetic pole position change standard deviation
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the position coordinates of the negative magnetic pole are recorded, and the standard deviation of N coordinate positions is calculated.
15. TT interval negative pole position change
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the displacement value of the negative magnetic pole at the adjacent moment is calculated, and the sum of the N-1 displacement values is calculated.
16. Standard deviation of variation of positive pole area value in TT interval
The method comprises the steps of uniformly dividing the whole T-wave band by N, respectively drawing N equal magnetic maps, dividing an area with the magnetic field intensity being more than or equal to 80% multiplied by the magnetic field intensity of a positive magnetic pole on the equal magnetic maps into positive magnetic pole areas, calculating the area of the area, and then taking the standard difference of the area change values of the N positive magnetic poles.
17. TT interval positive pole position variation standard deviation
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the position coordinates of the positive magnetic pole are recorded, and the standard deviation of N coordinate positions is calculated.
18. TT interval positive pole position change
The whole T wave band is evenly divided into N equal parts, N equal magnetic maps are respectively drawn, the displacement value of the positive magnetic pole at the adjacent moment is calculated, and the sum of the N-1 displacement values is calculated.
Step S2, discrete wavelet transform is carried out on T wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument to acquire low-frequency signal components, inverse transform is carried out on the low-frequency signal components to acquire reconstructed signals, and frequency domain characteristic parameters are extracted based on the low-frequency signal components and the reconstructed signals.
Specifically, db4 discrete wavelet transform is performed on T-wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument, a low-frequency signal component after the fourth-order transform is a4 by taking a fourth-order transform as an example, the a4 components are subjected to inverse transform, and a4 signal component is reconstructed to correspond to a 0-7.8 Hz frequency band of the original signal. Because the one-dimensional wavelet transform processes one channel magnetocardiogram signal at a time, 3 × N frequency domain parameters are obtained for the N-channel magnetocardiogram system.
In an embodiment of the present invention, each channel of the multi-channel magnetocardiograph includes the following three frequency domain characteristic parameters:
1. and the frequency domain signal energy represents the energy ratio of the low-frequency signal component to the original signal.
2. The frequency domain signal power represents the power of the reconstructed signal.
3. The frequency domain signal to original signal power ratio represents the power ratio of the reconstructed signal to the original signal.
And step S3, extracting information theory characteristic parameters based on the T wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument.
In an embodiment of the present invention, the number of the information theory characteristic parameters is N +2, where N is the number of channels of the multi-channel magnetocardiograph. Specifically, the information theory characteristic parameters include:
1. the TT interval Shannon entropy represents the corresponding Shannon entropy of the T wave band on each heart magnetic channel, and has N values.
2. The standard deviation of the TT interval Bonney index, which represents the standard deviation of the Bonney index of the T-wave band magnetocardiogram signal, has 1 value.
3. The TT interval singular value entropy represents the singular value index of the T wave band magnetocardiogram signal and has 1 value in total.
Therefore, extraction of myocardial ischemia characteristic parameters based on the multi-channel magnetocardiogram instrument is completed. In the clinic, the obtained myocardial ischemia characteristic parameters are generally applied by:
A. manual screening method
According to the normal range of the time domain characteristic parameters and the information theory characteristic parameters, the morbidity of myocardial ischemia can be directly predicted through manual screening.
B. Machine learning method
And the time domain characteristic parameters, the frequency domain characteristic parameters and the information theory characteristic parameters are used as the input of the machine learning method, so that the normal person and the myocardial ischemia patient are distinguished according to the output of the machine learning method.
As shown in fig. 2, in an embodiment, the myocardial ischemia feature parameter extracting system based on a multi-channel magnetocardiograph of the present invention includes a first extracting module 21, a second extracting module 22 and a third extracting module 23.
The first extraction module 21 is configured to obtain an isomagnetic map and a current density map of a T-wave band based on the magnetic field intensity and channel position obtained by each channel of the multi-channel magnetocardiogram apparatus; extracting time domain characteristic parameters based on the isomagnetic map and the current density map;
the second extraction module 22 is configured to perform discrete wavelet transform on T-wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument to acquire low-frequency signal components, and perform inverse transform on the low-frequency signal components to acquire reconstructed signals; extracting frequency domain feature parameters based on the low frequency signal component and the reconstructed signal;
the third extraction module 23 is configured to extract information theory characteristic parameters based on T-wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram apparatus.
The specific working principle of the first extraction module 21, the second extraction module 22 and the third extraction module 23 corresponds to each step in the myocardial ischemia characteristic parameter extraction method based on the multi-channel magnetocardiogram apparatus, and therefore, the detailed description thereof is omitted here.
It should be noted that the division of the modules of the above system is only a logical division, and the actual implementation may be wholly or partially integrated into one physical entity, or may be physically separated. And these modules can be realized in the form of software called by processing element; or may be implemented entirely in hardware; and part of the modules can be realized in the form of calling software by the processing element, and part of the modules can be realized in the form of hardware. For example, the x module may be a processing element that is set up separately, or may be implemented by being integrated in a chip of the apparatus, or may be stored in a memory of the apparatus in the form of program code, and the function of the x module may be called and executed by a processing element of the apparatus. Other modules are implemented similarly. In addition, all or part of the modules can be integrated together or can be independently realized. The processing element described herein may be an integrated circuit having signal processing capabilities. In implementation, each step of the above method or each module above may be implemented by an integrated logic circuit of hardware in a processor element or an instruction in the form of software.
For example, the above modules may be one or more integrated circuits configured to implement the above methods, such as: one or more Application Specific Integrated Circuits (ASICs), or one or more microprocessors (DSPs), or one or more Field Programmable Gate Arrays (FPGAs), etc. For another example, when one of the above modules is implemented in the form of a processing element scheduler code, the processing element may be a general-purpose processor, such as a Central Processing Unit (CPU) or other processor capable of calling program code. For another example, these modules may be integrated together and implemented in the form of a system-on-a-chip (SOC).
The storage medium of the invention stores a computer program which is executed by a processor to realize the myocardial ischemia characteristic parameter extraction method based on the multi-channel magnetocardiograph. Preferably, the storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
As shown in fig. 3, in one embodiment, the terminal of the present invention includes a processor 31 and a memory 32.
The memory 32 is used for storing computer programs.
Preferably, the memory 32 comprises: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
The processor 31 is connected to the memory 32 and is configured to execute the computer program stored in the memory 32, so that the terminal executes the method for extracting myocardial ischemia feature parameters based on the multi-channel magnetocardiogram apparatus.
Preferably, the processor 31 may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; the integrated circuit may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, or discrete hardware components.
As shown in fig. 4, in an embodiment, the myocardial ischemia feature parameter extracting system based on a multi-channel magnetocardiograph of the present invention includes the above-mentioned terminal 41 and the multi-channel magnetocardiograph 42. The multi-channel magnetocardiogram instrument 42 collects data for the terminal to extract myocardial ischemia characteristic parameters.
In one embodiment of the present invention, the multi-channel magnetocardiograph 42 includes 36 channels.
In conclusion, the myocardial ischemia characteristic parameter extraction method, the myocardial ischemia characteristic parameter extraction system, the storage medium and the terminal can extract time domain magnetocardiogram parameters, frequency domain parameters and information theory parameters, and have strong anti-interference performance and high accuracy; the sensitivity and the reliability of early diagnosis of myocardial ischemia are effectively improved, and the influence of signal noise and environmental interference on the diagnosis result is reduced. Therefore, the invention effectively overcomes various defects in the prior art and has high industrial utilization value.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.

Claims (8)

1. A myocardial ischemia characteristic parameter extraction method based on a multi-channel magnetocardiogram instrument is characterized by comprising the following steps:
acquiring an isomagnetic map and a current density map of a T wave band based on the magnetic field intensity and channel position acquired by each channel of the multi-channel magnetocardiogram instrument; extracting time domain characteristic parameters based on the isomagnetic map and the current density map;
performing discrete wavelet transformation on T-wave band magnetocardiogram signals acquired by a multi-channel magnetocardiogram instrument to acquire low-frequency signal components, and performing inverse transformation on the low-frequency signal components to acquire reconstructed signals; extracting frequency domain feature parameters based on the low frequency signal component and the reconstructed signal;
extracting information theory characteristic parameters based on T wave band magnetocardiogram signals acquired by a multi-channel magnetocardiogram instrument;
the characteristic parameters of the information theory comprise TT interval Shannon entropy, TT interval kini index and TT interval singular value entropy corresponding to each heart magnetic channel;
the time domain characteristic parameters comprise a maximum current angle of a T peak, a maximum current angle of a TT interval, a minimum current angle of a TT interval, a maximum current angle change value of a TT interval, a maximum current position change standard deviation of the TT interval, a maximum current position displacement sum of the TT interval, a magnetic field angle of a positive magnetic pole and a negative magnetic pole of the T peak, a positive magnetic pole and negative magnetic pole intensity ratio of the T peak, a maximum magnetic field angle of the TT interval, a minimum magnetic field angle of the TT interval, a magnetic field angle change value of the TT interval, a magnetic pole distance displacement change value of the TT interval, a negative magnetic pole area change standard deviation of the TT interval, a negative magnetic pole position change sum of the TT interval, a positive magnetic pole area change standard deviation of the TT.
2. The myocardial ischemia feature parameter extraction method based on the multi-channel magnetocardiograph according to claim 1, wherein the number of the time domain feature parameters is 18, the number of the frequency domain feature parameters is 3 × N, the number of the information theory feature parameters is N +2, where N is the number of channels of the multi-channel magnetocardiograph.
3. The myocardial ischemia feature parameter extraction method based on multi-channel magnetocardiograph according to claim 1, characterized in that, for each channel of the multi-channel magnetocardiograph, the following three frequency domain feature parameters are included: frequency domain signal energy, frequency domain signal power, and frequency domain signal to original signal power ratio.
4. A myocardial ischemia characteristic parameter extraction system based on a multi-channel magnetocardiogram instrument is characterized by comprising a first extraction module, a second extraction module and a third extraction module;
the first extraction module is used for acquiring an isomagnetic map and a current density map based on the magnetic field intensity and channel position acquired by each channel of the multi-channel magnetocardiogram instrument; extracting time domain characteristic parameters based on the T wave band isomagnetic map and the current density map;
the second extraction module is used for performing discrete wavelet transformation on T-wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument to acquire low-frequency signal components, and performing inverse transformation on the low-frequency signal components to acquire reconstructed signals; extracting frequency domain feature parameters based on the low frequency signal component and the reconstructed signal;
the third extraction module is used for extracting information theory characteristic parameters based on T wave band magnetocardiogram signals acquired by the multi-channel magnetocardiogram instrument;
the time domain characteristic parameters comprise a maximum current angle of a T peak, a maximum current angle of a TT interval, a minimum current angle of a TT interval, a maximum current angle change value of a TT interval, a maximum current position change standard deviation of the TT interval, a maximum current position displacement sum of the TT interval, a magnetic field angle of a positive magnetic pole and a negative magnetic pole of the T peak, a positive magnetic pole and negative magnetic pole intensity ratio of the T peak, a maximum magnetic field angle of the TT interval, a minimum magnetic field angle of the TT interval, a magnetic field angle change value of the TT interval, a magnetic pole distance displacement change value of the TT interval, a negative magnetic pole area change standard deviation of the TT interval, a negative magnetic pole position change sum of the TT interval, a positive magnetic pole area change standard deviation of the TT.
5. A storage medium having stored thereon a computer program, wherein the program, when executed by a processor, implements the method for myocardial ischemia feature parameter extraction based on multi-channel magnetocardiograph of any one of claims 1 to 3.
6. A terminal, comprising: a processor and a memory;
the memory is used for storing a computer program;
the processor is used for executing the computer program stored in the memory to enable the terminal to execute the myocardial ischemia characteristic parameter extraction method based on the multi-channel magnetocardiograph, which is described in any one of claims 1 to 3.
7. A myocardial ischemia feature parameter extraction system based on a multi-channel magnetocardiogram instrument, characterized by comprising the terminal of claim 6 and the multi-channel magnetocardiogram instrument.
8. The myocardial ischemia feature parameter extraction system based on multichannel magnetocardiograph according to claim 7, wherein the multichannel magnetocardiograph comprises 36 channels.
CN201810027163.7A 2018-01-11 2018-01-11 Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal Active CN108245152B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810027163.7A CN108245152B (en) 2018-01-11 2018-01-11 Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810027163.7A CN108245152B (en) 2018-01-11 2018-01-11 Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal

Publications (2)

Publication Number Publication Date
CN108245152A CN108245152A (en) 2018-07-06
CN108245152B true CN108245152B (en) 2021-06-18

Family

ID=62726246

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810027163.7A Active CN108245152B (en) 2018-01-11 2018-01-11 Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal

Country Status (1)

Country Link
CN (1) CN108245152B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110074774B (en) * 2019-04-28 2022-04-12 漫迪医疗仪器(上海)有限公司 Analysis method, system, medium and terminal for heart ventricular septal abnormality based on magnetocardiogram
CN113317793B (en) * 2021-06-11 2023-02-17 宁波大学 Magnetocardiogram high-frequency signal analysis method, storage medium, and electronic device
CN115985505B (en) * 2023-01-19 2023-12-12 北京未磁科技有限公司 Multidimensional fusion myocardial ischemia auxiliary diagnosis model and construction method thereof
CN116189902B (en) * 2023-01-19 2024-01-02 北京未磁科技有限公司 Myocardial ischemia prediction model based on magnetocardiogram video data and construction method thereof
CN117084684B (en) * 2023-10-19 2024-02-02 山东大学齐鲁医院 Characteristic parameter extraction method and system based on electrocardio current density map extension field
CN117100276B (en) * 2023-10-23 2024-01-12 山东大学齐鲁医院 Cardiac function detection system, computer storage medium and terminal
CN117113064B (en) * 2023-10-23 2024-02-02 杭州诺驰生命科学有限公司 Multi-dimensional magnetocardiogram characteristic parameter extraction method and system
CN117137492B (en) * 2023-11-01 2024-02-09 山东大学齐鲁医院 Coronary artery blood flow abnormality detection system, storage medium, and terminal
CN117297613B (en) * 2023-11-14 2024-04-05 杭州诺驰生命科学有限公司 Disease feature extraction method based on magnetocardiogram and cardiovascular disease prediction system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
UA100450C2 (en) * 2011-05-04 2012-12-25 Владимир Николаевич Сосницкий METHOD OF MEASUREMENT OF QT, QRS, ST-T-CARDIOCYCLE INTERVALS AND DEVICES FOR ITS PERFORMANCE
CN102512158B (en) * 2011-12-31 2013-09-04 济南汇医融工科技有限公司 Electrocardiosignal quality evaluation method and device based on high-dimensional fuzzy recognition
US8954140B2 (en) * 2012-09-27 2015-02-10 Samsung Electronics Co., Ltd. Method and system for determining QRS complexes in electrocardiogram signals
CN104586387B (en) * 2015-01-19 2017-02-01 秦皇岛市惠斯安普医学***股份有限公司 Method for extracting and fusing time, frequency and space domain multi-parameter electroencephalogram characters
CN107361760B (en) * 2017-07-12 2020-06-12 中国科学院上海微***与信息技术研究所 Magnetocardiogram instrument diagnosis system and magnetocardiogram instrument using the same

Also Published As

Publication number Publication date
CN108245152A (en) 2018-07-06

Similar Documents

Publication Publication Date Title
CN108245152B (en) Myocardial ischemia characteristic parameter extraction method and system, storage medium and terminal
Liu et al. Real-time multilead convolutional neural network for myocardial infarction detection
CN108577825B (en) Myocardial ischemia lesion position positioning method and system, storage medium and terminal
CN110074774B (en) Analysis method, system, medium and terminal for heart ventricular septal abnormality based on magnetocardiogram
Wu et al. Extracting deep features from short ECG signals for early atrial fibrillation detection
Xiang et al. ECG-based heartbeat classification using two-level convolutional neural network and RR interval difference
CN101919704B (en) Heart sound signal positioning and segmenting method
CN106137185A (en) A kind of epileptic chracter wave detecting method based on structure of transvers plate small echo
Zhang et al. Detection of atrial fibrillation from variable-duration ECG signal based on time-adaptive densely network and feature enhancement strategy
CN112932498A (en) T wave morphology classification system with strong generalization capability based on deep learning
Liu et al. A novel P-QRS-T wave localization method in ECG signals based on hybrid neural networks
Liu et al. An efficient neural network-based method for patient-specific information involved arrhythmia detection
Spicher et al. ECG delineation using a piecewise Gaussian derivative model with parameters estimated from scale-dependent algebraic expressions
CN115486849A (en) Electrocardiosignal quality evaluation method, device and equipment
Tariq et al. Automatic multimodal heart disease classification using phonocardiogram signal
Kishore et al. An innovative machine learning approach for classifying ECG signals in healthcare devices
CN114287950A (en) Heart disease computer-assisted classification method based on continuous coherence and Fourier transform
Kim et al. Lightweight Convolutional Neural Network for Real-Time Arrhythmia Classification on Low-Power Wearable Electrocardiograph
Farag A Matched Filter-Based Convolutional Neural Network (Cnn) for Inter-Patient Ecg Classification and Arrhythmia Detection at the Edge
Khan et al. Performance analysis of modified zero crossing counts method for heart arrhythmias detection and implementation in HDL
CN117100276B (en) Cardiac function detection system, computer storage medium and terminal
Brucal et al. Development of a 12-Lead ECG Signal Processing Algorithm Using NI LabVIEW® and NI ELVIS®
CN117137492B (en) Coronary artery blood flow abnormality detection system, storage medium, and terminal
CN117084684B (en) Characteristic parameter extraction method and system based on electrocardio current density map extension field
CN112818830B (en) Magnetocardiogram signal noise reduction method, system, medium and device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant