CN113288169A - Method, device and equipment for identifying waveform of electrocardiographic waveform signal - Google Patents

Method, device and equipment for identifying waveform of electrocardiographic waveform signal Download PDF

Info

Publication number
CN113288169A
CN113288169A CN202110579117.XA CN202110579117A CN113288169A CN 113288169 A CN113288169 A CN 113288169A CN 202110579117 A CN202110579117 A CN 202110579117A CN 113288169 A CN113288169 A CN 113288169A
Authority
CN
China
Prior art keywords
point
waveform signal
target
value
derivative
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
CN202110579117.XA
Other languages
Chinese (zh)
Other versions
CN113288169B (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.)
Neusoft Corp
Original Assignee
Neusoft Corp
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 Neusoft Corp filed Critical Neusoft Corp
Priority to CN202110579117.XA priority Critical patent/CN113288169B/en
Publication of CN113288169A publication Critical patent/CN113288169A/en
Application granted granted Critical
Publication of CN113288169B publication Critical patent/CN113288169B/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/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The embodiment of the application discloses a method, a device and equipment for identifying the waveform of an electrocardiographic waveform signal, wherein the method comprises the steps of firstly identifying an R point, a Q point and an S point of a complete heart beat and time values respectively corresponding to all the points from the electrocardiographic waveform signal by acquiring the electrocardiographic waveform signal comprising a plurality of complete heart beats; obtaining an average value of first time intervals between adjacent R points; further determining a starting point time value and an ending point time value of any complete heartbeat; and obtaining a target T-wave waveform signal and a target P-wave waveform signal based on the determined starting point time value, Q point time value, S point time value and end point time value, and further determining a T point and a P point. The processing mode of the electrocardiographic waveform signal is simple, and the identification result of the electrocardiographic waveform signal can be obtained quickly. And the identification of five characteristic points in the electrocardiographic waveform signal can be realized, and a relatively complete identification result of the electrocardiographic waveform signal is obtained.

Description

Method, device and equipment for identifying waveform of electrocardiographic waveform signal
Technical Field
The application relates to the field of data processing, in particular to a method, a device and equipment for identifying waveforms of an electrocardiographic waveform signal.
Background
The electrocardiographic waveform signal is a waveform signal which is generated by monitoring the beating process of the heart of a patient and represents the beating condition of the heart of the patient. Based on the electrocardiographic waveform signals, an analysis of the cardiovascular health of the patient may be performed.
The waveform of the electrocardiogram waveform signal is identified, and the characteristic points in the electrocardiogram waveform signal can be determined. At present, the waveform identification result of the electrocardiographic waveform signal is not accurate enough, all waveforms in the electrocardiographic waveform signal cannot be identified, and the obtained identification result is not complete enough. Therefore, how to realize complete and accurate waveform identification of the electrocardiographic waveform signal is an urgent problem to be solved.
Disclosure of Invention
In view of this, embodiments of the present application provide a method, an apparatus, and a device for identifying waveforms of an electrocardiographic waveform signal, which can completely identify different waveforms in the electrocardiographic waveform signal, so as to obtain a more accurate waveform identification result of the electrocardiographic waveform signal.
In order to solve the above problem, the technical solution provided by the embodiment of the present application is as follows:
a method of waveform identification of an electrocardiographic waveform signal, the method comprising:
acquiring an electrocardiographic waveform signal comprising a plurality of complete heart beats;
identifying an R point, a Q point and an S point of the complete heartbeat from the electrocardiographic waveform signal, and acquiring a time value of the R point, a time value of the Q point and a time value of the S point;
calculating the average value of the first time interval between the adjacent R points according to the time value of the R point;
determining a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat according to an average value of first time intervals between the adjacent R points, wherein the target complete heartbeat is any one of the complete heartbeats;
intercepting an electrocardiographic waveform signal between the starting point time value of the target complete heartbeat and the time value of the Q point of the target complete heartbeat from the electrocardiographic waveform signal as a target P-wave waveform signal;
intercepting an electrocardiographic waveform signal between the time value of the S point of the complete target heart beat and the time value of the end point of the complete target heart beat from the electrocardiographic waveform signal as a target T-wave waveform signal;
determining a P point from the target P wave waveform signal by deriving the target P wave waveform signal;
and determining a T point from the target T wave waveform signal by performing derivation on the target T wave waveform signal.
In a possible implementation manner, the determining a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat according to an average value of first time intervals between the adjacent R points includes:
multiplying the average value of the first time intervals among the R points by a first proportional value to obtain a first time interval, and multiplying the average value of the first time intervals among the R points by a second proportional value to obtain a second time interval;
and subtracting the first time interval from the time value of the R point of the complete heartbeat of the target to obtain a starting point time value of the complete heartbeat of the target, and adding the second time interval to the time value of the R point of the complete heartbeat of the target to obtain an ending point time value of the complete heartbeat of the target.
In one possible implementation, the determining a P point from the target P-wave waveform signal by deriving the target P-wave waveform signal includes:
performing derivation on the target P-wave waveform signal to obtain a derivative value of each sampling point of the target P-wave waveform signal;
searching a first derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming the first derivative value into a first derivative value sequence;
searching a second derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming a second derivative value sequence by the second derivative values;
acquiring the time value of a sampling point corresponding to each first derivative value in the first derivative value sequence and the time value of a sampling point corresponding to each second derivative value in the second derivative value sequence;
intercepting an electrocardiographic waveform signal between sampling points corresponding to first derivative values and sampling points corresponding to second derivative values adjacent to time values from the target P-wave waveform signal as a candidate P-point waveform signal;
and determining the sampling point with the maximum amplitude absolute value in the candidate P point waveform signal as the P point.
In one possible implementation, before determining a sampling point with the largest absolute value of amplitude in the candidate P-point waveform signal as a P-point, the method further includes:
calculating a second time interval between a starting point of a target candidate P-point waveform signal and an ending point of the target candidate P-point waveform signal, wherein the target candidate P-point waveform signal is any one of the candidate P-point waveform signals;
calculating a first product of a derivative value corresponding to a starting point of the target candidate P point waveform signal and a derivative value corresponding to an ending point of the target candidate P point waveform signal;
and if the second time interval is smaller than a first threshold value or the first product is larger than zero, rejecting the target candidate P point waveform signal.
In one possible implementation, the determining a T point from the target T-wave waveform signal by deriving the target T-wave waveform signal includes:
performing derivation on the target T-wave waveform signal to obtain a derivative value of each sampling point of the target T-wave waveform signal;
searching a third derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a third derivative value sequence by the third derivative values;
searching a fourth derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a fourth derivative value sequence by the fourth derivative values;
acquiring the time value of a sampling point corresponding to each third derivative value in the third derivative value sequence and the time value of a sampling point corresponding to each fourth derivative value in the fourth derivative value sequence;
intercepting an electrocardiographic waveform signal between sampling points corresponding to third derivative values and sampling points corresponding to fourth derivative values adjacent to time values from the target T-wave waveform signal as a candidate T-point waveform signal;
and determining the sampling point with the maximum amplitude absolute value in the candidate T-point waveform signal as the T point.
In one possible implementation, before determining a sampling point with the largest absolute value of amplitude in the candidate T-point waveform signal as a T-point, the method further includes:
calculating a third time interval between a starting point of a target candidate T-point waveform signal and an ending point of the target candidate T-point waveform signal, wherein the target candidate T-point waveform signal is any one of the candidate T-point waveform signals;
calculating a second product of a derivative value corresponding to a starting point of the target candidate T-point waveform signal and a derivative value corresponding to an ending point of the target candidate T-point waveform signal;
and if the third time interval is smaller than a second threshold value or the second product is larger than zero, rejecting the target candidate T-point waveform signal.
In one possible implementation, before identifying the R, Q, and S points of the complete heartbeat from the electrocardiographic waveform signal, the method further comprises:
and performing wavelet transformation and nonlinear filtering processing on the electrocardiographic waveform signal to obtain the electrocardiographic waveform signal again.
An apparatus for waveform identification of an electrocardiographic waveform signal, the apparatus comprising:
an acquisition unit for acquiring an electrocardiographic waveform signal including a plurality of complete heart beats;
the identification unit is used for identifying an R point, a Q point and an S point of the complete heartbeat from the electrocardiographic waveform signal and acquiring a time value of the R point, a time value of the Q point and a time value of the S point;
the first calculating unit is used for calculating the average value of the first time interval between the adjacent R points according to the time value of the R point;
a first determining unit, configured to determine a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat according to an average value of first time intervals between the adjacent R points, where the target complete heartbeat is any one of the complete heartbeats;
a second determining unit, configured to intercept, from the electrocardiographic waveform signal, an electrocardiographic waveform signal between a starting point time value of the target complete heartbeat and a time value of a Q point of the target complete heartbeat, as a target P-wave waveform signal;
a third determining unit, configured to intercept, from the electrocardiographic waveform signal, an electrocardiographic waveform signal between a time value of an S point of the target complete heartbeat and a time value of an end point of the target complete heartbeat, as a target T-wave waveform signal;
a P-point determining unit configured to determine a P-point from the target P-wave waveform signal by deriving the target P-wave waveform signal;
a T-point determining unit configured to determine a T-point from the target T-wave waveform signal by deriving the target T-wave waveform signal.
An apparatus for waveform identification of an electrocardiographic waveform signal, comprising: the device comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the waveform identification method of the electrocardiogram waveform signal.
A computer-readable storage medium having stored therein instructions which, when run on a terminal device, cause the terminal device to execute the method of waveform identification of an electrocardiographic waveform signal as described above.
Therefore, the embodiment of the application has the following beneficial effects:
the embodiment of the application provides a method, a device and equipment for identifying the waveform of an electrocardiographic waveform signal, wherein the method comprises the steps of firstly identifying an R point, a Q point and an S point of a complete heart beat and time values respectively corresponding to all the points from the electrocardiographic waveform signal by acquiring the electrocardiographic waveform signal comprising a plurality of complete heart beats; based on the time value of the R point, the time interval corresponding to a complete heart beat, namely the average value of the first time intervals between the adjacent R points can be obtained; further determining the starting point time value and the ending point time value of any complete heartbeat; based on the determined starting point time value and the Q point time value, a target P wave waveform signal can be obtained; based on the determined S point time value and the determined end point time value, a target T wave waveform signal can be obtained; and respectively carrying out derivation calculation on the target P wave waveform signal and the target T wave waveform signal to determine a P point and a T point. The processing mode of the electrocardiographic waveform signal is simple, and the identification result of the electrocardiographic waveform signal can be obtained quickly. And based on the R point, the Q point and the S point obtained by identification, the P point and the T point can be further determined, the identification of five characteristic points in the electrocardiogram waveform signal is realized, and a more complete and accurate identification result of the electrocardiogram waveform signal is obtained. Through the complete and accurate electrocardio waveform signals obtained by identification, doctors can conveniently diagnose diseases according to the electrocardio waveform signals.
Drawings
Fig. 1 is a schematic diagram of an example scenario provided in an embodiment of the present application;
fig. 2 is a flowchart of a waveform identification method for an electrocardiographic waveform signal according to an embodiment of the present application;
FIG. 3 is a schematic diagram of an electrocardiographic waveform signal including two complete heart beats according to an embodiment of the present application;
fig. 4 is a schematic structural diagram of a waveform identification apparatus for an electrocardiographic waveform signal according to an embodiment of the present application.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, embodiments accompanying the drawings are described in detail below.
In order to facilitate understanding of the technical solutions provided in the present application, the following description will be made on the background related to the present application.
After researching the traditional method for identifying the electrocardiographic waveform signal, the electrocardiographic waveform signal comprises five characteristic points which are respectively an R point, a Q point, an S point, a P point and a T point. By analyzing the waveform between the five characteristic points, the analysis of the health condition of the patient based on the electrocardiographic waveform signal can be realized. At present, most of analyses aiming at the electrocardio waveform signals are the analyses aiming at QRS waveforms, and the analyses of P waves and T waves are lacked, so that the waveform identification of the electrocardio waveform signals is incomplete. In addition, in analysis of an electrocardiographic waveform signal, a method such as differential filtering, an artificial neural network, a differential method, and template matching is generally used. The accuracy of the analysis method is low, the analysis process is complex, and the requirement on the computing power of waveform identification hardware equipment is high.
Based on the above, the embodiment of the application provides a method, a device and equipment for identifying the waveform of an electrocardiographic waveform signal. In order to facilitate understanding of the waveform identification method of an electrocardiographic waveform signal provided in the embodiments of the present application, the following description is made with reference to a scene example shown in fig. 1. Referring to fig. 1, the figure is a schematic diagram of a framework of an exemplary application scenario provided in an embodiment of the present application. The method can be applied to the terminal device 101.
In practical application, the terminal device 101 may obtain an electrocardiographic waveform signal including a plurality of complete heart beats, obtain, from the obtained electrocardiographic waveform signal, R, Q, and S points of the complete heart beats, and time values corresponding to the R, Q, and S points, respectively, and based on the time values corresponding to the R points, may calculate an average value of a first time interval between adjacent R points. The average value of the first time intervals between the R points can represent the time interval of the complete heartbeat, and the starting point time value and the ending point time value of the target complete heartbeat can be determined by using the first time intervals between the R points. Based on the starting point time value and the Q point time value, a target P wave waveform signal can be determined; based on the S-point time value and the end point time value, a target T-wave waveform signal may be determined. The point P is a certain point of the target P-wave waveform signal with the derivative of 0, and the point T is a certain point of the target T-wave waveform signal with the derivative of 0. By derivation of the target P-wave waveform signal and the target T-wave waveform signal, a P point and a T point can be calculated, and the electrocardio waveform signal can be completely and accurately identified.
Those skilled in the art will appreciate that the block diagram shown in fig. 1 is only one example in which embodiments of the present application may be implemented. The scope of applicability of the embodiments of the present application is not limited in any way by this framework.
Based on the above description, the waveform identification method of an electrocardiographic waveform signal provided by the present application will be described in detail below with reference to the accompanying drawings.
Referring to fig. 2, which is a flowchart illustrating a method for identifying a waveform of an electrocardiographic waveform signal according to an embodiment of the present application, as shown in fig. 2, the method may include S201-S208:
s201: an electrocardiographic waveform signal comprising a plurality of complete heart beats is acquired.
The electrocardiographic waveform signal is a waveform signal containing information related to the heart beat obtained by monitoring the heart beat.
Since the heart beats regularly, the resulting electrocardiographic waveform signal is periodic. The method comprises the steps of obtaining an electrocardiographic waveform signal comprising a plurality of complete heart beats, determining a complete heart beat based on the periodicity of the electrocardiographic waveform signal, and further realizing the identification of the electrocardiographic waveform signal of the complete heart beat.
Referring to fig. 3, a schematic diagram of an electrocardiographic waveform signal including two complete heart beats is provided according to an embodiment of the present application. The starting point of an electrocardiographic waveform signal of a complete heart beat is SP, and the ending point is EP. The normal electrocardiographic waveform signal of a complete heart beat is composed of a T wave, a QRS complex and a P wave. Wherein the P-wave represents an atrial depolarization process of the heart; the PR interval is a P-wave and a PR segment, and represents the beginning of atrial depolarization to ventricular depolarization; the QRS complex represents the depolarization process of the ventricles; ST-T represents the entire ventricular repolarization process; the QT interval represents the entire ventricular activity process. By determining the point P, the point Q, the point R, the point S and the point T, each stage in the electrocardiographic waveform signal of a complete heart beat can be further determined, so that the electrocardiographic waveform signal can be analyzed.
S202: and identifying the R point, the Q point and the S point of the complete heart beat from the electrocardiographic waveform signal, and acquiring the time value of the R point, the time value of the Q point and the time value of the S point.
The points R, Q, and S in the complete heartbeat in the electrocardiographic waveform signal are points where the turn is more pronounced in the electrocardiographic waveform signal. And identifying the R point, the Q point and the S point of the complete heart beat in the electrocardio waveform signal to obtain the R point, the Q point and the S point in the complete heart beat. For an electrocardiographic waveform signal, the abscissa represents the time-of-day value. Based on the identified R point, Q point, and S point, the time value of the corresponding R point, the time value of the corresponding Q point, and the time value of the corresponding S point can be obtained.
It should be noted that the embodiments of the present application do not limit the specific method for identifying the R point, the Q point, and the S point of the complete heartbeat in the electrocardiographic waveform signal. Specifically, the electrocardio waveform signal can be detected by a method such as a differential threshold method.
S203: and calculating the average value of the first time interval between the adjacent R points according to the time value of the R points.
The R point is more obvious in the electrocardiographic waveform signal, and is usually the point with the largest absolute value of amplitude in the electrocardiographic waveform signal of each complete heart beat. The time interval between adjacent R points may be used to determine the time interval of an entire heartbeat of the electrocardiographic waveform signal based on the periodicity of the electrocardiographic waveform signal.
And calculating to obtain the average value of the first time interval between the adjacent R points according to the determined time value of the R point. The accurate first time interval between the adjacent R points can be obtained by calculating the average value, and errors caused by the first time interval between partial adjacent R points can be avoided.
S204: and determining a starting point time value of the complete heartbeat of the target and an end point time value of the complete heartbeat of the target according to the average value of the first time interval between the adjacent R points, wherein the complete heartbeat of the target is any one of the complete heartbeats.
The average value of the first time interval between the adjacent R points can be used as the time interval of the electrocardiographic waveform signal of the complete heart beat, and any complete heart beat, namely the starting point time value and the ending point time value of the target complete heart beat, can be determined by using the average value of the first time interval between the adjacent R points. Based on the starting point time value and the ending point time of the complete heartbeat of the target, the range of the complete heartbeat of the target can be determined.
In a possible implementation manner, an embodiment of the present application provides a specific implementation manner for determining a start point time value of a target complete heartbeat and an end point time value of the target complete heartbeat according to an average value of first time intervals between adjacent R points, which is specifically referred to below.
S205: and intercepting the electrocardiographic waveform signal between the starting point time value of the complete heartbeat of the target and the time value of the Q point of the complete heartbeat of the target from the electrocardiographic waveform signal as a target P wave waveform signal.
The P-wave is the first band of a complete heartbeat. And determining the electrocardiographic waveform signal between the starting point of the target complete heartbeat and the Q point of the target complete heartbeat as a target P-wave waveform signal. And intercepting the electrocardiographic waveform signal between the starting point time value of the complete heartbeat of the target and the time value of the Q point of the complete heartbeat of the target from the electrocardiographic waveform signal as a target P wave waveform signal.
S206: and intercepting the electrocardiographic waveform signal between the time value of the S point of the complete heartbeat of the target and the time value of the end point of the complete heartbeat of the target from the electrocardiographic waveform signal as a target T-wave waveform signal.
The T wave is the last wave band of a complete heartbeat. And determining the electrocardiographic waveform signal between the point S of the complete heartbeat of the target and the end point of the complete heartbeat of the target as a target T-wave waveform signal. And intercepting the electrocardiographic waveform signal between the time value of the S point of the complete target heart beat and the time value of the end point of the complete target heart beat from the electrocardiographic waveform signal as a target T-wave waveform signal.
S207: and determining a P point from the target P wave waveform signal by differentiating the target P wave waveform signal.
The P-point is a point where the derivative in the target P-wave waveform signal is 0. And (4) carrying out derivation on the target P-wave waveform signal, and determining a P point in the target P-wave waveform signal by using a derivative obtained by derivation.
In a possible implementation manner, the embodiment of the present application provides a specific implementation manner of determining a P point from a target P-wave waveform signal by performing derivation on the target P-wave waveform signal, which is described in detail below.
S208: and determining a T point from the target T wave waveform signal by carrying out derivation on the target T wave waveform signal.
Similarly, the T point is also a point in the target T-wave waveform signal at which the derivative is 0. And (4) carrying out derivation on the target T wave waveform signal, and determining a T point in the target T wave waveform signal by using a derivative obtained by derivation.
In a possible implementation manner, the present application provides a specific implementation manner of determining a T point from a target T-wave waveform signal by performing derivation on the target T-wave waveform signal, which is described in detail below.
Based on the related contents of the above S201-S208, it can be known that by determining the R point, the Q point, and the S point that are more obvious in the complete heartbeat first, the accurate start point time value and the end point time value of the target complete heartbeat can be determined based on the accurate R point and further according to the average value of the first time interval between the adjacent R points. Based on the starting point time value, the Q point time value, the end point time value and the S point time value, the target complete heartbeat can be divided to obtain a P wave and a T wave, and complete analysis of all waveforms in the complete heartbeat is achieved. And according to the obtained P wave and T wave, the P point and the T point are determined by calculating the derivative, the P point and the T point can be accurately determined by simple calculation, and the requirement on the computing capability of hardware is reduced. Based on the results of the electrocardiographic waveform signal recognition for a complete heart beat, a doctor is helped to diagnose the symptoms of the patient.
The positions of the R points in the complete heartbeat are relatively fixed, and the starting point time value and the ending point time value of the complete heartbeat can be determined based on the determined average value of the first time intervals between the adjacent R points and the time value of the R points.
Correspondingly, an embodiment of the present application provides a specific implementation manner for determining a start point time value of a target complete heartbeat and an end point time value of the target complete heartbeat according to an average value of first time intervals between adjacent R points, which specifically includes the following two steps:
a1: and multiplying the average value of the first time intervals among the R points by a first proportional value to obtain a first time interval, and multiplying the average value of the first time intervals among the R points by a second proportional value to obtain a second time interval.
Referring to fig. 3, based on the position of the R point in the target complete heartbeat, the complete heartbeat may be divided into two stages by using the time value of the R point, which are the time interval from the start point time value to the time value of the R point and the time interval from the time value of the R point to the end point time value. Correspondingly, the first proportional value may be a proportional value of a time interval from the start point time value to the R point time value to a time interval of the target complete heartbeat, and the second proportional value may be a proportional value of a time interval from the R point time value to the end point time value to a time interval of the target complete heartbeat.
The first proportion value and the second proportion value can be preset according to the relevant data of the complete heartbeat. For example, the first proportional value may be set to 0.4 and the second proportional value may be set to 0.6.
Multiplying the average value of the first time interval between the R points by a first proportional value to obtain a first time interval; and multiplying the average value of the first time interval between the R points by a second proportional value to obtain a second time interval. The first time interval and the second time interval are two time intervals obtained by dividing the complete heartbeat of the target based on the R point.
A2: and subtracting the first time interval from the time value of the R point of the complete heartbeat of the target to obtain a starting point time value of the complete heartbeat of the target, and adding the second time interval to the time value of the R point of the complete heartbeat of the target to obtain an ending point time value of the complete heartbeat of the target.
And calculating to obtain the starting point time value of the complete heartbeat of the target by using the time value of the R point of the complete heartbeat of the target and the first time interval. And subtracting the first time interval from the time value of the R point of the complete heartbeat of the target to obtain the time value of the starting point of the complete heartbeat of the target.
And calculating to obtain the time value of the end point of the complete heartbeat of the target by using the time value of the R point of the complete heartbeat of the target and the second time interval. And adding the time value of the R point of the complete heartbeat of the target to the second time interval to obtain the time value of the end point of the complete heartbeat of the target.
In the embodiment of the application, the starting point time value and the ending point time value of the complete heartbeat can be determined according to the time value of the R points, the average value of the first time interval between the R points, the first proportion value and the second proportion value. Because the R point is obvious in the electrocardiogram waveform signal, the determined time value of the R point is accurate, and further the accurate initial point time value and the accurate end point time value of the complete heartbeat can be obtained based on the time value of the R point, so that the complete heartbeat can be accurately divided, and the complete heartbeat can be accurately analyzed subsequently.
In order to accurately determine the P point, a derivative may be performed on the target P-wave waveform signal, and a candidate P-point waveform signal is determined first, and then the P point is determined based on the candidate P-point waveform signal.
The embodiment of the present application provides a specific implementation manner for determining a point P from a target P-wave waveform signal by deriving the target P-wave waveform signal, which specifically includes the following six steps:
b1: and carrying out derivation on the target P wave waveform signal to obtain a derivative value of each sampling point of the target P wave waveform signal.
In order to determine the waveform change condition of each sampling point in the target P-wave waveform signal, derivation is performed on the target P-wave waveform signal to obtain a derivative value of each sampling point.
B2: and searching a first derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming the first derivative value into a first derivative value sequence.
The derivative value of the sampling point which is larger than the left and right adjacent sampling points shows that the slope of the sampling point is larger than that of the left and right adjacent sampling points, and the waveform change of the sampling point is obvious. The derivative values of the sample points may be taken as first derivative values and form a first derivative value sequence. The resulting first derivative value sequence may be represented as { Y'pmax-1,Y’pmax-2,Y’pmax-3,……,Y’pmax-nN is the number of the first derivative value.
For example, if the derivative values of the respective sampling points in the target P-wave waveform signal are 1, 2, 1, 3, 0, 2, 3, -1, -2, 3, and 3 can be used as the first derivative values to form a first derivative value sequence {2, 3, 3 }. It should be noted that, since the derivative values of the first sampling point and the last sampling point only have one adjacent sampling point, the first sampling point and the last sampling point may not be compared with the adjacent sampling point.
B3: and searching a second derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming a second derivative value sequence by the second derivative values.
The derivative value smaller than the left and right adjacent sampling points indicates that the slope of the sampling point is smaller than the slope of the left and right adjacent sampling points, and indicates that the waveform change of the sampling point is obvious. The derivative values of the sample points may be taken as second derivative values and form a second derivative value sequence. The resulting second derivative value sequence may be represented as { Y'pmin-1,Y’pmin-2,Y’pmin-3,……,Y’pmin-mAnd m is the number of the second derivative value.
Taking the derivative values of the sampling points in the target P-wave waveform signal as 1, 2, 1, 3, 0, 2, 3, -1 and-2 as examples, 1 and 0 can be taken as the second derivative values to form a second derivative value sequence {1,0 }. It should be noted that, since the derivative values of the first sampling point and the last sampling point only have one adjacent sampling point, the first sampling point and the last sampling point may not be compared with the adjacent sampling point. B4: and acquiring the time value of the sampling point corresponding to each first derivative value in the first derivative value sequence and the time value of the sampling point corresponding to each second derivative value in the second derivative value sequence.
And acquiring the time value of the sampling point corresponding to the first derivative value aiming at the sampling point corresponding to each first derivative value in the first derivative value sequence. And acquiring the time value of the sampling point corresponding to the second derivative value aiming at the sampling point corresponding to each second derivative value in the second derivative value sequence.
B5: and intercepting the electrocardiographic waveform signals between the sampling points corresponding to the first derivative values and the sampling points corresponding to the second derivative values adjacent to the time values from the target P-wave waveform signal as candidate P-point waveform signals.
In order to determine the waveform signal of the P point in the target P-wave waveform signal, the waveform signal of the P point, that is, the candidate P-point waveform signal, is determined. Since the P-point is a certain point of the P-wave with a derivative of 0, the candidate P-point waveform signal includes a peak.
If the sampling point corresponding to the first derivative value and the sampling point corresponding to the second derivative value are adjacent to each other in time, the electrocardiographic waveform signal between the sampling point corresponding to the first derivative value and the sampling point corresponding to the second derivative value includes a portion with wavy waveform, and may include a portion with wave crest. And taking the electrocardiographic waveform signal between the two sampling points as a candidate P-point waveform signal. For example, it can be expressed as { (X)pmax-1,Xpmin-1),(Xpmax-2,Xpmin-2),……,(Xpmax-x,Xpmin-x)}. Where x denotes the total number of candidate P-point waveform signals.
B6: and determining the sampling point with the maximum amplitude absolute value in the candidate P point waveform signal as the P point.
And determining the sampling point with the maximum amplitude absolute value as the P point from the determined candidate P point waveform signals.
Based on the above, the candidate P-point waveform signal is determined by the derivative value, and then the sampling point with the maximum amplitude absolute value is determined from the candidate P-point waveform signal as the P point, so that the accuracy of the determined P point can be improved by determining the candidate P-point waveform signal, and the obtained waveform identification result of the electrocardiogram waveform signal is more accurate.
In some possible cases, the interference of noise in the acquired electrocardiographic waveform signal is large, which may affect the determination of the sampling point with the maximum absolute value of the amplitude, resulting in the determination of the P point being not accurate enough.
Further, the embodiment of the present application provides another waveform identification method for an electrocardiographic waveform signal, before determining a sampling point with the largest absolute value of amplitude in a candidate P-point waveform signal as a P point, the method further includes the following three steps:
c1: a second time interval between a start point of the target candidate P-point waveform signal and an end point of the target candidate P-point waveform signal is calculated, and the target candidate P-point waveform signal is any one of the candidate P-point waveform signals.
The time interval occupied by the waveform fluctuation caused by the noise is short. Based on the time interval, the time interval of the target candidate P-point waveform signal can be calculated, and whether the target candidate P-point waveform signal is a noise signal or not can be judged.
A second time interval between the start point and the end point, that is, a time interval of the target candidate P-point waveform signal is calculated based on the start point and the end point of the target candidate P-point waveform signal. The target candidate P-point waveform signal is any one of the obtained candidate P-point waveform signals. The second time interval may be represented by the following equation:
Δti=Xpmin-i-Xpmax-i (1)
wherein i represents the number of target candidate P point waveform signals, and the value of i is a positive integer less than or equal to x.
C2: a first product of the derivative value corresponding to the start point of the target candidate P-point waveform signal and the derivative value corresponding to the end point of the target candidate P-point waveform signal is calculated.
Part of target candidate P-point waveform signals may be relatively flat transition bands and do not include peak bands, and further screening of the target candidate P-point waveform signals is required.
By calculating derivative values corresponding to the starting point and the ending point of the target candidate P point waveform signal, the waveform change of the target candidate P point waveform signal can be judged. A first product of the derivative value corresponding to the start point of the target candidate P-point waveform signal and the derivative value corresponding to the end point of the target candidate P-point waveform signal is calculated. The first product may be represented by:
εi=Y’pmin-i×Y’pmax-i (2)
c3: and if the second time interval is smaller than the first threshold value or the first product is larger than zero, rejecting the target candidate P point waveform signal.
The first threshold is the time interval of the waveform caused by noise. The first threshold may be derived from an average of time intervals of a waveform of the noise disturbance.
If the second time interval is smaller than the first threshold, it can be said that the time interval of the target candidate P-point waveform signal is short, possibly due to jitter caused by noise, and not a normal peak band.
If the first product is greater than zero, it indicates that no peak appears in the target candidate P-point waveform signal. And when the target candidate P point waveform signal meets the condition that the second time interval is smaller than the first threshold value or the first product is larger than zero, the target candidate P point waveform signal cannot meet the condition of determining the P point, and removing.
Specifically, the target candidate P-point waveform signal can be determined by the following formulas (3) and (4):
Figure BDA0003085367290000161
εi<0 (4)
wherein F represents the sampling frequency of the electrocardiographic waveform, and for example, the specific value may be 1000 hz.
Based on the above, by determining the second time interval and the first product of the target candidate P-point waveform signal, it can be determined whether the target candidate P-point waveform signal has a peak band, and the unsatisfactory target candidate P-point waveform signal is rejected. Therefore, more accurate target candidate P point waveform signals are obtained, and accurate P points are obtained.
Correspondingly, there may also be a plurality of peaks in the T-wave of the electrocardiographic waveform signal. In order to accurately determine the T point, the derivative may be first obtained on the target T-wave waveform signal to determine a candidate T-point waveform signal, and then the T point may be determined based on the determined T-point waveform signal.
The embodiment of the present application provides a method for determining a T point from a target T wave waveform signal by deriving the target T wave waveform signal, including:
d1: and carrying out derivation on the target T-wave waveform signal to obtain a derivative value of each sampling point of the target T-wave waveform signal.
In order to determine the waveform change condition of each sampling point in the target T-wave waveform signal, the target T-wave waveform signal is subjected to derivation to obtain a derivative value of each sampling point.
D2: and searching a third derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a third derivative value sequence by the third derivative values.
The derivative value of the sampling point which is larger than the left and right adjacent sampling points shows that the slope of the sampling point is larger than that of the left and right adjacent sampling points, and the waveform change of the sampling point is obvious. The derivative values of the sample points may be taken as third derivative values and form a third derivative value sequence. The resulting sequence of third derivative values may be represented as { Y'tmax-1,Y’tmax-2,Y’tmax-3,……,Y’tmax-aAnd b, wherein a is the number of the third derivative value.
For example, if the derivative values of the respective sampling points in the target T-wave waveform signal are 1, 2, 1, 3, 0, 2, 3, -1, -2, 3, and 3 can be used as the third derivative values to form a third derivative value sequence {2, 3, 3 }. It should be noted that, since the derivative values of the first sampling point and the last sampling point only have one adjacent sampling point, the first sampling point and the last sampling point may not be compared with the adjacent sampling point.
D3: and searching a fourth derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a fourth derivative value sequence by the fourth derivative values.
The derivative value smaller than the left and right adjacent sampling points indicates that the slope of the sampling point is smaller than the slope of the left and right adjacent sampling points, and indicates that the waveform change of the sampling point is obvious. The derivative values of the sample points may be taken as fourth derivative values and form a fourth derivative value sequence. The obtained fourth derivative value sequence can be expressed as { Y'tmin-1,Y’tmin-2,Y’tmin-3,……,Y’tmin-bB is the number of the fourth derivative value.
Taking the derivative values of the sampling points in the target T-wave waveform signal as 1, 2, 1, 3, 0, 2, 3, -1 and-2 as examples, 1 and 0 can be taken as the fourth derivative values to form a fourth derivative value sequence {1,0 }. It should be noted that, since the derivative values of the first sampling point and the last sampling point only have one adjacent sampling point, the first sampling point and the last sampling point may not be compared with the adjacent sampling point.
D4: and acquiring the time value of the sampling point corresponding to each third derivative value in the third derivative value sequence and the time value of the sampling point corresponding to each fourth derivative value in the fourth derivative value sequence.
And acquiring the time value of the sampling point corresponding to the third derivative value aiming at the sampling point corresponding to each third derivative value in the third derivative value sequence. And acquiring the time value of the sampling point corresponding to the fourth derivative value aiming at the sampling point corresponding to each fourth derivative value in the fourth derivative value sequence.
D5: and intercepting the electrocardiographic waveform signals between the sampling points corresponding to the third derivative values and the sampling points corresponding to the fourth derivative values adjacent to the time values from the target T-wave waveform signals as candidate T-point waveform signals.
In order to determine the waveform signal of the T point in the target T-wave waveform signal, the waveform signal of the T point, that is, the candidate T-point waveform signal, is determined. Since the T-point is a certain point where the derivative of the T-wave is 0, the candidate T-point waveform signal includes a peak.
Sampling corresponding to the third derivative valueIf the sampling point corresponding to the fourth derivative value is a sampling point adjacent to the time value, the electrocardiographic waveform signal between the sampling point corresponding to the third derivative value and the sampling point corresponding to the fourth derivative value includes a portion with waveform fluctuation, and possibly includes a portion with a peak. And taking the electrocardiographic waveform signal between the two sampling points as a candidate T-point waveform signal. For example, it can be expressed as { (X)tmax-1,Xtmin-1),(Xtmax-2,Xtmin-2),……,(Xtmax-y,Xtmin-y)}. Wherein y is the total number of candidate T-point waveform signals.
D6: and determining the sampling point with the maximum amplitude absolute value in the candidate T-point waveform signal as the T point.
And determining the sampling point with the maximum amplitude absolute value as the T point from the determined candidate T point waveform signals.
Based on the above, the candidate T-point waveform signal is determined by the derivative value, and then the sampling point with the maximum absolute value of the amplitude is determined from the candidate T-point waveform signal as the T point.
In some possible cases, the interference of noise in the acquired electrocardiographic waveform signal is large, which may affect the determination of the sampling point with the maximum absolute value of the amplitude, resulting in the determination of the T point being not accurate enough.
Further, the embodiment of the present application provides another waveform identification method for an electrocardiographic waveform signal, before determining a sampling point with the largest absolute value of amplitude in a candidate T-point waveform signal as a T point, the method further includes the following three steps:
e1: and calculating a third time interval between the starting point of the target candidate T-point waveform signal and the end point of the target candidate T-point waveform signal, wherein the target candidate T-point waveform signal is any one of the candidate T-point waveform signals.
The time interval occupied by the waveform fluctuation caused by the noise is short. Based on the time interval, the time interval of the target candidate T-point waveform signal can be calculated, and whether the target candidate T-point waveform signal is a noise signal or not can be judged.
And calculating a third time interval between the starting point and the end point according to the starting point and the end point of the target candidate T-point waveform signal, namely the time interval of the target candidate T-point waveform signal. The target candidate T-point waveform signal is any one of the obtained candidate T-point waveform signals. The third time interval may be represented by the following formula:
Δtj=Xtmin-j-Xtmax-j (5)
wherein j represents the number of target candidate T point waveform signals, and the value of j is a positive integer less than or equal to y.
E2: and calculating a second product of the derivative value corresponding to the starting point of the target candidate T-point waveform signal and the derivative value corresponding to the ending point of the target candidate T-point waveform signal.
Part of target candidate T-point waveform signals may be relatively gentle transition bands and do not include peak bands, and further screening of the target candidate T-point waveform signals is required.
By calculating derivative values corresponding to the starting point and the ending point of the target candidate T-point waveform signal, the waveform change of the target candidate T-point waveform signal can be judged. And calculating a second product of the derivative value corresponding to the starting point of the target candidate T-point waveform signal and the derivative value corresponding to the ending point of the target candidate T-point waveform signal. The second product may be represented by:
εj=Y’tmin-j×Y’tmax-j (6)
e3: and if the third time interval is smaller than a second threshold value or the second product is larger than zero, rejecting the target candidate T-point waveform signal.
The second threshold is the time interval of the waveform caused by the noise. The second threshold may be derived from an average of time intervals of a waveform of the noise disturbance.
If the third time interval is smaller than the second threshold, it can be said that the time interval of the target candidate T-point waveform signal is short, possibly due to jitter caused by noise.
And if the second product is larger than zero, the situation that no peak appears in the target candidate T point waveform signal is shown. And when the target candidate T point waveform signal meets the condition that the third time interval is smaller than a second threshold value or the second product is larger than zero, the target candidate T point waveform signal cannot meet the condition of determining the T point, and removing.
Specifically, the target candidate T-point waveform signal can be determined by the formulas (7) and (8):
Figure BDA0003085367290000201
εj<0 (8)
wherein F represents the sampling frequency of the electrocardiographic waveform, and for example, the specific value may be 1000 hz.
Based on the above, by determining the third time interval and the second product of the target candidate T-point waveform signal, it can be determined whether the target candidate T-point waveform signal has a peak band, and the target candidate T-point waveform signal that does not meet the requirement is rejected. Therefore, more accurate target candidate T point waveform signals are obtained, and accurate T points are obtained.
The detected electrocardiographic waveform signal may contain more noise, the characteristics of the electrocardiographic waveform signal are not obvious, and the electrocardiographic waveform signal can be preprocessed before the R point, the Q point and the S point of the complete heart beat of the electrocardiographic waveform signal.
The embodiment of the present application provides another waveform identification method for an electrocardiographic waveform signal, before identifying an R point, a Q point, and an S point of a complete heartbeat from the electrocardiographic waveform signal, the method further includes:
and performing wavelet transformation and nonlinear filtering processing on the electrocardiograph waveform signals to obtain the electrocardiograph waveform signals again.
In order to reduce the influence of noise in the electrocardiographic waveform signal, wavelet transformation may be performed on the electrocardiographic waveform signal first, and then noise reduction processing may be performed on the electrocardiographic waveform signal. And then carrying out nonlinear filtering processing on the electrocardiogram waveform signal so as to enhance the characteristics of the P wave, the QRS wave group and the T wave in the electrocardiogram waveform signal.
Based on the above, the electrocardiographic waveform signals obtained through wavelet transformation and nonlinear filtering processing are more accurate, processing of the electrocardiographic waveform signals is convenient to achieve, and more accurate identification results are obtained.
Based on the waveform identification method of the electrocardiographic waveform signal provided by the method embodiment, the embodiment of the application also provides a waveform identification device of the electrocardiographic waveform signal, and the waveform identification device of the electrocardiographic waveform signal will be described with reference to the accompanying drawings.
Referring to fig. 4, the figure is a schematic structural diagram of a waveform identification apparatus for an electrocardiographic waveform signal according to an embodiment of the present application. As shown in fig. 4, the waveform identifying apparatus for an electrocardiographic waveform signal includes:
an acquisition unit 401 configured to acquire an electrocardiographic waveform signal including a plurality of complete heart beats;
an identifying unit 402, configured to identify an R point, a Q point, and an S point of the complete heartbeat from the electrocardiographic waveform signal, and obtain a time value of the R point, a time value of the Q point, and a time value of the S point;
a first calculating unit 403, configured to calculate an average value of first time intervals between adjacent R points according to the time value of the R point;
a first determining unit 404, configured to determine, according to an average value of first time intervals between the adjacent R points, a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat, where the target complete heartbeat is any one of the complete heartbeats;
a second determining unit 405, configured to intercept, from the electrocardiographic waveform signal, an electrocardiographic waveform signal between a starting point time value of the target complete heartbeat and a time value of a Q point of the target complete heartbeat, as a target P-wave waveform signal;
a third determining unit 406, configured to intercept, from the electrocardiographic waveform signal, an electrocardiographic waveform signal between a time value of an S point of the target complete heartbeat and a time value of an end point of the target complete heartbeat, as a target T-wave waveform signal;
a P-point determining unit 407, configured to determine a P point from the target P-wave waveform signal by performing derivation on the target P-wave waveform signal;
a T-point determining unit 408, configured to determine a T-point from the target T-wave waveform signal by performing derivation on the target T-wave waveform signal.
In a possible implementation manner, the first determining unit 404 includes:
the first determining subunit is configured to multiply an average value of the first time intervals between the R points by a first proportional value to obtain a first time interval, and multiply an average value of the first time intervals between the R points by a second proportional value to obtain a second time interval;
and the second determining subunit is configured to subtract the first time interval from the time value of the R point of the complete heartbeat of the target to obtain a start point time value of the complete heartbeat of the target, and add the second time interval to the time value of the R point of the complete heartbeat of the target to obtain an end point time value of the complete heartbeat of the target.
In a possible implementation manner, the P point determining unit 407 includes:
the first derivation subunit is configured to perform derivation on the target P-wave waveform signal to obtain a derivative value of each sampling point of the target P-wave waveform signal;
the first searching subunit is used for searching a first derivative value which is greater than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming the first derivative value into a first derivative value sequence;
the second searching subunit is used for searching a second derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming a second derivative value sequence by the second derivative value;
the first obtaining subunit is configured to obtain a time value of a sampling point corresponding to each first derivative value in the first derivative value sequence and a time value of a sampling point corresponding to each second derivative value in the second derivative value sequence;
the third determining subunit is used for intercepting an electrocardiographic waveform signal between sampling points corresponding to the first derivative values and the second derivative values adjacent to the time value from the target P-wave waveform signal as a candidate P-point waveform signal;
and the P point determining subunit is used for determining the sampling point with the maximum amplitude absolute value in the candidate P point waveform signal as the P point.
In one possible implementation, the apparatus further includes:
a second calculation unit configured to calculate a second time interval between a start point of a target candidate P-point waveform signal and an end point of the target candidate P-point waveform signal, the target candidate P-point waveform signal being any one of the candidate P-point waveform signals;
a third calculating unit, configured to calculate a first product of a derivative value corresponding to a starting point of the target candidate P-point waveform signal and a derivative value corresponding to an ending point of the target candidate P-point waveform signal;
a first rejecting unit, configured to reject the target candidate P-point waveform signal if the second time interval is smaller than a first threshold or the first product is greater than zero.
In a possible implementation manner, the T point determining unit 408 includes:
the second derivation subunit is configured to perform derivation on the target T-wave waveform signal to obtain a derivative value of each sampling point of the target T-wave waveform signal;
the third searching subunit is used for searching a third derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a third derivative value sequence by the third derivative value;
the fourth searching subunit is used for searching a fourth derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal and forming the fourth derivative value into a fourth derivative value sequence;
the second obtaining subunit is configured to obtain a time value of a sampling point corresponding to each third derivative value in the third derivative value sequence and a time value of a sampling point corresponding to each fourth derivative value in the fourth derivative value sequence;
the fourth determining subunit is used for intercepting an electrocardiographic waveform signal between a sampling point corresponding to a third derivative value and a sampling point corresponding to a fourth derivative value, which are adjacent to each other in time value, from the target T-wave waveform signal as a candidate T-point waveform signal;
and the T point determining subunit is used for determining the sampling point with the maximum amplitude absolute value in the candidate T point waveform signal as the T point.
In one possible implementation, the apparatus further includes:
a fourth calculation unit configured to calculate a third time interval between a start point of a target candidate T-point waveform signal and an end point of the target candidate T-point waveform signal, the target candidate T-point waveform signal being any one of the candidate T-point waveform signals;
a fifth calculating unit, configured to calculate a second product of a derivative value corresponding to a starting point of the target candidate T-point waveform signal and a derivative value corresponding to an ending point of the target candidate T-point waveform signal;
and the second rejecting unit is used for rejecting the target candidate T point waveform signal if the third time interval is smaller than a second threshold value or the second product is larger than zero.
In one possible implementation, the apparatus further includes:
and the filtering unit is used for performing wavelet transformation and nonlinear filtering processing on the electrocardiographic waveform signal to obtain the electrocardiographic waveform signal again.
In addition, an embodiment of the present application further provides a waveform identification apparatus for an electrocardiographic waveform signal, including: a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the method for identifying waveforms of an electrocardiographic waveform signal according to any one of the above embodiments when executing the computer program.
In addition, an embodiment of the present application further provides a computer-readable storage medium, where instructions are stored, and when the instructions are executed on a terminal device, the terminal device is caused to execute the waveform identification method for an electrocardiographic waveform signal according to any one of the above embodiments.
The embodiment of the application provides a waveform identification device and equipment of an electrocardiographic waveform signal, wherein by acquiring the electrocardiographic waveform signal comprising a plurality of complete heart beats, an R point, a Q point and an S point of the complete heart beats and time values respectively corresponding to all the points are firstly identified from the electrocardiographic waveform signal; based on the time value of the R point, the time interval corresponding to a complete heart beat, namely the average value of the first time intervals between the adjacent R points can be obtained; further determining the starting point time value and the ending point time value of any complete heartbeat; based on the determined starting point time value and the Q point time value, a target P wave waveform signal can be obtained; based on the determined S point time value and the determined end point time value, a target T wave waveform signal can be obtained; and respectively carrying out derivation calculation on the target P wave waveform signal and the target T wave waveform signal to determine a P point and a T point. The processing mode of the electrocardiographic waveform signal is simple, and the identification result of the electrocardiographic waveform signal can be obtained quickly. And based on the R point, the Q point and the S point obtained by identification, the P point and the T point can be further determined, the identification of five characteristic points in the electrocardiogram waveform signal is realized, and a more complete and accurate identification result of the electrocardiogram waveform signal is obtained. Through the complete and accurate electrocardio waveform signals obtained by identification, doctors can conveniently diagnose diseases according to the electrocardio waveform signals.
It should be noted that, in the present specification, the embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments may be referred to each other. For the system or the device disclosed by the embodiment, the description is simple because the system or the device corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
It should be understood that in the present application, "at least one" means one or more, "a plurality" means two or more. "and/or" for describing an association relationship of associated objects, indicating that there may be three relationships, e.g., "P and/or B" may indicate: p only, B only and P and B are present simultaneously, wherein P, B may be singular or plural. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship. "at least one of the following" or similar expressions refer to any combination of these items, including any combination of single item(s) or plural items. For example, at least one (one) of P, b, or c, may represent: p, b, c, "P and b", "P and c", "b and c", or "P and b and c", wherein P, b, c may be single or plural.
It is further noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may reside in random access memory (RPM), memory, Read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (10)

1. A method of waveform identification of an electrocardiographic waveform signal, the method comprising:
acquiring an electrocardiographic waveform signal comprising a plurality of complete heart beats;
identifying an R point, a Q point and an S point of the complete heartbeat from the electrocardiographic waveform signal, and acquiring a time value of the R point, a time value of the Q point and a time value of the S point;
calculating the average value of the first time interval between the adjacent R points according to the time value of the R point;
determining a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat according to an average value of first time intervals between the adjacent R points, wherein the target complete heartbeat is any one of the complete heartbeats;
intercepting an electrocardiographic waveform signal between the starting point time value of the target complete heartbeat and the time value of the Q point of the target complete heartbeat from the electrocardiographic waveform signal as a target P-wave waveform signal;
intercepting an electrocardiographic waveform signal between the time value of the S point of the complete target heart beat and the time value of the end point of the complete target heart beat from the electrocardiographic waveform signal as a target T-wave waveform signal;
determining a P point from the target P wave waveform signal by deriving the target P wave waveform signal;
and determining a T point from the target T wave waveform signal by performing derivation on the target T wave waveform signal.
2. The method according to claim 1, wherein determining a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat according to an average value of the first time interval between the adjacent R points comprises:
multiplying the average value of the first time intervals among the R points by a first proportional value to obtain a first time interval, and multiplying the average value of the first time intervals among the R points by a second proportional value to obtain a second time interval;
and subtracting the first time interval from the time value of the R point of the complete heartbeat of the target to obtain a starting point time value of the complete heartbeat of the target, and adding the second time interval to the time value of the R point of the complete heartbeat of the target to obtain an ending point time value of the complete heartbeat of the target.
3. The method of claim 1, wherein determining a P-point from the target P-wave waveform signal by deriving the target P-wave waveform signal comprises:
performing derivation on the target P-wave waveform signal to obtain a derivative value of each sampling point of the target P-wave waveform signal;
searching a first derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming the first derivative value into a first derivative value sequence;
searching a second derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target P-wave waveform signal, and forming a second derivative value sequence by the second derivative values;
acquiring the time value of a sampling point corresponding to each first derivative value in the first derivative value sequence and the time value of a sampling point corresponding to each second derivative value in the second derivative value sequence;
intercepting an electrocardiographic waveform signal between sampling points corresponding to first derivative values and sampling points corresponding to second derivative values adjacent to time values from the target P-wave waveform signal as a candidate P-point waveform signal;
and determining the sampling point with the maximum amplitude absolute value in the candidate P point waveform signal as the P point.
4. The method according to claim 3, wherein before determining as the P point the sample point of the candidate P point waveform signal having the largest absolute value of amplitude, the method further comprises:
calculating a second time interval between a starting point of a target candidate P-point waveform signal and an ending point of the target candidate P-point waveform signal, wherein the target candidate P-point waveform signal is any one of the candidate P-point waveform signals;
calculating a first product of a derivative value corresponding to a starting point of the target candidate P point waveform signal and a derivative value corresponding to an ending point of the target candidate P point waveform signal;
and if the second time interval is smaller than a first threshold value or the first product is larger than zero, rejecting the target candidate P point waveform signal.
5. The method according to claim 1, wherein determining a T point from the target T-wave waveform signal by deriving the target T-wave waveform signal comprises:
performing derivation on the target T-wave waveform signal to obtain a derivative value of each sampling point of the target T-wave waveform signal;
searching a third derivative value which is larger than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a third derivative value sequence by the third derivative values;
searching a fourth derivative value smaller than the derivative values of the left and right adjacent sampling points from the derivative values of each sampling point of the target T-wave waveform signal, and forming a fourth derivative value sequence by the fourth derivative values;
acquiring the time value of a sampling point corresponding to each third derivative value in the third derivative value sequence and the time value of a sampling point corresponding to each fourth derivative value in the fourth derivative value sequence;
intercepting an electrocardiographic waveform signal between sampling points corresponding to third derivative values and sampling points corresponding to fourth derivative values adjacent to time values from the target T-wave waveform signal as a candidate T-point waveform signal;
and determining the sampling point with the maximum amplitude absolute value in the candidate T-point waveform signal as the T point.
6. The method according to claim 5, wherein before determining as T point the sample point of the candidate T point waveform signal having the largest absolute value of amplitude, the method further comprises:
calculating a third time interval between a starting point of a target candidate T-point waveform signal and an ending point of the target candidate T-point waveform signal, wherein the target candidate T-point waveform signal is any one of the candidate T-point waveform signals;
calculating a second product of a derivative value corresponding to a starting point of the target candidate T-point waveform signal and a derivative value corresponding to an ending point of the target candidate T-point waveform signal;
and if the third time interval is smaller than a second threshold value or the second product is larger than zero, rejecting the target candidate T-point waveform signal.
7. The method of any one of claims 1-6, wherein prior to identifying R, Q, and S points of the full heartbeat from the electrocardiographic waveform signal, the method further comprises:
and performing wavelet transformation and nonlinear filtering processing on the electrocardiographic waveform signal to obtain the electrocardiographic waveform signal again.
8. An apparatus for waveform identification of an electrocardiographic waveform signal, the apparatus comprising:
an acquisition unit for acquiring an electrocardiographic waveform signal including a plurality of complete heart beats;
the identification unit is used for identifying an R point, a Q point and an S point of the complete heartbeat from the electrocardiographic waveform signal and acquiring a time value of the R point, a time value of the Q point and a time value of the S point;
the first calculating unit is used for calculating the average value of the first time interval between the adjacent R points according to the time value of the R point;
a first determining unit, configured to determine a starting point time value of a target complete heartbeat and an ending point time value of the target complete heartbeat according to an average value of first time intervals between the adjacent R points, where the target complete heartbeat is any one of the complete heartbeats;
a second determining unit, configured to intercept, from the electrocardiographic waveform signal, an electrocardiographic waveform signal between a starting point time value of the target complete heartbeat and a time value of a Q point of the target complete heartbeat, as a target P-wave waveform signal;
a third determining unit, configured to intercept, from the electrocardiographic waveform signal, an electrocardiographic waveform signal between a time value of an S point of the target complete heartbeat and a time value of an end point of the target complete heartbeat, as a target T-wave waveform signal;
a P-point determining unit configured to determine a P-point from the target P-wave waveform signal by deriving the target P-wave waveform signal;
a T-point determining unit configured to determine a T-point from the target T-wave waveform signal by deriving the target T-wave waveform signal.
9. An apparatus for waveform identification of an electrocardiographic waveform signal, comprising: a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the method of waveform identification of an electrocardiographic waveform signal according to any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium having stored therein instructions which, when run on a terminal device, cause the terminal device to execute the method of waveform identification of an electrocardiographic waveform signal according to any one of claims 1 to 7.
CN202110579117.XA 2021-05-26 2021-05-26 Waveform identification method, device and equipment for electrocardiographic waveform signals Active CN113288169B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110579117.XA CN113288169B (en) 2021-05-26 2021-05-26 Waveform identification method, device and equipment for electrocardiographic waveform signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110579117.XA CN113288169B (en) 2021-05-26 2021-05-26 Waveform identification method, device and equipment for electrocardiographic waveform signals

Publications (2)

Publication Number Publication Date
CN113288169A true CN113288169A (en) 2021-08-24
CN113288169B CN113288169B (en) 2024-05-17

Family

ID=77325264

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110579117.XA Active CN113288169B (en) 2021-05-26 2021-05-26 Waveform identification method, device and equipment for electrocardiographic waveform signals

Country Status (1)

Country Link
CN (1) CN113288169B (en)

Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1057333A (en) * 1996-08-15 1998-03-03 Nec Corp P wave detection method in electrocardiogram analysis, electrocardiograph with analysis function and holter electrocardiogram analysis device
US5758654A (en) * 1996-08-29 1998-06-02 Harley Street Software Ltd. ECG P QRS T onset and peak detection method
CN1370503A (en) * 2001-02-20 2002-09-25 流图私人有限公司 Harmless monitoring utilizing biological impedance of breast and dynamic blood parameter of electrocardiogram
WO2010000009A1 (en) * 2008-07-02 2010-01-07 Cardanal Pty Ltd Improved detection of cardiac dysfunction
CN101799974A (en) * 2010-03-12 2010-08-11 上海交通大学 Electrocardio signal transmission method based on self-adaptive codebook
KR20120116213A (en) * 2011-04-12 2012-10-22 부산대학교 산학협력단 System for discriminating premature ventricular contraction, method for discriminating premature ventricular contraction and the recording media storing the program performing the said method
US8433395B1 (en) * 2009-11-03 2013-04-30 Vivaquant Llc Extraction of cardiac signal data
TW201505606A (en) * 2013-08-07 2015-02-16 Univ Nat Cheng Kung Electrocardiography signal extraction method
CN105877739A (en) * 2016-02-25 2016-08-24 姜坤 Clinical examination method of electrocardio intelligent analyzing system
CN107616791A (en) * 2017-07-26 2018-01-23 浙江大学 The lead electrocardio piecewise-linear reconstruction method of standard 12 based on adaptive electrocardiosignal region segmentation
CN108324268A (en) * 2018-02-26 2018-07-27 河南善仁医疗科技有限公司 A kind of analysis method of electrocardiogram caardiophonogram
CN108836303A (en) * 2018-03-23 2018-11-20 南京大学 ECG Gathering System and its acquisition method based on general-purpose platform
CN109009086A (en) * 2018-07-06 2018-12-18 四川斐讯信息技术有限公司 A kind of R wave of electrocardiosignal detection method and system
CN109009047A (en) * 2018-05-24 2018-12-18 东软集团股份有限公司 A kind of arrhythmia detection method, device and equipment
CN109171711A (en) * 2018-11-14 2019-01-11 中国计量大学 A kind of quick P wave detecting method based on extremum method
CN109691994A (en) * 2019-01-31 2019-04-30 英菲泰克(天津)科技有限公司 A kind of rhythm of the heart analysis method based on electrocardiogram
CN109820501A (en) * 2018-11-12 2019-05-31 浙江清华柔性电子技术研究院 A kind of recognition methods of R wave of electrocardiosignal, device, computer equipment
CN110327032A (en) * 2019-03-27 2019-10-15 苏州平稳芯跳医疗科技有限公司 It is a kind of singly to lead the accurate recognizer of electrocardiosignal PQRST wave joint
CN110680302A (en) * 2019-08-30 2020-01-14 西北大学 Automatic identification method for electrocardiosignal characteristic wave
CN111067511A (en) * 2019-11-21 2020-04-28 东软集团股份有限公司 Method, device and equipment for determining key points in electrocardiogram data
CN111134659A (en) * 2020-01-17 2020-05-12 乐普(北京)医疗器械股份有限公司 Method and device for detecting P wave and T wave in electrocardiosignal

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1057333A (en) * 1996-08-15 1998-03-03 Nec Corp P wave detection method in electrocardiogram analysis, electrocardiograph with analysis function and holter electrocardiogram analysis device
US5758654A (en) * 1996-08-29 1998-06-02 Harley Street Software Ltd. ECG P QRS T onset and peak detection method
CN1370503A (en) * 2001-02-20 2002-09-25 流图私人有限公司 Harmless monitoring utilizing biological impedance of breast and dynamic blood parameter of electrocardiogram
WO2010000009A1 (en) * 2008-07-02 2010-01-07 Cardanal Pty Ltd Improved detection of cardiac dysfunction
US8433395B1 (en) * 2009-11-03 2013-04-30 Vivaquant Llc Extraction of cardiac signal data
CN101799974A (en) * 2010-03-12 2010-08-11 上海交通大学 Electrocardio signal transmission method based on self-adaptive codebook
KR20120116213A (en) * 2011-04-12 2012-10-22 부산대학교 산학협력단 System for discriminating premature ventricular contraction, method for discriminating premature ventricular contraction and the recording media storing the program performing the said method
TW201505606A (en) * 2013-08-07 2015-02-16 Univ Nat Cheng Kung Electrocardiography signal extraction method
CN105877739A (en) * 2016-02-25 2016-08-24 姜坤 Clinical examination method of electrocardio intelligent analyzing system
CN107616791A (en) * 2017-07-26 2018-01-23 浙江大学 The lead electrocardio piecewise-linear reconstruction method of standard 12 based on adaptive electrocardiosignal region segmentation
CN108324268A (en) * 2018-02-26 2018-07-27 河南善仁医疗科技有限公司 A kind of analysis method of electrocardiogram caardiophonogram
CN108836303A (en) * 2018-03-23 2018-11-20 南京大学 ECG Gathering System and its acquisition method based on general-purpose platform
CN109009047A (en) * 2018-05-24 2018-12-18 东软集团股份有限公司 A kind of arrhythmia detection method, device and equipment
CN109009086A (en) * 2018-07-06 2018-12-18 四川斐讯信息技术有限公司 A kind of R wave of electrocardiosignal detection method and system
CN109820501A (en) * 2018-11-12 2019-05-31 浙江清华柔性电子技术研究院 A kind of recognition methods of R wave of electrocardiosignal, device, computer equipment
CN109171711A (en) * 2018-11-14 2019-01-11 中国计量大学 A kind of quick P wave detecting method based on extremum method
CN109691994A (en) * 2019-01-31 2019-04-30 英菲泰克(天津)科技有限公司 A kind of rhythm of the heart analysis method based on electrocardiogram
CN110327032A (en) * 2019-03-27 2019-10-15 苏州平稳芯跳医疗科技有限公司 It is a kind of singly to lead the accurate recognizer of electrocardiosignal PQRST wave joint
CN110680302A (en) * 2019-08-30 2020-01-14 西北大学 Automatic identification method for electrocardiosignal characteristic wave
CN111067511A (en) * 2019-11-21 2020-04-28 东软集团股份有限公司 Method, device and equipment for determining key points in electrocardiogram data
CN111134659A (en) * 2020-01-17 2020-05-12 乐普(北京)医疗器械股份有限公司 Method and device for detecting P wave and T wave in electrocardiosignal

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
毛玲;张国敏;孙即祥;: "基于位置估计与识别后处理的心电信号P波检测", 信号处理, no. 06 *
毛玲;张国敏;孙即祥;: "心电图ST段形态分析方法研究", 信号处理, no. 09 *

Also Published As

Publication number Publication date
CN113288169B (en) 2024-05-17

Similar Documents

Publication Publication Date Title
Sayadi et al. A model-based Bayesian framework for ECG beat segmentation
Dohare et al. An efficient new method for the detection of QRS in electrocardiogram
Li et al. A low-complexity data-adaptive approach for premature ventricular contraction recognition
Kathirvel et al. An efficient R-peak detection based on new nonlinear transformation and first-order Gaussian differentiator
Krasteva et al. QRS template matching for recognition of ventricular ectopic beats
Mehta et al. K-means algorithm for the detection and delineation of QRS-complexes in Electrocardiogram
CN108968941B (en) Arrhythmia detection method, device and terminal
Rahul et al. A novel and lightweight P, QRS, and T peaks detector using adaptive thresholding and template waveform
Last et al. Multi-component based cross correlation beat detection in electrocardiogram analysis
CN109124620B (en) Atrial fibrillation detection method, device and equipment
Kalidas et al. Cardiac arrhythmia classification using multi-modal signal analysis
Islam et al. A morphology alignment method for resampled heartbeat signals
Malik et al. An adaptive QRS detection algorithm for ultra-long-term ECG recordings
Pangerc et al. Robust detection of heart beats in multimodal records using slope-and peak-sensitive band-pass filters
Kuzilek et al. Electrocardiogram beat detection enhancement using Independent Component Analysis
Rao et al. P-and T-wave delineation in ECG signals using parametric mixture Gaussian and dynamic programming
Bruun et al. Automatic atrial fibrillation detection: A novel approach using discrete wavelet transform and heart rate variability
Gholam-Hosseini et al. Detection and extraction of the ECG signal parameters
Khaustov et al. An open-source standard T-wave alternans detector for benchmarking
Colak Preprocessing effects in time–frequency distributions and spectral analysis of heart rate variability
CN111643070B (en) T wave starting point determining method and device, storage medium and electronic equipment
Chen et al. Multi-window detection for P-wave in electrocardiograms based on bilateral accumulative area
CN113288169B (en) Waveform identification method, device and equipment for electrocardiographic waveform signals
Tun et al. Analysis of computer aided identification system for ECG characteristic points
Yoon et al. Automatic detection of congestive heart failure and atrial fibrillation with short RR interval time series

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