US-12616384-B2 - Cardiac diastolic function assessment method, device and system

US12616384B2US 12616384 B2US12616384 B2US 12616384B2US-12616384-B2

Abstract

Disclosed is a cardiac diastolic function assessment method applicable to the field of cardiac monitoring. The method comprises: acquiring vibration information of the thoracic cavity body surface of an object in a noninvasive manner; preprocessing the vibration information to generate hemodynamic-related information; determining a target wave group on the basis of the hemodynamic-related information; determining the highest peak on the target wave group, determining a rising edge amplitude before the highest peak as a first characteristic value, and determining, as a second characteristic value, an amplitude between the highest peak and the subsequent lowest valley; and generating an indicating parameter on the basis of the first characteristic value and the second characteristic value, and assessing a cardiac diastolic function of the object on the basis of the indicating parameter.

Inventors

Pengbo Liu
Shaochun Zhuang
Lingjun Zeng
ZHENGPEI CHU

Assignees

Shenzhen Darma Technology Co., Ltd.

Dates

Publication Date: 20260505
Application Date: 20190520

Claims (13)

1 . A diastolic function assessment method, performed by one or more processors executing one or more computer programs stored in a memory, comprising steps of: non-invasively and continuously acquiring vibration information on a body surface corresponding to a thoracic cavity of a subject in a supine position through one or more fiber-optic sensors which are connected to the one or more processors; wherein the one or more fiber-optic sensors are configured to be placed under the subject's right shoulder and around the subject's right shoulder scapula, a sensing area of the one or more fiber-optic sensors is at least 20 square centimeters and covers the body surface area of the right shoulder scapula; optical fibers are distributed in the sensing area of the one or more fiber-optic sensors; the one or more fiber-optic sensors are sensitive to changes in vibration displacement; the vibration information contains breathing signals, hemodynamic signals, body motion signals and noise signals; a waveform of the vibration information includes breathing envelopes generated by the breathing signals; the hemodynamic signals, the body motion signals and the noise signals are superposed on the breathing envelopes; a horizontal axis of the waveform represents time, and a vertical axis represents normalized vibration information which is dimensionless; preprocessing the vibration information to generate hemodynamic related information; comprising: filtering the vibration information below 2 Hz to remove breathing signals and body motion signals and filtering the vibration information above 45 Hz to remove the noise signals, thereby generating the hemodynamic related information of 2-45 Hz; performing a first-order differential processing on the hemodynamic related information to generate a first-order differential information with a frequency band of 2-45 Hz, performing a second-order differential processing on the hemodynamic related information to generate a second-order differential information with a frequency band of 2-45 Hz; generating, by performing energy integration on the hemodynamic related information, vibration energy information comprising a first energy envelope and a second energy envelope which represent energy accumulation during a systolic process and an early diastole of the subject's heart; synchronizing the hemodynamic related information, or the first-order differential information, or the second-order differential information, and the vibration energy information on the same time axis, and performing heartbeat segmentation; determining, in one cardiac cycle, a highest peak of the hemodynamic related information, a highest peak of the first-order differential information, or a highest peak of the second-order differential information, wherein the highest peak of the hemodynamic related information or the first-order differential information or the second-order differential information represents a shock caused by blood flowing into an aortic arch after aortic ejection; determining a target time window; wherein the first energy envelope contains the highest peak of the hemodynamic related information, or the highest peak of the first-order differential information, or the highest peak of the second-order differential information, while the second energy envelope does not have the highest peak of the hemodynamic related information, or the highest peak of the first-order differential information, or the highest peak of the second-order differential information; determine a time duration of the second energy envelope as the target time window; determining, wave clusters within the target time window of the second-order differential information, as a target wave group; determining a highest peak on the target wave group of the second-order differential information; determining a rising edge amplitude before the highest peak on the target wave group as a first characteristic value, and determining an amplitude between the same highest peak on the target wave group and a subsequent lowest valley on the target wave group as a second characteristic value; generating an indicating parameter for assessing a diastolic function of the subject, comprising: determining a ratio of the second characteristic value to the first characteristic value as the indicating parameter; and determining an elevated filling pressure state if the indicating parameter is greater than a threshold; wherein the threshold is 1.405 or depends on a specific people group; and outputting, through an output device, a result of the assessing.
2 . The method of claim 1 , wherein the hemodynamic related information is: data in one cardiac cycle; or data that is superimposed and averaged in a unit of cardiac cycle within a preset time period.
3 . The method of claim 1 , wherein the step of performing heartbeat segmentation comprises: performing the highest peak or a lowest valley search on the hemodynamic related information, or the first-order differential information, or the second-order differential information at a search interval between 0.6 seconds and 1 second; and performing the heartbeat segmentation based on the repetitive highest peaks or the repetitive lowest valleys on the hemodynamic related information or the first-order differential information or the second-order differential information.
4 . The method of claim 1 , wherein the step of performing heartbeat segmentation comprises: acquiring electrocardiography (ECG) data through an electrocardiogram sensor when acquiring the vibration information; and performing heartbeat segmentation on the hemodynamic related information or the first-order differential information or the second-order differential information based on the ECG data.
5 . The method of claim 1 , wherein the step of filtering the vibration information uses one or more of low-pass filtering, band-pass filtering, Infinite Impulse Response (IIR) filtering, Finite Impulse Response (FIR) filtering, wavelet filtering, zero-phase bidirectional filtering, polynomial smoothing filtering, integral transformation, and differential transformation, to filter the vibration information at least once to generate the hemodynamic related information.
6 . A diastolic function assessment method, performed by one or more processors executing one or more computer programs stored in a memory, comprising steps of: non-invasively and continuously acquiring the vibration information on the body surface corresponding to a thoracic cavity of a subject in a supine position through the one or more fiber-optic sensors configured to be placed under the subject's right shoulder and around a right shoulder scapula, a sensing area of the one or more fiber-optic sensors is at least 20 square centimeters and covers the body surface area of the right shoulder scapula of the subject; optical fibers are distributed in the sensing area of the one or more fiber-optic sensors; preprocessing the vibration information to generate hemodynamic related information with a frequency band of 2-45 Hz; comprising: filtering, noise removal and signal scaling; performing a first-order differential processing on the hemodynamic related information to generate a first-order differential information with a frequency band of 2-45 Hz, or, performing a second-order differential processing on the hemodynamic related information to generate a second-order differential information with a frequency band of 2-45 Hz; generating, by performing energy integration on the hemodynamic related information, vibration energy information comprising a first energy envelope and a second energy envelope which represent energy accumulation during a systolic process and an early diastole of the subject's heart; synchronizing the hemodynamic related information, or the first-order differential information, or the second-order differential information, and the vibration energy information on the same time axis, and performing heartbeat segmentation; determining, in one cardiac cycle, a highest peak of the hemodynamic related information, or a highest peak of the first-order differential information, or a highest peak of the second-order differential information, wherein the highest peak of the hemodynamic related information or the first-order differential information or the second-order differential information represents a shock caused by blood flowing into an aortic arch after aortic ejection; determining a target time window; wherein the first energy envelope contains the highest peak of the hemodynamic related information, or the highest peak of the first-order differential information, or the highest peak of the second-order differential information, while the second energy envelope does not have the highest peak of the hemodynamic related information, or the highest peak of the first-order differential information, or the highest peak of the second-order differential information; determine a time duration of the second energy envelope as the target time window; and determining wave clusters within the target time window of the first-order differential information as a target wave group; determining a highest peak on the target wave group of the first-order differential information; determining a rising edge amplitude before the highest peak on the target wave group as a first characteristic value, and determining an amplitude between the same highest peak on the target wave group and a subsequent lowest valley on the target wave group as a second characteristic value; generating an indicating parameter for assessing a diastolic function of the subject, comprising: determining a ratio of the second characteristic value to the first characteristic value as the indicating parameter; and determining an elevated filling pressure state if the indicating parameter is greater than a threshold; wherein the threshold is 1.294 or depends on a specific people group; and outputting, through an output device, a result of the assessing.
7 . The diastolic function assessment method of claim 6 , wherein the step of performing heartbeat segmentation comprises: performing the highest peak or a lowest valley search on the hemodynamic related information, or the first-order differential information, or the second-order differential information at a search interval between 0.6 seconds and 1 second; and performing the heartbeat segmentation based on the repetitive highest peaks or the repetitive lowest valleys on the hemodynamic related information or the first-order differential information or the second-order differential information.
8 . The diastolic function assessment method of claim 6 , wherein the step of performing heartbeat segmentation comprises: acquiring electrocardiography (ECG) data through an electrocardiogram sensor when acquiring the vibration information; and performing heartbeat segmentation on the hemodynamic related information or the first-order differential information or the second-order differential information based on the ECG data.
9 . The diastolic function assessment method of claim 6 , wherein the step of filtering uses one or more of low-pass filtering, band-pass filtering, Infinite Impulse Response (IIR) filtering, Finite Impulse Response (FIR) filtering, wavelet filtering, zero-phase bidirectional filtering, polynomial smoothing filtering, integral transformation, and differential transformation, to filter the vibration information at least once to generate the hemodynamic related information.
10 . A diastolic function assessment system based on machine learning, comprising: one or more processors, which are programmed to perform the steps of: receiving vibration information on a body surface corresponding to a thoracic cavity of a subject in a supine position through one or more fiber-optic sensors which are connected to the one or more processors as input information for training; analyzing the input information for training to establish an assessment model by machine learning, and; receiving the vibration information on the body surface corresponding to the subject's thoracic cavity through the one or more fiber-optic sensors; and performing an assessment to the subject's diastolic function by the assessment model; wherein the one or more fiber-optic sensors are configured to be placed under the subject's right shoulder and around the subject's right shoulder scapula, a sensing area of the one or more fiber-optic sensors is at least 20 square centimeters and covers the body surface area of the right shoulder scapula; optical fibers are distributed in the sensing area of the one or more fiber-optic sensors; the one or more fiber-optic sensors are sensitive to changes in vibration displacement; the vibration information contains breathing signals, hemodynamic signals, body motion signals and noise signals; a waveform of the vibration information includes breathing envelopes generated by the breathing signals; the hemodynamic signals, the body motion signals and the noise signals are superposed on the breathing envelopes; a horizontal axis of the waveform represents time, and a vertical axis represents normalized vibration information which is dimensionless; wherein the assessment model performs steps of: preprocessing the vibration information to generate hemodynamic related information; comprising: filtering the vibration information below 2 Hz to remove breathing signals and body motion signals and filtering the vibration information above 45 Hz to remove the noise signals, thereby generating the hemodynamic related information of 2-45 Hz; performing a first-order differential processing on the hemodynamic related information to generate a first-order differential information, or, performing a second-order differential processing on the hemodynamic related information to generate a second-order differential information; generating, by performing energy integration on the hemodynamic related information, vibration energy information comprising a first energy envelope and a second energy envelope which represent energy accumulation during a systolic process and an early diastole of the subject's heart; synchronizing the hemodynamic related information, or the first-order differential information, or the second-order differential information, and the vibration energy information on the same time axis, and performing heartbeat segmentation; determining, in one cardiac cycle, a highest peak of the hemodynamic related information, or a highest peak of the first-order differential information, or a highest peak of the second-order differential information, wherein the highest peak of the hemodynamic related information or the first-order differential information or the second-order differential information represents a shock caused by blood flowing into an aortic arch after aortic ejection; determining a target time window; wherein the first energy envelope contains the highest peak of the hemodynamic related information, or the highest peak of the first-order differential information, or the highest peak of the second-order differential information, while the second energy envelope does not have the highest peak of the hemodynamic related information, or the highest peak of the first-order differential information, or the highest peak of the second-order differential information; determine a time duration of the second energy envelope as the target time window; and determining, wave clusters within the target time window of the first-order differential information or the second-order differential information, as a target wave group; determining a highest peak on the target wave group of the first-order differential information or the second-order differential information; determining a rising edge amplitude before the highest peak on the same target wave group as a first characteristic value and determining an amplitude between the same highest peak on the same target wave group and a subsequent lowest valley on the same target wave group as a second characteristic value; and generating an indicating parameter for assessing a diastolic function of the subject, comprising: determining a ratio of the second characteristic value to the first characteristic value as the indicating parameter; and determining an elevated filling pressure state if the indicating parameter is greater than a threshold; wherein the threshold is 1.294 when determining the first characteristic value and the second characteristic value based on the first-order differential information; and the threshold is 1.405 when determining the first characteristic value and the second characteristic value based on the second-order differential information; or the threshold depends on a specific people group.
11 . The diastolic function assessment system based on machine learning of claim 10 , wherein the step of performing heartbeat segmentation comprises: performing the highest peak or a lowest valley search on the hemodynamic related information, or the first-order differential information, or the second-order differential information at a search interval between 0.6 seconds and 1 second; and performing the heartbeat segmentation based on the repetitive highest peaks or the repetitive lowest valleys on the hemodynamic related information or the first-order differential information or the second-order differential information.
12 . The diastolic function assessment system based on machine learning of claim 10 , wherein the step of performing heartbeat segmentation comprises: acquiring electrocardiography (ECG) data through an electrocardiogram sensor when acquiring the vibration information; and performing heartbeat segmentation on the hemodynamic related information or the first-order differential information or the second-order differential information based on the ECG data.
13 . The diastolic function assessment system based on machine learning of claim 10 , wherein the step of filtering the vibration information uses one or more of low-pass filtering, band-pass filtering, Infinite Impulse Response (IIR) filtering, Finite Impulse Response (FIR) filtering, wavelet filtering, zero-phase bidirectional filtering, polynomial smoothing filtering, integral transformation, and differential transformation, to filter the vibration information at least once to generate the hemodynamic related information.

Description

CROSS REFERENCE TO RELATED APPLICATIONS The present application is a 35 U.S.C. § 371 National Phase conversion of International (PCT) Patent Application No. PCT/CN2019/087642, filed on May 20, 2019, the disclosure of which is incorporated by reference herein. The PCT International Patent Application was filed and published in Chinese. FIELD OF THE INVENTION The present invention relates to the field of cardiac monitoring, and particularly relates to a non-invasive diastolic function assessment method, device and system. BACKGROUND OF THE INVENTION Heart failure (abbreviated as HF) is a clinical syndrome with multiple etiologies and pathogenesis. With the aging of the population and an increasing survival rate of patients with acute myocardial infarction, the number of patients with chronic heart failure is increasing rapidly. Patients with heart failure suffer from a chronic state to an acute worsening state, and suffer from an accompanied deterioration of diastolic function, such as an elevated filling pressure. Elevated filling pressure will cause the heart's function to enter a rapid vicious circle, but the patient itself will not feel the symptoms until the filling pressure continues to rise for about 20 days and need to be admitted to the hospital urgently; while at this time, the impairment of the heart function is caused and is irreversible. When the patient is identified in an elevated filling pressure status, timely intervention is required to avoid further deterioration. This has become the consensus of clinicians. At present, there are implantable products used to evaluate the diastolic function, but the cost is relatively high, and if it is only used for monitoring, patients are less likely to accept. Therefore, a more friendly and more convenient product is needed for monitoring the diastolic function. SUMMARY OF THE INVENTION Technical Problem An object of the present invention is to provide a method, device, system, and computer-readable storage medium for accessing a cardiac diastolic function of a subject. Solutions to the Problem Technical Solutions In a first aspect, the present invention provides a cardiac diastolic function assessment method, comprising steps of: acquiring vibration information on a body surface corresponding to a subject's thoracic cavity in a noninvasive manner;preprocessing the vibration information to generate hemodynamic related information;determining a target wave group based on the hemodynamic related information;determining the highest peak on the target wave group; determining a rising edge amplitude before the highest peak as a first characteristic value; and determining an amplitude between the highest peak and the subsequent lowest valley as a second characteristic value; andgenerating an indicating parameter based on the first characteristic value and the second characteristic value; and assessing a diastolic function of the subject based on the indicating parameter; In a second aspect, the present invention provides a computer-readable storage medium having computer programs stored thereon, which when being executed by a processor, cause the processor to perform the steps of the above-mentioned cardiac diastolic function assessment method. In a third aspect, the present invention provides a diastolic function assessment device, comprising: one or more processors; a memory; and one or more computer programs, wherein the one or more computer programs are stored in the memory, and configured to be executed by the one or more processors; and the one or more processors execute the one or more computer programs to perform the steps of the above-mentioned diastolic function assessment method. In a fourth aspect, the present invention provides a cardiac diastolic function assessment system, comprising: one or more vibration sensors for acquiring vibration information on a body surface corresponding to a subject's thoracic cavity surface; andthe diastolic function assessment device, as described above, connected to the one or more vibration sensors. Advantages of the Preset Invention Advantages The method of the present invention monitors the diastolic function by acquiring the vibration information of the subject without intruding his body, it is a passively measuring, and can realize continuous monitoring. The subject only needs to lie on the measuring device to perform the measurement, and no need for professional assistance. The method has the advantages of high measurement accuracy and simple operation, can improve the comfort of the tester, and can be applied to scenes such as hospitals and homes. The diastolic function assessment system provided in the present invention can evaluate the diastolic function of the subject, and then prompt a warning in advance when deterioration appear, so as to help the subject avoid deterioration. BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a flowchart of a diastolic function assessment method in accordance with a first embo