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

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

V

avg

=

V

max

⁢

(

1

+

p

⁢

⁢

W

q

)

2

,

where V avg is the average velocity, V max is the maximum velocity, p and q are the patient-specific mapping parameters, and

W

=

R

⁢

f

⁢

⁢

ρ

μ

,

where f is a frequency of the pulsatile flow, R is a cross-sectional area of the cross-section of the coronary artery, ρ is the density of blood, and μ the dynamic viscosity of blood.

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.