WO2022248355A1 - Noninvasive measurement of left ventricular compliance - Google Patents
Noninvasive measurement of left ventricular compliance Download PDFInfo
- Publication number
- WO2022248355A1 WO2022248355A1 PCT/EP2022/063723 EP2022063723W WO2022248355A1 WO 2022248355 A1 WO2022248355 A1 WO 2022248355A1 EP 2022063723 W EP2022063723 W EP 2022063723W WO 2022248355 A1 WO2022248355 A1 WO 2022248355A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- left ventricular
- temporal variation
- heart
- volume
- processing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
- A61B8/0883—Clinical applications for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/04—Measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/486—Diagnostic techniques involving arbitrary m-mode
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
- A61B8/065—Measuring blood flow to determine blood output from the heart
Definitions
- the present invention relates to the field of heart monitoring, and in particular, to the field of measuring left ventricular compliance of a subject’s heart.
- Left ventricular compliance is a measurement that can be used of an indicator of likelihood of diastolic heart failure.
- left ventricular compliance is defined as the differential of a volume-pressure curve (dV/dP) within the left ventricle of the patient’s heart during a diastolic phase of the patient’s heart.
- dV/dP volume-pressure curve
- myocardial stiffness/compliance proposes to measure myocardial stiffness/compliance through direct elastographic measurement.
- myocardial compliance and ventricular compliance are related, they are not the same; as myocardial compliance is an intrinsic local tissue property whereas ventricular compliance can depend on other factors such as chamber shape.
- a computer-implemented method for performing non-invasive estimation of left ventricular compliance for a subject is provided.
- the computer-implemented method comprises: obtaining ultrasound data of the subject’s heart, wherein the ultrasound data was captured over a time period comprising at least a portion of a cardiac cycle of the heart and comprises velocity data along a direction of mitral inflow; determining, by processing the velocity data, a velocity component along the direction of mitral inflow of the heart over the time period; determining a temporal variation of the left ventricular pressure from the determined velocity component by processing the velocity component using a physical model; determining, by processing the obtained ultrasound data, a temporal variation of the left ventricle volume; and estimating the left ventricular compliance by processing the determined temporal variation of the left ventricle volume and the determined temporal variation of the left ventricular pressure.
- the present disclosure therefore provides an approach for estimating (a value or function of) left ventricular compliance by processing only ultrasound data.
- a non-invasive procedure is proposed for predicting left ventricular compliance.
- a non-invasive procedure reduces a likelihood of adverse effects (e.g. reduces a risk of infection), improves a speed of obtaining a measure of left ventricular compliance (e.g. as no surgical intervention is required) and/or can be performed by fewer or less experienced clinicians than an invasive mechanism.
- the present disclosure effectively proposes to determine a time differential (“temporal variation”) of left ventricular pressure and a time differential (“temporal variation”) of left ventricular volume. These two features are then processed to predict a left ventricular compliance.
- a temporal variation of left ventricular pressure and a temporal variation of left ventricular volume can both be obtained through appropriate processing of ultrasound data (of a subject’s heart). Moreover, it is recognized that these two features can be processed to calculate left ventricular compliance.
- an underlying recognition of the present disclosure is that an approach for determining left ventricular compliance can be (re)cast or formulated into steps that can be executed by appropriate processing of ultrasound data (alone).
- the step of estimating the left ventricular compliance may comprise determining the ratio between the temporal variation of the left ventricle volume and the temporal variation of the left ventricular pressure.
- the left ventricular compliance may be estimated by dividing the temporal variation of the left ventricle volume by the temporal variation of the left ventricular pressure.
- the step of determining the temporal variation of the left ventricular pressure may comprise processing the velocity component using the Navier-Stokes equation to determine the temporal variation of the left ventricular pressures.
- the Navier-Stokes equation can thereby act as the physical model, and provide a mechanism for deriving a temporal variation of left ventricular pressure from one or more velocity components of velocity data of the ultrasound data.
- the step of determining the temporal variation of the left ventricle volume may comprise determining the temporal variation of the left ventricle volume by processing the velocity component. It is recognized that a temporal variation of left ventricular volume is responsive to (i.e. changes in response to changes of) a velocity of blood entering the left ventricle (i.e. via the mitral valve). This means it is possible to calculate, determine or predict the temporal variation of the left ventricle volume based on a velocity component of blood flow through the mitral valve.
- the step of determining the temporal variation of the left ventricle volume may comprise multiplying the velocity component by a surface area value, the surface area value representing a surface area size of the mitral annulus of the subject’s heart.
- the computer-implemented method further comprises a step of determining the surface area value, said step comprising: determining, by processing the ultrasound data, a first volume of the left ventricle, of the subject’s heart wherein the first volume is a volume of the left ventricle of the subject’s heart at a first point during the portion of the cardiac cycle of the subject’s heart; determining, by processing the ultrasound data, a second volume of the left ventricle, of the subject’s heart wherein the second volume is a volume of the left ventricle of the subject’s heart at a second, later point during the portion of the cardiac cycle of the subject’s heart; determining the difference between the second volume and the first volume; determining the integral of the velocity component between the first point and the second point during the portion of the cardiac cycle of the subject’s heart; and determining the ratio between the difference between the second volume and the first volume and the integral, to thereby calculate the surface area value.
- the computer-implemented method is adapted wherein: the ultrasound data was captured over a time period comprising at least a portion of a cardiac cycle of the heart that includes an end-systolic phase of the cardiac cycle and an end-diastolic phase of the cardiac cycle, the first point is a point in the end-systolic phase of the cardiac cycle; and the second point is a point in the end-diastolic phase of the cardiac cycle.
- the surface area value may be predetermined.
- the surface area value may be a value that represents a (global) average size for a surface area of the mitral valve.
- the ultrasound data is ultrasound data that was captured over a time period comprising at least a portion of a cardiac cycle of the heart that includes an end-diastolic phase of the cardiac cycle;
- the temporal variation of the left ventricle volume is a temporal variation of the left ventricle volume during the end-diastolic phase of the cardiac cycle;
- the temporal variation of the left ventricular pressure is a temporal variation of the left ventricular pressure during the end-diastolic phase of the cardiac cycle.
- the computer-implemented method may be adapted wherein: the step of determining a temporal variation of the left ventricle volume comprises processing the ultrasound data to determine a first set of measurements, each measurement in the first set of measurements representing the temporal variation of the left ventricle volume at a different, respective time point during the time period over which the ultrasound data was captured; the step of determining a temporal variation of the left ventricular pressure comprises processing the at least one velocity component using a physical model to determine a second set of measurements, each measurement in the second set of measurements corresponding to a measurement of the first set of measurements and representing the temporal variation of the left ventricular pressure at the same time point as the corresponding measurement of the first set of measurements; and the step of estimating the left ventricular compliance comprises processing the first set of measurements and the second set of measurements to estimate the left ventricular compliance.
- the step of estimating the left ventricular compliance comprises: multiplying each measurement of the first set of measurements by the corresponding measurement in the second set of measurements, to produce a multiplied set of measurements; averaging the multiplied set of measurements to produce an averaged multiplied measurement; squaring each measurement of the second set of measurements to produce a squared set of measurements; averaging the squared set of measurements to produce an averaged squared measurement; and estimating the left ventricular compliance by determining a ratio between the averaged multiplied measurement and the averaged squared measurement.
- the computer program product may be a (non-transitory) computer storage medium.
- the processing arrangement is configured to: obtain ultrasound data of the subject’s heart, wherein the ultrasound data was captured over a time period comprising at least a portion of a cardiac cycle of the heart and comprises velocity along a direction of mitral inflow; determine, by processing the velocity, a velocity component along the direction of mitral inflow of the heart over the time period; determine, by processing the obtained ultrasound data, a temporal variation of the left ventricle volume; determine a temporal variation of the left ventricular pressure from the determined velocity component by processing the velocity component using a physical model; and estimate the left ventricular compliance by processing the determined temporal variation of the left ventricle volume and the determined temporal variation of the left ventricular pressure.
- an imaging system comprising: the processing arrangement previously described; and an ultrasound system configured to capture ultrasound data of the subject’s heart and provide the ultrasound data to the processing arrangement.
- processing system comprising: the processing arrangement or the imaging system; and an output interface configured to provide a user-perceptible output responsive to the left ventricular compliance estimated by the processing arrangement.
- Fig. 1 is a flowchart illustrating a method according to an embodiment
- Fig. 2 provides a visual representation of functions derived from ultrasound data
- Fig. 3 is a graphical illustration of a function of left ventricular pressure
- Fig. 4 is a graphical illustration of a function of the time derivative of left ventricular volume
- Fig. 5 illustrates a processing arrangement
- the invention provides a mechanism for performing a (non-invasive) estimate of left ventricular compliance.
- a temporal variation of left ventricular volume and a temporal variation of left ventricular pressure are determined by processing ultrasound data of a subject. These two temporal variations are then processed to estimate left ventricular compliance.
- Embodiments are based on the realization that a measure for left ventricular compliance of a subject can be derived from measurements that can be estimated from ultrasound data of the subject’s heart.
- appropriate processing of ultrasound data facilitates determination of a measurement of left ventricular compliance.
- Fig. 1 is a flowchart illustrating a (computer-implemented) method 100 according to an embodiment.
- the illustrated method 100 is non-invasive, e.g. as it makes use of ultrasound data alone.
- the method generates a predicted left ventricular compliance, which may be in the form of a function (e.g. a function of time) and/or one or more discrete values/measurements of a left ventricular compliance.
- the method 100 comprises a step 110 of obtaining ultrasound data 101 of the subject’s heart.
- the ultrasound data is data that has been captured over a time period comprising at least a portion of a cardiac cycle of the heart.
- the ultrasound data may be obtained, for example, directly from an ultrasound imaging device or from a memory (e.g. that stores data generated by the ultrasound imaging device).
- the ultrasound data comprises velocity data along a direction of mitral inflow.
- Velocity data or “blood velocity data”, is any ultrasound data that can be processed to identify a velocity of blood flow in at least the direction of mitral inflow (i.e. blood flow through the mitral valve).
- Suitable examples of velocity data includes (color) M-mode acquisition data, i.e. a (color) M-mode trace, and/or other Doppler ultrasound data.
- the method 100 performs a step 120 of processing the velocity data to determine one or more velocity components along the direction of mitral inflow of the heart over the time period.
- Suitable velocity components include an inertial component and a convective component.
- a velocity component may be a velocity component as a function of (at least) time.
- velocity data may provide a blood velocity as a function of time and distance, i.e. a blood velocity function v r (r t), where v r represents a measure of blood velocity, r represents a distance in a direction along the mitral valve from the left atrium to the left ventricle and t represents time.
- v r represents a measure of blood velocity
- r represents a distance in a direction along the mitral valve from the left atrium to the left ventricle and t represents time.
- the blood velocity function may be a standard or conventional output of an (color) M-mode trace of the heart of the subject, as would be readily appreciated by the skilled person.
- Processing the velocity data, in step 120, to determine one or more velocity components may comprise calculating an internal component (— p — -), or temporal derivative component, and/or a d convective component (—OIL. or distance derivate component, of the velocity data p is density, and dr is assumed to be known or a predetermined value.
- r represent different positions along the direction of mitral inflow.
- Different anatomical elements may be associated with different values of r, representing the position of the anatomical element along the direction of mitral inflow.
- An r value for a particular anatomical element (such as the mitral valve) can be identified, e.g. by segmenting the ultrasound data to identify the anatomical element and identifying the corresponding r value for the position of the ultrasound data. This can be achieved because the positional relationship between positions in the ultrasound data and corresponding position along the r scale is predetermined or derivable (as the ultrasound data will contain a representation of the direction along the mitral valve).
- the method 100 then moves to a step 130 of determining a temporal variation of the left ventricular pressure (dR/dt) from the determined velocity component (in step 120) by processing the velocity component using a physical model.
- a suitable physical model for use in step 130 is provided by the Navier Stokes equation as follows: 1
- the Navier Stokes equation outputs a distance variation of the pressure - (r, t) as a function of distance and time, i.e. a pressure gradient.
- This pressure gradient may then be integrated (from the position of the atrium to the position of the ventricle on the r-scale) to obtain the left ventricular pressure as a function of time (i.e. to obtain P(t)).
- the position of the atrium and the ventricle on the r- scale may be readily identifiable through appropriate processing of the velocity data and/or ultrasound data.
- Obtaining the temporal variation of the left ventricle function can then be carried out by differentiating the left ventricular pressure as a function of time with respect to time, i.e. to calculate dV/dt.
- the method 100 also performs a step 140 of processing the ultrasound data to determine a temporal variation of the left ventricle volume (i.e. fV/ft). This can be performed, for instance, by processing image data (e.g. 2D or 3D pixel/voxel data) to determine a difference or a change in volume at a particular point in time.
- image data e.g. 2D or 3D pixel/voxel data
- the ultrasound data may comprise pixel/voxel data as a function of time.
- This pixel/voxel data could be segmented, e.g. using an automated segmentation approach, to determine a volume of the left ventricle as a function of time. This function could then be differentiated to extract a temporal variation of the left ventricle volume as a function of time.
- each frame of the pixel/voxel data may be segmented to identify a region representing the left ventricle, i.e. the predicted bounds of the left ventricle.
- the volume of the left ventricle can be determined based on the predicted bounds of the left ventricle for each frame.
- These discrete volume values can be converted into a function, e.g. a best fit function, to thereby produce the predicted volume of the left ventricle as a function of time. This produced function may then be differentiated to generate the temporal variation of the left ventricle volume as a function of time.
- v r (t) a function of blood velocity as a function of time
- v T (r, t) the value for r is selected to identify a predicted/relative position of the mitral annulus in the mitral valve, e.g. is a predetermined value.
- the value for r may be selected based on a known position of the mitral annulus with respect to r, e.g. based on an assessment of the ultrasound data and/or velocity data.
- an r value may be selected that provides the maximum value v r (r, t) (for any t). This example reflects an expectation that the maximum blood velocity will occur at the mitral annulus (e.g. where the cross-sectional area is smallest).
- determining the temporal variation of the left ventricle volume may comprise multiplying the velocity component by a surface area value S, the surface area value representing a surface area size of the mitral annulus of the subject’s heart.
- the value for S could be assumed to be constant, e.g. a predetermined or pre-estimated value.
- the value for S could be an average value (e.g. a global average or an average from amongst similar subjects).
- it may be possible to calibrate for S e.g. calculate a subject-specific value for S. This can be performed by looking at the end-systolic and end-diastolic volume, e.g. as measured from voxel data of the heart obtained at the end of the systolic phase and at the end of the diastolic phase.
- the method may comprise a step of determining the surface area value S.
- This step may comprise: determining, by processing the ultrasound data, a first volume of the left ventricle, of the subject’s heart wherein the first volume is a volume of the left ventricle of the subject’s heart at a first point during the portion of the cardiac cycle of the subject’s heart (e.g. at the end of the systolic phase); determining, by processing the ultrasound data, a second volume of the left ventricle, of the subject’s heart wherein the second volume is a volume of the left ventricle of the subject’s heart at a second, later point during the portion of the cardiac cycle of the subject’s heart (e.g.
- One advantage of processing velocity data using such approaches is the synchronous measurements of volume and pressure from data obtained during the same cardiac cycle, and where the same data is used (i.e. data sampled at the same time points).
- step 150 of estimating the left ventricular compliance by processing the determined temporal variation of the left ventricle volume and the determined temporal variation of the left ventricular pressure.
- step 150 may comprise determining one or more values of left ventricular compliance by processing the determined temporal variation of the left ventricular volume and the temporal variation of the left ventricular pressure.
- step 150 comprises determining the ratio between the temporal variation of the left ventricle volume and the temporal variation of the left ventricular pressure.
- a specific value or sample for left ventricular compliance could be calculated for each of a plurality of different points in time (e.g. by selecting an appropriate time value for t).
- an average left ventricular compliance could be calculated by averaging multiple values for left ventricular compliance that have been calculated using equation (2) for different values of t.
- the values of t may be selected so that they all fall within a predetermined phase of the motion of the subject’s heart, e.g. all fall within the diastolic phase of the subject’s heart movement and/or all within a single cardiac cycle of the subject’s heart.
- the values of t may be evenly distributed.
- the step of estimating the left ventricular compliance may comprise using a first set of measurements/samples (of the temporal variation of the left ventricle volume) and a second set of measurements/samples (of the temporal variation of the left ventricle pressure).
- Each measurement in the first/second set of measurements represents a measurement at a particular point in time (e.g. for a different value of t).
- Each measurement may represent a measurement at a particular point in time t, e.g. within a predetermined phase of the motion of the subject’s heart.
- Each measurement in the first set of measurement corresponds to a measurement in the second set of measurements, e.g. measurement representing a same point in time.
- the total number of measurements in the first set is identical to the total number of measurements in the second set.
- the first and second sets of measurements may be associated with a particular period of time, e.g. part of a single cardiac cycle of the subject, such as part of the end-diastolic phase of the subject.
- all measurements in the first/second set are measurements obtained for values of t that fall within this particular period of time.
- the values of t (for the measurements of the first/second set of measurements) may be selected so that they all fall within a predetermined phase of the motion of the subject’s heart, e.g. all fall within the diastolic phase of the subject’s heart movement and/or all within a single cardiac cycle of the subject’s heart.
- the values of t may be evenly distributed.
- a value for the left ventricular compliance k may then be calculated using the following equation: where E(f(t)) is an averaging function, which is calculated by the following equation where f(t) represents some arbitrary function of t (e.g. (t) ⁇ (t) or
- Equation (3) the value for E (— (t) —(£)) in equation (3) is calculated by determining, for each measurement in the first set of measurements, the product of said measurement and the corresponding measurement of the second set of measurements. The products produced in this step are then summed, and the summed products are divided by the total number of products (or the total number of measurements in the first/second set of measurements).
- a value for the left ventricular compliance k could be calculated using the following equation: It has previously been explained how the first/second set of measurements may be measurements representing a particular period of time (e.g. are all within a predetermined phase of the subject’s heart, such as the diastolic phase or the end of the diastolic phase).
- the time period is selected so that the denominator (of equation (3) or (5), or any other equation for calculating left ventricular compliance k that makes use of a ratio) is large.
- 3P , 2 measurements should contain a time ti for which the value of (t is greater than the value of
- the time period represented by the first/second set of measurements may be centered around this point in time ti, e.g. to represent a predetermined fraction of the cardiac cycle.
- a plurality of values for the left ventricular compliance k could be generated for different versions of the first and second sets of measurement.
- Each version of the first and second sets of measurements may correspond to sets of measurements associated with different time period (e.g. different parts of the motion of the subject’s heart).
- time periods associated with different versions of the first/second sets of measurements may represent a time period of a different cardiac cycle of a subject. In other examples, time periods associated with different versions of first/second sets of measurements may be a time period of a sliding window. Other approaches will be apparent to the skilled person.
- the method 100 may further comprise a step 160 of controlling a user-perceptible output responsive to the left ventricular compliance estimated by the processing arrangement.
- step 160 may comprise controlling the user-perceptible output to provide a visual representation of the (value(s) of the) left ventricular compliance calculated in step 150. This may be in the form of a numerical indicator (e.g. a visual representation of a number), a graphical indicator, a color indicator (e.g. different colors representing different magnitudes of the left ventricular compliance) and so on.
- step 160 may comprise controlling the user-perceptible output to provide a visual, audible and/or haptic output that is responsive to the left ventricular compliance. For instance, if a value of the left ventricular compliance breaches some predetermined threshold, a visual, audible and/or haptic alert may be generated and output.
- step 160 may comprise controlling the user-perceptible output to provide a user-perceptible output of the left ventricular compliance (e.g. a measure).
- Step 160 is not essential.
- the left ventricular compliance could act as a valuable input or feature for other (automated) method of analyzing characteristics of the subject’s heart.
- the velocity data comprises color M-mode ultrasound data at the mitral inflow during a diastolic phase of the cardiac cycle of the heart.
- This data can be obtained using 4-chamber view or 3-chamber view.
- the Doppler beam is preferably aligned with the main direction of the mitral flow.
- Color M-mode ultrasound data ensures adequate temporal sampling to resolve the temporal derivative component in the Navier-Stokes equation.
- the calculated temporal variation in left ventricular pressure (dP/dt) based on a temporal variation in left atrium pressure (dP a /dt).
- the temporal variation in left atrium pressure can be calculated using the ultrasound data. In some examples, this can be obtained from a left atrium pressure monitor, e.g. positioned in the subject.
- a pulmonary capillary wedge pressure can be used to represent the left atrium pressure, as the two are correlated with or proportional to one another. Thus, the left atrium pressure may be obtained from a pulmonary capillary wedge pressure monitor/sensor.
- Use of pulmonary capillary wedge pressure is less invasive, as it can be acquired with venous access, reducing a risk to the subject.
- Fig. 2 provides various representations of partly processed velocity M-mode ultrasound data.
- the E diastolic filling waves 201 and the A diastolic filling waves 202 can be identified.
- Each visual representation represents information captured over a period of time that encompasses three and a half heart cycles.
- a first visual representation 210 represents velocity M-mode at the mitral valve after performing some post-processing including smoothing and dealiasing. This effectively represents an example of velocity data (of ultrasound data obtained in step 110 of Fig. 1).
- the first visual representation 210 effectively provides a visual representation of the blood velocity function v r (r, t), where the x-axis represents time t, the y-axis represents distance r (in a direction along the mitral valve), and the saturation/grey scale value represents the (blood) velocity v r .
- the second visual representation 220 represents a convective component of the velocity M-mode ultrasound data.
- the saturation/greyscale value instead represents the
- convective component —pv r —y ⁇ (rather than the blood velocity).
- the convective component is calculated through appropriate processing of the blood velocity function.
- the third visual representation 230 represents an inertial component of the velocity M- mode ultrasound data.
- the saturation/greyscale value instead represents the inertial component (— p ⁇ ), which is calculated through appropriate processing of the blood velocity function.
- the fourth visual representation 240 represent a pressure gradient, i.e. the output of the Navier-Stokes equation, which is determined using the convective component and the inertial component. As previously explained, the pressure gradient can be integrated to get an estimate of left ventricular pressure.
- Fig. 3 provides a representation of the integrated pressure gradient, i.e. a function of left ventricular pressure over time.
- the y-axis represents a relative pressure and the x-axis represent time. Approaches for generating the function of left ventricular pressure have been previously described.
- Fig. 4 provides a representation of blood velocity at the mitral valve. This may effectively be a slice of the first visual representation 210, i.e. for a value of r that identifies a location of the mitral valve.
- the blood velocity at the mitral valve can be used to calculate or represent the temporal variation of the left ventricular volume ⁇ W fdt).
- each step of any herein illustrated flow chart may represent a different action performed by a processing system, and may be performed by a respective module of the processing system.
- Embodiments may therefore make use of a processing system.
- the processing arrangement can be implemented in numerous ways, with software and/or hardware, to perform the various functions required.
- a processor is one example of a processing arrangement which employs one or more microprocessors that may be programmed using software (e.g., microcode) to perform the required functions.
- a processing arrangement may however be implemented with or without employing a processor, and also may be implemented as a combination of dedicated hardware to perform some functions and a processor (e.g., one or more programmed microprocessors and associated circuitry) to perform other functions.
- Examples of processing arrangement components that may be employed in various embodiments of the present disclosure include, but are not limited to, conventional microprocessors, application specific integrated circuits (ASICs), and field-programmable gate arrays (FPGAs).
- a processor or processing arrangement may be associated with one or more storage media such as volatile and non-volatile computer memory such as RAM, PROM, EPROM, and EEPROM.
- the storage media may be encoded with one or more programs that, when executed on one or more processors and/or processing systems, perform the required functions.
- Various storage media may be fixed within a processor or processing arrangement or may be transportable, such that the one or more programs stored thereon can be loaded into a processor or processing system.
- Fig. 5 illustrates a processing system 500 comprising an imaging system 510 and an output interface 520 according to an embodiment.
- the imaging system is itself an embodiment of the invention.
- the imaging system 510 comprises a processing arrangement 511 and an ultrasound system 512.
- the processing arrangement 511 itself an embodiment of the invention, is configured to carry out any previously described method. Suitable examples of a processing arrangement have been previously described.
- the ultrasound system 512 configured to capture ultrasound data of the subject’s heart and provide the ultrasound data to the processing arrangement.
- the output interface 520 is configured to provide a user-perceptible output responsive to the left ventricular compliance estimated by the processing arrangement.
- the operation of the output interface may be controlled by the processing arrangement 511, e.g. by carrying out step 160 described with reference to Fig. 1.
- a computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.
- a suitable medium such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Hematology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP22733861.3A EP4346617B1 (en) | 2021-05-28 | 2022-05-20 | Noninvasive measurement of left ventricular compliance |
| JP2023572803A JP2024519394A (ja) | 2021-05-28 | 2022-05-20 | 左室コンプライアンスの非侵襲的測定 |
| CN202280038484.1A CN117412712A (zh) | 2021-05-28 | 2022-05-20 | 对左心室顺应性的无创测量 |
| US18/564,652 US20240260927A1 (en) | 2021-05-28 | 2022-05-20 | Noninvasive measurement of left ventricular compliance |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP21290035.1 | 2021-05-28 | ||
| EP21290035.1A EP4094696A1 (en) | 2021-05-28 | 2021-05-28 | Noninvasive measurement of left ventricular compliance |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2022248355A1 true WO2022248355A1 (en) | 2022-12-01 |
Family
ID=76695667
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/EP2022/063723 Ceased WO2022248355A1 (en) | 2021-05-28 | 2022-05-20 | Noninvasive measurement of left ventricular compliance |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20240260927A1 (https=) |
| EP (2) | EP4094696A1 (https=) |
| JP (1) | JP2024519394A (https=) |
| CN (1) | CN117412712A (https=) |
| WO (1) | WO2022248355A1 (https=) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN120636716B (zh) * | 2025-08-15 | 2025-11-25 | 国家心血管病中心 | 一种基于医学影像的左室顺应性评估方法及系统 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6780155B2 (en) * | 2001-12-18 | 2004-08-24 | Koninklijke Philips Electronics | Method and system for ultrasound blood flow imaging and volume flow calculations |
| US7204798B2 (en) * | 2003-01-24 | 2007-04-17 | Proteus Biomedical, Inc. | Methods and systems for measuring cardiac parameters |
| US8211024B2 (en) * | 2005-06-06 | 2012-07-03 | Siemens Medical Solutions Usa, Inc. | Medical ultrasound pressure gradient measurement |
| US20120172727A1 (en) * | 2010-12-30 | 2012-07-05 | Boston Scientific Scimed, Inc. | Imaging system |
| WO2013051943A1 (en) * | 2011-10-03 | 2013-04-11 | Surf Technology As | Nonlinear imaging with dual band pulse complexes |
| US20130245441A1 (en) * | 2012-03-13 | 2013-09-19 | Siemens Medical Solutions Usa, Inc. | Pressure-Volume with Medical Diagnostic Ultrasound Imaging |
| WO2013147262A1 (ja) * | 2012-03-30 | 2013-10-03 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、画像処理装置及び画像処理方法 |
| US12102415B2 (en) * | 2019-11-07 | 2024-10-01 | The Curators Of The University Of Missouri | Systems and methods to improve management and monitoring of cardiovascular disease |
-
2021
- 2021-05-28 EP EP21290035.1A patent/EP4094696A1/en not_active Withdrawn
-
2022
- 2022-05-20 JP JP2023572803A patent/JP2024519394A/ja active Pending
- 2022-05-20 EP EP22733861.3A patent/EP4346617B1/en active Active
- 2022-05-20 US US18/564,652 patent/US20240260927A1/en active Pending
- 2022-05-20 WO PCT/EP2022/063723 patent/WO2022248355A1/en not_active Ceased
- 2022-05-20 CN CN202280038484.1A patent/CN117412712A/zh active Pending
Non-Patent Citations (2)
| Title |
|---|
| DEMARIA A N ET AL: "Mitral valve early diastolic closing velocity in the echocardiogram: Relation to sequential diastolic flow and ventricular compliance", AMERICAN JOURNAL OF CARDIOLOGY, CAHNERS PUBLISHING CO., NEWTON, MA, US, vol. 37, no. 5, 1 April 1976 (1976-04-01), pages 693 - 700, XP026333936, ISSN: 0002-9149, [retrieved on 19760401], DOI: 10.1016/0002-9149(76)90362-3 * |
| HANOWELL ET AL: "Transesophageal echocardiography in the perioperative assessment of intravascular volume", SEMINARS IN ANESTHESIA, SAUNDERS, CO, NEW YORK, NY, US, vol. 17, no. 3, 1 September 1998 (1998-09-01), pages 252 - 266, XP005452778, ISSN: 0277-0326, DOI: 10.1016/S0277-0326(98)80023-8 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117412712A (zh) | 2024-01-16 |
| US20240260927A1 (en) | 2024-08-08 |
| EP4346617B1 (en) | 2024-12-04 |
| EP4094696A1 (en) | 2022-11-30 |
| EP4346617A1 (en) | 2024-04-10 |
| JP2024519394A (ja) | 2024-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11576639B2 (en) | Method, system, device and medium for determining a blood flow velocity in a vessel | |
| CN112074866B (zh) | 4d mr图像数据中的流分析 | |
| KR101939778B1 (ko) | 필요 혈류량 결정 방법 및 장치, 혈류 영상 생성 방법 및 장치, 심근 관류 영상 처리 방법 및 장치 | |
| US10299862B2 (en) | Three-dimensional quantitative heart hemodynamics in medical imaging | |
| US9999361B2 (en) | Method and system for non-invasive assessment of coronary artery disease | |
| JP4489770B2 (ja) | 動く臓器を描出した画像シリーズ内の最適画像を識別する方法及びシステム | |
| US12115014B2 (en) | Most relevant x-ray image selection for hemodynamic simulation | |
| EP3606437B1 (en) | Method and apparatus for physiological functional parameter determination | |
| KR101334064B1 (ko) | 혈관영상에서의 움직임 추정 방법 및 그 장치 | |
| CN110637322B (zh) | 用于数字化图像中的边缘检测的系统、方法以及计算机可读存储介质 | |
| US10438364B2 (en) | System and method for registration of perfusion images | |
| CN111275755A (zh) | 基于人工智能的二尖瓣瓣口面积检测方法、系统和设备 | |
| Medvedofsky et al. | Feasibility of left ventricular global longitudinal strain measurements from contrast-enhanced echocardiographic images | |
| US20240260927A1 (en) | Noninvasive measurement of left ventricular compliance | |
| Wick et al. | Detection of cardiac quiescence from B-mode echocardiography using a correlation-based frame-to-frame deviation measure | |
| EP4280965B1 (en) | Blood velocity determination | |
| KR20180008134A (ko) | 분획 혈류 예비력 예측 방법 | |
| JP2023552331A5 (https=) | ||
| JPWO2018159708A1 (ja) | 血流解析装置および方法並びにプログラム | |
| EP4329590B1 (en) | System and method for estimating value for target cardiac parameter | |
| Turco et al. | Fully automated quantification of left and right ventricular volumes throughout the cardiac cycle from magnetic resonance imaging | |
| Bakker et al. | Preclinical validation of the advection diffusion flow estimation method using computational patient specific coronary tree phantoms | |
| Moladoust et al. | Radial strain assessment of the interventricular septum wall by a new technique in healthy subjects | |
| Herway et al. | Left ventricular systolic function | |
| Vardoulis et al. | Articles in PresS. Am J Physiol Heart Circ Physiol (April 19, 2013). doi: 10.1152/ajpheart. 00963.2012 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22733861 Country of ref document: EP Kind code of ref document: A1 |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 2023572803 Country of ref document: JP |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 202280038484.1 Country of ref document: CN |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 2022733861 Country of ref document: EP |
|
| ENP | Entry into the national phase |
Ref document number: 2022733861 Country of ref document: EP Effective date: 20240102 |