Electrocardiosignal detection effect evaluation method and system

Abstract

The invention relates to an electrocardiosignal detection effect evaluation method and system, and belongs to the technical field of electrocardiograph monitoring. The method comprises the steps that electrocardiosignals with R-wave classification marks are obtained, differential conversion is conducted on the electrocardiosignals, and electrocardiograph differential signals are obtained; constructing a standard R wave differential waveform corresponding to each type of R waves based on the electrocardiogram differential signal; calculating the similarity ratio of each R-wave segment in the electrocardiogram differential signal to the standard R-wave differential waveform of the corresponding type; and obtaining a detection effect evaluation result of each type of R waves based on the similarity ratio corresponding to each type of R waves. The problem that an objective and effective detection effect evaluation method is lacked in the prior art is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of electrocardiogram monitoring, and in particular to a method and system for evaluating detection effect of electrocardiogram signals. Background Art

[0002] The ECG waveform extracted from the body surface under dynamic conditions is different from the ECG under static conditions. It is affected by external interference, human body movement, and lead wire connection. The obtained waveform is full of various noises, electromyographic drift, and motion artifacts, which brings great obstacles to computer automatic recognition.

[0003] Generally speaking, the accuracy of the detection results of long-term dynamic ECG data is not high and cannot meet the level of medical application. ECG analysis technicians need to intervene to correct the waveforms misjudged by the computer to avoid subsequent diagnostic errors. In the whole process, the recognition of QRS complex is the basis for subsequent intervention. The accuracy of QRS complex recognition directly affects the efficiency and workload of subsequent technician intervention. At present, the accuracy of existing detection methods has not reached the thousandth level based on the MIT ECG database. The evaluation of detection effect mainly relies on manual judgment, which means that the 24-hour manual intervention of technicians greatly increases the workload of ECG technicians. There is a lack of objective and effective detection effect evaluation methods. Summary of the invention

[0004] In view of the above analysis, the embodiments of the present invention aim to provide a method and system for evaluating detection effect of an electrocardiogram signal, so as to solve the problem of lack of an objective and effective detection effect evaluation method in the prior art.

[0005] On the one hand, an embodiment of the present invention provides a method for evaluating detection effect of an electrocardiogram signal, comprising the following steps:

[0006] Acquire an electrocardiogram signal with an R wave classification mark, and perform differential conversion on the electrocardiogram signal to obtain an electrocardiogram differential signal;

[0007] Based on the ECG differential signal, a standard R wave differential waveform corresponding to each type of R wave is constructed;

[0008] Calculate the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform;

[0009] The detection effect evaluation result of each type of R wave is obtained based on the similarity ratio corresponding to each type of R wave.

[0010] Based on the further improvement of the above method, the standard R wave differential waveform corresponding to each type of R wave is constructed based on the ECG differential signal, including:

[0011] For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal;

[0012] A bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window;

[0013] The standard vertebral R wave differential waveform corresponding to this type of R wave is obtained according to the central differential waveform and the boundary differential waveform of the standard R wave corresponding to this type of R wave.

[0014] Based on the further improvement of the above method, a bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window, including:

[0015] Determine the left initial position corresponding to the left sliding window and the right initial position corresponding to the right sliding window respectively according to the width of the central differential waveform of the standard R wave;

[0016] Using a left sliding window to slide leftward from the left initial position in all such R wave segments, and determining the left boundary differential waveform of the standard R wave corresponding to the such R wave according to the differential mean of the R wave in each left sliding window;

[0017] A right sliding window is used to slide rightward from the right initial position in all such R wave segments, and the right boundary differential waveform of the standard R wave corresponding to this type of R wave is determined according to the differential mean of the R wave in each right sliding window.

[0018] Based on the further improvement of the above method, for each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal, including:

[0019] For each signal point of the central differential waveform of the standard R wave, the differential mean of all sampling points corresponding to the signal point in this type of R wave is calculated as the value of the signal point; and the central differential waveform of the standard R wave is obtained.

[0020] Based on the further improvement of the above method, a left sliding window is used to slide leftward from the left initial position in all such R wave segments, and the left boundary differential waveform of the standard R wave corresponding to the such R wave is determined according to the differential mean of the R wave in each left sliding window, including:

[0021] S2221, taking each left sliding window at the left initial position as the current sliding window;

[0022] Calculate the differential mean of each sampling point in the current sliding window as the value of the corresponding signal point in the differential waveform of the left boundary of the standard R wave; the mean of each sampling point in the current sliding window constitutes a first set, and the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set; slide the current sliding window to the left by the first step length as the current sliding window;

[0023] S2222, calculating the differential mean of each sampling point in the current sliding window as the value of the corresponding signal point in the differential waveform of the left boundary of the standard R wave; the mean of each sampling point in the current sliding window constitutes a first set, and the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set;

[0024] S2223. If the positive and negative signs corresponding to the current sliding window are different from those corresponding to the previous sliding window, stop sliding the window to obtain the left boundary differential waveform of the standard R wave; otherwise, slide the current sliding window to the left by the first step length as the current sliding window, and return to step S2222.

[0025] Based on the further improvement of the above method, the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set, including:

[0026] If the ratio of the number of values ​​with positive signs in the first set to the total number of values ​​in the first set exceeds a first threshold, the current sliding window corresponds to a positive sign, otherwise the current sliding window corresponds to a negative sign.

[0027] Based on the further improvement of the above method, the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform is calculated, including:

[0028] Perform integral reduction on each type of standard R wave differential waveform to obtain a standard R wave reduction waveform;

[0029] For each type of standard R wave, extract each R wave segment corresponding to the type in the ECG differential signal, perform integral restoration, and use the restored R wave segment as the target R wave;

[0030] The ratio of each standard R-wave restoration waveform to the corresponding signal point of each target R-wave of the corresponding type is calculated to construct a similarity ratio set.

[0031] Based on the further improvement of the above method, the detection effect evaluation results of each type of R wave are obtained based on the similarity ratio corresponding to each type of R wave, including:

[0032] It is determined whether the similarity ratio corresponding to each type of R wave satisfies the normal distribution. If so, the detection result of this type of R wave is accurate, otherwise, the detection result is inaccurate.

[0033] On the other hand, an embodiment of the present invention provides a detection effect evaluation system for an electrocardiogram signal, comprising the following modules:

[0034] An ECG signal conversion module is used to obtain an ECG signal with an R wave classification mark, and perform differential conversion on the ECG signal to obtain an ECG differential signal;

[0035] A standard R wave construction module is used to construct a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal;

[0036] A similarity calculation module is used to calculate the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform;

[0037] The evaluation module is used to obtain the detection effect evaluation result of each type of R wave based on the similarity ratio corresponding to each type of R wave.

[0038] Based on the further improvement of the above system, the standard R wave construction module constructs a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal, including:

[0039] For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal;

[0040] A bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window;

[0041] The standard vertebral R wave differential waveform corresponding to this type of R wave is obtained according to the central differential waveform and the boundary differential waveform of the standard R wave corresponding to this type of R wave.

[0042] The present invention performs differential conversion on the ECG signal, regards the ECG differential waveform as the myocardial polarization intensity at a specific position of myocardial conduction, thereby eliminating baseline interference, calculates the standard R wave differential waveform segment corresponding to each type according to the converted ECG signal, calculates the similarity of each R wave in the ECG signal segment with respect to the standard R wave based on the standard R wave differential waveform, and obtains the ECG signal detection effect evaluation result based on the similarity result, thereby weakening the ECG QRS detection method based on the ECG polarization waveform, and using its differential waveform instead, which can be more conducive to the effective evaluation of the R wave detection effect, realizes the objective and effective evaluation of the ECG detection effect, and helps to effectively improve and enhance the automation level of ECG analysis.

[0043] In the present invention, the above-mentioned technical solutions can also be combined with each other to achieve more preferred combination solutions. Other features and advantages of the present invention will be described in the subsequent description, and some advantages can become obvious from the description, or can be understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained through the contents particularly pointed out in the description and the drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] The accompanying drawings are only used for the purpose of illustrating specific embodiments and are not to be considered as limiting the present invention. In the entire drawings, the same reference symbols represent the same components;

[0045] Figure 1 It is a flow chart of a detection effect evaluation method of an electrocardiogram signal according to an embodiment of the present invention;

[0046] Figure 2 The block diagram of the detection effect evaluation system of the electrocardiogram signal according to the embodiment of the present invention. DETAILED DESCRIPTION

[0047] The preferred embodiments of the present invention are described in detail below in conjunction with the accompanying drawings, wherein the accompanying drawings constitute a part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, but are not used to limit the scope of the present invention.

[0048] A specific embodiment of the present invention discloses a method for evaluating the detection effect of an electrocardiogram signal. Figure 1 As shown, the following steps are included:

[0049] S1, obtaining an electrocardiogram signal with an R wave classification mark, performing differential conversion on the electrocardiogram signal, and obtaining an electrocardiogram differential signal;

[0050] S2, constructing a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal;

[0051] S3, calculating the similarity ratio between each R wave segment in the ECG differential signal and the standard R wave differential waveform of the corresponding type;

[0052] S4. Obtain the detection effect evaluation result of each type of R wave based on the similarity ratio corresponding to each type of R wave.

[0053] The present invention performs differential conversion on the ECG signal, regards the ECG differential waveform as the myocardial polarization intensity at a specific position of myocardial conduction, thereby eliminating baseline interference, calculates the standard R wave differential waveform segment corresponding to each type according to the converted ECG signal, calculates the similarity of each R wave in the ECG signal segment with respect to the standard R wave based on the standard R wave differential waveform, and obtains the ECG signal detection effect evaluation result based on the similarity result, thereby weakening the ECG QRS detection method based on the ECG polarization waveform, and using its differential waveform instead, which can be more conducive to the effective evaluation of the R wave detection effect, realizes the objective and effective evaluation of the ECG detection effect, and helps to effectively improve and enhance the automation level of ECG analysis.

[0054] At present, in the process of QRS complex identification, the main focus is on processing the original ECG waveform, using filtering and variable threshold judgment methods, with the emphasis on strengthening the distinction of the R wave. The most representative one is the Pan&Tomskins algorithm, which uses a digital bandpass filter (5-15Hz), squares after differentiation, and then performs a moving window integration (150ms). The result of this is to highlight the mutation characteristics of the R wave and reduce the false positive probability caused by the tall T wave; the results are then screened and tested by dynamic thresholds to calibrate the R wave position for subsequent waveform clustering and classification.

[0055] During implementation, the ECG signal segment with R wave classification mark obtained in step S1 can be an ECG signal with R wave type marked after Pan & Tomskins algorithm. It should be noted that the ECG signal with R wave classification mark marks the center position of the R wave, generally the peak point, and the type of the R wave mark point. For example, the R wave type includes three types: normal (N), supraventricular premature beat (S), and ventricular premature beat (V).

[0056] In order to eliminate noise interference, the acquired ECG signal is differentially converted, and the differential waveform of the ECG is regarded as the myocardial polarization intensity at a specific position of myocardial conduction.

[0057] For example, the original ECG signal is represented by X = {x0, x1, x2, ..., x N The first signal point d1 after conversion is calculated by (x1-x0)×f, and the second signal point d2 after conversion is calculated by (x2-x1)×f. And so on, the ECG differential signal D={d1,d2,…,d N}, where N represents the total number of signal points in the ECG signal segment. Since the first signal point x0 in the original ECG signal has no differential value, it is deleted and the original ECG signal is represented by X = {x1, x2, …, x N}. Where f represents the sampling frequency.

[0058] By differentiating the waveform at each point, the dynamic information of all original ECG data is retained, the AC component of the ECG in the low-frequency band is emphasized, the DC component and baseline drift component are weakened, which is conducive to highlighting the ECG characteristics.

[0059] A standard R wave differential waveform corresponding to each type of R wave is constructed based on the ECG differential signal. Specifically, step S2 includes:

[0060] S21. For each type of R wave, obtain a central differential waveform of a standard R wave corresponding to the type of R wave according to a differential mean value of the type of R wave in the electrocardiographic differential signal;

[0061] S22, using a bidirectional sliding window to synchronously slide all such R wave segments, and obtaining a boundary differential waveform of a standard R wave corresponding to such R wave according to the differential mean value of the R wave in each sliding window;

[0062] S23. Obtain a standard vertebral R wave differential waveform corresponding to this type of R wave according to the central differential waveform and the boundary differential waveform of the standard R wave corresponding to this type of R wave.

[0063] The construction of the standard R wave differential waveform includes the construction of the central part waveform and the construction of the boundary waveform.

[0064] For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal, including:

[0065] S211. For each signal point of the central differential waveform of the standard R wave, calculate the differential mean of all sampling points corresponding to the signal point in this type of R wave as the value of the signal point; and obtain the central differential waveform of the standard R wave.

[0066] Since the original ECG signal X and the ECG differential signal D are of equal length and correspond one to one, the R wave marker in the original signal can find the corresponding position in the ECG differential signal.

[0067] For example, to construct the central differential waveform of the standard R wave, it is necessary to first determine the width of the central differential waveform of the standard R wave. According to the beating law of the heart, the R wave is generally greater than 60ms. If the sampling interval is 1ms during implementation, the central waveform width can be set to 60 sampling intervals, with the R wave peak as the center, and the left and right widths are 30 sampling intervals respectively, that is, there are 30 signal points on the left and right.

[0068] For each signal point of the central differential waveform of the standard R wave, the differential mean of all sampling points corresponding to the signal point in this type of R wave is calculated as the value of the signal point in the standard R wave.

[0069] For example, the first signal point on the left side of the center position of the central differential waveform of the normal standard R wave is extracted, and the differential value corresponding to the first sampling point on the left side of each normal R wave mark point in the ECG differential signal is calculated, and the average value is calculated to obtain the value of the first signal point on the left side of the center position of the central differential waveform of the standard R wave. Similarly, after constructing the value of each signal point in the central wave differential shape, the central waveform is obtained, which is expressed as S i ={s i-k ,…,s i-1 ,s i1 ,…,s ik}. Among them, s i-1 is the first signal point on the left side of the center position of the central differential waveform of the i-th type of standard R wave, s i1is the first signal point to the right of the center position of the central differential waveform of the i-th type of standard R wave, and 2k is the width of the central differential waveform.

[0070] After the central differential waveform is constructed, it is necessary to construct the boundary waveform of the standard R wave differential waveform.

[0071] Specifically, a bidirectional sliding window is used to synchronously slide all such R wave segments, and a boundary differential waveform of a standard R wave corresponding to such R wave is obtained according to the differential mean value of the R wave in each sliding window, including:

[0072] S221, determining a left initial position corresponding to the left sliding window and a right initial position corresponding to the right sliding window respectively according to the width of the central differential waveform of the standard R wave;

[0073] S222, using a left sliding window to slide leftward from the left initial position in all such R wave segments, and determining a left boundary differential waveform of a standard R wave corresponding to such R wave according to a differential mean value of the R wave in each left sliding window;

[0074] S223, using a right sliding window to slide rightward from the right initial position in all such R wave segments, and determining the right boundary differential waveform of the standard R wave corresponding to this type of R wave according to the differential mean of the R wave in each right sliding window.

[0075] First, the initial position of the left sliding window is determined based on half the width of the standard R wave center differential waveform. For example, the width of the center waveform is 60 signal points, that is, 30 signal points on the left and right of the peak. The left initial position is the position of the 31st sampling point on the left side of the R mark point, and the right initial position is the position of the 31st sampling point on the right side of the R mark point.

[0076] The process of determining the left boundary differential waveform and the right boundary differential waveform of the standard R wave is the same, the difference lies in the sliding direction, and the following is explained by taking the construction of the left boundary differential waveform as an example.

[0077] When using a sliding window, you must first determine the window size and sliding step size. During implementation, you can determine the size based on computational efficiency and accuracy. For example, if the sliding window size is determined to be 10 sampling intervals, each sliding window corresponds to 10 sampling points, and the sliding step size is 1 sampling interval, then one sampling point is slid to the left each time.

[0078] Specifically, in step S222, the left boundary differential waveform of the standard R wave corresponding to the type of R wave is determined in the following manner:

[0079] S2221, taking each left sliding window at the left initial position as the current sliding window; calculating the differential mean of each sampling point in the current sliding window as the value of the signal point corresponding to the left boundary differential waveform of the standard R wave; the mean of each sampling point in the current sliding window constitutes a first set, and the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set; sliding the current sliding window to the left by the first step length as the current sliding window;

[0080] S2222, calculating the differential mean of each sampling point in the current sliding window as the value of the corresponding signal point in the differential waveform of the left boundary of the standard R wave; the mean of each sampling point in the current sliding window constitutes a first set, and the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set;

[0081] S2223. If the positive and negative signs corresponding to the current sliding window are different from those corresponding to the previous sliding window, stop sliding the window to obtain the left boundary differential waveform of the standard R wave; otherwise, slide the current sliding window to the left by the first step length as the current sliding window, and return to step S2222.

[0082] For example, the left initial position is the 31st sampling point on the left side of the R wave marker, and the sequence D = {d1, d2, …, d N} at a certain position R d There is an R wave marker (the R marker point may be between two sampling points or at a certain sampling point), and the type is normal. The first sampling point to the left of the marker point is d i , then the left sliding window is initially positioned at the marked point R d The 31st to 40th sampling points on the left (counting from right to left), the data in the sliding window is d i-30 ~d i-39 .

[0083] Extract the data of the 31st to 40th sampling points on the left side of each normal R wave mark point, and calculate the mean of each sampling point, that is, calculate the mean of the values ​​of the 31st sampling point on the left side of the normal R wave peak as the value of the 31st signal point on the left side of the standard R wave peak, calculate the mean of the values ​​of the 32nd sampling point on the left side of the normal R wave peak as the value of the 32nd signal point on the left side of the standard R wave peak, and so on. These means constitute the first set.

[0084] According to the signs of the means in the first set, the sign corresponding to the current window is determined. If the ratio of the number of values ​​with positive signs in the first set to the total number of values ​​in the first set exceeds the first threshold, the current sliding window corresponds to a positive sign, otherwise the current sliding window corresponds to a negative sign. During implementation, the first threshold can be set to 2 / 3.

[0085] After obtaining the sign of the first sliding window, the current sliding window is slid to the left by the first step as the current sliding window; at this time, the data covered by the sliding window is the data of the 32nd to 41st sampling points on the left side of the R wave mark point, and the mean value at each moment is calculated according to the above process as the data of the 32nd to 41st signal points on the left side of the standard R wave peak (at this time, the 32nd to 40th sampling points have been calculated, so only the mean value of the 41st sampling point needs to be calculated), and the means of the 32nd to 41st sampling points constitute a new first set, and the sign corresponding to the current window is determined according to the sign of these means in the new first set. If the sign corresponding to the current sliding window is different from the sign corresponding to the previous sliding window, the sliding window is stopped, and the left boundary waveform of the standard R wave is obtained, that is, the value of the 31st to 41st signal points on the left side of the standard R wave is obtained. If the sign corresponding to the current sliding window is the same as the sign corresponding to the previous sliding window, the current sliding window is slid to the left by the first step as the current sliding window, and the step S2222 is returned.

[0086] The same process is used to obtain the right boundary differential waveform of the standard R wave. It should be noted that at this time, the sliding window starts to slide to the right from the 31st sampling point on the right side of the R wave mark point.

[0087] The center differential waveform, left boundary differential waveform, and right boundary differential waveform of the standard R wave are obtained, that is, the standard R wave differential waveform corresponding to the R wave type is obtained. The constructed standard R wave differential waveform excludes the influence of DC component and baseline drift; only the medium and low frequency AC component is retained, which is the dynamic information of the ECG waveform. It can make detection judgment more accurate.

[0088] After obtaining the standard R wave differential waveform, all ECG R waves of this type are similarly calculated with reference to the standard R wave. Specifically, step S3 calculates the similarity ratio between each R wave segment in the ECG differential signal and the standard R wave differential waveform of its corresponding type, including:

[0089] S31, performing integral reduction on each type of standard R wave differential waveform to obtain a standard R wave reduction waveform;

[0090] S32, for each type of standard R wave, extract each R wave segment corresponding to the type from the ECG differential signal, perform integral restoration, and use the restored R wave segment as the target R wave;

[0091] S33, calculating the ratio of the corresponding signal point of each standard R wave restoration waveform and each target R wave of the corresponding type to construct a similarity ratio set.

[0092] Since the standard R wave is obtained based on the differential waveform, it is necessary to restore it, that is, to convert the standard R wave in the form of integration. The standard R wave differential waveform is represented by S i ={s i-m ,…,si-1 ,s i1 ,…,s in}, where s i-1 is the first signal point on the left side of the center position of the ith type standard R wave differential waveform, s i1 is the first signal point on the right side of the center position of the i-th standard R-wave differential waveform, s i-m is the left boundary signal point of the i-th standard R-wave differential waveform, s in is the right boundary signal point of the i-th standard R wave differential waveform. The total number of signal points in the standard R wave differential waveform is m+n.

[0093] Integral restoration is to use the following formula to calculate the corresponding integral restored value for each signal point in the standard R wave differential waveform:

[0094]

[0095] Among them, sr ip represents the value of the pth signal point in the i-th standard R-wave restoration waveform, s iq represents the value of the qth signal point in the i-th standard R-wave differential waveform, and r0 represents the preset baseline value. p∈[1,m+n]. During implementation, r0 takes a positive value so that each sr ip The standard R wave restoration waveform is expressed as SR i ={sr i-m ,…,sr i-1 ,sr i1 ,...,sr in For the standard R wave of the i-th type, extract each R wave segment corresponding to the type in the ECG differential signal. That is, find an R wave marker of this type, extract its left m sampling points and right n sampling points to form an R wave segment of this type, and use the same method as above to perform integral reduction to obtain the target waveform, which is expressed as {xr i-m ,…,xr i-1 ,xr i1 ,...,xr in}.

[0096] There are two differences between the standard R wave restoration waveform and the target R wave and the usual ECG wave: they are inverted and the difference in the baseline starting point is removed, thus avoiding the influence of baseline drift.

[0097] Construct the corresponding similarity ratio set as If there are M i-th type R waves, M similarity ratio sets will be constructed.

[0098] The data set in the ratio set is the shape ratio of the R waveform. The waveform of the R wave is affected by many random factors, including electrophysiology, muscle movement, external electromagnetic interference, and the characteristics of the sampling circuit itself. These random factors create a ratio set that obeys or approximately obeys a normal distribution.

[0099] Therefore, calculate the mean μ and variance σ of all values ​​in the M similarity ratio sets 2 , according to μ and variance σ 2 Evaluate the QRS wave detection effect.

[0100] Specifically, based on the similarity ratio corresponding to each type of R wave, the detection effect evaluation result of each type of R wave is obtained, including:

[0101] It is determined whether the similarity ratio corresponding to each type of R wave satisfies the normal distribution. If so, the detection result of this type of R wave is accurate, otherwise, the detection result is inaccurate.

[0102] The detection effect is determined by whether all the values ​​in the M similarity ratio sets meet the normal distribution. If they do, the detection result of the i-th category is accurate, otherwise the detection result is inaccurate. Assuming that there are Nm numbers in the M similarity ratio sets, if the following three conditions are met, they meet the normal distribution:

[0103] (1) The ratio of the number of all values ​​falling within the range [μ-σ, μ+σ] to the total number Nm is approximately equal to the second threshold;

[0104] (2) The ratio of the number of all values ​​falling within the range [μ-2σ,μ+2σ] to the total number Nm is approximately equal to the third threshold;

[0105] (3) The ratio of the number of all values ​​falling within the range of [μ-3σ,μ+3σ] to the total number Nm is approximately equal to the fourth threshold.

[0106] In implementation, the second threshold, the third threshold and the fourth threshold may be set to 68%, 95% and 99.7% respectively. When the ratio falls within the range of 5% above and below the threshold, it may be considered to be approximately equal.

[0107] The above is called the "68-95-99.7" rule. If the normal distribution is satisfied, the smaller the variance, the more accurate the classification of the class; if the variance is too large, the classification is rough.

[0108] For each target R wave, the similarity of each target R wave is obtained according to the ratio of the number of elements in the corresponding similarity ratio set that fall within the confidence interval to the number of elements in the similarity ratio set.

[0109] A target R wave and a standard R wave have the following similarity ratio sets:

[0110] {1.44, 1.22, 1.09, 1.33, 1.15, 1.29, 1.15, 1.27, 1.27, 1.13, 1.16, 1.15, 1.15, 1.19, 1.15, 1.3, 1.41, 1.2, 1.24, 1.22, 1.04, 1.25, 1.21, 1.01, 1.14, 1.25, 1.19, 1.2, 1.38, 1.21, 1.19, 1.09, 1.36, 1.17, 1.35, 1.1, 1.11, 1. 19, 1.24, 1.29, 1.11, 1.18, 1.09, 1.08, 1.28, 1.36, 1.18, 1.19, 1.21, 1.11, 1.15, 1.27, 1.07, 1.15, 0.92, 1.24, 1.12, 1.24, 1.18, 1.32}; if μ=1.2, σ=0.05, the corresponding confidence interval is expressed as [μ-zσ, μ+zσ], where z represents the confidence scale corresponding to the confidence level, for example, the corresponding value of 94.5% confidence level is 2. Then the number of ratios in the similarity ratio set that fall within the confidence interval is 45, so the similarity value of the target R wave is 45 / 60=75%.

[0111] The accuracy of detection can be improved by eliminating target R waves with small similarity. Eliminating target R waves with small similarity can reduce the variance of the normal distribution, which is convenient for the subsequent use of ECG signals.

[0112] A specific embodiment of the present invention discloses a detection effect evaluation system for an electrocardiogram signal, such as Figure 2 As shown, it includes the following modules:

[0113] An ECG signal conversion module is used to obtain an ECG signal with an R wave classification mark, and perform differential conversion on the ECG signal to obtain an ECG differential signal;

[0114] A standard R wave construction module is used to construct a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal;

[0115] A similarity calculation module is used to calculate the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform;

[0116] The evaluation module is used to obtain the detection effect evaluation result of each type of R wave based on the similarity ratio corresponding to each type of R wave.

[0117] Preferably, the standard R wave construction module constructs a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal, including:

[0118] For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal;

[0119] A bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window;

[0120] The standard vertebral R wave differential waveform corresponding to this type of R wave is obtained according to the central differential waveform and the boundary differential waveform of the standard R wave corresponding to this type of R wave.

[0121] The above method embodiments and system embodiments are based on the same principle, and their related parts can be used for reference, and can achieve the same technical effect. The specific implementation process refers to the above embodiment, which will not be repeated here.

[0122] Those skilled in the art will appreciate that all or part of the processes of the above-mentioned embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, wherein the computer-readable storage medium is a disk, an optical disk, a read-only storage memory, or a random access memory, etc.

[0123] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by any technician familiar with the technical field within the technical scope disclosed by the present invention should be covered within the protection scope of the present invention.

Claims

1. A method for evaluating detection effect of an electrocardiogram signal, characterized in that: The following steps are involved: Acquire an electrocardiogram signal with an R wave classification mark, and perform differential conversion on the electrocardiogram signal to obtain an electrocardiogram differential signal; Based on the ECG differential signal, a standard R wave differential waveform corresponding to each type of R wave is constructed; Calculate the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform; The detection effect evaluation result of each type of R wave is obtained based on the similarity ratio corresponding to each type of R wave.

2. The detection effect evaluation method of the electrocardiogram signal according to claim 1, characterized in that: Based on the ECG differential signal, a standard R wave differential waveform corresponding to each type of R wave is constructed, including: For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal; A bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window; The standard vertebral R wave differential waveform corresponding to this type of R wave is obtained according to the central differential waveform and the boundary differential waveform of the standard R wave corresponding to this type of R wave.

3. The detection effect evaluation method of the electrocardiogram signal according to claim 2, characterized in that: A bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window, including: Determine the left initial position corresponding to the left sliding window and the right initial position corresponding to the right sliding window respectively according to the width of the central differential waveform of the standard R wave; Using a left sliding window to slide leftward from the left initial position in all such R wave segments, and determining the left boundary differential waveform of the standard R wave corresponding to the such R wave according to the differential mean of the R wave in each left sliding window; A right sliding window is used to slide rightward from the right initial position in all such R wave segments, and the right boundary differential waveform of the standard R wave corresponding to this type of R wave is determined according to the differential mean of the R wave in each right sliding window.

4. The detection effect evaluation method of the electrocardiogram signal according to claim 2, characterized in that: For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal, including: For each signal point of the central differential waveform of the standard R wave, the differential mean of all sampling points corresponding to the signal point in this type of R wave is calculated as the value of the signal point; and the central differential waveform of the standard R wave is obtained.

5. The method for evaluating detection effect of an electrocardiogram signal according to claim 3, characterized in that: Using a left sliding window to slide leftward from the left initial position in all such R wave segments, determining the left boundary differential waveform of the standard R wave corresponding to the such R wave according to the differential mean of the R wave in each left sliding window, including: S2221, taking each left sliding window at the left initial position as the current sliding window; Calculate the differential mean of each sampling point in the current sliding window as the value of the corresponding signal point in the differential waveform of the left boundary of the standard R wave; the mean of each sampling point in the current sliding window constitutes a first set, and the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set; slide the current sliding window to the left by the first step length as the current sliding window; S2222, calculating the differential mean of each sampling point in the current sliding window as the value of the corresponding signal point in the differential waveform of the left boundary of the standard R wave; the mean of each sampling point in the current sliding window constitutes a first set, and the positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set; S2223. If the positive and negative signs corresponding to the current sliding window are different from those corresponding to the previous sliding window, stop sliding the window to obtain the left boundary differential waveform of the standard R wave; otherwise, slide the current sliding window to the left by the first step length as the current sliding window, and return to step S2222.

6. The detection effect evaluation method of the electrocardiogram signal according to claim 5, characterized in that: The positive and negative signs corresponding to the current sliding window are obtained according to the positive and negative signs of each value in the first set, including: If the ratio of the number of values ​​with positive signs in the first set to the total number of values ​​in the first set exceeds a first threshold, the current sliding window corresponds to a positive sign, otherwise the current sliding window corresponds to a negative sign.

7. The detection effect evaluation method of an electrocardiogram signal according to claim 1, characterized in that: Calculate the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform, including: Perform integral reduction on each type of standard R wave differential waveform to obtain a standard R wave reduction waveform; For each type of standard R wave, extract each R wave segment corresponding to the type in the ECG differential signal, perform integral restoration, and use the restored R wave segment as the target R wave; The ratio of each standard R-wave restoration waveform to the corresponding signal point of each target R-wave of the corresponding type is calculated to construct a similarity ratio set.

8. The method for evaluating detection effect of an electrocardiogram signal according to claim 1, characterized in that: Based on the similarity ratio corresponding to each type of R wave, the detection effect evaluation result of each type of R wave is obtained, including: It is determined whether the similarity ratio corresponding to each type of R wave satisfies the normal distribution. If so, the detection result of this type of R wave is accurate, otherwise, the detection result is inaccurate.

9. A detection effect evaluation system for an electrocardiogram signal, characterized in that: Includes the following modules: An ECG signal conversion module is used to obtain an ECG signal with an R wave classification mark, and perform differential conversion on the ECG signal to obtain an ECG differential signal; A standard R wave construction module is used to construct a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal; A similarity calculation module is used to calculate the similarity ratio between each R wave segment in the ECG differential signal and its corresponding type of standard R wave differential waveform; The evaluation module is used to obtain the detection effect evaluation result of each type of R wave based on the similarity ratio corresponding to each type of R wave.

10. The detection effect evaluation system of the electrocardiogram signal according to claim 9, characterized in that: The standard R wave construction module constructs a standard R wave differential waveform corresponding to each type of R wave based on the ECG differential signal, including: For each type of R wave, the central differential waveform of the standard R wave corresponding to the type of R wave is obtained according to the differential mean value of the type of R wave in the ECG differential signal; A bidirectional sliding window is used to synchronously slide all such R wave segments, and the boundary differential waveform of the standard R wave corresponding to this type of R wave is obtained according to the differential mean value of the R wave in each sliding window; The standard vertebral R wave differential waveform corresponding to this type of R wave is obtained according to the central differential waveform and the boundary differential waveform of the standard R wave corresponding to this type of R wave.

Images