Method and System for Multi-Scale Anatomical and Functional Modeling of Coronary Circulation

• US 20130132054A1
• Filed: 11/09/2012
• Published: 05/23/2013
• Est. Priority Date: 11/10/2011
• Status: Active Grant

• Current Assignee: Siemens Healthineers AG (Siemens AG)
• Original Assignee: Siemens Healthcare GMBH (Siemens AG)
• Inventors: Sharma, Puneet, Georgescu, Bogdan, Mihalef, Viorel, Kamen, Ali, Comaniciu, Dorin, Itu, Lucian Mihai
• Granted Patent: US 10,162,932 B2
• Time in Patent Office: Days
• US Class Current: 703/9
• CPC Class Codes: G16B 5/00 ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
• 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/30 for calculating health indices; for individual health risk assessment
• CPC Class Codes: G16H 50/50 for simulation or modelling of medical disorders
• CPC Class Codes: G16Z 99/00 Subject matter not provided for in other main groups of this subclass

Abstract

A method and system for multi-scale anatomical and functional modeling of coronary circulation is disclosed. A patient-specific anatomical model of coronary arteries and the heart is generated from medical image data of a patient. A multi-scale functional model of coronary circulation is generated based on the patient-specific anatomical model. Blood flow is simulated in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation. Hemodynamic quantities, such as fractional flow reserve (FFR), are computed to determine a functional assessment of the stenosis, and virtual intervention simulations are performed using the multi-scale function model of coronary circulation for decision support and intervention planning.

Claims

1. A method, comprising:

generating a patient-specific anatomical model of coronary arteries and the heart from medical image data of a patient;

generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model; and

simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation.

2. The method of claim 1, wherein generating a patient-specific anatomical model of coronary arteries and the heart from medical image data of a patient comprises:

generating a 4D geometric model of the coronary arteries from 4D medical image data; and

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

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

segmenting the coronary arteries are segmented in each of a plurality of frames of the 4D medical image data; and

generating a geometric surface model for the segmented coronary arteries in each of the plurality of frames of the 4D medical image data.

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

extracting individual models of each of a plurality of heart components in each of a plurality of frames of the 4D medical image data; and

integrating the individual models for the plurality of heart components in each of the plurality of frames of the 4D medical image data by establishing mesh point correspondence between the individual models.

5. The method of claim 1, wherein generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model comprises:

generating a 3D computation model for each of one or more stenosis regions in the coronary arteries;

generating 1D computation models for non-stenosis regions of the coronary arteries and the aorta; and

representing microvasculature vessels using 0D lumped models.

6. The method of claim 5, wherein the 3D computation model for each stenosis region is a rigid wall 3D model and 0D interface models between the 3D computation model for each stenosis region and the 1D computation models for non-stenosis regions of the coronary arteries adjacent to each stenosis region concentrate compliance of the stenosis region.

7. The method of claim 5, wherein generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model further comprises:

generating a structured tree model for the vascular tree of the patient.

8. The method of claim 5, wherein generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model further comprises:

generating a reduced order model of the heart from full-order anatomical and hemodynamic models of the heart.

9. The method of claim 8, wherein generating a reduced order model of the heart from full-order anatomical and hemodynamic models of the heart comprises:

estimating motion and mechanical parameters of one or more heart components based on the anatomical and hemodynamic models of the heart; and

determining boundary conditions for computational fluid dynamic simulations based on the motion and mechanical parameters of the one or more heart components.

10. The method of claim 1, wherein simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation comprises:

simulating blood flow in the at least one stenosis region using the multi-scale function model of coronary circulation based on boundary conditions determined from the anatomical model of the coronary arteries and the heart.

11. The method of claim 5, wherein simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation comprises:

performing computational fluid dynamics (CFD) simulations in the 3D computation model for each stenosis region and the 1D computation models; and

coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models.

12. The method of claim 11, wherein coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models comprises:

deriving an inflow boundary condition of a system tree model by coupling a 1D computation model representing the aorta to the left ventricle of a heart model.

13. The method of claim 11, wherein coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models comprises:

imposing boundary conditions representing an influence of heart contractions on 1D computation models of epicardial coronary vessels using a 3D strain map extracted from the medical image data.

14. The method of claim 11, wherein coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models comprises:

determining extracellular pressure applied to the 1D computation models of coronary vessels based on locations of the coronary vessels using the 0D lumped models.

15. The method of claim 11, wherein coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models comprises:

coupling the 1D computation models to the 0D lumped models via wall shear stress terms.

16. The method of claim 11, wherein coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models comprises:

coupling the 3D computation model to adjacent 1D computation models using 0D interface models.

17. The method of claim 1, further comprising:

calculating a hemodynamic quantity to determine a functional significance of the at least one stenosis region based on the simulated blood flow through the at least one stenosis region.

18. The method of claim 17, wherein calculating a hemodynamic quantity to determine a functional significance of the at least one stenosis region based on the simulated blood flow through the at least one stenosis region comprises:

calculating a fractional flow reserve (FFR) of the at least one stenosis region based on the computation blood flow through the at least one stenosis region.

19. The method of claim 1, further comprising:

simulating a virtual intervention in at least one stenosis region using the multi-scale function model of coronary circulation.

20. The method of claim 19, wherein simulating a virtual intervention in at least one stenosis region using the multi-scale function model of coronary circulation comprises:

simulating a balloon inflation by virtually reducing an obstruction from the at least one stenosis region in the multi-scale function model of coronary circulation and re-simulating the blood flow through the at least one stenosis region.

21. The method of claim 19, wherein simulating a virtual intervention in at least one stenosis region using the multi-scale function model of coronary circulation comprises:

simulating a stent implantation by introducing a virtual stent model to at least one stenosis region in the multi-scale function model of coronary circulation and re-simulating the blood flow through the at least one stenosis region.

22. The method of claim 19, wherein simulating a virtual intervention in at least one stenosis region using the multi-scale function model of coronary circulation comprises:

simulating a coronary artery bypass graft (CABG) by adding a bypass vessel adjacent to the at least one stenosis region in the multi-scale function model of coronary circulation and re-simulating the blood flow through the at least one stenosis region.

23. An apparatus comprising:

means for generating a patient-specific anatomical model of coronary arteries and the heart from medical image data of a patient;

means for generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model; and

means for simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation.

24. The apparatus of claim 23, wherein the means for generating a patient-specific anatomical model of coronary arteries and the heart from medical image data of a patient comprises:

means for generating a 4D geometric model of the coronary arteries from 4D medical image data; and

means for generating a 4D anatomical model of the heart from the 4D medical image data.

25. The apparatus of claim 23, wherein the means for generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model comprises:

means for generating a 3D computation model for each of one or more stenosis regions in the coronary arteries;

means for generating 1D computation models for non-stenosis regions of the coronary arteries and the aorta; and

means for representing microvasculature vessels using 0D lumped models.

26. The apparatus of claim 25, wherein the 3D computation model for each stenosis region is a rigid wall 3D model and 0D interface models between the 3D computation model for each stenosis region and the 1D computation models for non-stenosis regions of the coronary arteries adjacent to each stenosis region concentrate compliance of the stenosis region.

27. The apparatus of claim 25, wherein the means for generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model further comprises:

means for generating a reduced order model of the heart from full-order anatomical and hemodynamic models of the heart.

28. The apparatus of claim 23, wherein the means for simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation comprises:

means for simulating blood flow in the at least one stenosis region using the multi-scale function model of coronary circulation based on boundary conditions determined from the anatomical model of the coronary arteries and the heart.

29. The apparatus of claim 25, wherein the means for simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation comprises:

means for performing computational fluid dynamics (CFD) simulations in the 3D computation model for each stenosis region and the 1D computation models; and

means for coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models.

30. The apparatus of claim 23, further comprising:

means for calculating a hemodynamic quantity to determine a functional significance of the at least one stenosis region based on the simulated blood flow through the at least one stenosis region.

31. The apparatus of claim 23, further comprising:

simulating a virtual intervention in at least one stenosis region using the multi-scale function model of coronary circulation.

32. A non-transitory computer readable medium storing computer program instructions, which when operated on a processor cause the processor to perform operations comprising:

generating a patient-specific anatomical model of coronary arteries and the heart from medical image data of a patient;

generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model; and

simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation.

33. The non-transitory computer readable medium of claim 32, wherein generating a patient-specific anatomical model of coronary arteries and the heart from medical image data of a patient comprises:

generating a 4D geometric model of the coronary arteries from 4D medical image data; and

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

34. The non-transitory computer readable medium of claim 32, wherein generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model comprises:

generating a 3D computation model for each of one or more stenosis regions in the coronary arteries;

generating 1D computation models for non-stenosis regions of the coronary arteries and the aorta; and

representing microvasculature vessels using 0D lumped models.

35. The non-transitory computer readable medium of claim 34, wherein the 3D computation model for each stenosis region is a rigid wall 3D model and 0D interface models between the 3D computation model for each stenosis region and the 1D computation models for non-stenosis regions of the coronary arteries adjacent to each stenosis region concentrate compliance of the stenosis region.

36. The non-transitory computer readable medium of claim 34, wherein generating a multi-scale functional model of coronary circulation based on the patient-specific anatomical model further comprises:

generating a reduced order model of the heart from full-order anatomical and hemodynamic models of the heart.

37. The non-transitory computer readable medium of claim 32, wherein simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation comprises:

simulating blood flow in the at least one stenosis region using the multi-scale function model of coronary circulation based on boundary conditions determined from the anatomical model of the coronary arteries and the heart.

38. The non-transitory computer readable medium of claim 34, wherein simulating blood flow in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation comprises:

performing computational fluid dynamics (CFD) simulations in the 3D computation model for each stenosis region and the 1D computation models; and

coupling the 3D computation model for each stenosis region, the 1D computation models, and the 0D lumped models.

39. The non-transitory computer readable medium of claim 32, the operations further comprising:

calculating a hemodynamic quantity to determine a functional significance of the at least one stenosis region based on the simulated blood flow through the at least one stenosis region.

40. The non-transitory computer readable medium of claim 32, the operations further comprising:

simulating a virtual intervention in at least one stenosis region using the multi-scale function model of coronary circulation.