Methods for the non-invasive determination of heart and pulmonary pressures

• US 10,390,713 B2
• Filed: 03/16/2014
• Issued: 08/27/2019
• Est. Priority Date: 03/15/2013
• Status: Active Grant

• Current Assignee: Ghassan S. Kassab, Yunlong Huo
• Original Assignee: Ghassan S. Kassab, Yunlong Huo
• Inventors: Kassab, Ghassan S., Huo, Yunlong
• Primary Examiner(s): Moher, Amanda Lauritzen
• Application Number:
  US14/215,027
• Publication Number:
  US 20140275961A1
• Time in Patent Office: 1,990 Days
• CPC Class Codes: A61B 5/02133 by using induced vibration of the blood vessel
• CPC Class Codes: A61B 5/022 by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
• CPC Class Codes: A61B 5/0263 using NMR
• CPC Class Codes: A61B 5/0265 using electromagnetic means, e.g. electromagnetic flowmeter
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: A61B 6/503 for diagnosis of the heart
• CPC Class Codes: A61B 6/504 for diagnosis of blood vessels, e.g. by angiography
• CPC Class Codes: A61B 8/04 Measuring blood pressure
• CPC Class Codes: A61B 8/06 Measuring blood flow measuring volume flow in general G01F, e.g. G01F1/66, G01F1/72; measuring speed of fluids in general G01P5/00
• CPC Class Codes: A61B 8/0883 for diagnosis of the heart
• CPC Class Codes: A61B 8/485 involving measuring strain or elastic properties
• CPC Class Codes: A61B 8/5223 for extracting a diagnostic or physiological parameter from medical diagnostic data for algorithms to analyse biomedical images G06T7/0012
• CPC Class Codes: G01R 33/20 involving magnetic resonance medical aspects A61B5/055; magnetic resonance gyrometers G01C19/60
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G16H 50/30 for calculating health indices; for individual health risk assessment

Abstract

The disclosure of the present application provides methods for the noninvasive determination of cardiac and pulmonary pressures based off of patient-specific medical images and one or more computations. In one exemplary embodiment, the method comprises noninvasively generating a patient-specific medical image of at least a portion of a heart, determining a free vibration measurement from the left ventricular free wall based on the patient-specific image, and performing at least one computation using the free vibration measurement to noninvasively determine a cardiac or pulmonary pressure.

Claims

1. A method for noninvasively determining one or more pressure measurements comprising the steps of:

generating a patient-specific medical image of at least one portion of a heart of a patient noninvasively;

determining, from the medical image, a measured position of a left ventricular free wall of the heart;

determining from the patient-specific medical image at least two measurements from a portion of the heart, one of the at least two measurements comprising a free vibration measurement determined from the left ventricular free wall of the heart shown in the medical image, and one of the at least two measurements relating to the measured position of the left ventricular free wall;

calculating at least one targeted pressure within one or more portions of at least one of the heart and an adjacent vessel using one or more computations and at least one of the at least two measurements, the at least one targeted pressure comprising a first targeted pressure comprising a left ventricular end-diastolic pressure and at least one of the one or more computations for determining the left ventricular end-diastolic pressure by iterative solution that numerically optimizes a calculated position of the left ventricular free wall to match the measured position of the left ventricular free wall;

calculating a blood flow velocity and determining a pressure gradient from the blood flow velocity, the pressure gradient determined from magnetic resonance imaging blood flow velocity measurements through a right ventricle of the heart;

noninvasively determining venous pressure; and

using the venous pressure as reference for at least one of determining and validating the pressure gradient;

wherein at least one of the one or more computations determines an end-diastolic Young'"'"'s modulus value of the left ventricular free wall using the free vibration measurement;

wherein at least one of the one or more computations determines patient-specific material parameters in diastolic strain energy function of the heart, where through numerical optimization a Young'"'"'s modulus value from strain energy function equals the value of the end-diastolic Young'"'"'s modulus value of the left ventricular free wall; and

wherein the preceeding steps are performed using a computer, the computer comprising a processor capable of processing obtained data.

2. The method of claim 1, wherein the targeted pressure is selected from the group consisting of a left ventricle end-diastolic pressure, a right ventricle pressure, and a pulmonary pressure.

3. The method of claim 1, wherein the patient-specific medical image of at least a portion of the heart is generated through magnetic resonance imaging or echocardiography.

4. The method of claim 1, wherein:

at least one of the at least one measurement determined from the patient-specific medical image relates to the position of the interventricular septum; and

at least one of the computations used in the calculating at least one targeted pressure determines right ventricular end-diastolic pressure by iterative solution such that a predicted end-diastolic interventricular septum position matches the at least one measurement related to the position of the interventricular septum.

5. The method of claim 1, wherein determining a pressure gradient further comprises calculating the pressure gradient using a Navier-Stokes equation.

6. The method of claim 1, further comprising the step of determining left ventricular systolic pressure by noninvasively measuring aortic pressure;

left ventricular systolic pressure being at aortic pressure.

7. The method of claim 6, further comprising:

determining, from the medical image, a position of a left ventricle systolic free wall of the heart, wherein;

at least one of the at least two measurements determined from the patient-specific medical image relates to the position of the left ventricle systolic free wall;

at least one of the at least one targeted pressure comprises a third targeted pressure comprising right ventricular systolic pressure; and

at least one of the one or more computations used in the calculating the third targeted pressure determines patient-specific material parameters of the heart using a systolic second Piola-Kirchhoff stress value and numerical optimization such that a predicted left ventricle systolic free wall position equates with the at least one measurement related to the position of the left ventricle systolic free wall.

8. The method of claim7, further comprising determining, from the medical image, a position of a systolic interventricular septum of the heart, wherein:

at least one of the at least two measurements determined from the patient-specific medical image relates to the position of a systolic interventricular septum; and

at least one of the one or more computations used in the calculating the at least one targeted pressure determines right ventricular systolic pressure by iterative solution such that a calculated systolic interventricular septum position matches the at least one measurement related to the position of the systolic interventricular septum.

9. The method of claim 8, further comprising the step of calculating a fourth targeted pressure within one or more portions of at least one of the heart and an adjacent vessel, wherein the third targeted pressure comprises pulmonary pressure derived from the right ventricular systolic pressure and the pressure gradient.

10. The method of claim 6, further comprising the steps of:

assessing a pressure difference across an interventricular septum of the heart in view of the left ventricular systolic pressure to determine right ventricular systolic pressure.

11. The method of claim 10, wherein the aortic pressure is obtained through blood pressure cuff measurements.

12. The method of claim 1, further comprising the step of:

calculating a second targeted pressure of the at least one targeted pressure within one or more portions of at least one of the heart and an adjacent vessel using at least one of the one or more computations and at least one of the at least two measurements.

13. The method of claim 12, wherein the second targeted pressure comprises left ventricular systolic pressure and the step of calculating the second targeted pressure comprises calculating left ventricular systolic pressure by noninvasively measuring aortic pressure of the patient, wherein left ventricular systolic pressure is at aortic pressure.

14. A method for noninvasively determining one or more pressure measurements comprising the steps of:

generating a patient-specific medical image of at least one portion of a heart of a patient noninvasively;

determining, from the medical image, a first free vibration measurement from a left ventricle systolic free wall of the heart;

determining, from the medical image, a position of the left ventricle systolic free wall of the heart and a position of a systolic interventricular septum of the heart;

determining a second free vibration measurement related to the position of the systolic interventricular septum from the medical image;

calculating left ventricular systolic pressure by noninvasively measuring aortic pressure of the patient, wherein left ventricular systolic pressure is at aortic pressure;

calculating a blood flow velocity and determining a pressure gradient from the blood flow velocity, the pressure gradient determined from magnetic resonance imaging blood flow velocity measurements through a right ventricle of the heart;

calculating an end-diastolic Young'"'"'s modulus value of the left ventricular free wall using the second free vibration measurement;

determining a first patient-specific material parameter of the heart in diastolic strain energy function using a first numerical optimization such that a Young'"'"'s modulus value from strain energy function matches the end-diastolic Young'"'"'s modulus value of the left ventricular free wall;

determining a second patient-specific material parameter of the heart using a systolic second Piola-Kirchhoff stress value and a second numerical optimization such that a left ventricle systolic free wall position matches the position of the left ventricle systolic free wall of the heart determined from the medical image;

determining a right ventricular systolic pressure by iterative solution such that a predicted systolic interventricular septum position matches the position of the systolic interventricular septum determined from the medical image;

noninvasively determining venous pressure;

using the venous pressure as reference for validating the pressure gradient; and

wherein the preceeding steps are performed using a computer, the computer comprising a processor capable of processing obtained data.