System for the determination of vessel geometry and flow characteristics

• US 7,738,626 B2
• Filed: 01/31/2006
• Issued: 06/15/2010
• Est. Priority Date: 02/04/2005
• Status: Active Grant

• Current Assignee: Koninklijke Philips N.V.
• Original Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Inventors: Weese, Juergen, Bredno, Joerg, Groth, Alexandra
• Primary Examiner(s): Glick; Edward J
• Assistant Examiner(s): Taningco; Alexander H
• Application Number:
  US11/815,113
• Publication Number:
  US 20080192887A1
• Time in Patent Office: 1,596 Days
• Field of Search: 378/4, 378/41, 378/98.12
• US Class Current: 378/41
• CPC Class Codes: A61B 6/481 involving the use of contrast agents
• CPC Class Codes: A61B 6/504 for diagnosis of blood vessels, e.g. by angiography
• CPC Class Codes: G06T 11/003 Reconstruction from projections, e.g. tomography
• CPC Class Codes: G06T 2207/10112 Digital tomosynthesis [DTS]
• CPC Class Codes: G06T 2207/10116 X-ray image
• CPC Class Codes: G06T 2207/30104 Vascular flow; Blood flow; Perfusion
• CPC Class Codes: G06T 2211/404 Angiography
• CPC Class Codes: G06T 7/0012 Biomedical image inspection

Abstract

The invention relates to a method and a system for the simultaneous reconstruction of the three-dimensional vessel geometry and the flow characteristics in a vessel system. According to one realization of the method, vessel segments (41) of a parametric models are fitted to differently oriented X-ray projections (P₁, P_k, P_N) of the vessel system that are generated during the passage of a bolus of contrast agent, wherein the fitting takes the imaged contrast agent dynamics and physiological a priori knowledge into account. In an alternative embodiment, the vessel geometry is reconstructed progressively along each vessel, wherein a new segment of a vessel is added based on the continuity of the vessel direction, vessel radius and reconstructed contrast agent dynamics.

Claims

1. Investigation system for the determination of vessel geometry and flow characteristics in a vessel system, comprising:

an X-ray device for generating differently oriented projections of a vessel system during a passage of a bolus of contrast agent therethrough, the projections including data representative of contrast agent dynamics within the vessel system;

a data processing unit adapted to receive data representative of at least some of the differently oriented projections from the X-ray device and adapted to reconstruct a three-dimensional model of a geometry and flow characteristics of the vessel system from received data generated in a single temporal sequence with the use of predetermined modeled contrast agent dynamics.

2. The system according to claim 1, wherein the data processing unit comprises a parametric model of both the vessel geometry and the contrast agent dynamics.

3. The system according to claim 2, wherein the data processing unit is adapted to reconstruct an instance of the model that fits both to the sequence of X-ray projections (P₁, P_k, P_N) and to plausibility rules for vessel geometry and contrast agent dynamics.

4. The system according to claim 2, wherein the parametric model comprises elements corresponding to local vessel segments and local contrast agent dynamics, wherein each of said elements is associated with a local plausibility score value.

5. The system according to claim 4, wherein the data processing unit is adapted to reconstruct an instance of the model by determination of an optimal path with respect to the local plausibility score values.

6. The system according to claim 3, wherein the data processing unit is adapted to reconstruct the flow characteristics and the vessel geometry in three dimensions progressively along the vessels.

7. The system according to claim 1, wherein the data processing unit is adapted to predict a function (c(t, {right arrow over (p)})), wherein c represents an amount of contrast agent at time t at a point {right arrow over (p)}, the function c yielding a time-intensity curve of the vessel structure.

8. The system according to claim 1, further comprising an input device, and wherein the data processing unit is adapted to start the reconstruction of the vessel geometry from user-defined seed points.

9. The system according to claim 1, further comprising a display unit for displaying the reconstructed vessel geometry and/or the reconstructed flow characteristics.

10. A method for the determination of the vessel geometry and the flow characteristics in a vessel system, comprising:

a) injecting a bolus of contrast agent into a vessel system;

b) generating a temporal sequence of X-ray projections (P₁, P_k, P_N) of the vessel system during a passage of said bolus from different directions;

c) reconstructing a three-dimensional vessel geometry and flow characteristics of at least a portion of the vessel system with the use of projections generated in a single temporal sequence and predetermined model contrast agent dynamics.

11. The method recited in claim 10, wherein the reconstructing comprises:

extracting at least one set of image features from the received data;

using an optimization algorithm, comparing a predetermined geometric and flow-dynamic model of a vascular subsystem with the extracted image features; and

determining from the comparison a best fit therebetween, the best fit comprising the reconstructed three-dimensional model.

12. The method recited in claim 11, wherein the image features comprise a position, a local orientation, and a local diameter of dark, elongated structures in the received data, and the extracting comprises using an eigenvalue analysis of a Hessian matrix extracted at a plurality of scales.

13. The method recited in claim 11, wherein the image features comprise at least one of a time of bolus arrival, a length of the bolus, and a dispersion of a front of the bolus, and the extracting comprises fitting the image features to a time-intensity curve at a plurality of positions calculated from the received data.

14. The method recited in claim 13, wherein the optimization algorithm using and determining comprise:

dividing the vessel system into a plurality of tubular vessel segments;

fitting parameters extracted from the image features with the predetermined model for each tubular vessel segment;

calculating a plausibility score for the predetermined model, the image features, and the time-intensity curve;

inserting the plausibility score into a three-dimensional map;

calculating a minimal cost path through the three-dimensional map; and

reconstructing a section of the vessel system from the minimal cost path.

15. The method recited in claim 14, wherein the minimal-cost-path calculating comprises using a front propagation algorithm.

16. The method recited in claim 13, wherein the optimization algorithm using and determining comprise using an iterative tracking algorithm adapted to add successive tubular segments, starting from a defined point, a direction of the adding determined from a calculated highest plausibility score.

17. The method recited in claim 16, wherein the plausibility score is calculated as a weighted combination of a consistency of a hypothetical direction of a next tubular segment in three dimensions as compared with all planar directions extracted as image features.

18. The method recited in claim 16, wherein the plausibility score is calculated from the time-intensity curve of a preceding tubular segment in comparison with a time-intensity curve extracted for a next tubular segment in a hypothetical direction.

19. A record carrier on which a computer program for the determination of the vessel geometry and the flow characteristics in a vessel system is stored, said program being adapted to execute a method according to claim 10.

20. A system for determining vessel geometry and flow characteristics in a vessel system, comprising:

an input for receiving X-ray image data representative of a single temporal sequence of differently oriented projections of a vessel system during a passage of a bolus of contrast agent therethrough, the projections including data representative of contrast agent dynamics within the vessel system;

a data processing unit adapted to reconstruct a three-dimensional model of a geometry and flow characteristics of the vessel system from the received data with the use of predetermined modeled contrast agent dynamics; and

an output for outputting the reconstructed three-dimensional model.

21. The system recited in claim 20, wherein the data processing unit comprises a plurality of modules comprising:

a first module comprising an algorithm for extracting at least one set of image features from the received data;

a second module comprising a predetermined geometric and flow-dynamic model of a vascular subsystem; and

a third module comprising an optimization algorithm for comparing the predetermined model of the second module with the extracted image features from the first module and determining from the comparison a best fit therebetween, the best fit comprising the reconstructed three-dimensional model.

22. The system recited in claim 21, wherein the image features comprise a position, a local orientation, and a local diameter of dark, elongated structures in the received data, and the first module algorithm comprises an eigenvalue analysis of a Hessian matrix extracted at a plurality of scales.

23. The system recited in claim 21, wherein the image features comprise at least one of a time of bolus arrival, a length of the bolus, and a dispersion of a front of the bolus, and the first module algorithm comprises fitting the image features to a time-intensity curve at a plurality of positions calculated from the received data.

24. The system recited in claim 23, wherein the third module comprises code segments adapted to:

divide the vessel system into a plurality of tubular vessel segments;

fit parameters extracted from the image features with the predetermined model of the second module for each tubular vessel segment;

calculate a plausibility score for the predetermined model, the image features, and the time-intensity curve;

insert the plausibility score into a three-dimensional map;

calculate a minimal cost path through the three-dimensional map; and

reconstruct a section of the vessel system from the minimal cost path.

25. The system recited in claim 24, wherein the minimal-cost-path calculation comprises using a front propagation algorithm.

26. The system recited in claim 23, wherein the third module comprises an iterative tracking algorithm adapted to add successive tubular segments, starting from a defined point, a direction of the adding determined from a calculated highest plausibility score.

27. The system recited in claim 26, wherein the plausibility score is calculated as a weighted combination of a consistency of a hypothetical direction of a next tubular segment in three dimensions as compared with all planar directions extracted as image features.

28. The system recited in claim 26, wherein the plausibility score is computed from the time-intensity curve of a preceding tubular segment in comparison with a time-intensity curve extracted for a next tubular segment in a hypothetical direction.