A method for analyzing static electrocardiogram

Abstract

A method for analyzing a static electrocardiogram includes: establishing a database; grouping the database; determining a grouping threshold; processing a static electrocardiogram to be determined; and determining the degree of coronary artery obstruction.

Description

Technical field:

[0001] The present invention relates to a method for analyzing an electrocardiogram, and in particular to a method for analyzing a static electrocardiogram. Background technology:

[0002] Generally speaking, a static electrocardiogram can detect whether the heart rate is too slow or too fast, and whether there is persistent arrhythmia. However, some people, even if the degree of coronary artery blockage exceeds 70%, will not have chest pain or chest tightness when resting, but will only have chest pain or chest tightness when exercising. Static electrocardiogram may not be able to detect this paroxysmal heart disease, so an exercise electrocardiogram must be used to perform a stress test on the heart.

[0003] Under normal circumstances, the heart rate will increase and blood pressure will rise after exercise. If the heart rate and blood pressure of the subject slow down or decrease, it indicates that the heart may be seriously lacking oxygen. Although exercise electrocardiogram is non-invasive, it is dangerous for people with potential cardiovascular disease. Once the coronary artery is blocked, exercise will increase the heart load and oxygen consumption, which may easily induce acute myocardial infarction.

[0004] In addition, some physical conditions are not suitable for exercise ECG examinations, such as chest pain caused by unstable myocardial infarction, osteoporosis problems, being too old to run, or people who seldom exercise and find it difficult to even brisk walk or run for 5 or 6 minutes.

[0005] Although exercise electrocardiogram can understand whether the heart is lacking oxygen and whether the blood vessels are blocked, it cannot tell the severity of the coronary artery blockage. This part requires cardiac computed tomography (CT) examination to measure the percentage of blood vessel blockage. Summary of the invention:

[0006] The present invention provides a method for analyzing static electrocardiograms, comprising: establishing a database comprising a plurality of static electrocardiogram data, each of which comprises at least a static electrocardiogram waveform and a plurality of corresponding components, wherein the total number of components of one of the components is i, and at least comprises the three earliest essential modal functions generated in time sequence, which are C1, C2, and C3 in sequence, wherein |Ci|=√(∫(Ci·Ci)dt), and Ctotal=√(∑|Ci|·|Ci|); grouping the static electrocardiogram data, wherein the static electrocardiogram data lacking a corresponding cardiac computed tomography result forms a control group, wherein the static electrocardiogram data comprises a corresponding cardiac computed tomography result, and the result comprises a coronary artery obstruction degree greater than a blockage threshold, forming a first priority group, wherein another corresponding cardiac computed tomography result, and the result comprises another coronary artery obstruction degree less than a blockage threshold. The static electrocardiogram data is used to form a second priority group; for each static electrocardiogram data in the database, a corresponding activity index is calculated = (W1·|C1|+W2·|C2|+W3·|C3|) / |Ctotal|, where W1, W2, and W3 are weight parameters assigned to C1, C2, and C3 respectively, and then all activity pointers in the same group after grouping are summed up and averaged to obtain a control group grouping gate value VM1 and a first priority group grouping gate value VM2. VM2, a second priority group threshold value VM3, wherein W1, W2, and W3 are selected according to the following principles: (1) (W1+W2+W3) is less than 1; (2) W1≧W2≧W3; (3) VM2 and VM3 must be statistically significantly different from VM1 respectively; a static electrocardiogram to be determined is obtained, and an activity index to be determined is calculated; and the activity index to be determined is compared with VM1, VM2, and VM3 to determine the corresponding degree of coronary artery blockage. Description of the drawings:

[0007] Figure 1 It is a static electrocardiogram and its decomposed components.

[0008] Figure 2 The figure is a flow chart of the method for analyzing a resting electrocardiogram according to the present invention.

[0009] Figure 3 FIG. 4 is a schematic diagram showing the relative relationship between group gate values ​​according to an embodiment of the present invention.

[0010] Figure 4 FIG. 4 is a schematic diagram showing the relative relationship between group gate values ​​according to another embodiment of the present invention.

[0011] Reference numerals:

[0012] 10. Resting ECG

[0013] 11 - 18 Component 1 - 8

[0014] 19 Remainder

[0015] 21 Establish a database

[0016] 22 Database grouping

[0017] 23 Determine the grouping threshold

[0018] 24 Process the static electrocardiogram to be judged

[0019] 25 Judge the degree of coronary artery blockage Specific implementation method:

[0020] Using a computer for data processing, in one embodiment, Hilbert - Huang Transformation (HHT) can use the empirical mode decomposition method (EMD) to decompose the original signal into multiple intrinsic mode functions (IMFs), and then use the intrinsic mode functions as the expansion basis. Other commonly used time - frequency analyses also include the short - time Fourier transform (STFT), wavelet transform (WT), etc., the Wigner - Ville distribution (WVD), and the Choi - Williams distribution (CWD).

[0021] Reference Figure 1 , generally speaking, for the static electrocardiogram obtained by an electrocardiogram measuring instrument, after removing the noise, the static electrocardiogram x(t)10 can be decomposed by EMD into the 1st component (C1)11, the 2nd component (C2)12, the 3rd component (C3)13... the 8th component (C8)18, and the remainder (r)19. In terms of the decomposition time series, the 1st component (C1)11 is generated first, followed by the 2nd component (C2)12, and so on for other components. When doing signal processing, the first 8 components 11 - 18 can be used to represent the original static electrocardiogram 10, or maybe only the first 6 components are taken.

[0022] In one embodiment, define the energy index of each component ∣Ci∣ = √(∫(Ci·Ci)dt), and the integration range, for example, is taken as 4 seconds after the signal is stable, or even 10 seconds if it is a bit longer. And the energy index of the static electrocardiogram Ctotal = √(∑∣Ci∣·∣Ci∣).

[0023] Reference Figure 2, which is one embodiment of the present invention, is a method for analyzing a static electrocardiogram, which can determine whether there is coronary artery obstruction and the severity of the obstruction. The method comprises: step 21, establishing a database; step 22, grouping the database; step 23, determining the grouping threshold; step 24, processing the static electrocardiogram to be determined; step 25, determining the degree of coronary artery obstruction.

[0024] In steps 21 and 22, static state ECG data are collected to establish a database and group them. Some people further underwent dynamic ECG examinations, perhaps because the doctor believed that there were no symptoms indicating the need for further examinations, or the examinees did not undergo cardiac computed tomography (CT) because of the heavy cost burden or time constraints. Such static state ECG data are classified as the control group, assuming that their coronary arteries are not blocked, or the blockage is extremely mild, and the coronary arteries can be considered to be in a state that does not require special medical treatment.

[0025] With an electrocardiogram accompanying the CT results, the data on the blockage of the coronary arteries can be known. In one embodiment, the degree of blockage is 50%. If it is greater than 50%, it is an emergency and requires invasive treatment, such as stenting. This type of static electrocardiogram is classified as the first priority group. If it is less than 50%, it can be treated with drugs and is classified as the second priority group. In this way, according to the degree of coronary artery blockage, the database contains three groups of data, one group is a control group that does not require special medical treatment, one group is a first priority group that requires emergency treatment, and one group is a second priority group that does not require emergency treatment.

[0026] In one embodiment, in the database, electrocardiograms without accompanying CT results are classified as a control group and do not require special medical treatment; on the other hand, electrocardiograms with accompanying CT results are classified as a priority group and may require special medical treatment.

[0027] In one embodiment, each static ECG in the database is decomposed into the 1st to 8th components using EMD, and the corresponding energy indexes |C1|, |C2|, |C3|…|C8|, and Ctotal are calculated. Therefore, each static ECG data in the database, in addition to the original static ECG, will at least have the calculated energy index. And each static ECG data of the first priority group and the second priority group will also have coronary artery occlusion data.

[0028] In one embodiment, each static ECG data in the database will have at least the calculated energy index in addition to the original static ECG data, and each static ECG data in the priority group will also have the data of coronary artery occlusion.

[0029] In step 23, the grouping threshold is determined. The activity index corresponding to each static electrocardiogram data is calculated as (W1·|C1|+W2·|C2|+W3·|C3|) / |Ctotal|, where W1, W2, and W3 are weight parameters assigned to |C1|, |C2|, and |C3|, respectively. The activity pointers of all static electrocardiogram data in the same group are summed up and averaged to obtain the grouping threshold value VM1 of the control group, the grouping threshold value VM2 of the first priority group, and the grouping threshold value VM3 of the second priority group.

[0030] In one embodiment, the group gate value of the control group is VM1, and the group gate value of the priority group is VM2.

[0031] The selection of weight parameters W1, W2, and W3 is based on the following principles: (1) (W1+W2+W3) is less than 1; (2) W1≧W2≧W3; (3) VM2 and VM3 should be statistically significantly different from VM1, that is, the statistical correlation should be less than 0.05. In one embodiment, W1 is 0.3-0.4, W2 is 0.3-0.4, and W3 is 0.1-0.2. Figure 3 In one embodiment, VM1 is 22.30%, VM2 is 30.75%, VM3 is 26.72%, the statistical correlation between VM3 and VM1 is less than 0.001, and the statistical correlation between VM2 and VM1 is less than 0.008.

[0032] refer to Figure 4 In one embodiment, the grouping threshold value VM1 of the control group is 22.30%, the grouping threshold value VM2 of the priority group is 27.96%, and the statistical correlation between VM2 and VM1 is less than 0.001.

[0033] After determining each grouping gate value and each weight parameter, each static electrocardiogram data containing coronary artery blockage data in the database also has a corresponding activity index, so that a corresponding relationship between the activity index and the coronary artery blockage data can be established.

[0034] Step 24, processing the static electrocardiogram to be determined. In one embodiment, a static electrocardiogram to be determined is obtained from outside the database, and then decomposed into 8 components by EMD, and then the activity index to be determined is calculated using these 8 components.

[0035] Step 25, determine the degree of coronary artery obstruction. Compare the activity pointer to be determined with the group gate value of each group in the database. In one embodiment, less than 22.30% is a condition that does not require special medical treatment, between 22.30% and 26.72% is a condition that does not require emergency medical treatment, and greater than 26.72% is a condition that may require emergency medical treatment.

[0036] In one embodiment, a rate less than 22.30% is a condition that does not require special medical treatment, and a rate greater than 22.30% is a condition that may require medical treatment.

[0037] In one embodiment, the established correspondence between the activity index and the coronary artery blockage data may be referred to, and the percentage of coronary artery blockage may be obtained from the activity index to be determined.

Claims

1. A method for analyzing a static electrocardiogram, comprising: Establishing a database includes collecting multiple static electrocardiogram data, each of which includes at least one static electrocardiogram waveform and corresponding multiple components, characterized in that: One of the components has a total number of components i and includes at least the three essential modal functions that are generated earliest in time sequence, which are C1, C2, and C3 in order, wherein |Ci|=√(∫(Ci·Ci)dt), and Ctotal=√(∑|Ci|·|Ci|); The static electrocardiogram data are grouped, wherein each piece of the static electrocardiogram data lacking a corresponding cardiac computed tomography result is grouped to form a control group, wherein each piece of the static electrocardiogram data containing a corresponding cardiac computed tomography result and the result containing a corresponding coronary artery occlusion data greater than a occlusion threshold is grouped to form a first priority group, wherein each piece of the static electrocardiogram data containing another corresponding cardiac computed tomography result and the result containing another corresponding coronary artery occlusion data less than a occlusion threshold is grouped to form a second priority group; For each static electrocardiogram data in the database, a corresponding activity index = (W1·|C1|+W2·|C2|+W3·|C3|) / |Ctotal| is calculated, where W1, W2, and W3 are weight parameters assigned to C1, C2, and C3, respectively. Then, all activity indexes in the same group after grouping are summed up and averaged to obtain a control group gate value VM1, a first priority group gate value VM2, and a second priority group gate value VM3, respectively. The selection of W1, W2, and W3 is based on the following principles: (1) (W1+W2+W3) is less than 1; (2) W1≧W2≧W3; (3) VM2 and VM3 must be statistically significantly different from VM1; Obtaining a static electrocardiogram to be determined and calculating an activity index to be determined; and The activity index to be determined is compared with VM1, VM2, and VM3 to determine the degree of blockage of a corresponding coronary artery.

2. The method according to claim 1, characterized in that The occlusion threshold of the coronary artery is 50%.

3. The method according to claim 1, characterized in that The essential mode function is generated using an empirical mode decomposition method.

4. The method according to claim 1, characterized in that: The corresponding coronary artery obstruction data contained in each piece of the static electrocardiogram data in the first priority group and the second priority group is expressed as a percentage.

5. The method according to claim 4, characterized in that In the first priority group and the second priority group, a corresponding relationship between an activity index and coronary artery obstruction data is established.

6. A method for analyzing a static electrocardiogram, characterized in that: Include: Establishing a database, comprising collecting a plurality of static electrocardiogram data, each of which comprises at least a static electrocardiogram waveform and a plurality of corresponding components, wherein a total number of components of the components is i, and at least comprises the three essential modal functions generated earliest in time sequence, which are C1, C2, and C3 in sequence, wherein |Ci|=√(∫(Ci·Ci)dt), and Ctotal=√(∑|Ci|·|Ci|); The static electrocardiogram data are grouped into groups, wherein each piece of the static electrocardiogram data lacking a corresponding cardiac computed tomography result is grouped into a control group, and each piece of the static electrocardiogram data containing a corresponding cardiac computed tomography result and the result containing a piece of coronary artery occlusion data is grouped into a priority group; For each static electrocardiogram data in the database, a corresponding activity index = (W1·|C1|+W2·|C2|+W3·|C3|) / |Ctotal| is calculated, where W1, W2, and W3 are weight parameters assigned to C1, C2, and C3, respectively. Then, all activity indexes in the same group after grouping are summed up and averaged to obtain a control group grouping threshold value VM1 and a priority group grouping threshold value VM2, respectively. The selection of W1, W2, and W3 is based on the following principles: (1) (W1+W2+W3) is less than 1; (2) W1≧W2≧W3; (3) VM2 and VM1 must be statistically significantly different; Obtaining a static electrocardiogram to be determined and calculating an activity index to be determined; and The activity index to be determined is compared with VM1 and VM2 to determine a corresponding degree of coronary artery obstruction.

7. The method according to claim 6, characterized in that The essential mode function is generated using an empirical mode decomposition method.

Images