Method and system for non-invasive assessment of coronary artery disease

• US 9,119,540 B2
• Filed: 09/07/2011
• Issued: 09/01/2015
• Est. Priority Date: 09/16/2010
• Status: Active Grant

• Current Assignee: Siemens Healthineers AG (Siemens AG)
• Original Assignee: Siemens AG
• Inventors: Sharma, Puneet, Georgescu, Bogdan, Mihalef, Viorel, Chen, Terrence, Comaniciu, Dorin
• Primary Examiner(s): Shah, Kamini S
• Assistant Examiner(s): Hann, Jay B
• Application Number:
  US13/226,779
• Publication Number:
  US 20120072190A1
• Time in Patent Office: 1,455 Days
• Field of Search: None
• US Class Current: 1/1
• CPC Class Codes: A61B 2576/023 for the heart
• CPC Class Codes: A61B 5/026 Measuring blood flow A61B3/1233, A61B3/1241 take precedence
• CPC Class Codes: G06T 2207/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/10096 Dynamic contrast-enhanced magnetic resonance imaging [DCE-MRI]
• CPC Class Codes: G06T 2207/10116 X-ray image
• CPC Class Codes: G06T 2207/30048 Heart; Cardiac
• CPC Class Codes: G06T 2207/30104 Vascular flow; Blood flow; Perfusion
• CPC Class Codes: G06T 7/0016 involving temporal comparison
• CPC Class Codes: G16H 20/17 delivered via infusion or injection
• CPC Class Codes: G16H 30/20 for handling medical images, e.g. DICOM, HL7 or PACS
• CPC Class Codes: G16H 30/40 for processing medical images, e.g. editing
• CPC Class Codes: G16H 50/50 for simulation or modelling of medical disorders

Abstract

In order to assess coronary artery disease from medical image data, an anatomical model of a coronary artery is generated from the medical image data. A velocity of blood in the coronary artery is estimated based on a spatio-temporal representation of contrast agent propagation in the medical image data. Blood flow is simulated in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition.

Claims

1. A method for assessment of coronary artery disease based on medical image data comprising:

generating an anatomical model of a coronary artery from medical image data;

estimating a velocity of blood in the coronary artery based on a spatio-temporal representation of contrast agent propagation over a length of the coronary artery in the medical image data, wherein estimating a velocity of blood in the coronary artery comprises;

estimating a maximum velocity of blood on a cross-section of the coronary artery based on the spatio-temporal representation of contrast agent propagation over the length of the coronary artery in the medical image data, and

mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters; and

simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition.

2. The method of claim 1, wherein the medical image data comprises 4D medical image data, and the step of generating an anatomical model of a coronary artery from medical image data comprises:

generating a 4D anatomical model of the coronary artery from 4D medical image data.

3. The method of claim 2, wherein the step of generating a 4D anatomical model of the coronary artery from 4D medical image data comprises:

segmenting a coronary artery of interest in each of a plurality of frames of the 4D medical image data; and

generating a 3D model of the segmented coronary artery of interest in each of the plurality of frames of the 4D medical image data.

4. The method of claim 1, wherein the step of estimating a maximum velocity of blood on a cross-section of the coronary artery based on a spatio-temporal representation of contrast agent propagation over a length of the coronary artery in the medical image data comprises:

generating the spatio-temporal representation of contrast agent propagation over the length of the coronary artery based on time-intensity curves extracted from the medical image data; and

estimating a slope of the spatio-temporal representation of contrast agent propagation for each point in the medical image data over the length of the coronary artery.

5. The method of claim 4, wherein the step of generating the spatio-temporal representation of contrast agent propagation over the length of the coronary artery based on time-intensity curves extracted from the medical image data comprises:

mapping an inverted intensity magnitude of points in the medical image data over the length of the coronary artery over time.

6. The method of claim 1, wherein the step of simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition comprises:

solving Navier-Stokes equations constrained by the anatomical model of the coronary artery with the estimated velocity of the blood as an inlet velocity boundary condition.

7. The method of claim 1, wherein the step of simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition comprises:

simulating pressure and velocity values over the coronary artery over time.

8. The method of claim 1, further comprising:

calculating one or more hemodynamic parameters based on the CFD simulation.

9. The method of claim 8, wherein the one or more hemodynamic parameters comprise at least one of flow rate, pressure drop, and wall-shear stress.

10. The method of claim 8, wherein the one or more hemodynamic parameters comprises at least one of coronary flow reserve (CFR) and fractional flow reserve (FFR).

11. The method of claim 1, wherein the medical image data comprises a sequence of computed tomography (CT) data.

12. The method of claim 1, wherein the medical image data comprises a sequence of X-ray angiography data.

13. The method of claim 1, wherein the medical image data comprises a sequence of rotational angiography data.

14. The method of claim 1, wherein mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters comprises:

mapping the maximum velocity on the cross-section of the coronary artery to the average velocity representing the spatial average of the velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters determined using a series of CFD simulations having various prescribed flow values and various simulated maximum and average velocity values.

15. The method of claim 1, wherein mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters comprises:

calculating the average velocity representing the spatial average of the velocity profile over the cross-section of the coronary artery as

16. The method of claim 15, wherein the patient-specific parameters p and q are determined using a series of CFD simulations having various prescribed flow values, various values for one or more of W, f, R, ρ

, or μ

, and various simulated maximum and average velocity values.

17. An apparatus for assessment of coronary artery disease based on medical image data comprising:

means for generating an anatomical model of a coronary artery from medical image data;

means for estimating a velocity of blood in the coronary artery based on a spatio-temporal representation of contrast agent propagation over a length of the coronary artery in the medical image data, wherein the means for estimating a velocity of blood in the coronary artery comprises;

means for estimating a maximum velocity of blood on a cross-section of the coronary artery based on the spatio-temporal representation of contrast agent propagation over the length of the coronary artery in the medical image data,

means for mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters; and

means for simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition.

18. The apparatus of claim 17, wherein the means for estimating a maximum velocity of blood on a cross-section of the coronary artery based on a spatio-temporal representation of contrast agent propagation over a length of the coronary artery in the medical image data comprises:

means for generating the spatio-temporal representation of contrast agent propagation over the length of the coronary artery based on time-intensity curves extracted from the medical image data; and

means for estimating a slope of the spatio-temporal representation of contrast agent propagation for each point in the medical image data over the length of the coronary artery.

19. The apparatus of claim 18, wherein the means for generating the spatio-temporal representation of contrast agent propagation over the length of the coronary artery based on time-intensity curves extracted from the medical image data comprises:

means for mapping an inverted intensity magnitude of points in the medical image data over the length of the coronary artery over time.

20. The apparatus of claim 17, wherein the means for simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition comprises:

means for solving Navier-Stokes equations constrained by the anatomical model of the coronary artery with the estimated velocity of the blood as an inlet velocity boundary condition.

21. The apparatus of claim 17, further comprising:

means for calculating one or more hemodynamic parameters based on the CFD simulation.

22. The apparatus of claim 17, wherein the means for mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters comprises:

means for mapping the maximum velocity on the cross-section of the coronary artery to the average velocity representing a spatial average of a velocity profile over the cross section of the coronary artery using patient-specific mapping parameters determined using a series of CFD simulations having various prescribed flow values and various simulated maximum and average velocity values.

23. A non-transitory computer readable medium encoded with computer executable instructions for assessment of coronary artery disease based on medical image data, the computer executable instructions defining steps comprising:

generating an anatomical model of a coronary artery from medical image data;

estimating a velocity of blood in the coronary artery based on a spatio-temporal representation of contrast agent propagation over a length of the coronary artery in the medical image data, wherein estimating a velocity of blood in the coronary artery comprises;

estimating a maximum velocity of blood on a cross-section of the coronary artery based on the spatio-temporal representation of contrast agent propagation over the length of the coronary artery in the medical image data, and

mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters; and

simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition.

24. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of estimating a maximum velocity of blood on a cross-section of the coronary artery based on a spatio-temporal representation of contrast agent propagation over a length of the coronary artery in the medical image data comprise computer executable instructions defining the steps of:

generating the spatio-temporal representation of contrast agent propagation over the length of the coronary artery based on time-intensity curves extracted from the medical image data; and

estimating a slope of the spatio-temporal representation of contrast agent propagation for each point in the medical image data over the length of the coronary artery.

25. The non-transitory computer readable medium of claim 24, wherein the computer executable instructions defining the step of generating the spatio-temporal representation of contrast agent propagation over a length of the coronary artery based on time-intensity curves extracted from the medical image data comprise computer executable instructions defining the step of:

mapping an inverted intensity magnitude of points in the medical image data over the length of the coronary artery over time.

26. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of simulating blood flow in the anatomical model of the coronary artery using a computational fluid dynamics (CFD) simulation using the estimated velocity of the blood in the coronary artery as a boundary condition comprise computer executable instructions defining the step of:

solving Navier-Stokes equations constrained by the anatomical model of the coronary artery with the estimated velocity of the blood as an inlet velocity boundary condition.

27. The non-transitory computer readable medium of claim 23, further comprising computer executable instructions defining the step of:

calculating one or more hemodynamic parameters based on the CFD simulation.

28. The non-transitory computer readable medium of claim 23, wherein mapping the maximum velocity of the blood on the cross-section of the coronary artery to an average velocity representing a spatial average of a velocity profile over the cross-section of the coronary artery using patient-specific mapping parameters comprises:

mapping the maximum velocity on the cross-section of the coronary artery to the average velocity representing a spatial average of a velocity profile over the cross section of the coronary artery using patient-specific mapping parameters determined using a series of CFD simulations having various prescribed flow values and various simulated maximum and average velocity values.