Method, device, system and storage medium for evaluating heart health status

Abstract

The present disclosure relates to the technical field of medical instruments, in particular to a heart health state evaluation method, processing device, system and computer readable storage medium, by respectively acquiring exercise electrocardiogram data, resting electrocardiogram data and conventional electrocardiogram data output by electrocardiogram leads and analyzing the same to determine the possibility of existence of myocardial ischemia, coronary artery lesion, myocardial bridge and coronary microcirculation disorder, qualitatively evaluating the heart health state according to the determination result, and in the case of determining the possibility of existence of myocardial ischemia and / or coronary artery lesion, quantitatively evaluating the myocardial ischemia and / or coronary artery lesion, and finally evaluating the heart health state according to the qualitative evaluation result and the quantitative evaluation result. The technical scheme provided by the embodiments of the present application can more accurately qualitatively and quantitatively analyze the heart health state by combining exercise electrocardiogram data, resting electrocardiogram data and conventional ECG data.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of medical instruments, in particular to a heart health state evaluation method, processing device, system and computer readable storage medium. BACKGROUND

[0002] At present, there is a scheme for evaluating the heart health state based on electrocardiogram (ECG). Chinese patent application with publication number CN114742114A discloses a high-frequency QRS waveform curve analysis method, which evaluates the degree of myocardial ischemia under coronary stenosis by analyzing the electrocardiogram signal under exercise load, but this method only considers the myocardial ischemia caused by coronary stenosis and cannot more accurately evaluate the heart health state. SUMMARY

[0003] The purpose of the embodiments of the present application is to provide a heart health state evaluation method, processing device, system and computer readable storage medium, which can more accurately qualitatively and quantitatively analyze the heart health state

[0004] To achieve the above-mentioned purpose, the first aspect of the present application provides a heart health state evaluation method, comprising:

[0005] Respectively acquiring exercise electrocardiogram data, resting electrocardiogram data and conventional electrocardiogram data output through at least one electrocardiogram lead;

[0006] Determining the possibility of myocardial ischemia according to the exercise electrocardiogram data, the resting electrocardiogram data and the conventional electrocardiogram data;

[0007] Determining the possibility of coronary artery lesion according to the exercise electrocardiogram data;

[0008] Determining the possibility of myocardial bridge according to the exercise electrocardiogram data;

[0009] Determining the possibility of coronary microcirculation disorder according to the exercise electrocardiogram data and the resting electrocardiogram data;

[0010] Qualitatively evaluating the heart health state according to the determination results of the myocardial ischemia possibility, the coronary artery lesion possibility, the myocardial bridge possibility and the coronary microcirculation disorder possibility;

[0011] In the case of determining the possibility of myocardial ischemia and / or the possibility of coronary artery lesion, quantitatively evaluating the myocardial ischemia and / or the coronary artery lesion; and

[0012] Evaluating the heart health state according to the qualitative evaluation result and the quantitative evaluation result.

[0013] In the embodiments of the present application, the possibility of existence of myocardial ischemia is determined according to exercise electrocardiogram data, resting electrocardiogram data and conventional electrocardiogram data, and the method comprises:

[0014] A first lead positive index corresponding to each electrocardiogram lead is obtained according to first high-frequency QRS complex data obtained from the exercise electrocardiogram data;

[0015] A second lead positive index corresponding to each electrocardiogram lead is obtained according to second high-frequency QRS complex data obtained from the resting electrocardiogram data;

[0016] A third lead positive index corresponding to each electrocardiogram lead is obtained according to the conventional electrocardiogram data;

[0017] In the case that any one of the first lead positive index, the second lead positive index and the third lead positive index indicates positivity, it is determined that the possibility of existence of myocardial ischemia.

[0018] In the embodiments of the present application, in the case that only one of the first lead positive index, the second lead positive index and the third lead positive index indicates positivity, the possibility of existence of myocardial ischemia is low;

[0019] In the case that two of the first lead positive index, the second lead positive index and the third lead positive index indicate positivity, the possibility of existence of myocardial ischemia is high;

[0020] In the case that all of the first lead positive index, the second lead positive index and the third lead positive index indicate positivity, the possibility of existence of myocardial ischemia is the highest.

[0021] In the embodiments of the present application, the possibility of existence of coronary artery lesion is determined according to exercise electrocardiogram data, and the method comprises:

[0022] A first point sequence corresponding to each electrocardiogram lead is obtained according to first high-frequency QRS complex data corresponding to each electrocardiogram lead obtained from the exercise electrocardiogram data, and the first point sequence comprises a plurality of sampling points arranged in time sequence;

[0023] For the first point sequence corresponding to any electrocardiogram lead, in an iteration loop, any sampling point in the plurality of sampling points in the first point sequence is taken as a current reference sampling point, and sampling points in a first preset time interval from the current reference sampling point are iterated to determine whether there are a first sampling point and a second sampling point satisfying a preset condition in the first preset time interval, wherein the first sampling point is earlier than the second sampling point in time sequence;

[0024] In a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has the first waveform feature indicating the coronary artery lesion, so that it is determined that the coronary artery lesion possibility exists;

[0025] In a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition do not exist, the current traversal loop ends, and other sampling points in the first point sequence except the current reference sampling point are determined as new reference sampling points;

[0026] The next traversal loop is performed for the new reference sampling points, until it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, or all the sampling points to be traversed in the first point sequence are traversed.

[0027] In the embodiments of the present application, the determination of the coronary artery lesion possibility according to the exercise electrocardiogram data further includes:

[0028] In a case where the first sampling point and the second sampling point satisfying the preset condition do not exist after all the sampling points to be traversed in the first point sequence are traversed, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead does not have the first waveform feature indicating the coronary artery lesion;

[0029] In a case where the first high-frequency QRS complex data corresponding to all the electrocardiogram leads does not have the first waveform feature indicating the coronary artery lesion, it is determined that the coronary artery lesion possibility does not exist.

[0030] In the embodiments of the present application, the preset condition includes:

[0031] The amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value, and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value; or

[0032] The amplitude drop relative value between the first sampling point and the second sampling point reaches the first preset threshold value, or the voltage drop degree between the first sampling point and the second sampling point reaches the second preset threshold value.

[0033] In the embodiments of the present application, the determination of the coronary artery lesion possibility according to the exercise electrocardiogram data further includes:

[0034] The candidate first high-frequency QRS complex data is determined from the first high-frequency QRS complex data corresponding to each electrocardiogram lead, and the candidate first high-frequency QRS complex data does not have the first waveform feature indicating the coronary artery lesion;

[0035] The second point sequence corresponding to the candidate first high-frequency QRS complex data is obtained, and the second point sequence includes a plurality of sampling points arranged in time sequence;

[0036] For the second point sequence corresponding to any candidate first high-frequency QRS complex data, in a traversal loop, any sampling point in the plurality of sampling points in the second point sequence is taken as a current reference sampling point, and sampling points in a second preset time interval from the current reference sampling point are traversed to determine whether there are third sampling points and fourth sampling points satisfying another preset condition in the second preset time interval, wherein the third sampling point is earlier than the fourth sampling point in time sequence;

[0037] In a case where it is determined that there are third sampling points and fourth sampling points satisfying another preset condition, it is determined that the arbitrary candidate first high-frequency QRS complex data has a second waveform feature indicating a myocardial bridge, thereby determining the possibility of a myocardial bridge;

[0038] In a case where it is determined that there are no third sampling points and fourth sampling points satisfying another preset condition, the current traversal loop ends, and other sampling points in the second point sequence except the current reference sampling point are determined as new reference sampling points;

[0039] The next traversal loop is performed for the new reference sampling points, until it is determined that there are third sampling points and fourth sampling points satisfying another preset condition, or all sampling points in the second point sequence that need to be traversed are traversed.

[0040] In the embodiments of the present application, the determination of the possibility of a myocardial bridge according to the exercise electrocardiogram data comprises:

[0041] In a case where all sampling points in the second point sequence that need to be traversed are traversed and there are no third sampling points and fourth sampling points satisfying another preset condition, it is determined that the corresponding candidate first high-frequency QRS complex data does not have a second waveform feature indicating a myocardial bridge;

[0042] In a case where all candidate first high-frequency QRS complex data do not have a second waveform feature indicating a myocardial bridge, it is determined that there is no myocardial bridge.

[0043] In the embodiments of the present application, another preset condition comprises:

[0044] The amplitude drop relative value between the third sampling point and the fourth sampling point reaches a third preset threshold value, and the voltage drop degree between the third sampling point and the fourth sampling point reaches a fourth preset threshold value, and the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point shows a continuous downward trend; or

[0045] The amplitude drop relative value between the third sampling point and the fourth sampling point reaches a third preset threshold value or the voltage drop degree between the third sampling point and the fourth sampling point reaches a fourth preset threshold value, and the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point shows a continuous downward trend.

[0046] In the embodiments of the present application, the determining the possibility of the existence of the myocardial bridge according to the exercise electrocardiogram data comprises:

[0047] determining candidate first high-frequency QRS complex data from the first high-frequency QRS complex data corresponding to each electrocardiogram lead, the candidate first high-frequency QRS complex data being free of the first waveform feature indicating the coronary artery lesion;

[0048] For any candidate first high-frequency QRS complex data, at least one sampling point in the candidate first high-frequency QRS complex data located in a preset time period is determined as a reference sampling point, the preset time period being a time period from the moment of the end of exercise in the load exercise electrocardiogram detection process to a moment after the end of exercise;

[0049] determining the to-be-traversed sampling points according to the determined reference sampling points, wherein the to-be-traversed sampling points are sampling points in a preset time interval before the reference sampling points in the candidate first high-frequency QRS complex data, and the preset time interval includes a plurality of continuous subintervals;

[0050] traversing the sampling points in the plurality of subintervals to determine whether there is a target sampling point in each subinterval that meets another preset condition corresponding to the subinterval;

[0051] In the case where there is a target sampling point in any one of the plurality of subintervals that meets the corresponding preset condition, it is determined that there is a possibility of the existence of the myocardial bridge.

[0052] In the embodiments of the present application, the determining the possibility of the existence of the coronary microcirculation disorder according to the exercise electrocardiogram data and the resting electrocardiogram data comprises:

[0053] determining whether there is target first high-frequency QRS complex data in the first high-frequency QRS complex data corresponding to the first lead positive indicator indicating the positivity;

[0054] In the case where it is determined that there is the target first high-frequency QRS complex data, it is determined that there is a possibility of the existence of the coronary microcirculation disorder.

[0055] In the embodiments of the present application, in the case where it is determined that there is the target first high-frequency QRS complex data, the possibility of the existence of the coronary microcirculation disorder is determined, comprising:

[0056] determining a target electrocardiogram lead corresponding to the target first high-frequency QRS complex data;

[0057] determining whether the second lead positive indicator corresponding to the target electrocardiogram lead indicates the positivity;

[0058] In the case where the second lead positive indicator corresponding to the target electrocardiogram lead indicates the positivity, it is determined that there is a possibility of the existence of the coronary microcirculation disorder.

[0059] In the embodiments of the present application, determining the possibility of existence of coronary microcirculation disorder according to the exercise ECG data and the resting ECG data further comprises:

[0060] determining a first maximum root mean square voltage according to the root mean square voltages in the first high-frequency QRS complex data corresponding to all ECG leads and a second maximum voltage according to the voltages in the second high-frequency QRS complex data corresponding to all ECG leads;

[0061] determining the number of ECG leads corresponding to the positive indication of the second lead positive indicator;

[0062] In the case that the number of ECG leads is less than the preset value and at least the first maximum root mean square voltage is greater than the first voltage threshold, it is determined that there is no coronary microcirculation disorder.

[0063] In the embodiments of the present application, the quantitative evaluation of myocardial ischemia comprises:

[0064] obtaining a first lead number according to the first lead positive indicator and a second lead number according to the second lead positive indicator, wherein the first lead number is the number of ECG leads corresponding to the positive or critical indication of the first lead positive indicator, and the second lead number is the number of ECG leads corresponding to the positive or critical indication of the second lead positive indicator;

[0065] analyzing the resting ECG data to determine a target high-frequency morphology index;

[0066] normalizing the first lead number, the second lead number and the target high-frequency morphology index to obtain a first ischemia index, a second ischemia index and a third ischemia index, respectively;

[0067] obtaining a target ischemia index according to the first ischemia index, the second ischemia index and the third ischemia index.

[0068] In the embodiments of the present application, the target ischemia index is the product of the first ischemia index, the second ischemia index and the third ischemia index, or the minimum one of the first ischemia index, the second ischemia index and the third ischemia index is taken as the target ischemia index.

[0069] The second aspect of the present application provides a processing device, comprising a processor and a memory, the memory storing instructions, the processor being configured to call and execute the instructions from the memory to implement the above-mentioned heart health status evaluation method.

[0070] The third aspect of the present application provides a heart health status evaluation system, comprising:

[0071] An electrocardiosignal collection device includes at least one electrode forming at least one electrocardiogram lead for collecting electrocardiosignals of a subject; and

[0072] The processing device described above.

[0073] The fourth aspect of the present application provides a machine readable storage medium, which stores instructions for causing a machine to execute the above-mentioned heart health status evaluation method.

[0074] The technical solution provided by the embodiments of the present application can more accurately qualitatively and quantitatively analyze myocardial ischemia by combining exercise electrocardiosignal data, resting electrocardiosignal data and conventional ECG data. In addition, compared with only focusing on the steep drop wave in the first three minutes of exercise, the technical solution can improve the accuracy of coronary artery lesion qualification by judging whether there is a steep drop wave in the high-frequency QRS complex data in the whole load exercise detection process, thereby improving the recognition sensitivity of coronary artery lesions without significantly reducing or even not reducing the specificity, and providing a qualitative identification method for myocardial bridge and coronary microcirculation disorder. Therefore, the evaluation results of myocardial ischemia qualitative and quantitative analysis, coronary artery lesion qualitative and quantitative analysis, myocardial bridge qualitative analysis and coronary microcirculation disorder evaluation can be comprehensively evaluated to more comprehensively and accurately evaluate the heart health status.

[0075] Other features and advantages of the embodiments of the present application will be described in detail in the following specific embodiments. BRIEF DESCRIPTION OF DRAWINGS

[0076] The accompanying drawings are included to provide a further understanding of the embodiments of the present application, and constitute a part of the specification, and are used together with the following specific embodiments to explain the embodiments of the present application, but do not constitute a limitation on the embodiments of the present application. In the drawings:

[0077] Figure 1 An example flowchart of a heart health status evaluation method according to an embodiment of the present application is schematically shown.

[0078] Figure 2A An exercise high-frequency QRS waveform curve is schematically shown, in which a steep drop waveform appears after 3 minutes of exercise.

[0079] Figure 2B An example flowchart of a method for determining the possibility of existence of coronary artery lesions according to an embodiment of the present application is schematically shown.

[0080] Figure 3 An example flowchart of a method for determining the possibility of existence of coronary artery lesions according to another embodiment of the present application is schematically shown.

[0081] Figure 4An example flowchart of a method for determining the likelihood of the presence of a coronary artery lesion according to another embodiment of the present application is schematically shown.

[0082] Figure 5 An example flowchart of a method for determining the likelihood of the presence of a myocardial bridge according to an embodiment of the present application is schematically shown.

[0083] Figure 6 An example flowchart of a method for determining the likelihood of the presence of a myocardial bridge according to an embodiment of the present application is schematically shown.

[0084] Figure 7 An example flowchart of a method for determining the likelihood of the presence of a myocardial bridge according to another embodiment of the present application is schematically shown.

[0085] Figure 8 An example flowchart of a method for determining the likelihood of the presence of a myocardial bridge according to another embodiment of the present application is schematically shown.

[0086] Figure 9 An example flowchart of a method for determining the likelihood of the presence of a myocardial bridge according to yet another embodiment of the present application is schematically shown.

[0087] Figure 10 An example flowchart of a method for determining the likelihood of the presence of a myocardial bridge according to yet another embodiment of the present application is schematically shown.

[0088] Figure 11 An example flowchart of a method for determining the likelihood of the presence of a coronary microcirculatory disorder according to an embodiment of the present application is schematically shown.

[0089] Figure 12 An example block diagram of a heart health status assessment system according to an embodiment of the present application is schematically shown. DETAILED DESCRIPTION

[0090] The specific embodiments of the present application will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are merely exemplary and are not intended to limit the present application.

[0091] It should be noted that the acquisition, transmission, storage, use, processing, etc. of data in the technical solutions of the present application comply with relevant provisions of laws and regulations. In the embodiments of the present application, some industry existing solutions, components, models, etc. may be mentioned, which should be considered as exemplary, and the purpose is only to illustrate the feasibility of the implementation of the technical solutions of the present application, but it does not mean that the applicant has or will necessarily use the solutions.

[0092] The heart health state evaluation method provided in the embodiments of the present application can be applied to a terminal, a server, an interactive system including a terminal and a server, and is implemented through the interaction of the terminal and the server, which is not specifically limited herein. The terminal can be, but is not limited to, various personal computers, notebook computers, smart phones, tablet computers, electrocardiogram monitoring devices, and portable wearable devices, and the server can be implemented by an independent server or a server cluster composed of multiple servers.

[0093] The "exercise electrocardiogram data" or equivalent terms mentioned in the embodiments of the present application can refer to electrocardiogram data collected by a subject during a load exercise electrocardiogram detection process. Load exercise electrocardiogram detection is a method of increasing cardiac load through a certain amount of exercise to collect electrocardiogram data of the subject, and analyzing the heart health status of the subject based on the collected electrocardiogram data. It is widely used in the detection of heart disease and cardiovascular disease. The exercise electrocardiogram data can include multiple QRS complexes reflecting changes in left and right ventricular depolarization potential and time. Each QRS complex is a set of Q waves, R waves, and S waves in an electrocardiogram. Based on the QRS complex in the exercise electrocardiogram data, a corresponding high-frequency QRS waveform curve can be analyzed. The high-frequency QRS waveform curve can be a curve obtained by connecting multiple sampling points in time sequence. The high-frequency QRS waveform curve, also known as the high-frequency QRS time intensity curve, can be used to represent the trend of the root mean square voltage of the high-frequency component of the QRS complex of the subject with respect to time during the entire load exercise electrocardiogram detection process, that is, to reflect the energy change trend during the entire load exercise electrocardiogram detection process. The high-frequency QRS waveform curve is presented through a high-frequency QRS waveform graph. In the high-frequency QRS waveform graph, the horizontal coordinate can be time, corresponding to the detection time of the load exercise electrocardiogram detection process, and the unit can be min (minute). The vertical coordinate can be the root mean square voltage (RMS voltage), which can also be understood as intensity or amplitude, and the unit can be uV (microvolt).

[0094] In an embodiment, the load exercise electrocardiogram detection process includes multiple stages, specifically including three stages of a resting stage, an exercise stage, and a recovery stage in sequence, and the exercise electrocardiogram data includes electrocardiogram data of each stage. It can be understood that the division of stages is not limited to this, and can be divided according to actual conditions. In the embodiments of the present application, the pre-exercise stage is located in the resting stage, the exercise stage is located in the exercise stage, and the post-exercise stage is located in the recovery stage.

[0095] The "resting electrocardio data" or its equivalent terms mentioned in the embodiments of the present application can refer to electrocardio data collected in a resting electrocardio detection process. The subject is in a resting state in the resting electrocardio detection process. The resting electrocardio data can include a plurality of QRS complexes. The QRS complexes in the resting electrocardio data are sequentially subjected to alignment, averaging, and high-frequency filtering processing to obtain high-frequency QRS complex data, or the QRS complexes in the resting electrocardio data are sequentially subjected to high-frequency filtering, alignment, and averaging processing to obtain high-frequency QRS complex data, or high-frequency electrocardio data is extracted from the resting electrocardio data by analysis, and the QRS complexes in the high-frequency electrocardio data are sequentially subjected to alignment and averaging processing to obtain high-frequency QRS complex data, which is not specifically limited here. The high-frequency QRS complex data corresponds to the high-frequency QRS envelope curve, and the corresponding high-frequency QRS envelope curve can be formed based on the high-frequency QRS complex data. As can be seen, the corresponding high-frequency QRS envelope curve can be obtained by data processing of the resting electrocardio data.

[0096] The "routine ECG data" or its equivalent terms mentioned in the embodiments of the present application can refer to ECG data that can be obtained from a routine electrocardiogram (ECG) signal. In the embodiments of the present application, the routine ECG data can be obtained by analyzing and processing the resting electrocardio data, for example, by low-pass filtering the resting electrocardio data to filter out the high-frequency components and retain the low-frequency effective electrocardio components, and the routine ECG data is obtained after the filtering processing.

[0097] The heart health status evaluation method according to the embodiments of the present application respectively performs myocardial ischemia qualitative evaluation, coronary artery lesion (which can include, for example, coronary stenosis and coronary sclerosis) qualitative evaluation, myocardial bridge qualitative evaluation, and coronary microcirculation disorder qualitative evaluation on the subject based on the acquired exercise electrocardio data, resting electrocardio data, and routine ECG data, respectively, comprehensively qualitatively evaluates the heart health status based on the qualitative evaluation results of the above-mentioned items, and quantitatively evaluates the heart health status based on at least one of the myocardial ischemia quantitative evaluation, coronary artery lesion quantitative evaluation, myocardial bridge quantitative evaluation, and coronary microcirculation disorder quantitative evaluation on the basis of the comprehensive qualitative evaluation of the heart health status.

[0098] Figure 1 An example flowchart of the heart health status evaluation method according to the embodiments of the present application is schematically shown. Specifically, as shown in Figure 1 The heart health status evaluation method can include the following steps.

[0099] In step S101, exercise electrocardio data, resting electrocardio data, and routine electrocardiogram (ECG) data output by at least one electrocardiogram lead are respectively acquired.

[0100] Specifically, in the process of collecting electrocardiogram, at least one electrode patch can be used to collect signals from the subject. For example, taking 10 electrode patches as an example, the electrode patches can be distributed on the chest and limbs of the human body to form 12 electrocardiogram leads (such as V1, V2, V3, V4, V5, V6, I, II, III, aVL, aVF and aVR), corresponding to 12 groups of electrocardiogram data. It can be understood that 10 electrode patches are only used as an example and do not specifically limit the number of electrode patches, and more or fewer electrode patches can be used according to actual needs. The electrocardiogram signal collection of the subject in the process of load exercise can obtain exercise electrocardiogram data, and the electrocardiogram signal collection of the subject in the resting state can obtain resting electrocardiogram data and also obtain conventional ECG data.

[0101] In step S102, the possibility of myocardial ischemia is determined according to the exercise electrocardiogram data, the resting electrocardiogram data and the conventional ECG data.

[0102] Specifically, in the embodiments of the present application, the exercise electrocardiogram data can be analyzed and processed to obtain high-frequency QRS complex data, which can be specifically obtained by analyzing the high-frequency components of the QRS complex in the exercise electrocardiogram data. Specifically, the exercise electrocardiogram data includes ECG (electrocardiogram) corresponding to each heartbeat of the subject during the entire load exercise electrocardiogram detection process, and the ECG includes QRS complex. The exercise electrocardiogram data is divided into a plurality of electrocardiogram data subsets by a window function according to time sequence and a preset moving step length, and each electrocardiogram data subset includes ECG corresponding to multiple heartbeats. For each electrocardiogram data subset, the ECG or QRS complex corresponding to the multiple heartbeats included therein is sequentially aligned, averaged and band-pass filtered to obtain a corresponding high-frequency QRS complex (high-frequency band of QRS complex), and the root mean square voltage of the high-frequency QRS complex is calculated to obtain a corresponding root mean square voltage as the root mean square voltage corresponding to the electrocardiogram data subset. Thus, a time-intensity data point set arranged in time sequence can be obtained, each data point corresponding to a time and a root mean square voltage (or intensity or amplitude), and the high-frequency QRS complex data can include the data point set (for convenience, it can be referred to as first high-frequency QRS complex data below). Based on the time-intensity data point set, a corresponding high-frequency QRS time-intensity curve can be obtained, in other words, the visualization of the time-intensity data point set is the aforementioned high-frequency QRS time-intensity curve.

[0103] It can be understood that the window length of the window function and the preset moving step length can be customized according to actual needs, for example, the window length is set to 10 seconds, and the preset moving step length is set to 10 seconds or one heartbeat period, which refers to the time interval between adjacent two heartbeats, which is not specifically limited herein. According to time sequence refers to the order of the collection time of the signal / the detection time of the load exercise electrocardiogram detection process.

[0104] The first high frequency QRS complex data can be used to obtain a lead positive indicator corresponding to each ECG lead (for ease of distinction, hereinafter referred to as a first lead positive indicator). The first lead positive indicator can be obtained using methods disclosed in the prior art. For example, in an example, the amplitude drop relative value and the amplitude absolute value can be calculated from the first high frequency QRS complex data to form the first lead positive indicator. A set function can be used to process the first high frequency QRS complex data to calculate the amplitude drop relative value and the amplitude absolute value between the two reference points with the most rapid drop in RMS voltage in the high frequency QRS waveform within a period of time before exercise, during exercise, and within a period of time after exercise. The amplitude absolute value is obtained by subtracting the RMS voltage of the two reference points, and the amplitude drop relative value is calculated as a percentage based on the amplitude absolute value and the RMS voltage value of the reference point with the highest RMS voltage. The amplitude drop relative value and the amplitude absolute value are mainly used to evaluate the blood flow change of the human heart.

[0105] The resting ECG data can be analyzed to obtain high frequency QRS complex data (for ease of distinction, hereinafter referred to as second high frequency QRS complex data). The specific operation of analyzing the resting ECG data to obtain the corresponding high frequency QRS complex data is described above and will not be repeated here. Thus, the envelope curve formed when the second high frequency QRS complex data is visualized is the aforementioned high frequency QRS envelope curve. The second high frequency QRS complex data can be used to obtain a lead positive indicator corresponding to each ECG lead (for ease of distinction, hereinafter referred to as a second lead positive indicator). For example, in an example, the second high frequency QRS envelope curve based on the second high frequency QRS complex data corresponding to each ECG lead can be analyzed to obtain a first total area of each amplitude reduction region on the high frequency QRS envelope curve and a second total area below the high frequency QRS envelope curve, and the ratio of the first total area to the second total area is taken as the high frequency morphology index corresponding to the resting lead. The corresponding lead positive indicator is determined according to the age of the subject and the high frequency morphology index corresponding to each resting lead. For example, if the age is ≥ 50 years old, the second lead positive indicator is determined to be positive when the high frequency morphology index is ≥ 8%, and if the age is < 50 years old, the second lead positive indicator is determined to be positive when the high frequency morphology index is ≥ 15%.

[0106] The ST segment of the conventional ECG data can be analyzed to determine the positive indicator of the corresponding ECG lead (for ease of distinction, hereinafter referred to as a third lead positive indicator). For example, if the ST segment in the conventional ECG data has a horizontal depression or a downward slope depression ≥ 0.1 mV (millivolt) and a duration > 2 min (minutes), the corresponding third lead positive indicator is determined to be positive.

[0107] If the first lead positive indicator, the second lead positive indicator, and the third lead positive indicator all indicate positive, it can be determined that there is the highest possibility (e.g., a probability of 100%) of myocardial ischemia. In addition, if any one of the first lead positive indicator, the second lead positive indicator, and the third lead positive indicator indicates positive, it indicates that there is a possibility of myocardial ischemia, i.e., a lower possibility (e.g., a probability less than 30%). If any two of the first lead positive indicator, the second lead positive indicator, and the third lead positive indicator indicate positive, it indicates that there is a high suspicion of myocardial ischemia, i.e., a higher possibility (e.g., a probability more than 70%).

[0108] In step S103, the possibility of the existence of a coronary artery lesion is determined according to the exercise electrocardiogram data.

[0109] The coronary artery lesion may, for example, be a coronary artery macrovascular lesion, and specifically can include coronary stenosis and coronary sclerosis. When the possibility of the existence of coronary stenosis or coronary sclerosis is determined according to the exercise electrocardiogram data, it can be considered that there is a possibility of the existence of a coronary artery lesion.

[0110] Specifically, in an embodiment of the present application, the method disclosed in the Chinese patent application with the publication number CN114742114A, previously applied for by the applicant of the present application, can be used to determine the possibility of the existence of coronary stenosis. The application discloses a high-frequency QRS waveform curve analysis method, which includes obtaining a high-frequency QRS waveform curve (first high-frequency QRS waveform curve) corresponding to the exercise electrocardiogram data; selecting the first high-frequency QRS waveform curve in a first time period as a first reference waveform curve; selecting a point with the minimum root mean square voltage on the first reference waveform curve as a first reference point, and a point with the maximum root mean square voltage earlier in time than the first reference point as a second reference point; determining a first amplitude drop relative value according to the root mean square voltages of the first reference point and the second reference point; if the first amplitude drop relative value is greater than or equal to a first preset threshold value, it is determined that the first high-frequency QRS waveform curve satisfies a preset condition, i.e., the high-frequency QRS waveform curve has a waveform feature indicating coronary stenosis, and thus the possibility of the existence of coronary stenosis can be judged or determined.

[0111] The high-frequency QRS waveform curve analysis method disclosed in the above embodiments mainly focuses on whether the first amplitude drop relative value is greater than or equal to the first preset threshold in the QRS waveform curve segment corresponding to the first 3 minutes in the motion in the high-frequency QRS waveform curve, that is, mainly judges whether there is a steep drop waveform in the high-frequency QRS waveform curve in the first 3 minutes in the motion, and if there is, it is determined that there is coronary stenosis. However, the actual situation may be that the position (time point) where the first amplitude drop relative value greater than or equal to the first preset threshold appears in the high-frequency QRS waveform curve is located after 3 minutes, that is, the high-frequency QRS waveform curve appears a steep drop waveform after 3 minutes in the motion (as shown in Figure 2A The method of the above embodiments may ignore this situation and fail to determine the existence of coronary stenosis. In addition, this method also ignores the identification of the possibility of coronary arteriosclerosis.

[0112] In view of this, the inventors of the present application have further studied and proposed an innovative method for determining the possibility of coronary artery disease, which can improve the accuracy of determining the possibility of coronary stenosis and also determine the possibility of coronary arteriosclerosis compared with the method provided in the above embodiments disclosed before.

[0113] In the preferred embodiments of the present application, the general inventive concept of determining whether there is a coronary artery disease can include obtaining first high-frequency QRS wave group data from exercise electrocardiogram data. The first high-frequency QRS wave group data can include a set of data points arranged in time sequence (i.e., chronological order), where each data point corresponds to a time point (time stamp) and an amplitude value (or intensity value), which can be, for example, an RMS voltage value. The first high-frequency QRS waveform curve can be essentially a curve obtained by connecting these data points in time sequence, or a curve obtained by removing minor fluctuations through curve smoothing processing.

[0114] A point sequence for each electrocardiogram lead can be obtained from the first high-frequency QRS wave group data, which can include a plurality of data points sampled or selected from a plurality of data point sets of the first high-frequency QRS wave group data at a set sampling period (sampling interval) or sampling frequency. Each sampling point in the point sequence is traversed to determine whether the first high-frequency QRS wave group data has a waveform feature indicative of a coronary artery disease, such as a steep drop wave, and if so, it can be determined that there is a possibility of a coronary artery disease.

[0115] Figure 2B An example flowchart of a method for determining the possibility of a coronary artery disease according to an embodiment of the present application is schematically shown. As Figure 2B shown, specifically, in this embodiment, step S103 can include the following steps.

[0116] In step S1031, a point sequence (may be referred to as a first high-frequency QRS time-intensity point sequence) corresponding to each electrocardiogram lead is obtained according to the first high-frequency QRS wave group data corresponding to each electrocardiogram lead. The point sequence can include a plurality of sampling points arranged in time sequence.

[0117] In step S1032, for the point sequence corresponding to any electrocardiogram lead, any sampling point in the plurality of sampling points in the point sequence is taken as a reference sampling point, and the sampling points in a preset time interval from the reference sampling point are traversed to determine whether there is a first sampling point and a second sampling point that meet a preset condition in the preset time interval. The sampling points in the preset time interval include the reference sampling point. In an example, the preset condition can be that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold (hereinafter referred to as the first preset condition). In another example, the preset condition can be that the amplitude drop relative value between the first sampling point and the second sampling point reaches the first preset threshold or the voltage drop degree between the first sampling point and the second sampling point reaches the second preset threshold (hereinafter referred to as the second preset condition).

[0118] The amplitude drop relative value may, for example, be: (RMS voltage of the first sampling point - RMS voltage of the second sampling point) / RMS voltage of the first sampling point. In an example, the voltage drop degree can be a voltage drop absolute value, which can be the RMS voltage of the first sampling point minus (-) the RMS voltage of the second sampling point. In an alternative example, the voltage drop degree can be the ratio of the voltage drop absolute value to the maximum voltage value, i.e. (RMS voltage of the first sampling point - RMS voltage of the second sampling point) / maximum voltage value, which can be understood as the maximum power, which can be used to reflect the maximum heart pumping function of the subject, i.e. the peak value of the first high-frequency QRS wave group data (or the first high-frequency QRS waveform curve) will not exceed the maximum voltage value. The maximum voltage value can be determined according to the maximum value of the root mean square voltage in the high-frequency QRS wave group data corresponding to all electrocardiogram leads (for example, the maximum RMS voltage is taken as an even integer). In this example, the second preset threshold can be in the form of a percentage as the first preset threshold.

[0119] In step S1033, in the case where the preset condition is met, it is determined that there is a possibility of coronary artery lesion. If the first sampling point and the second sampling point that meet the preset condition appear during the traversal, it indicates that the first high-frequency QRS wave group data has a waveform feature indicating coronary artery lesion, such as a steep drop wave, in which case it can be determined that there is a possibility of coronary artery lesion.

[0120] In step S1034, if the preset condition is not met, the step S1032 is repeated after replacing the reference sampling point. In this step, if all the sampling points in the preset time interval are traversed and no first sampling point and second sampling point meeting the preset condition are found, the current traversal loop ends, and the reference sampling point can be replaced, and the next traversal loop is entered. For example, the adjacent sampling point of the current reference sampling point is taken as the replaced (new) reference sampling point, and the sampling points in the preset time interval from the replaced reference sampling point are traversed based on the replaced reference sampling point to determine whether there are first sampling point and second sampling point meeting the preset condition. If there are first sampling point and second sampling point meeting the preset condition, it indicates that the first high-frequency QRS complex data has the waveform feature indicating the coronary artery lesion, and the traversal ends. If the current traversal loop ends and no first sampling point and second sampling point meeting the preset condition are found, the next traversal loop is entered, and the process is repeated until the first sampling point and second sampling point meeting the preset condition are found, or all the sampling points in the point sequence to be traversed are traversed. If all the sampling points to be traversed are traversed and the preset condition is still not met, it is determined that there is no waveform feature indicating the coronary artery lesion.

[0121] The traversal manner or rule of the sampling points can be various. Figure 3 An example flowchart of a method for determining the presence of a coronary artery lesion according to an embodiment of the present application is schematically shown. As shown in the figure, the method comprises the following steps. Figure 3As shown, in this embodiment, for the first high-frequency QRS complex data corresponding to each electrocardiogram lead, the first sampling point in the point sequence in time sequence can be taken as a reference sampling point (which can be referred to as an initial reference point), the next sampling point of the initial reference point in time sequence is taken as a traversal starting point to perform traversal in time sequence, and the amplitude drop relative value and the voltage drop degree between the initial reference point (the first sampling point in this embodiment) and the current traversal sampling point (the second sampling point in this embodiment) are calculated. It is judged whether the initial reference point and the current traversal sampling point satisfy a preset condition. If the two sampling points satisfy the preset condition, it can be determined that the first high-frequency QRS complex data corresponding to the lead exists a waveform feature indicating coronary artery lesion, for example, a steep drop wave, so that it can be determined that there is a possibility of coronary artery lesion, and the traversal is stopped. If the preset condition is not satisfied, it is judged whether the time interval between the current traversal sampling point and the initial reference point exceeds a preset time length, and whether the current traversal sampling point is the last sampling point in the point sequence. If the time interval between the current traversal sampling point and the initial reference point does not exceed the preset time length and the current traversal sampling point is not the last sampling point, the next sampling point of the current traversal sampling point is switched to, the next sampling point is taken as a new current traversal sampling point, and the same operation is performed on the new current traversal sampling point. Specifically, the amplitude drop relative value and the voltage drop degree between the initial reference point and the new current traversal sampling point are calculated. It is judged whether the initial reference point and the new current traversal sampling point satisfy the preset condition. If the preset condition is satisfied, it can be determined that the first high-frequency QRS complex data corresponding to the lead exists a steep drop wave, that is, it can be determined that there is a possibility of coronary artery lesion, and the traversal is stopped. If the preset condition is not satisfied, the above operation can be repeated to traverse the subsequent sampling points, and if the preset condition is still not satisfied, the traversal can be continued until the time interval between the current traversal sampling point and the initial reference point reaches the preset time length (the preset time interval in this embodiment) or the current traversal sampling point is the last sampling point in the point sequence. If the time interval between the current traversal sampling point and the initial reference point reaches the preset time length or the current traversal sampling point is the last sampling point in the point sequence, the next sampling point of the initial reference point is switched to, the next sampling point is taken as a reference sampling point, and the next traversal loop is performed in the same manner as the above traversal operation. If the preset condition is still not satisfied, the traversal loop can be continuously performed until the preset condition is satisfied or until all the sampling points except the last sampling point in the point sequence are taken as reference sampling points to complete the traversal. That is, when the previous sampling point of the last sampling point is taken as a reference sampling point to complete the traversal, and the preset condition is still not satisfied, it is indicated that there is no waveform feature indicating coronary artery lesion, for example, a steep drop wave, for the first high-frequency QRS complex data corresponding to the electrocardiogram lead.

[0122] In this embodiment, the preset time length can be, for example, 1 minute to 3 minutes. If the preset condition is that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value (the first preset condition), the value range of the first preset threshold value can be 30% to 40%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the second preset threshold value can be 1 uV to 2 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value can be 10% to 20%.

[0123] If the preset condition is that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value, or, the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value (the first preset condition), the value range of the first preset threshold value can be 50% to 60%. In this case, the value of the second preset threshold value is associated with the preset time length. In general, the value of the second preset threshold value is positively correlated with the preset time length, that is, the greater the preset time length, the greater the second preset threshold value. For example, taking the value range of the preset time length as 1 minute to 3 minutes, if the preset time length is 1 minute, in the case where the voltage drop degree is the voltage drop absolute value, the value range of the second preset threshold value is 1 uV to 2 uV, and in the case where the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 10% to 20%. If the preset time length is 2 minutes, in the case where the voltage drop degree is the voltage drop absolute value, the value range of the second preset threshold value is 2 uV to 3 uV, and in the case where the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 20% to 30%. If the preset time length is 3 minutes, in the case where the voltage drop degree is the voltage drop absolute value, the value range of the second preset threshold value is 3 uV to 4 uV, and in the case where the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 30% to 40%. Those skilled in the art can understand that the specific values of the above preset threshold values are exemplary, and are intended to illustrate that the second preset threshold value is associated with the preset time length. For a smaller preset time length, if a larger second preset threshold value is used, the steep drop wave will be misjudged as a non-steep drop wave. For a larger preset time length, if a smaller second preset threshold value is used, the non-steep drop wave will be misjudged as a steep drop wave, thereby reducing the recognition accuracy. The specific values can be adjusted according to actual applications.

[0124] The first high-frequency QRS complex data corresponding to all electrocardiogram leads can be processed as above. If the first high-frequency QRS complex data corresponding to at least one electrocardiogram lead has waveform characteristics indicative of coronary artery lesion, such as a steep drop wave, it can be determined that there is a possibility of coronary artery lesion. Conversely, if the first high-frequency QRS complex data corresponding to all electrocardiogram leads does not have waveform characteristics indicative of coronary artery lesion, it can be determined that there is no coronary artery lesion.

[0125] Although the above embodiment processes the multiple sampling points in the point sequence in the time sequence from early to late, those skilled in the art can understand that the above method is also applicable to the multiple sampling points in the time sequence from late to early, in which case, the amplitude drop relative value and the voltage drop degree between the current traversal sampling point and the reference sampling point are calculated (in this case, the current traversal sampling point is earlier in time than the reference sampling point).

[0126] Figure 4 An example flowchart of a method for determining the possibility of coronary artery lesion according to another embodiment of the present application is schematically shown. As shown in FIG. 6, the method comprises the following steps. Figure 4As shown, in this embodiment, a window function can be preset, and the window of the window function includes a start point and an end point, and the distance between the start point and the end point is a window length (in this embodiment, the window length is the preset time interval described above). The window length can be, for example, 1 minute to 3 minutes. For the first high-frequency QRS complex data corresponding to each electrocardiogram lead, the point sequence of the first high-frequency QRS complex data is traversed using the window function. Specifically, the earliest first sampling point in the point sequence in time sequence can be taken as the start point of the window, and the sampling point with the maximum RMS voltage among the multiple sampling points in the point sequence included in the window (the first sampling point in this embodiment) and the sampling point with the minimum RMS voltage among the sampling points after the first sampling point in time sequence (the second sampling point in this embodiment) can be determined. The amplitude drop relative value between the first sampling point and the second sampling point and the voltage drop degree are calculated. It is determined whether the first sampling point and the second sampling point satisfy a preset condition. If the two sampling points satisfy the preset condition, it can be determined that the first high-frequency QRS complex data corresponding to the lead exists a sharp drop wave, that is, it can be determined that there is a possibility of coronary artery lesion, and the traversal is stopped. If the preset condition is not satisfied, the window of the window function is slid backward by a preset step, and the above operation is repeatedly performed. If there is still no situation satisfying the preset condition, the window is slid backward by the preset step again until the preset condition is satisfied or until the end point of the window reaches or exceeds the last sampling point in the point sequence. That is, when the end point of the window reaches or exceeds the last sampling point after the window is slid by the preset step, and there is still no situation satisfying the preset condition, it means that there is no waveform feature indicating coronary artery lesion, such as a sharp drop wave, for the first high-frequency QRS complex data of the electrocardiogram lead. In this embodiment, the preset step can be, for example, an integer multiple of the time interval (for example, a sampling period) between adjacent two sampling points, and is preferably the time interval between adjacent two sampling points, for example, 10 seconds. Of course, those skilled in the art can understand that it is also feasible to set the preset step to other values.

[0127] In this embodiment, if the preset condition is that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold (the first preset condition), the value range of the first preset threshold can be 30% to 40%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the second preset threshold can be 1 uV to 2 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold can be 10% to 20%.

[0128] If the preset condition is that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value, or the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value (the first preset condition), the value range of the first preset threshold value can be 50% to 60%. In this case, the value of the second preset threshold value is associated with the window length. In general, the value of the second preset threshold value is positively correlated with the window length, that is, the longer the window length, the larger the second preset threshold value. For example, taking the value range of the window length as 1 minute to 3 minutes as an example, if the window length is 1 minute, the value range of the second preset threshold value is 1 uV to 2 uV when the voltage drop degree is the voltage drop absolute value, and the value range of the second preset threshold value is 10% to 20% when the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value. If the window length is 2 minutes, the value range of the second preset threshold value is 2 uV to 3 uV when the voltage drop degree is the voltage drop absolute value, and the value range of the second preset threshold value is 20% to 30% when the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value. If the window length is 3 minutes, the value range of the second preset threshold value is 3 uV to 4 uV when the voltage drop degree is the voltage drop absolute value, and the value range of the second preset threshold value is 30% to 40% when the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value. Those skilled in the art can understand that the specific values of the above preset threshold values are exemplary, and the purpose is to explain that the second preset threshold value is associated with the window length. For a shorter window length, if a larger second preset threshold value is used, the steep drop wave will be misjudged as a non-steep drop wave. For a longer window length, if a smaller second preset threshold value is used, the non-steep drop wave will be misjudged as a steep drop wave, thereby reducing the recognition accuracy. The specific value can be adjusted according to actual application.

[0129] The first high-frequency QRS complex data corresponding to all electrocardiogram leads can be processed as above. If the first high-frequency QRS complex data corresponding to at least one electrocardiogram lead has a waveform feature indicating coronary artery lesion, such as a steep drop wave, it can be determined that there is a possibility of coronary artery lesion. Conversely, if the first high-frequency QRS complex data corresponding to all electrocardiogram leads does not have a waveform feature indicating coronary artery lesion, it can be determined that there is no coronary artery lesion.

[0130] Although the window function used in the above embodiment starts from the earliest sampling point in the point sequence, those skilled in the art can understand that it can also start from the latest sampling point in the point sequence.

[0131] Figure 5An example flowchart illustrating a method for determining the likelihood of coronary artery disease according to another embodiment of this application is shown. Figure 5 As shown, in this embodiment, a window function can be preset. The window of this window function includes a start point and an end point, and the distance between the start point and the end point is the window length (in this embodiment, the window length is the preset time interval mentioned above). The window length can range from, for example, 1 minute to 3 minutes. For the first high-frequency QRS complex data corresponding to each ECG lead, the start point of the window can be aligned with the earliest sampling point in the time sequence (i.e., these two points coincide in time). Using this first sampling point as a reference point, each sampling point within the window is traversed sequentially. Specifically, the relative value of amplitude decrease and the degree of voltage decrease between the reference point and the currently traversed sampling point within the window can be calculated. It is determined whether the reference point and the currently traversed sampling point meet preset conditions. If these two sampling points meet the preset conditions, it can be determined that the first high-frequency QRS complex data corresponding to the lead has waveform characteristics indicating coronary artery disease, such as a steep drop wave, and the possibility of coronary artery disease can be determined, and the traversal stops. If the preset conditions are not met, the next sampling point within the window is traversed. If the reference point and the currently traversed sampling point meet the preset conditions, the traversal stops; otherwise, it continues until the last sampling point within the traversal window. If the preset conditions are still not met after all sampling points within the window have been traversed, the traversal loop ends, and the window of the window function can be slid backward (time-wise backward) by a preset step size, and the above operation is repeated (the next traversal loop). If the preset conditions are still not met after the traversal loop ends, the window can be slid backward by the preset step size again until the preset conditions are met, or until the end of the window reaches or exceeds the last sampling point in the point sequence. That is, if the preset conditions are still not met after the window has slid by the preset step size and its end reaches or exceeds the last sampling point, it indicates that there are no waveform characteristics indicating coronary artery lesions, such as steep descent waves, for the first high-frequency QRS complex data of that ECG lead. In this embodiment, the preset step size can be, for example, an integer multiple of the time interval (e.g., sampling period) between two adjacent sampling points, preferably, the time interval between two adjacent sampling points, such as 10 seconds. Of course, those skilled in the field can understand that setting the preset step size to other values ​​is also feasible.

[0132] In this embodiment, if the preset condition is that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value (the first preset condition), the value range of the first preset threshold value can be 30% to 40%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the second preset threshold value can be 1 uV to 2 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value can be 10% to 20%.

[0133] If the preset condition is that the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value, or, the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value (the first preset condition), the value range of the first preset threshold value can be 50% to 60%. In this case, the value of the second preset threshold value is associated with the window length. In general, the value of the second preset threshold value is positively correlated with the window length, that is, the longer the window length, the larger the second preset threshold value. For example, taking the value range of the window length as 1 minute to 3 minutes as an example, if the window length is 1 minute, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the second preset threshold value is 1 uV to 2 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 10% to 20%. If the window length is 2 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the second preset threshold value is 2 uV to 3 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 20% to 30%. If the window length is 3 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the second preset threshold value is 3 uV to 4 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 30% to 40%. Those skilled in the art can understand that the specific values of the above preset threshold values are exemplary, and the purpose is to explain that the second preset threshold value is associated with the window length. For a shorter window length, if a larger second preset threshold value is used, the steep drop wave will be misjudged as a non-steep drop wave. For a longer window length, if a smaller second preset threshold value is used, the non-steep drop wave will be misjudged as a steep drop wave, thereby reducing the recognition accuracy. The specific values can be adjusted according to actual application.

[0134] The first high-frequency QRS complex data corresponding to all electrocardiogram leads can be processed as above. If the first high-frequency QRS complex data corresponding to at least one electrocardiogram lead has a waveform feature indicative of coronary artery lesion, such as a steep drop wave, it can be determined that there is a possibility of coronary artery lesion. Conversely, if the first high-frequency QRS complex data corresponding to all electrocardiogram leads does not have a waveform feature indicative of coronary artery lesion, it can be determined that there is no coronary artery lesion.

[0135] Although the window function adopted in the above embodiment starts from the earliest sampling point in time sequence in the point sequence, it can be understood by those skilled in the art that it can also start from the latest sampling point in time sequence in the point sequence.

[0136] In the above embodiment, if it is determined in the traversal process that there are two sampling points that satisfy the preset condition, but the first high-frequency QRS complex data (or the first high-frequency QRS waveform curve) can have a significant fluctuation amplitude between the two sampling points, the cause of this situation can be coronary spasm or other factors such as microcirculation disorder. In order to exclude this interference factor, in the preferred embodiment of the present application, after the first sampling point and the second sampling point that satisfy the preset condition are determined, the first-order difference of the RMS voltage / amplitude / value / intensity of the sampling points between the first sampling point and the second sampling point is calculated to determine a third sampling point corresponding to a minimum value (minimum point) and at least one fourth sampling point corresponding to a maximum value (maximum point). Here, the minimum point refers to a point where the first-order difference result changes from negative to positive, and the maximum point refers to a point where the first-order difference result changes from positive to negative. For any one of the at least one fourth sampling point, the absolute value of the RMS voltage difference between the fourth sampling point and the third sampling point (referred to as the first RMS voltage difference) is calculated, and the RMS voltage difference between the first sampling point and the second sampling point (referred to as the second RMS voltage difference) is calculated. If the ratio of the first RMS voltage difference to the second RMS voltage difference is greater than or equal to a third preset threshold, it is determined that there is a significant fluctuation amplitude, it is determined that the two sampling points (the first sampling point and the second sampling point) of the first high-frequency QRS complex data do not have a waveform feature indicative of coronary artery lesion (it is considered that the first sampling point and the second sampling point do not satisfy the preset condition), and the traversal continues according to the case where the preset condition is not satisfied in the above embodiment. If the ratio is less than the third preset threshold, the significant fluctuation amplitude can be excluded, the first sampling point and the second sampling point satisfy the preset condition, and it can be determined that there is a steep drop wave. In the case where the number of fourth sampling points is multiple, it is preferred that the ratio obtained for each fourth sampling point is less than the third preset threshold, and the significant fluctuation amplitude can be excluded. In an example, the third preset threshold can be, for example, 40%, 50%, etc. Using this way to exclude the interference of significant fluctuation amplitude on the judgment of whether there is a steep drop wave can further improve the evaluation specificity.

[0137] Compared with the scheme for judging whether there is coronary stenosis disclosed in the prior patent application, the coronary artery lesion judgment method provided in the above embodiments can significantly improve the sensitivity without significantly reducing or even not reducing the specificity.

[0138] In further embodiments of the present application, the type of coronary artery lesion, i.e., coronary stenosis or coronary sclerosis, can be further determined. In this embodiment, step S103 can further include the following steps.

[0139] In step S1035, a first lead positive indicator corresponding to the first high-frequency QRS complex data with the waveform feature indicative of coronary artery lesion is determined. The determination method of the first lead positive indicator is as described above, which will not be repeated here.

[0140] In step S1035, in the case where the first lead positive indicator indicates positive, the type of coronary artery lesion is determined to be coronary stenosis. Where the critical case (i.e., the state between positive and negative) is considered, the case where the first lead positive indicator indicates positive includes the case where the first lead positive indicator indicates positive or critical.

[0141] In step S1036, in the case where the first lead positive indicator indicates negative, the type of coronary artery lesion is determined to be coronary sclerosis.

[0142] In step S104, the possibility of the presence of myocardial bridge is determined according to the exercise electrocardiogram data.

[0143] Figure 6 An example flowchart of a method for determining the possibility of the presence of myocardial bridge according to an embodiment of the present application is schematically shown. As shown, specifically, in this embodiment, step S104 can include the following steps. Figure 6

[0144] In step S1041, candidate first high-frequency QRS complex data without the waveform feature indicative of coronary artery lesion, such as the sharp drop wave, is determined from the first high-frequency QRS complex data corresponding to each electrocardiogram lead. In analyzing the first high-frequency QRS complex data to determine whether there is myocardial bridge, the first high-frequency QRS complex data with the waveform feature indicative of coronary artery lesion can be excluded first. That is, if the first high-frequency QRS complex data has the waveform feature indicative of coronary artery lesion, such as the sharp drop wave, the first high-frequency QRS complex data will not be used to determine whether there is myocardial bridge. The determination method of whether the first high-frequency QRS complex data has the waveform feature indicative of coronary artery lesion is as described in the embodiments of step S103, which will not be repeated here.

[0145] ​In step S1042, a point sequence (may be referred to as a first high-frequency QRS time-intensity point sequence) of each candidate first high-frequency QRS complex data is obtained, which can include a plurality of sampling points arranged in time sequence. The plurality of sampling points can be a plurality of data points sampled or selected from a plurality of data point sets of the candidate first high-frequency QRS complex data at a set sampling period (sampling interval) or sampling frequency.

[0146] In step S1043, for the point sequence of any candidate first high-frequency QRS complex data, any sampling point in the plurality of sampling points of the point sequence is taken as a reference sampling point, and the sampling points within a preset time interval from the reference sampling point are traversed to determine whether there are an earlier first sampling point and a later second sampling point in the preset time interval that satisfy a preset condition.

[0147] In an example case, the preset condition can be that the amplitude drop relative value between the first sampling point and the second sampling point reaches a fourth preset threshold value, and the voltage drop degree between the first sampling point and the second sampling point reaches a fifth preset threshold value, and the root mean square voltage of each sampling point between the first sampling point and the second sampling point is in a continuous downward trend (hereinafter referred to as a third preset condition). In another example case, the preset condition can be that the amplitude drop relative value between the first sampling point and the second sampling point reaches the fourth preset threshold value or the voltage drop degree between the first sampling point and the second sampling point reaches the fifth preset threshold value, and the root mean square voltage of each sampling point between the first sampling point and the second sampling point is in a continuous downward trend (hereinafter referred to as a fourth preset condition).

[0148] The amplitude drop relative value may, for example, be: (first sampling point root mean square voltage (RMS voltage) - second sampling point RMS voltage) / first sampling point RMS voltage. In an example, the voltage drop degree can be a voltage drop absolute value, which can be the first sampling point root mean square voltage (RMS voltage) minus (-) the second sampling point RMS voltage. In an alternative example, the voltage drop degree can be the ratio of the voltage drop absolute value to the maximum voltage value, i.e. (first sampling point root mean square voltage (RMS voltage) - second sampling point RMS voltage) / maximum voltage value, which can be understood as the maximum power, which can be used to reflect the maximum heart pumping function of the subject, i.e. the peak value of the candidate first high-frequency QRS complex data will not exceed the maximum voltage value, and the maximum voltage value is determined according to the maximum value of the root mean square voltage in the high-frequency QRS complex data of all electrocardiogram leads. In this example, the fifth preset threshold value can be in the form of a percentage as the fourth preset threshold value.

[0149] In step S1044, if the preset condition is met, it is determined that there is a possibility of myocardial bridge. During the traversal, if the first sampling point and the second sampling point meeting the preset condition appear, it means that the candidate first high-frequency QRS complex data has waveform characteristics indicating myocardial bridge, such as a slow descent wave. In this case, it can be determined that there is a possibility of myocardial bridge. Here, the slow descent wave has a more gentle downward trend than the steep descent wave in the displayed waveform.

[0150] In step S1045, if the preset condition is not met, the reference sampling point is replaced and step S1043 is repeated. In this step, if all the sampling points in the preset time interval are traversed and no first sampling point and second sampling point meeting the preset condition appear, the current traversal loop ends, and the reference sampling point is replaced. The next traversal loop is entered, for example, the adjacent sampling point of the current reference sampling point is taken as the replaced (new) reference sampling point. The sampling points in the preset time interval from the replaced reference sampling point are traversed based on the replaced reference sampling point as the reference to determine whether there are first sampling points and second sampling points meeting the preset condition. If there are first sampling points and second sampling points meeting the preset condition, it means that the candidate first high-frequency QRS complex data line has waveform characteristics indicating myocardial bridge, and the traversal ends. If the current traversal loop ends and no first sampling point and second sampling point meeting the preset condition appear, the next traversal loop is entered, and this process continues until a first sampling point and a second sampling point meeting the preset condition appear, or all the sampling points to be traversed in the point sequence are traversed. If all the sampling points to be traversed are traversed and the preset condition is still not met, it is determined that there is no waveform characteristic indicating myocardial bridge.

[0151] There can be various ways or rules for traversing the sampling points. Figure 7 An example flowchart of a method for determining the possibility of myocardial bridge according to an embodiment of the present application is schematically shown. As shown in FIG. 6, the method comprises the following steps. Figure 7As shown, in this embodiment, for each candidate first high-frequency QRS complex data (i.e., the complex data does not have waveform features indicative of coronary artery lesions), the first sampling point in time sequence in its point sequence can be taken as a reference sampling point (which can be referred to as an initial reference point), the next sampling point in time sequence of the initial reference point can be taken as a traversal starting point to perform traversal in time sequence, and the amplitude drop relative value and voltage drop degree between the initial reference point (the first sampling point in this embodiment) and the current traversal sampling point (the second sampling point in this embodiment) can be calculated. It is determined whether the initial reference point and the current traversal sampling point satisfy a preset condition. If the two sampling points satisfy the preset condition, it can be determined that the candidate first high-frequency QRS complex data has waveform features indicative of myocardial bridge, such as a slow descent wave, thereby determining the possibility of myocardial bridge, and stopping traversal. If the preset condition is not satisfied, it is determined whether the time interval between the current traversal sampling point and the initial reference point exceeds a preset time length, and whether the current traversal sampling point is the last sampling point in the point sequence. If the time interval between the current traversal sampling point and the initial reference point does not exceed the preset time length and the current traversal sampling point is not the last sampling point, the next sampling point of the current traversal sampling point is switched to, the next sampling point is taken as a new current traversal sampling point, and the same operation is performed on the new current traversal sampling point. Specifically, the amplitude drop relative value and voltage drop degree between the initial reference point and the new current traversal sampling point are calculated. It is determined whether the initial reference point and the new current traversal sampling point satisfy the preset condition. If the preset condition is satisfied, it can be determined that the candidate first high-frequency QRS complex data has a slow descent wave, i.e., the possibility of myocardial bridge is determined, and the traversal is stopped. If the preset condition is not satisfied, the above operation can be repeated to traverse the subsequent sampling points, and if the preset condition is still not satisfied, the traversal can be continued until the time interval between the current traversal sampling point and the initial reference point reaches the preset time length (the preset time interval in this embodiment) or the current traversal sampling point is the last sampling point in the point sequence. If the time interval between the current traversal sampling point and the initial reference point reaches the preset time length or the current traversal sampling point is the last sampling point in the point sequence, the next sampling point of the initial reference point is switched to, the next sampling point is taken as a reference sampling point, and the next traversal cycle is performed in the same manner as the above traversal operation. If the preset condition is still not satisfied, the traversal cycle can be continuously performed until the preset condition is satisfied or until all sampling points except the last sampling point in the point sequence are taken as reference sampling points to complete traversal. That is, when the previous sampling point of the last sampling point is taken as a reference sampling point to complete traversal, and the preset condition is still not satisfied, it is indicated that the candidate first high-frequency QRS complex data does not have waveform features indicative of myocardial bridge, such as a slow descent wave.

[0152] In this embodiment, the preset time length can be, for example, 4 to 6 minutes. If the preset condition includes that the amplitude drop relative value between the first sampling point and the second sampling point reaches the fourth preset threshold value and the voltage drop degree between the first sampling point and the second sampling point reaches the fifth preset threshold value (the third preset condition), the value range of the fourth preset threshold value can be 30% to 40%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold value can be 1 uV to 2 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold value can be 10% to 20%.

[0153] If the preset condition includes that the amplitude drop relative value between the first sampling point and the second sampling point reaches the fourth preset threshold value, or, the voltage drop degree between the first sampling point and the second sampling point reaches the fifth preset threshold value (the fourth preset condition), the value range of the fourth preset threshold value can be 50% to 60%. In this case, the value of the fifth preset threshold value is associated with the preset time length. In general, the value of the fifth preset threshold value is positively correlated with the preset time length, that is, the greater the preset time length, the greater the fifth preset threshold value. For example, taking the value range of the preset time length as 4 to 6 minutes, if the preset time length is 4 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the fifth preset threshold value is 1 uV to 2 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 10% to 20%. If the preset time length is 5 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the fifth preset threshold value is 2 uV to 3 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 20% to 30%. If the preset time length is 6 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the fifth preset threshold value is 3 uV to 4 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 30% to 40%. Those skilled in the art can understand that the specific values of the above preset threshold values are exemplary, and are intended to illustrate that the fifth preset threshold value is associated with the preset time length. For a smaller preset time length, if a larger fifth preset threshold value is used, the slow decline wave will be misjudged as a non-slow decline wave. For a larger preset time length, if a smaller fifth preset threshold value is used, the non-slow decline wave (such as a flat wave) will be misjudged as a slow decline wave, thereby reducing the recognition accuracy. Similarly, the preset time length is positively correlated with the fourth preset threshold value, so as to further improve the recognition accuracy. The specific values can be adjusted according to actual applications.

[0154] All candidate first high frequency QRS complex data can be processed as above, if at least one candidate first high frequency QRS complex data has waveform features indicative of myocardial bridge, such as a slow descent wave, it can be determined that there is a possibility of myocardial bridge. Conversely, if all candidate first high frequency QRS complex data do not have waveform features indicative of myocardial bridge, it can be determined that there is no myocardial bridge.

[0155] Although the above embodiment is to process multiple sampling points in the point sequence in time sequence from early to late, those skilled in the art can understand that the above method is also applicable to multiple sampling points in time sequence from late to early, in which case, the amplitude drop relative value and voltage drop degree between the current traversal sampling point and the reference sampling point are calculated (in this case, the current traversal sampling point is earlier in time than the reference sampling point).

[0156] Figure 8 An example flow chart of a method for determining the possibility of myocardial bridge according to another embodiment of the present application is schematically shown. As shown in the figure, the method comprises the following steps: Figure 8As shown, in this embodiment, a window function can be preset, the window of the window function including a start point and an end point, the distance between the start point and the end point being a window length (in this embodiment, the window length is the preset time interval described above). The window length can be, for example, 4 minutes to 6 minutes. For each candidate first high-frequency QRS complex data, the point sequence of the candidate first high-frequency QRS complex data is traversed using the window function. Specifically, the earliest first sampling point in the point sequence in time sequence can be taken as the start point of the window, and the sampling point with the maximum RMS voltage among the multiple sampling points in the point sequence included in the window (the first sampling point in this embodiment) and the sampling point with the minimum RMS voltage among the sampling points after the first sampling point in time sequence (the second sampling point in this embodiment) can be determined. The amplitude drop relative value between the first sampling point and the second sampling point and the voltage drop degree are calculated. It is determined whether the first sampling point and the second sampling point satisfy a preset condition. If the two sampling points satisfy the preset condition, it can be determined that the candidate first high-frequency QRS complex data has a slow descent wave, i.e., it can be determined that there is a possibility of myocardial bridge, and the traversal is stopped. If the preset condition is not satisfied, the window of the window function can be slid backward by a preset step, and the above operation can be repeatedly performed. If the situation where the preset condition is still not satisfied still does not appear, the window can be slid backward by the preset step again until the preset condition is satisfied, or until the end point of the window reaches or exceeds the last sampling point in the point sequence. That is, when the end point of the window reaches or exceeds the last sampling point after the window is slid by the preset step, and the situation where the preset condition is satisfied still does not appear, it means that there is no wave form feature indicating myocardial bridge, such as a slow descent wave, for the candidate first high-frequency QRS complex data. In this embodiment, the preset step can be, for example, an integer multiple of the time interval (e.g., a sampling period) between adjacent two sampling points, and is preferably the time interval between adjacent two sampling points, e.g., 10 seconds. Of course, those skilled in the art can understand that it is also feasible to set the preset step to other values.

[0157] In this embodiment, if the preset condition includes that the amplitude drop relative value between the first sampling point and the second sampling point reaches a fourth preset threshold and the voltage drop degree between the first sampling point and the second sampling point reaches a fifth preset threshold (a third preset condition), the value range of the fourth preset threshold can be 30% to 40%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold can be 1 uV to 2 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold can be 10% to 20%.

[0158] If the preset condition includes that the amplitude drop relative value between the first sampling point and the second sampling point reaches a fourth preset threshold value, or, the voltage drop degree between the first sampling point and the second sampling point reaches a fifth preset threshold value (the fourth preset condition), the value range of the fourth preset threshold value can be 50% to 60%. In this case, the value of the fifth preset threshold value is associated with the window length. In general, the value of the fifth preset threshold value is positively correlated with the window length, that is, the larger the window length, the larger the fifth preset threshold value. For example, taking the value range of the window length as 4 minutes to 6 minutes as an example, if the window length is 4 minutes, the value range of the fifth preset threshold value is 1 uV to 2 uV in the case of the voltage drop degree being the voltage drop absolute value, and the value range of the fifth preset threshold value is 10% to 20% in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value. If the window length is 5 minutes, the value range of the fifth preset threshold value is 2 uV to 3 uV in the case of the voltage drop degree being the voltage drop absolute value, and the value range of the second preset threshold value is 20% to 30% in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value. If the window length is 6 minutes, the value range of the fifth preset threshold value is 3 uV to 4 uV in the case of the voltage drop degree being the voltage drop absolute value, and the value range of the second preset threshold value is 30% to 40% in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value. Those skilled in the art can understand that the specific values of the above preset threshold values are exemplary, and are used to explain that the fifth preset threshold value is associated with the window length. For a smaller window length, if a larger fifth preset threshold value is used, the slow drop wave will be misjudged as a non-slow drop wave. For a larger window length, if a smaller fifth preset threshold value is used, the non-slow drop wave (such as a flat wave) will be misjudged as a slow drop wave, thereby reducing the recognition accuracy. Similarly, the window length is positively correlated with the fourth preset threshold value, so as to further improve the recognition accuracy. The specific values can be adjusted according to actual application.

[0159] The above processing can be performed on all candidate first high-frequency QRS complex data. If at least one candidate first high-frequency QRS complex data has a waveform feature indicating a myocardial bridge, such as a slow drop wave, it can be determined that there is a possibility of a myocardial bridge. Conversely, if all candidate first high-frequency QRS complex data do not have a waveform feature indicating a myocardial bridge, it can be determined that there is no myocardial bridge.

[0160] Although the window function used in the above embodiments starts from the earliest sampling point in the point sequence, those skilled in the art can understand that it can also start from the latest sampling point in the point sequence.

[0161] Figure 9An example flowchart illustrating a method for determining the likelihood of the presence of myocardial bridging according to another embodiment of this application is shown. Figure 9 As shown, in this embodiment, a window function can be preset. The window of this window function includes a start point and an end point, and the distance between the start point and the end point is the window length (in this embodiment, the window length is the preset time interval mentioned above). The value range of the window length can be, for example, 4 minutes to 6 minutes. For each candidate first high-frequency QRS complex data, the start point of the window can be aligned with the earliest sampling point in the point sequence in terms of time (i.e., these two points coincide in time). Using this first sampling point as a reference point, each sampling point within the window is traversed sequentially. Specifically, the relative value of amplitude decrease and the degree of voltage decrease between the reference point and the currently traversed sampling point within the window can be calculated. It is determined whether the reference point and the currently traversed sampling point meet preset conditions. If these two sampling points meet the preset conditions, it can be determined that the candidate first high-frequency QRS complex data has waveform characteristics indicating myocardial bridging, such as a slow-falling wave, and the possibility of myocardial bridging can be determined, and the traversal stops. If the preset conditions are not met, the next sampling point within the window is traversed. If the reference point and the currently traversed sampling point meet the preset conditions, the traversal stops; otherwise, it continues until the last sampling point within the traversal window. If the preset conditions are still not met after all sampling points within the window have been traversed, the traversal loop ends, and the window of the window function can be slid backward (time-wise backward) by a preset step size, and the above operation (the next traversal loop) is repeated. If no situation meets the preset conditions after the traversal loop ends, the window can be slid backward by the preset step size again until the preset conditions are met, or until the end of the window reaches or exceeds the last sampling point in the point sequence. That is, if no situation meets the preset conditions after the window has slid by a preset step size and its end reaches or exceeds the last sampling point, it indicates that there are no waveform features indicative of myocardial bridging, such as a slow-falling wave, for the candidate first high-frequency QRS complex data. In this embodiment, the preset step size can be, for example, an integer multiple of the time interval (e.g., sampling period) between two adjacent sampling points, preferably, the time interval between two adjacent sampling points, such as 10 seconds. Of course, those skilled in the art will understand that setting the preset step size to other values ​​is also feasible.

[0162] In this embodiment, if the preset condition includes that the amplitude drop relative value between the first sampling point and the second sampling point reaches the fourth preset threshold value and the voltage drop degree between the first sampling point and the second sampling point reaches the fifth preset threshold value (the third preset condition), the value range of the fourth preset threshold value can be 30% to 40%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold value can be 1 uV to 2 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold value can be 10% to 20%.

[0163] If the preset condition includes that the amplitude drop relative value between the first sampling point and the second sampling point reaches the fourth preset threshold value or the voltage drop degree between the first sampling point and the second sampling point reaches the fifth preset threshold value (the fourth preset condition), the value range of the fourth preset threshold value can be 50% to 60%. In this case, the value of the fifth preset threshold value is associated with the window length. In general, the value of the fifth preset threshold value is positively correlated with the window length, that is, the larger the window length, the larger the fifth preset threshold value. For example, taking the value range of the window length as 4 minutes to 6 minutes, if the window length is 4 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the fifth preset threshold value is 1 uV to 2 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 10% to 20%. If the window length is 5 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the fifth preset threshold value is 2 uV to 3 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 20% to 30%. If the window length is 6 minutes, in the case of the voltage drop degree being the voltage drop absolute value, the value range of the fifth preset threshold value is 3 uV to 4 uV, and in the case of the voltage drop degree being the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the second preset threshold value is 30% to 40%. Those skilled in the art can understand that the specific values of the above preset threshold values are exemplary, and are used to illustrate that the fifth preset threshold value is associated with the window length. For a smaller window length, if a larger fifth preset threshold value is used, the slow decline wave will be misjudged as a non-slow decline wave. For a larger window length, if a smaller fifth preset threshold value is used, the non-slow decline wave (such as a flat wave) will be misjudged as a slow decline wave, thereby reducing the recognition accuracy. Similarly, the window length is positively correlated with the fourth preset threshold value, so as to further improve the recognition accuracy. The specific values can be adjusted according to actual application.

[0164] The above processing can be performed on all candidate first high-frequency QRS complex data. If at least one candidate first high-frequency QRS complex data has a waveform feature indicative of a myocardial bridge, such as a slow descent wave, it can be determined that there is a possibility of a myocardial bridge. Conversely, if all candidate first high-frequency QRS complex data do not have a waveform feature indicative of a myocardial bridge, it can be determined that there is no myocardial bridge.

[0165] Although the window function used in the above embodiment starts from the earliest sampling point in the point sequence in terms of time sequence, it can be understood by those skilled in the art that it can also start from the latest sampling point in the point sequence in terms of time sequence.

[0166] In the above embodiment, the value of the preset time interval (e.g., preset time length, window length) for myocardial bridge determination is larger than the value of the preset time interval (e.g., preset time length, window length) for coronary artery lesion determination, the fourth preset threshold is greater than or equal to the first preset threshold, and the fifth preset threshold is greater than or equal to the second preset threshold. In addition, the time interval between two adjacent sampling points in the plurality of sampling points for myocardial bridge determination can be larger than the time interval between two adjacent sampling points in the plurality of sampling points for coronary artery lesion determination. For example, the time interval between two adjacent sampling points in the plurality of sampling points for myocardial bridge determination can be, for example, 12 seconds, and the time interval between two adjacent sampling points in the plurality of sampling points for coronary artery lesion determination can be, for example, 10 seconds.

[0167] In the preferred embodiment of the present application, in order to improve the accuracy of myocardial bridge determination, after obtaining the candidate first high-frequency QRS complex data, the candidate first high-frequency QRS complex data can be smoothed to smooth out the points with slight fluctuations. For example, the root mean square voltage / intensity / amplitude corresponding to a plurality of data points can be curve-smoothed in time sequence to obtain the candidate first high-frequency QRS complex data after smoothing.

[0168] In the preferred embodiment of the present application, the root mean square voltage of each sampling point between the first sampling point and the second sampling point in the preset condition can be determined in the following manner. Increasing the determination of this condition can further improve the accuracy of myocardial bridge determination or recognition.

[0169] Specifically, if it is determined that the candidate first high-frequency QRS complex data has a first sampling point and a second sampling point that satisfy the following conditions, i.e.,

[0170] the amplitude drop relative value between the first sampling point and the second sampling point reaches the fourth preset threshold and the voltage drop degree between the first sampling point and the second sampling point reaches the fifth preset threshold (the case of the third preset condition); or

[0171] The amplitude drop between the first sampling point and the second sampling point reaches a fourth preset threshold, or the voltage drop between the first sampling point and the second sampling point reaches a fifth preset threshold (the fourth preset condition). The first-order difference of the RMS voltage / amplitude / intensity of the sampling points in the candidate first high-frequency QRS complex data within the current preset time interval (for example, a preset time length, a window length) can be calculated to obtain a difference sequence, and the number of consecutive values greater than or equal to zero in the difference sequence is determined. If the number is less than or equal to a number threshold, for example, 2 (that is, the flat wave or rising wave can be excluded), it can be determined that the RMS voltage of each sampling point between the first sampling point and the second sampling point shows a continuous downward trend. At this time, it can be determined that the candidate first high-frequency QRS complex data has a waveform feature indicating a myocardial bridge, for example, a slow drop wave, and it can be determined that there is a possibility of a myocardial bridge.

[0172] In the preferred embodiment of the present application, in order to further improve the accuracy of myocardial bridge judgment or identification, for example, to further exclude the possibility of a steep drop wave in the candidate first high-frequency QRS complex data, after obtaining the difference sequence, all difference sequence groups are determined, each difference sequence group includes a number of consecutive difference values, and the sum of the number of difference values in each sequence group is calculated. For example, assuming that the difference sequence includes N difference values, the number of consecutive difference values is n, N and n are natural numbers, and n < N. The first sequence group can include the first difference value to the n-th difference value, the second sequence group can include the second difference value to the n+1-th difference value, and so on. If the number of consecutive values greater than or equal to zero in the difference sequence is less than or equal to 2, and the sum of the difference values of each difference sequence group is greater than or equal to a predetermined value, it can be determined that the candidate first high-frequency QRS complex data has a waveform feature indicating a myocardial bridge, for example, a slow drop wave, and it can be determined that there is a possibility of a myocardial bridge. The sum of the consecutive difference values can further exclude the steep drop, that is, exclude the case that the candidate first high-frequency QRS complex data has a steep drop wave. In an example, the number of consecutive difference values in each sequence group can be, for example, 5, and the predetermined value can be, for example, -15. Of course, those skilled in the art can understand that as long as the purpose of excluding the steep drop wave is achieved, other numerical values can be set according to actual application.

[0173] In theory, the RMS voltage drop of the first high-frequency QRS complex data caused by the myocardial bridge will continue until the exercise ends or even a period of time (for example, 20 seconds) after the exercise ends, and the duration is at least 3 minutes, and is a continuous slow drop. In alternative or additional embodiments of the present application, a technical solution can be designed according to this phenomenon to determine whether there is a possibility of a myocardial bridge.

[0174] Figure 10An example flowchart of a method for judging the possibility of existence of myocardial bridge according to another embodiment of the present application is schematically shown. As shown in the figure, in this embodiment, candidate first high-frequency QRS complex data can be determined from the first high-frequency QRS complex data corresponding to each electrocardiogram lead, which does not have waveform features indicative of coronary artery lesions, such as a sharp drop wave. Figure 10

[0175] For any candidate first high-frequency QRS complex data, at least one sampling point in the candidate first high-frequency QRS complex data within a preset time period is determined as a reference sampling point. In an example, the preset time period can be a time period from the time when the exercise ends in the exercise electrocardiogram detection process to a time after the exercise ends, for example, from the time when the exercise ends to 20-30 seconds after the exercise ends. In this embodiment, the reference sampling point in step S1043 is determined from the sampling points within a period of time after the exercise ends.

[0176] According to the determined reference sampling point, the sampling points to be traversed are determined, wherein the sampling points to be traversed are the sampling points within a preset time interval before the reference sampling point in the candidate first high-frequency QRS complex data. In an example, the preset time interval is within the exercise phase of the load exercise electrocardiogram detection process. The preset time interval can include a plurality of continuous sub-intervals (i.e., adjacent to each other but not overlapping). For example, taking the reference sampling point as the exercise end point and the time point of the exercise end point as 9 minutes (counted from the start of the load exercise electrocardiogram detection process), the preset time interval can be 3-6 minutes before the reference sampling point, and the preset time interval can include a first sub-interval, a second sub-interval and a third sub-interval. The first sub-interval can be 3-4 minutes before the reference sampling point (not including 3 minutes but including 4 minutes), and taking the start of the load exercise electrocardiogram detection process as the zero point, the first sub-interval can be represented as: [5 minutes, 6 minutes); the second sub-interval can be 4-5 minutes before the reference sampling point, which can also be represented as: [4 minutes, 5 minutes); and the third sub-interval can be 5-6 minutes before the reference sampling point, which can also be represented as: [3 minutes, 4 minutes). In the embodiments of the present application, if the number of at least one reference sampling point includes multiple, each reference sampling point can have a corresponding preset time interval. That is, the time position (timestamp) of the reference sampling point can be different, and the (time) position of the corresponding preset time interval can also be different, and the (time) position of each sub-interval of the corresponding preset time interval can also be different.

[0177] The sampling points in the plurality of sub-intervals can be traversed to determine whether there is a target sampling point in each sub-interval that meets the preset condition corresponding to the sub-interval. In the embodiments of the present application, the preset conditions corresponding to different sub-intervals can be different.​

[0178] In an example case, the preset condition can be that the root mean square voltage of each sampling point between the target sampling point and the reference sampling point is in a continuous decreasing trend, and the amplitude drop relative value between the target sampling point and the reference sampling point reaches a fourth preset threshold and the voltage drop degree between the target sampling point and the reference sampling point reaches a fifth preset threshold (a third preset condition). In another example case, the preset condition can be that the root mean square voltage of each sampling point between the target sampling point and the reference sampling point is in a continuous decreasing trend, and the amplitude drop relative value between the target sampling point and the reference sampling point reaches the fourth preset threshold or the voltage drop degree between the target sampling point and the reference sampling point reaches the fifth preset threshold (a fourth preset condition).

[0179] In the embodiments of the present application, the preset conditions corresponding to different sub-intervals can be different. Specifically, the fourth preset threshold and the fifth preset threshold corresponding to different sub-intervals can be different. The farther the sub-interval is from the reference sampling point in time, the greater the fourth preset threshold and the fifth preset threshold are.

[0180] For example, if the preset condition includes that the amplitude drop relative value between the target sampling point and the reference sampling point reaches the fourth preset threshold and the voltage drop degree between the target sampling point and the reference sampling point reaches the fifth preset threshold (the third preset condition), the value range of the fourth preset threshold can be 30% to 40% for the first sub-interval. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold can be 2 uV to 3 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold can be 30% to 40%. For the second sub-interval, the value range of the fourth preset threshold can be 40% to 50%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold can be 3 uV to 4 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold can be 40% to 50%. For the third sub-interval, the value range of the fourth preset threshold can be 50% to 60%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold can be 4 uV to 5 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold can be 50% to 60%.

[0181] If the preset condition includes that the amplitude drop relative value between the target sampling point and the reference sampling point reaches a fourth preset threshold value and the voltage drop degree between the target sampling point and the reference sampling point reaches a fifth preset threshold value (the third preset condition), for the first sub-interval, the value range of the fourth preset threshold value can be 40% to 50%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold value can be 3 uV to 4 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold value can be 40% to 50%. For the second sub-interval, the value range of the fourth preset threshold value can be 50% to 60%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold value can be 4 uV to 5 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold value can be 50% to 60%. For the third sub-interval, the value range of the fourth preset threshold value can be 60% to 70%. In this case, if the voltage drop degree is the voltage drop absolute value, the value range of the fifth preset threshold value can be 5 uV to 6 uV. If the voltage drop degree is the ratio of the voltage drop absolute value to the maximum voltage value, the value range of the fifth preset threshold value can be 60% to 70%.

[0182] During the traversal, if there is a target sampling point satisfying the corresponding preset condition in any one of the plurality of sub-intervals, it can be determined that there is a possibility of myocardial bridge. The traversal manner can include sequentially traversing the to-be-traversed sampling points in the reverse order of time sequence, or performing parallel traversal for the plurality of sub-intervals. For the case where the at least one reference sampling point includes a plurality of reference sampling points, if no target reference point satisfying the preset condition is found after completing the traversal for a certain reference sampling point, the traversal operation is switched to the next reference sampling point to be re-performed. If no target reference point satisfying the preset condition is found after the traversal is performed for all the reference sampling points, it can be determined that the candidate first high-frequency QRS complex data does not have the waveform feature indicating myocardial bridge.

[0183] All candidate first high-frequency QRS complex data can be processed as above, and if at least one candidate first high-frequency QRS complex data has a waveform feature indicating myocardial bridge, such as a slow descent wave, it can be determined that there is a possibility of myocardial bridge. Conversely, if all candidate first high-frequency QRS complex data do not have the waveform feature indicating myocardial bridge, it can be determined that there is no myocardial bridge.

[0184] The determination method of the root mean square voltage of each sampling point between the target sampling point and the reference sampling point in the preset condition in this embodiment can be the same as that in the previous embodiments, which will not be described here.

[0185] In step S105, the possibility of existence of the coronary microcirculation disorder is determined according to the exercise electrocardiogram data and the resting electrocardiogram data.

[0186] Figure 11 An example flowchart of a method for determining the possibility of existence of the coronary microcirculation disorder according to an embodiment of the present application is schematically shown. As shown, specifically, in this embodiment, step S105 can include the following steps. Figure 11

[0187] The first high-frequency QRS complex data can be obtained according to the exercise electrocardiogram data, and the second high-frequency QRS complex data can be obtained according to the resting electrocardiogram data.

[0188] In step S1051, the first high-frequency QRS complex data in which the lead positivity index indicates positive is determined from all the first high-frequency QRS complex data.

[0189] In step S1052, it is determined whether there is target first high-frequency QRS complex data in the first high-frequency QRS complex data in which the lead positivity index indicates positive, the target first high-frequency QRS complex data not having a waveform feature (which can be referred to as a first waveform feature here) indicative of coronary artery lesion and / or a waveform feature (which can be referred to as a second waveform feature here) indicative of myocardial bridge. When analyzing the first high-frequency QRS complex data to determine whether there is the coronary microcirculation disorder, the first high-frequency QRS complex data having the first waveform feature (such as a steep drop wave) and / or the second waveform feature (such as a slow drop wave) can be excluded first. That is, if the first high-frequency QRS complex data has the first waveform feature and / or the second waveform feature, the first high-frequency QRS complex data is no longer used to determine whether there is the coronary microcirculation disorder. The method for determining whether the first high-frequency QRS complex data has the waveform feature indicative of the coronary artery lesion is as described in the embodiment regarding step S103, and the method for determining whether the first high-frequency QRS complex data has the waveform feature indicative of the myocardial bridge is as described in the embodiment regarding step S104, which will not be repeated here.

[0190] In step S1053, the possibility of existence of the coronary microcirculation disorder is determined in the case where the target first high-frequency QRS complex data is determined to exist.

[0191] In order to further improve the accuracy of determining the possibility of existence of the coronary microcirculation disorder, it can be further determined whether the lead positivity index of the target second high-frequency QRS complex data corresponding to the target electrocardiogram lead associated with the target first high-frequency QRS complex data indicates positive.

[0192] ​Specifically, when a target first high-frequency QRS complex data is determined, a target electrocardiogram lead associated with the target first high-frequency QRS complex data can be determined, i.e., the target first high-frequency QRS complex data is output through the target electrocardiogram lead. The target second high-frequency QRS complex data output through the target electrocardiogram lead can be analyzed to determine whether the lead positivity index of the target second high-frequency QRS complex data indicates positivity.

[0193] In the case where the lead positivity index of the target second high-frequency QRS complex data indicates positivity, it is determined that there is a possibility of coronary microcirculation disorder. For example, if the target electrocardiogram lead corresponding to the target first high-frequency QRS complex data is the V1 lead, and the lead positivity index corresponding to the second high-frequency QRS complex data corresponding to the V1 lead in the resting state indicates positivity, it is determined that there is a possibility of coronary microcirculation disorder.

[0194] In the embodiments of the present application, step S105 can further include determining that there is no coronary microcirculation disorder. Specifically, in step S1054, a first maximum RMS voltage is determined according to the maximum value of the RMS voltage in the first high-frequency QRS complex data corresponding to all electrocardiogram leads, and a second maximum voltage is determined according to the maximum value of the voltage in the second high-frequency QRS complex data corresponding to all electrocardiogram leads.

[0195] In an example, the first maximum RMS voltage can be the maximum voltage value in the above embodiments, and for single exercise load detection, the maximum value of the RMS voltage in the first high-frequency QRS complex data corresponding to all electrocardiogram leads can be rounded up to an even integer to obtain a voltage value as the first maximum RMS voltage. In another example, the maximum value of the RMS voltage in the first high-frequency QRS complex data corresponding to all electrocardiogram leads can also be taken as the first maximum RMS voltage.

[0196] Similarly, the second maximum voltage can round up the maximum value of the voltage / amplitude in the second high-frequency QRS complex data (or the second high-frequency QRS envelope curve) corresponding to all electrocardiogram leads to an even integer to obtain a voltage value as the second maximum voltage. In another example, the maximum value of the voltage in the second high-frequency QRS complex data corresponding to all electrocardiogram leads can also be taken as the second maximum voltage.

[0197] In step S1055, the number of positive leads in at least one electrocardiogram lead corresponding to the second QRS complex data whose second lead positivity index indicates positivity is determined according to the second high-frequency QRS complex data;

[0198] In step S1056, the number of positive leads is compared with a preset value, the first maximum RMS voltage is compared with a first voltage threshold, and the second maximum voltage is compared with a second voltage threshold;

[0199] In step S1057, if the number of positive leads is less than the preset value, and at least the first maximum RMS voltage is greater than the first voltage threshold, it indicates that the heart-lung function of the subject is strong, and it can be determined that there is no coronary microcirculation disorder. At least the first maximum RMS voltage greater than the first voltage threshold can include that the first maximum RMS voltage is greater than the first voltage threshold, and the second maximum voltage is greater than the second voltage threshold, or the first maximum RMS voltage is greater than the first voltage threshold, but the second maximum voltage is less than or equal to the second voltage threshold. In an example, the preset value can be 3, for example, the first voltage threshold can be 8uV, and the second voltage threshold can be 14uV. However, those skilled in the art can understand that the preset value and the voltage threshold can be set to other values according to actual application.

[0200] In step S106, the heart health status is qualitatively evaluated according to the determination results of the myocardial ischemia possibility, the coronary artery lesion possibility, the myocardial bridge possibility, and the coronary microcirculation disorder possibility. Specifically, the qualitative evaluation result can include, for example, myocardial ischemia, coronary artery lesion combined with myocardial ischemia, myocardial bridge combined with myocardial ischemia, microcirculation disorder combined with myocardial ischemia, coronary artery lesion combined with myocardial ischemia combined with myocardial bridge, only coronary artery lesion, etc.

[0201] In step S107, if it is determined that there is a myocardial ischemia possibility and / or a coronary artery lesion possibility, the myocardial ischemia and / or the coronary artery lesion is quantitatively evaluated.

[0202] In step S108, the heart health status is evaluated according to the qualitative evaluation result and the quantitative evaluation result.

[0203] It should be noted that although the above embodiment describes the quantitative evaluation of the myocardial ischemia and / or coronary artery lesion possibility in step S107, it can be understood that the quantitative evaluation of the myocardial ischemia and / or coronary artery lesion can be performed after determining the existence of the possibility. Taking myocardial ischemia as an example, if it is determined that there is a myocardial ischemia possibility in step S102, the myocardial ischemia possibility can be quantitatively evaluated. This processing mode is also applicable to coronary artery lesion.

[0204] The quantitative evaluation methods of myocardial ischemia and coronary artery lesion will be described below.

[0205] Quantitative assessment of myocardial ischemia

[0206] In a case where the first lead positive indicator, the second lead positive indicator, and the third lead positive indicator all indicate positive, it can be determined that myocardial ischemia exists. The myocardial ischemia can be quantitatively evaluated.

[0207] In the embodiments of the present application, the method for quantitatively evaluating myocardial ischemia can include the following steps.

[0208] In step S1021, the first lead number and the second lead number are obtained according to the first lead positive indicator and the second lead positive indicator, respectively.

[0209] Specifically, the first lead number can be obtained according to the first lead positive indicator, which indicates the number of leads corresponding to the first lead positive indicator indicating positive. The second lead number can be obtained according to the second lead positive indicator, which indicates the number of leads corresponding to the second lead positive indicator indicating positive.

[0210] In an embodiment of the present application, the first lead number not only includes the number of leads indicated as positive by the first lead positive indicator, but also includes the number of leads indicated as critical by the first lead positive indicator. Similarly, the second lead number not only includes the number of leads indicated as positive by the second lead positive indicator, but also includes the number of leads indicated as critical by the second lead positive indicator. The determination method of the lead positive indicator can adopt a method known to those skilled in the art. The lead positive indicator can include positive, critical, and negative, which will not be described here.

[0211] In step S1022, the resting electrocardio data is analyzed to determine a target high-frequency morphology index.

[0212] Specifically, the resting electrocardio data can be analyzed to obtain the high-frequency morphology index of each lead. The maximum value is selected from the high-frequency morphology index and determined as the target high-frequency morphology index (i.e., the maximum index value). The high-frequency morphology index of the resting electrocardio signal is a key parameter in high-frequency electrocardiogram (HFECG) for evaluating myocardial ischemia. In particular, in the resting HFQRS technology, the high-frequency component (usually more than 100 Hz) of the QRS complex is analyzed to quantify the abnormal conduction of myocardial cells. This index can reflect the slight changes in myocardial electrical activity and has high sensitivity to early myocardial ischemia.

[0213] In step S1023, the first lead number, the second lead number, and the target high-frequency morphology index are normalized to obtain a first ischemia index, a second ischemia index, and a third ischemia index, respectively.

[0214] Specifically, it is known to those skilled in the art that the first lead number, the second lead number and the target high frequency morphology index can reflect the degree of myocardial ischemia. The greater the first lead number, the more serious the myocardial ischemia, the greater the second lead number, the more serious the myocardial ischemia, and the higher the target high frequency morphology index, the more serious the myocardial ischemia. Normalizing the first lead number, the second lead number and the target high frequency morphology index can map the first lead number, the second lead number and the target high frequency morphology index to a value in the range of [0, 1], i.e. an ischemia index. In the embodiments of the present application, the correspondence between different values in [0, 1] and the first lead number, the second lead number or the target high frequency morphology index can be set according to different grades of myocardial ischemia severity. In an example, taking the use of a 10-electrode patch, i.e. 12 ECG leads, as an example, it is assumed that the severity of myocardial ischemia can be divided into four grades of "mild", "moderate", "heavy" and "severe". If the first lead number is zero, it means that there is no myocardial ischemia, and the corresponding first ischemia index can be 1. If the first lead number is greater than or equal to 1 but less than 3, the corresponding myocardial ischemia severity is "mild", and the value of the corresponding first ischemia index can be 0.8. If the first lead number is greater than or equal to 3 but less than 5, the corresponding myocardial ischemia severity is "moderate", and the value of the corresponding first ischemia index can be 0.6. If the first lead number is greater than or equal to 5 but less than 7, the corresponding myocardial ischemia severity is "heavy", and the value of the corresponding first ischemia index can be 0.4. If the first lead number is greater than or equal to 7, the corresponding myocardial ischemia severity is "severe", and the value of the corresponding first ischemia index can be 0.2. In this example, the more serious the myocardial ischemia, the smaller the value of the first ischemia index.

[0215] The corresponding relationship between the second lead number and the value of the second ischemia index can also be processed similarly. In an example, still taking the 10-electrode patch with 12 ECG leads as an example, assuming that the severity of myocardial ischemia can be divided into four levels of “mild”, “moderate”, “heavy”, and “severe”, if the second lead number is zero, it means that there is no myocardial ischemia, and the corresponding second ischemia index can be 1. If the second lead number is greater than or equal to 1 but less than 4, the corresponding severity of myocardial ischemia is “mild”, and the value of the corresponding second ischemia index can be 0.8. If the second lead number is greater than or equal to 4 but less than 7, the corresponding severity of myocardial ischemia is “moderate”, and the value of the corresponding second ischemia index can be 0.6. If the second lead number is greater than or equal to 7 but less than 10, the corresponding severity of myocardial ischemia is “heavy”, and the value of the corresponding second ischemia index can be 0.4. If the second lead number is greater than or equal to 10, the corresponding severity of myocardial ischemia is “severe”, and the value of the corresponding second ischemia index can be 0.2. In this example, the more severe the myocardial ischemia, the smaller the value of the second ischemia index.

[0216] The corresponding relationship between the target high-frequency morphology and the value of the third ischemia index can also be processed similarly. In an example, still taking the 10-electrode patch with 12 ECG leads as an example, assuming that the severity of myocardial ischemia can be divided into four levels of “mild”, “moderate”, “heavy”, and “severe”, if the value of the target high-frequency morphology is less than 15%, it means that there is no myocardial ischemia, and the corresponding third ischemia index can be 1. If the value of the target high-frequency morphology is greater than or equal to 15% but less than 25%, the corresponding severity of myocardial ischemia is “mild”, and the value of the corresponding third ischemia index can be 0.8. If the value of the target high-frequency morphology is greater than or equal to 25% but less than 35%, the corresponding severity of myocardial ischemia is “moderate”, and the value of the corresponding third ischemia index can be 0.6. If the value of the target high-frequency morphology is greater than or equal to 35% but less than 45%, the corresponding severity of myocardial ischemia is “heavy”, and the value of the corresponding third ischemia index can be 0.4. If the value of the target high-frequency morphology is greater than or equal to 45%, the corresponding severity of myocardial ischemia is “severe”, and the value of the corresponding third ischemia index can be 0.2. In this example, the more severe the myocardial ischemia, the smaller the value of the third ischemia index.

[0217] The above division of the severity of myocardial ischemia, the corresponding relationship between different severity levels and the first lead number, the second lead number, and the target high-frequency morphology index, and the corresponding relationship between the values of different first lead numbers, second lead numbers, and target high-frequency morphology indexes and the values of the first ischemia index, the second ischemia index, and the third ischemia index are exemplary, and those skilled in the art can understand that different corresponding relationships can be set according to actual needs.

[0218] In step S1024, the target ischemia index is obtained based on the first ischemia index, the second ischemia index, and the third ischemia index.

[0219] Specifically, in one example, the product of the first ischemia index, the second ischemia index, and the third ischemia index can be used as the target ischemia index reflecting the severity of myocardial ischemia. In another example, the different influences of the first, second, and third ischemia indices on reflecting the severity of myocardial ischemia can be considered, and influence factors (i.e., weights) can be assigned to the first, second, and third ischemia indices respectively. In this example, the target ischemia index can be the sum of the first, second, and third ischemia indices multiplied by their respective influence factors. The calculation method of the target ischemia index can be as shown in the following formula (1):

[0220] 1+b I2+c I3 formula (1)

[0221] Among them, I target Let I1 be the target ischemia index, a be the influencing factor of the first ischemia index, I2 be the second ischemia index, b be the influencing factor of the second ischemia index, I3 be the third ischemia index, and c be the influencing factor of the third ischemia index. In one example, a, b, and c can satisfy a + b + c = 1.

[0222] In another example, the smallest of the first, second, and third ischemia indices can be used as the target ischemia index. The core logic of this example is: based on the number of leads in the first lead, the number of leads in the second lead, and the target high-frequency morphology index, a degree of myocardial ischemia can be assessed separately, and then the most severe of the three (highest grade) is selected as the final degree of myocardial ischemia.

[0223] Quantitative assessment of coronary artery lesions

[0224] As described above, when the possibility of coronary artery disease is determined, the maximum relative amplitude decrease and the maximum voltage decrease are determined from the relative amplitude decrease and voltage decrease (e.g., the absolute voltage decrease, or the ratio of the absolute voltage decrease to the maximum voltage value) for each electrocardiogram lead. The coronary artery disease is then quantitatively assessed based on the maximum relative amplitude decrease and the maximum voltage decrease.

[0225] For example, the coronary artery lesion can be set multiple risk levels (for example, D1, D2, D3, etc.), each risk corresponds to a maximum amplitude drop relative value interval and a maximum voltage drop degree interval, and the respective corresponding risk level is determined according to the interval in which the maximum amplitude drop relative value and the maximum voltage drop degree are respectively located. If the risk level corresponding to the maximum amplitude drop relative value is different from the risk level corresponding to the maximum voltage drop degree, the more serious risk level of the two risk levels is taken as the finally determined risk level.

[0226] In an embodiment of the present application, a processor is provided, configured to perform the method described in each of the above embodiments.

[0227] In an embodiment of the present application, a machine readable storage medium is provided, and the machine readable storage medium has instructions stored thereon for causing a machine to perform the method described in each of the above embodiments.

[0228] In an embodiment of the present application, a processing device is provided, comprising a processor and a memory, and the memory has instructions stored therein, and the processor is configured to call and execute the instructions from the memory to implement the method described in each of the above embodiments.

[0229] Examples of the processor can include, but are not limited to, a single-chip microcomputer, a microprocessor, a digital signal processor (DSP), a field programmable gate array (FPGA), an application-specific integrated circuit (ASIC), etc.

[0230] The memory can include a non-persistent memory in a computer readable medium, a random access memory (RAM) and / or a non-volatile memory such as a read-only memory (ROM) or a flash memory (flash RAM), and the memory includes at least one memory chip.

[0231] Figure 12 An example block diagram of a heart health status evaluation system according to an embodiment of the present application is schematically shown. As shown in the figure, Figure 12 In an embodiment of the present application, a heart health status evaluation system is provided, which can include:

[0232] The electrocardiosignal acquisition device 100 includes at least one electrode forming at least one electrocardiogram lead for acquiring the electrocardiosignal of the subject; and

[0233] The processing device 200 described above.

[0234] The ECG signal acquisition device 100 can also include a display configured to display various data and information, including but not limited to, for example, subject information, such as age, height, weight, etc., and ECG signal / QRS complex data / QRS waveform curve, etc. acquired through ECG leads.

[0235] The processing device 200 can be configured to perform various methods described in the above embodiments.

[0236] In particular, in an embodiment, a method performed by a computer for assessing a heart health status of a subject can include:

[0237] acquiring exercise ECG signals of the subject acquired through at least one ECG lead during a load exercise ECG detection process, and resting ECG signals of the subject acquired through at least one ECG lead during a resting ECG detection process;

[0238] analyzing the exercise ECG signals to obtain exercise ECG data, and analyzing the resting ECG signals to obtain resting ECG data and routine ECG data;

[0239] determining a likelihood of existence of myocardial ischemia according to the exercise ECG data, the resting ECG data, and the routine ECG data;

[0240] determining a likelihood of existence of coronary artery lesion according to the exercise ECG data;

[0241] determining a likelihood of existence of myocardial bridge according to the exercise ECG data;

[0242] determining a likelihood of existence of coronary microcirculation disorder according to the exercise ECG data and the resting ECG data;

[0243] qualitatively assessing the heart health status according to the determination results of the likelihood of existence of myocardial ischemia, the likelihood of existence of coronary artery lesion, the likelihood of existence of myocardial bridge, and the likelihood of existence of coronary microcirculation disorder;

[0244] in a case where the likelihood of existence of myocardial ischemia and / or the likelihood of existence of coronary artery lesion is determined, quantitatively assessing the myocardial ischemia and / or the coronary artery lesion; and

[0245] assessing the heart health status according to the qualitative assessment result and the quantitative assessment result.

[0246] In an embodiment, the determining of the likelihood of existence of myocardial ischemia according to the exercise ECG data, the resting ECG data, and the routine ECG data includes:

[0247] obtaining a first lead positive indicator corresponding to each ECG lead according to first high-frequency QRS complex data obtained from the exercise ECG data;

[0248] a second lead positive indicator corresponding to each electrocardiogram lead is obtained according to second high frequency QRS complex data obtained from the resting electrocardiogram data;

[0249] a third lead positive indicator corresponding to each electrocardiogram lead is obtained according to conventional electrocardiogram data;

[0250] in a case where any one of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicates positivity, it is determined that there is a possibility of myocardial ischemia.

[0251] in a case where only one of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicates positivity, the possibility of myocardial ischemia is low in an embodiment;

[0252] in a case where two of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicate positivity, the possibility of myocardial ischemia is high;

[0253] in a case where all of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicate positivity, the possibility of myocardial ischemia is the highest.

[0254] in an embodiment, determining the possibility of coronary artery lesion according to exercise electrocardiogram data comprises:

[0255] a first point sequence corresponding to each electrocardiogram lead is obtained according to first high frequency QRS complex data corresponding to each electrocardiogram lead obtained from exercise electrocardiogram data, the first point sequence comprising a plurality of sampling points arranged in time sequence;

[0256] for the first point sequence corresponding to any electrocardiogram lead, in a traversal loop, any sampling point in the plurality of sampling points in the first point sequence is taken as a current reference sampling point, and sampling points within a first preset time interval from the current reference sampling point are traversed to determine whether there are first and second sampling points satisfying a preset condition within the first preset time interval, wherein the first sampling point is earlier than the second sampling point in time sequence, and the sampling points within the first preset time interval include the current reference sampling point;

[0257] in a case where it is determined that there are first and second sampling points satisfying the preset condition, it is determined that the first high frequency QRS complex data corresponding to the any electrocardiogram lead has a first waveform feature indicating coronary artery lesion, thereby determining the possibility of coronary artery lesion;

[0258] in a case where it is determined that there are no first and second sampling points satisfying the preset condition, the current traversal loop ends, and other sampling points in the first point sequence except the current reference sampling point are determined as new reference sampling points;

[0259] The next traversal cycle is performed for a new reference sampling point until it is determined that there is a first sampling point and a second sampling point satisfying the preset condition, or all sampling points in the first point sequence that need to be traversed are traversed.

[0260] In an embodiment, determining the possibility of existence of the coronary artery lesion according to the motion electrocardiogram data further comprises:

[0261] In a case where there is no first sampling point and second sampling point satisfying the preset condition after all sampling points in the first point sequence that need to be traversed are traversed, it is determined that there is no first waveform feature indicating the coronary artery lesion in the corresponding first high-frequency QRS complex data.

[0262] In a case where there is no first waveform feature indicating the coronary artery lesion in the first high-frequency QRS complex data corresponding to all electrocardiogram leads, it is determined that there is no possibility of existence of the coronary artery lesion.

[0263] In an embodiment, the preset condition comprises:

[0264] The amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value, and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value; or

[0265] The amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold value, or the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold value.

[0266] In an embodiment, the first preset time interval includes a preset time length. For the first point sequence corresponding to any electrocardiogram lead, in a traversal cycle, any sampling point in the plurality of sampling points in the first point sequence is taken as a current reference sampling point, and the sampling points in the first preset time interval from the current reference sampling point are traversed to determine whether there is a first sampling point and a second sampling point satisfying the preset condition in the first preset time interval, comprising:

[0267] The first sampling point in the first point sequence that is the earliest in time sequence is taken as the current reference sampling point, and the next sampling point in time sequence from the current reference sampling point is traversed to determine whether the current traversal sampling point and the current reference sampling point satisfy the preset condition.

[0268] In a case where the current traversal sampling point and the current reference sampling point do not satisfy the preset condition, it is determined whether the time interval between the current traversal sampling point and the current reference sampling point exceeds the preset time length.

[0269] In a case where it is determined that the time interval between the current traversal sampling point and the current reference sampling point exceeds the preset time length, it is determined that the traversal cycle ends.

[0270] In an embodiment, in a case where it is determined that there exist the first sampling point and the second sampling point satisfying the preset condition, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has the first waveform feature indicating the coronary artery lesion, thereby determining that there exists the coronary artery lesion possibility, comprising:

[0271] In a case where it is determined that there exist the current iteration sampling point and the current reference sampling point satisfying the preset condition during the iteration loop, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has the first waveform feature indicating the coronary artery lesion, thereby determining that there exists the coronary artery lesion possibility.

[0272] In an embodiment, after the current iteration loop ends, other sampling points in the first point sequence except the current reference sampling point are determined as new reference sampling points, comprising:

[0273] The next sampling point of the current reference sampling point in time sequence is determined as the new reference sampling point.

[0274] In an embodiment, the next iteration loop is performed for the new reference sampling point, and the iteration loop is repeated until it is determined that there exist the first sampling point and the second sampling point satisfying the preset condition, or all the sampling points to be iterated in the first point sequence are iterated, comprising:

[0275] The next iteration loop is started from the next sampling point of the new reference sampling point in time sequence, and the iteration loop is repeated until it is determined that there exist the current iteration sampling point and the current reference sampling point satisfying the preset condition, or the current iteration sampling point is the last sampling point in the first point sequence.

[0276] In an embodiment, the first preset time interval includes a window length of a window function. For the first point sequence corresponding to the arbitrary electrocardiogram lead, in an iteration loop, any sampling point in the first point sequence is taken as a current reference sampling point, and the sampling points in the first preset time interval from the current reference sampling point are iterated to determine whether there exist the first sampling point and the second sampling point satisfying the preset condition in the first preset time interval, comprising:

[0277] The first sampling point with the maximum root mean square (RMS) voltage and the second sampling point with the minimum RMS voltage after the first sampling point in time sequence in the current window are determined by taking the first sampling point in the first point sequence in time sequence as a starting point of the window and iterating the sampling points in the first point sequence contained in the window.

[0278] It is determined whether the first sampling point and the second sampling point satisfy the preset condition.

[0279] In an embodiment, in a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has the first waveform feature indicating the coronary artery lesion, thereby determining that the coronary artery lesion possibility exists, comprising:

[0280] In a case where it is determined that the first sampling point and the second sampling point in the current window satisfy the preset condition, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has the first waveform feature indicating the coronary artery lesion, thereby determining that the coronary artery lesion possibility exists.

[0281] In an embodiment, in a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition do not exist, the current traversal loop ends, and other sampling points in the first point sequence except the current reference sampling point are determined as new reference sampling points, comprising:

[0282] In a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition do not exist in the current window, the current traversal loop ends, and the window is moved backward by the preset step in time sequence.

[0283] In an embodiment, a next traversal loop is performed for the new reference sampling point until it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, or all the sampling points to be traversed in the first point sequence are traversed, comprising:

[0284] A next traversal loop is performed for the sampling points in the window moved by the preset step, and the traversal loop is repeated until it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, or the end point of the current window reaches at least the last sampling point in the first point sequence.

[0285] In an embodiment, the first preset time interval includes a window length of a window function. For the first point sequence corresponding to the arbitrary electrocardiogram lead, in a traversal loop, any sampling point in the plurality of sampling points in the first point sequence is taken as a current reference sampling point, and the sampling points in the first preset time interval from the current reference sampling point are traversed to determine whether the first sampling point and the second sampling point satisfying the preset condition exist in the first preset time interval, comprising:

[0286] The first sampling point in the first point sequence that is earliest in time sequence is taken as a starting point of the window, and the sampling points in the first point sequence contained in the window are traversed to determine whether the first sampling point and the second sampling point satisfying the preset condition exist in the current window;

[0287] In an embodiment, in a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has a first waveform feature indicating coronary artery lesion, thereby determining that the coronary artery lesion possibility exists, comprising:

[0288] In a case where it is determined that the first sampling point and the second sampling point in the current window satisfy the preset condition, it is determined that the first high-frequency QRS complex data corresponding to the arbitrary electrocardiogram lead has a first waveform feature indicating coronary artery lesion, thereby determining that the coronary artery lesion possibility exists.

[0289] In an embodiment, in a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition do not exist, the current traversal loop ends, and other sampling points in the first point sequence except the current reference sampling point are determined as new reference sampling points, comprising:

[0290] In a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition do not exist in the current window, the current traversal loop ends, and the window is moved backward by a preset step in time sequence.

[0291] In an embodiment, a next traversal loop is performed for the new reference sampling point, and the traversal loop is repeated until it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, or all the sampling points to be traversed in the first point sequence are traversed, comprising:

[0292] A next traversal loop is performed on the sampling points in the window moved by the preset step, and the traversal loop is repeated until it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, or the end point of the current window reaches at least the last sampling point in the first point sequence.

[0293] In an embodiment, the determining the coronary artery lesion possibility from the motion electrocardiogram data further comprises:

[0294] In a case where it is determined that the first sampling point and the second sampling point satisfying the preset condition exist, a first-order difference of mean square voltage of the sampling points between the first sampling point and the second sampling point in the first point sequence is determined to determine a minimum value sampling point corresponding to a minimum value and at least one maximum value sampling point corresponding to a maximum value;

[0295] A first root mean square voltage difference value between the maximum value sampling point and the minimum value sampling point is calculated;

[0296] A second root mean square voltage difference value between the first sampling point and the second sampling point is calculated;

[0297] In a case where a ratio of the first root mean square voltage difference value and the second root mean square voltage difference value is greater than or equal to a third threshold value, the first sampling point and the second sampling point are re-determined as not satisfying the preset condition.

[0298] In an embodiment, determining the likelihood of existence of the coronary artery lesion according to the exercise electrocardiogram data further comprises:

[0299] determining a first lead positive indicator corresponding to the first high-frequency QRS complex data of the first waveform feature indicative of the coronary artery lesion;

[0300] in a case where the first lead positive indicator indicates positive or critical, determining that the type of the coronary artery lesion is coronary stenosis;

[0301] in a case where the first lead positive indicator indicates negative, determining that the type of the coronary artery lesion is coronary sclerosis.

[0302] In an embodiment, determining the likelihood of existence of the myocardial bridge according to the exercise electrocardiogram data comprises:

[0303] determining candidate first high-frequency QRS complex data from the first high-frequency QRS complex data corresponding to each electrocardiogram lead, the candidate first high-frequency QRS complex data being free of the first waveform feature indicative of the coronary artery lesion;

[0304] obtaining a second point sequence corresponding to the candidate first high-frequency QRS complex data, the second point sequence comprising a plurality of sampling points arranged in time sequence;

[0305] for the second point sequence corresponding to any candidate first high-frequency QRS complex data, in a traversal loop, taking any sampling point in the plurality of sampling points in the second point sequence as a current reference sampling point, and traversing the sampling points within a second preset time interval from the current reference sampling point to determine whether there are third and fourth sampling points satisfying another preset condition within the second preset time interval, wherein the third sampling point is earlier in time sequence than the fourth sampling point;

[0306] in a case where it is determined that there are third and fourth sampling points satisfying the another preset condition, determining that the any candidate first high-frequency QRS complex data has the second waveform feature indicative of the myocardial bridge, thereby determining the likelihood of existence of the myocardial bridge;

[0307] in a case where it is determined that there are no third and fourth sampling points satisfying the another preset condition, ending the current traversal loop, and determining other sampling points in the second point sequence except the current reference sampling point as new reference sampling points;

[0308] performing the next traversal loop for the new reference sampling points until it is determined that there are third and fourth sampling points satisfying the another preset condition, or all the sampling points in the second point sequence that need to be traversed are traversed.

[0309] In an embodiment, the second preset time interval includes a preset time length. For the second point sequence corresponding to any candidate first high-frequency QRS complex data, in an iteration loop, any sampling point in the second point sequence is taken as a current reference sampling point, and sampling points within the second preset time interval from the current reference sampling point are iterated to determine whether there are third and fourth sampling points satisfying another preset condition within the second preset time interval, including:

[0310] The first sampling point in the second point sequence in time sequence is taken as the current reference sampling point, and the next sampling point in time sequence from the current reference sampling point is taken as a current iteration sampling point to determine whether the current iteration sampling point and the current reference sampling point satisfy the another preset condition;

[0311] In a case where the current iteration sampling point and the current reference sampling point do not satisfy the another preset condition, it is determined whether a time interval between the current iteration sampling point and the current reference sampling point exceeds the preset time length;

[0312] In a case where it is determined that the time interval between the current iteration sampling point and the current reference sampling point exceeds the preset time length, it is determined that the iteration loop ends.

[0313] In an embodiment, in a case where it is determined that there are third and fourth sampling points satisfying the another preset condition, it is determined that the any candidate first high-frequency QRS complex data has the second waveform feature indicating a myocardial bridge, thereby determining the possibility of the myocardial bridge, including:

[0314] In a case where it is determined that there are the current iteration sampling point and the current reference sampling point satisfying the preset condition during the iteration loop, it is determined that the any candidate first high-frequency QRS complex data has the second waveform feature indicating the myocardial bridge, thereby determining the possibility of the myocardial bridge.

[0315] In an embodiment, after the current iteration loop ends, other sampling points in the second point sequence except the current reference sampling point are determined as new reference sampling points, including:

[0316] The next sampling point in time sequence from the current reference sampling point is determined as the new reference sampling point.

[0317] In an embodiment, a next iteration loop is performed for the new reference sampling point until it is determined that there are third and fourth sampling points satisfying the another preset condition, or all sampling points in the second point sequence to be iterated are iterated, including:

[0318] The next iteration loop is started from the next sampling point in time sequence from the new reference sampling point, and the iteration loop is repeated until it is determined that there are the current iteration sampling point and the current reference sampling point satisfying the preset condition, or the current iteration sampling point is the last sampling point in the second point sequence.

[0319] In an embodiment, the second preset time interval comprises a window length of a window function. For a second point sequence corresponding to any candidate first high frequency QRS complex data, in a traversal loop, any sampling point in the second point sequence is taken as a current reference sampling point, and sampling points within the second preset time interval from the current reference sampling point are traversed to determine whether there are third and fourth sampling points satisfying another preset condition in the second preset time interval, wherein the third sampling point is earlier than the fourth sampling point in time sequence, comprising:

[0320] Taking the first sampling point in the second point sequence as the starting point of the window, the sampling points in the second point sequence contained in the window are traversed to determine the third sampling point with the maximum root mean square (RMS) voltage and the fourth sampling point with the minimum RMS voltage after the third sampling point in the current window in time sequence;

[0321] Determining whether the third and fourth sampling points satisfy another preset condition.

[0322] In an embodiment, in a case where it is determined that there are third and fourth sampling points satisfying another preset condition, it is determined that the arbitrary candidate first high frequency QRS complex data has a second waveform feature indicating a myocardial bridge, thereby determining the possibility of a myocardial bridge, comprising:

[0323] In a case where it is determined that the third and fourth sampling points in the current window satisfy another preset condition, it is determined that the arbitrary candidate first high frequency QRS complex data has a second waveform feature indicating a myocardial bridge, thereby determining the possibility of a myocardial bridge.

[0324] In an embodiment, in a case where it is determined that there are no third and fourth sampling points satisfying another preset condition, the current traversal loop ends, and other sampling points in the second point sequence except the current reference sampling point are determined as new reference sampling points, comprising:

[0325] In a case where it is determined that there are no third and fourth sampling points satisfying another preset condition in the current window, the current traversal loop ends, and the window is moved backward by a preset step in time sequence.

[0326] In an embodiment, the new reference sampling points perform a next traversal loop until it is determined that there are third and fourth sampling points satisfying another preset condition, or all sampling points in the second point sequence to be traversed are traversed, comprising:

[0327] The next traversal loop is performed on the sampling points in the window after the preset step length is moved, and the traversal loop is repeated until it is determined that there are third sampling points and fourth sampling points that meet another preset condition, or the end point of the current window reaches at least the last sampling point in the second point sequence.

[0328] In an embodiment, the second preset time interval includes a window length of a window function. For a second point sequence corresponding to any candidate first high-frequency QRS complex data, in a traversal loop, any sampling point in the second point sequence is taken as a current reference sampling point, and sampling points in the second preset time interval from the current reference sampling point are traversed to determine whether there are third sampling points and fourth sampling points that meet another preset condition in the second preset time interval, including:

[0329] The first sampling point in the second point sequence that is earliest in time sequence is taken as a starting point of a window, and the sampling points in the second point sequence contained in the window are traversed to determine whether there are third sampling points and fourth sampling points that meet another preset condition in the current window.

[0330] In an embodiment, in a case where it is determined that there are third sampling points and fourth sampling points that meet another preset condition, it is determined that the second waveform feature indicating a myocardial bridge exists in the arbitrary candidate first high-frequency QRS complex data, thereby determining the possibility of the myocardial bridge, including:

[0331] In a case where it is determined that the third sampling points and the fourth sampling points in the current window meet another preset condition, it is determined that the second waveform feature indicating a myocardial bridge exists in the arbitrary candidate first high-frequency QRS complex data, thereby determining the possibility of the myocardial bridge.

[0332] In an embodiment, in a case where it is determined that there are no third sampling points and fourth sampling points that meet another preset condition, the current traversal loop ends, and other sampling points in the second point sequence except the current reference sampling point are determined as new reference sampling points, including:

[0333] In a case where it is determined that there are no third sampling points and fourth sampling points that meet another preset condition in the current window, the current traversal loop ends, and the window is moved backward by a preset step length in time sequence.

[0334] In an embodiment, a next traversal loop is performed for the new reference sampling point until it is determined that there are third sampling points and fourth sampling points that meet another preset condition, or all the sampling points in the second point sequence that need to be traversed are traversed, including:

[0335] The next traversal loop is performed on the sampling points in the window after the preset step size is moved, and the traversal loop is repeated until it is determined that there are third sampling points and fourth sampling points that meet another preset condition, or the end point of the current window reaches at least the last sampling point in the second point sequence.

[0336] In an embodiment, determining the possibility of the myocardial bridge according to the motion electrocardiogram data further comprises:

[0337] In a case where no third sampling point and fourth sampling point that meet another preset condition are found after all the sampling points in the second point sequence are traversed, it is determined that the corresponding candidate first high-frequency QRS complex data does not have the second waveform feature indicating the myocardial bridge.

[0338] In a case where all the candidate first high-frequency QRS complex data do not have the second waveform feature indicating the myocardial bridge, it is determined that there is no myocardial bridge.

[0339] In the embodiments of the present application, the other preset condition comprises:

[0340] The amplitude drop relative value between the third sampling point and the fourth sampling point reaches a fourth preset threshold, and the voltage drop degree between the third sampling point and the fourth sampling point reaches a fifth preset threshold, and the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point shows a continuous downward trend; or

[0341] The amplitude drop relative value between the third sampling point and the fourth sampling point reaches a fourth preset threshold, or the voltage drop degree between the third sampling point and the fourth sampling point reaches a fifth preset threshold, and the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point shows a continuous downward trend.

[0342] In an embodiment, the length of the second preset time interval is greater than the length of the first preset time interval.

[0343] In an embodiment, the fourth preset threshold is greater than the first preset threshold, and the fifth preset threshold is greater than the second preset threshold.

[0344] In an embodiment, in a case where it is determined that the amplitude drop relative value between the third sampling point and the fourth sampling point reaches the fourth preset threshold and / or the voltage drop degree between the third sampling point and the fourth sampling point reaches the fifth preset threshold, a first-order difference of the root mean square voltage of the sampling points in the second point sequence corresponding to the candidate first high-frequency QRS complex data in the current second preset time interval is calculated to obtain a difference sequence.

[0345] The number of consecutive values greater than or equal to zero in the difference sequence is determined.

[0346] If the number is less than or equal to a first number threshold, it is determined that the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point is in a sustained downward trend.

[0347] In an embodiment, all difference sequence groups are determined according to the difference sequence, each difference sequence group including a number of consecutive difference values;

[0348] The sum of the difference values of each difference sequence group is calculated.

[0349] If the number of consecutive values greater than or equal to zero in the difference sequence is less than or equal to a second number threshold, and the sum of the difference values of each difference sequence group is greater than or equal to a predetermined value, it is determined that the candidate first high-frequency QRS complex data has a waveform feature indicative of a myocardial bridge.

[0350] In an embodiment, determining the likelihood of the presence of a myocardial bridge according to the exercise electrocardiogram data includes:

[0351] Determining candidate first high-frequency QRS complex data from the corresponding first high-frequency QRS complex data of each electrocardiogram lead, the candidate first high-frequency QRS complex data not having a first waveform feature indicative of a coronary artery lesion;

[0352] Obtaining a corresponding second point sequence from the candidate first high-frequency QRS complex data, the second point sequence including a plurality of sampling points arranged in time sequence;

[0353] For any candidate first high-frequency QRS complex data, determining a second point sequence corresponding thereto, and determining at least one sampling point in a preset time period in the second point sequence as a reference sampling point;

[0354] Determining a to-be-traversed sampling point in the second point sequence according to the determined reference sampling point, wherein the to-be-traversed sampling point is a sampling point in a third preset time interval before the reference sampling point;

[0355] Traversing the to-be-traversed sampling point to determine whether a target sampling point exists, wherein the target sampling point and the reference sampling point satisfy another preset condition;

[0356] In a case where it is determined that the target sampling point exists, it is determined that the candidate first high-frequency QRS complex data has a second waveform feature indicative of a myocardial bridge, and thus the likelihood of the presence of a myocardial bridge is determined.

[0357] In an embodiment, the preset time period is a time period from the time when the exercise ends to a time after the exercise ends.

[0358] In an embodiment, the third preset time interval is located in the exercise phase of the load exercise electrocardiogram detection process.

[0359] In an embodiment, the third preset time interval includes a plurality of continuous sub-intervals, and different sub-intervals correspond to different preset conditions.

[0360] In an embodiment, determining the possibility of existence of the coronary microcirculation disorder according to the exercise electrocardiogram data and the resting electrocardiogram data includes:

[0361] determining whether the target first high-frequency QRS complex data exists in the first high-frequency QRS complex data corresponding to the first lead positive indicator indicating positive;

[0362] In a case where the target first high-frequency QRS complex data exists, determining the possibility of existence of the coronary microcirculation disorder.

[0363] In an embodiment, in a case where the target first high-frequency QRS complex data exists, determining the possibility of existence of the coronary microcirculation disorder includes:

[0364] determining a target electrocardiogram lead corresponding to the target first high-frequency QRS complex data;

[0365] determining whether the second lead positive indicator corresponding to the target electrocardiogram lead indicates positive;

[0366] In a case where the second lead positive indicator corresponding to the target electrocardiogram lead indicates positive, determining the possibility of existence of the coronary microcirculation disorder.

[0367] In an embodiment, determining the possibility of existence of the coronary microcirculation disorder according to the exercise electrocardiogram data and the resting electrocardiogram data further includes:

[0368] determining a first maximum root mean square voltage in the first high-frequency QRS complex data corresponding to the electrocardiogram lead and a second maximum voltage in the second high-frequency QRS complex data corresponding to the electrocardiogram lead;

[0369] determining a number of electrocardiogram leads corresponding to the second lead positive indicator indicating positive;

[0370] In a case where the number is less than a preset value, and at least the first maximum root mean square voltage is greater than a first voltage threshold, determining that the coronary microcirculation disorder does not exist.

[0371] In an embodiment, quantitatively evaluating the myocardial ischemia includes:

[0372] obtaining a first lead number according to the first lead positive indicator, and obtaining a second lead number according to the second lead positive indicator, wherein the first lead number is a number of electrocardiogram leads corresponding to the first lead positive indicator indicating positive or critical, and the second lead number is a number of electrocardiogram leads corresponding to the second lead positive indicator indicating positive or critical.

[0373] analyze the resting electrocardio data to determine a target high frequency morphology index;

[0374] normalize the first lead number, the second lead number, and the target high frequency morphology index to obtain a first ischemia index, a second ischemia index, and a third ischemia index, respectively;

[0375] obtain a target ischemia index according to the first ischemia index, the second ischemia index, and the third ischemia index.

[0376] In an embodiment, the target ischemia index is a product of the first ischemia index, the second ischemia index, and the third ischemia index.

[0377] The technical solution provided by the embodiments of the present application can more accurately qualitatively and quantitatively analyze the myocardial ischemia by combining the exercise electrocardio data, the resting electrocardio data, and the conventional ECG data. In addition, compared with only focusing on the steep drop wave in the first three minutes of exercise, the technical solution can improve the accuracy of coronary artery lesion qualification by judging whether the steep drop wave exists in the high frequency QRS complex data in the whole load exercise detection process, thereby improving the recognition sensitivity of the coronary artery lesion without significantly reducing or even not reducing the specificity, and providing a qualitative identification method for myocardial bridge and coronary microcirculation disorder. Therefore, the evaluation results of myocardial ischemia qualification and quantification, coronary artery lesion qualification and quantification, myocardial bridge qualification, and coronary microcirculation disorder can be comprehensively and accurately evaluated to more comprehensively and accurately evaluate the heart health status.

[0378] Those skilled in the art will understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can adopt a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can adopt a computer program product in the form of a computer program product embodied on one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program code.

[0379] The present application is described with reference to flowcharts and / or block diagrams according to the methods, devices (systems), and computer program products of the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to produce a machine, so that the instructions executed by the computer or other programmable data processing devices produce a device that implements the flowcharts and / or block diagrams. Figure 1 one flow or multiple flows and / or blocks Figure 1an apparatus to perform each block or blocks of the flow or flows and / or steps of the functions specified in the block or blocks.

[0380] These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a Figure 1 a flow or flows and / or steps of the functions specified in the block or blocks. Figure 1 an apparatus to perform each block or blocks of the flow or flows and / or steps of the functions specified in the block or blocks.

[0381] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions that are executed on the computer or other programmable apparatus provide steps for implementing the flow or flows and / or steps of the functions specified in the block or blocks. Figure 1 a flow or flows and / or steps of the functions specified in the block or blocks. Figure 1 an apparatus to perform each block or blocks of the flow or flows and / or steps of the functions specified in the block or blocks.

[0382] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0383] The memory can include non-persistent memory and / or volatile memory, such as random access memory (RAM) about which the computer stores information about an operating system, application software, and / or the like. Memory is an example of computer readable media.

[0384] Computer readable media includes permanent and non-permanent, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read only memory (ROM), electrically erasable programmable read only memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital versatile discs (DVDs) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible to a computing device. According to the definition herein, computer readable media does not include transitory media, such as modulated data signals and carrier waves.

[0385] It should also be noted that the terms "comprising", "comprises" or other variations 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 can also include other elements not expressly listed or inherent to such process, method, article, or apparatus.

[0386] The above embodiments are only used to illustrate the present application, but not to limit it. Instead of the above, various modifications and changes can be made to the application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the application shall fall into the scope of the claims of the application.

Claims

1. A cardiac health status assessment device, characterized by, A device configured to implement a cardiac health status assessment method, the cardiac health status assessment method comprising: obtaining exercise electrocardiogram data, resting electrocardiogram data, and routine electrocardiogram data, respectively, output by at least one electrocardiogram lead; determining a likelihood of presence of myocardial ischemia based on the exercise electrocardiogram data, the resting electrocardiogram data, and the routine electrocardiogram data; determining a likelihood of presence of coronary artery lesion based on the exercise electrocardiogram data; determining a likelihood of presence of myocardial bridge based on the exercise electrocardiogram data; determining a likelihood of presence of coronary microcirculation disorder based on the exercise electrocardiogram data and the resting electrocardiogram data; qualitatively assessing the cardiac health status based on the determination results of the likelihood of presence of myocardial ischemia, the likelihood of presence of coronary artery lesion, the likelihood of presence of myocardial bridge, and the likelihood of presence of coronary microcirculation disorder; in a case where the likelihood of presence of myocardial ischemia and / or the likelihood of presence of coronary artery lesion is determined, quantitatively assessing the myocardial ischemia and / or the coronary artery lesion; and assessing the cardiac health status based on the qualitative assessment result and the quantitative assessment result; wherein the determining the likelihood of presence of coronary artery lesion based on the exercise electrocardiogram data comprises: obtaining, from the first high-frequency QRS complex data corresponding to each electrocardiogram lead derived from the exercise electrocardiogram data, a first point sequence corresponding to each electrocardiogram lead, respectively, the first point sequence comprising a plurality of sampling points arranged in time sequence; for the first point sequence corresponding to any electrocardiogram lead, in an iteration loop, taking any sampling point in the plurality of sampling points in the first point sequence as a current reference sampling point, and traversing the sampling points within a first preset time interval from the current reference sampling point to determine whether there are first and second sampling points satisfying a preset condition within the first preset time interval, wherein the first sampling point is earlier in time sequence than the second sampling point; in a case where the first and second sampling points satisfying the preset condition are determined, determining that the first high-frequency QRS complex data corresponding to the any electrocardiogram lead has a first waveform feature indicative of coronary artery lesion, thereby determining the likelihood of presence of coronary artery lesion; in a case where the first and second sampling points satisfying the preset condition are not determined, ending the current iteration loop, and determining other sampling points in the first point sequence except the current reference sampling point as new reference sampling points; performing the next iteration loop for the new reference sampling points until the first and second sampling points satisfying the preset condition are determined, or all the sampling points in the first point sequence that need to be traversed are traversed; in a case where no first and second sampling points satisfying the preset condition are found after all the sampling points in the first point sequence that need to be traversed are traversed, determining that the first high-frequency QRS complex data corresponding to the any electrocardiogram lead does not have the first waveform feature indicative of coronary artery lesion; in a case where the first high-frequency QRS complex data corresponding to all the electrocardiogram leads do not have the first waveform feature indicative of coronary artery lesion, determining the likelihood of absence of coronary artery lesion; the preset condition comprises: the amplitude drop relative value between the first sampling point and the second sampling point reaches a first preset threshold and the voltage drop degree between the first sampling point and the second sampling point reaches a second preset threshold; or the amplitude drop relative value between the first sampling point and the second sampling point reaches the first preset threshold or the voltage drop degree between the first sampling point and the second sampling point reaches the second preset threshold.

2. The cardiac health status evaluation apparatus according to claim 1, characterized by The determining the possibility of myocardial ischemia according to the exercise ECG data, the resting ECG data and the conventional ECG data comprises: obtaining a first lead positive indicator corresponding to each ECG lead according to first high-frequency QRS complex data obtained from the exercise ECG data; obtaining a second lead positive indicator corresponding to each ECG lead according to second high-frequency QRS complex data obtained from the resting ECG data; obtaining a third lead positive indicator corresponding to each ECG lead according to the conventional ECG data; in a case where any one of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicates positivity, determining the possibility of myocardial ischemia.

3. The heart health status evaluation device according to claim 2, wherein, in a case where only one of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicates positivity, the possibility of myocardial ischemia is low; in a case where two of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicate positivity, the possibility of myocardial ischemia is high; in a case where all of the first lead positive indicator, the second lead positive indicator and the third lead positive indicator indicate positivity, the possibility of myocardial ischemia is the highest.

4. The cardiac health status evaluation apparatus according to claim 1, characterized by The determining the possibility of myocardial bridge according to the exercise ECG data comprises: determining candidate first high-frequency QRS complex data from the first high-frequency QRS complex data corresponding to each ECG lead, the candidate first high-frequency QRS complex data being free of a first waveform feature indicative of coronary artery lesion; obtaining a second point sequence corresponding to the candidate first high-frequency QRS complex data according to the candidate first high-frequency QRS complex data, the second point sequence comprising a plurality of sampling points arranged in time sequence; for the second point sequence corresponding to any candidate first high-frequency QRS complex data, taking any one of the plurality of sampling points in the second point sequence as a current reference sampling point in an iteration loop, and iterating through the sampling points within a second preset time interval from the current reference sampling point to determine whether there is a third sampling point and a fourth sampling point satisfying another preset condition within the second preset time interval, wherein the third sampling point is earlier in time sequence than the fourth sampling point; in a case where it is determined that there is a third sampling point and a fourth sampling point satisfying another preset condition, determining that the any candidate first high-frequency QRS complex data has a second waveform feature indicative of myocardial bridge, thereby determining the possibility of myocardial bridge. In a case where it is determined that there are no third sampling point and fourth sampling point satisfying another preset condition, the current traversal loop ends, and other sampling points in the second point sequence except the current reference sampling point are determined as new reference sampling points; The next traversal loop is performed for the new reference sampling points until it is determined that there are third sampling point and fourth sampling point satisfying another preset condition, or all sampling points in the second point sequence that need to be traversed are traversed.

5. The cardiac health status evaluation apparatus according to claim 4, characterized by The determining the possibility of the existence of the myocardial bridge according to the motion electrocardio data comprises: In a case where all sampling points in the second point sequence that need to be traversed are traversed and there are no third sampling point and fourth sampling point satisfying another preset condition, it is determined that the corresponding candidate first high-frequency QRS complex data does not have the second waveform feature indicating the myocardial bridge; In a case where all candidate first high-frequency QRS complex data do not have the second waveform feature indicating the myocardial bridge, it is determined that there is no myocardial bridge.

6. The cardiac health status evaluation apparatus according to claim 4, characterized by The another preset condition comprises: The amplitude drop relative value between the third sampling point and the fourth sampling point reaches a third preset threshold value, the voltage drop degree between the third sampling point and the fourth sampling point reaches a fourth preset threshold value, and the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point shows a continuous downward trend; or The amplitude drop relative value between the third sampling point and the fourth sampling point reaches a third preset threshold value or the voltage drop degree between the third sampling point and the fourth sampling point reaches a fourth preset threshold value, and the root mean square voltage of each sampling point between the third sampling point and the fourth sampling point shows a continuous downward trend.

7. The cardiac health status evaluation apparatus according to claim 3, characterized by The determining the possibility of the existence of the coronary microcirculation disorder according to the motion electrocardio data and the resting electrocardio data comprises: Determining whether there is target first high-frequency QRS complex data in the first high-frequency QRS complex data corresponding to the first lead positive indicator indicating positive, the target first high-frequency QRS complex data not having the first waveform feature and the second waveform feature; In a case where it is determined that there is target first high-frequency QRS complex data, it is determined that there is a possibility of the existence of the coronary microcirculation disorder.

8. The cardiac health status evaluation apparatus according to claim 7, characterized by The determining the possibility of the existence of the coronary microcirculation disorder in a case where it is determined that there is target first high-frequency QRS complex data comprises: Determining a target electrocardiogram lead corresponding to the target first high-frequency QRS complex data; Determining whether the second lead positive indicator corresponding to the target electrocardiogram lead indicates positive; In a case where the second lead positive indicator corresponding to the target electrocardiogram lead indicates positive, it is determined that there is a possibility of the existence of the coronary microcirculation disorder.

9. The cardiac health status evaluation apparatus according to claim 3, characterized by The determining the possibility of the existence of the coronary microcirculation disorder according to the motion electrocardio data and the resting electrocardio data further comprises: Determining a first maximum root mean square voltage according to the root mean square voltage in the first high-frequency QRS complex data corresponding to all electrocardiogram leads and a second maximum voltage according to the voltage in the second high-frequency QRS complex data corresponding to all electrocardiogram leads; Determining the number of electrocardiogram leads corresponding to the second lead positive indicator indicating positive; In a case where the number of electrocardiogram leads is less than a preset value and at least the first maximum root mean square voltage is greater than a first voltage threshold value, it is determined that there is no coronary microcirculation disorder.

10. The cardiac health status evaluation apparatus according to claim 9, characterized by The quantitative assessment of myocardial ischemia comprises: obtaining a first lead number according to the first lead positive indicator, and obtaining a second lead number according to the second lead positive indicator, wherein the first lead number is the number of electrocardiogram leads corresponding to positive or critical indicated by the first lead positive indicator, and the second lead number is the number of electrocardiogram leads corresponding to positive or critical indicated by the second lead positive indicator; analyzing the resting electrocardio data to determine a target high frequency morphology index; normalizing the first lead number, the second lead number and the target high frequency morphology index to obtain a first ischemia index, a second ischemia index and a third ischemia index respectively; obtaining a target ischemia index according to the first ischemia index, the second ischemia index and the third ischemia index.

11. The cardiac health status evaluation apparatus according to claim 10, characterized by The target ischemia index is the product of the first ischemia index, the second ischemia index and the third ischemia index, or the minimum one of the first ischemia index, the second ischemia index and the third ischemia index is taken as the target ischemia index.

12. A cardiac health status assessment system, characterized by, comprise: an electrocardio signal acquisition device comprising at least one electrode forming at least one electrocardiogram lead for acquiring electrocardio signals of a subject; and a heart health status evaluation apparatus according to any one of claims 1 to 11. The machine readable storage medium has instructions stored thereon for causing a machine to perform the heart health status evaluation method according to any one of claims 1 to 11.

13. A machine-readable storage medium, characterized in that, ​

Images