Image processing method, system and apparatus for processing an image representing tubular structure and for constructing a path related to said structure

• US 6,496,188 B1
• Filed: 01/04/2000
• Issued: 12/17/2002
• Est. Priority Date: 01/04/1999
• Status: Expired due to Term

• Current Assignee: Koninklijke Philips N.V.
• Original Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Inventors: Deschamps, Thomas, Cohen, Laurent, Makram-Ebeid, Shérif
• Primary Examiner(s): Vo, Cliff N.
• Application Number:
  US09/477,221
• Time in Patent Office: 1,078 Days
• Field of Search: 345/418, 345/419, 345/420, 345/423, 345/424, 345/427
• US Class Current: 345/419
• CPC Class Codes: G06T 2207/10068 Endoscopic image
• CPC Class Codes: G06T 2207/30172 Centreline of tubular or elongated structure
• CPC Class Codes: G06T 7/13 Edge detection
• CPC Class Codes: G06V 10/457 by analysing connectivity, e.g. edge linking, connected component analysis or slices
• CPC Class Codes: G06V 2201/03 Recognition of patterns in medical or anatomical images

Abstract

The invention relates to an image processing method for processing an image representing a tubular structure having walls, comprising steps for determining a flight path (EP) inside this tubular structure between a first and a second predetermined end points, said flight path being both the shortest path between said end points and the farthest from the structure walls (NB1). Said steps may comprise locating the structure wall points, determining a surface at a predetermined constant distance from said wall points, inside the structure, for forming a central region (CR), and determining, in said central region, the shortest path (EP) between the first and second end points.

Claims

1. An image processing method for processing an image representing a tubular structure having walls, comprising the steps of:

determining a first end point and a second end point for a flight path inside the tubular structure; and

determining the flight path inside the tubular structure between the first and second end points which is the shortest path between the first and second end points at a maximum constant distance from the walls of the tubular structure irrespective of the shape of the tubular structure.

2. A method as claimed in claim 1, wherein the step of determining the flight path comprises the steps of:

determining the location of the structure wall points,

determining a surface at a predetermined constant distance from the wall points to thereby form a central region inside the structure, and

determining the shortest path between the first and second end points in the central region whereby the shortest path thus-determined constitutes the flight path.

3. A method as claimed in claim 2, wherein the step of determining the location of the structure wall points comprises the steps of:

imposing substantially low gray level values to internal voxels of the structure, and

propagating a first front from a first initial front reduced to the first end point, based on gray level values, which determines a first terminal front formed of the structure wall points when the second end point is reached.

4. A method as claimed in claim 3, wherein the step of determining a surface at a predetermined distance from the wall points comprises the steps of:

propagating a second front from a second initial front which is the first terminal front, based on distance values, towards the center of the structure, which determines a second terminal front whose voxels are equidistant to the second initial front, which also determines the central region inside the second terminal front.

5. A method as claimed in claim 4, wherein the step of determining the shortest path between the first and second end points in the central region comprises the step of:

propagating a third front inside the central region from the first end point to the second end point, thus determining the shortest path which is the maximum distance from the structure walls.

6. A method as claimed in claim 5, wherein the step of propagating the third front is performed based on a distance map calculated from the propagation of the second front which is subjected to a threshold and then inverted.

7. A method as claimed in claim 5, further comprising the step of performing the propagation of the first and second fronts along surfaces of minimal energy integrated between the initial fronts and the terminal fronts according to the Eikonal Equation.

8. A method as claimed in claim 5, further comprising the steps of:

basing the propagation of the first front on gray level values and voxel distances,

basing the propagation of the second front on the distance from the current front to the second initial front in the propagation of the second front, and

basing the propagation of the third front on the propagation time between the first and second end points in the central region.

9. A system comprising a programmed computer of a workstation or a special purpose processor having circuit means, which are arranged to perform the method as claimed in claim 1, and comprising:

means to display images processed according to said method, and

virtual camera means to perform virtual endoscopy along the flight path constructed using said method.

10. An imaging apparatus comprising:

means to acquire 3-D medical images of tubular structures; and

a system as claimed in claim 9 to perform virtual endoscopy in the tubular structures.

11. A method as claimed in claim 3, further comprising the step of performing the propagation of the first front along surfaces of minimal energy integrated between the first initial front and the first terminal front according to the Eikonal Equation.

12. A method as claimed in claim 11, further comprising the step of basing the propagation of the first front on gray level values and voxel distances.

13. A method as claimed in claim 4, further comprising the step of performing the propagation of the second front along surfaces of minimal energy integrated between the initial fronts and the terminal fronts according to the Eikonal Equation.

14. A method as claimed in claim 13, further comprising the step of:

basing the propagation of the second front on the distance from the current front to the second initial front in the propagation of the second front.

15. A method as claimed in claim 4, further comprising the step of performing the propagation of the first and second fronts along surfaces of minimal energy integrated between the initial fronts and the terminal fronts according to the Eikonal Equation.

16. A method as claimed in claim 15, further comprising the steps of:

basing the propagation of the first front on gray level values and voxel distances, and

basing the propagation of the second front on the distance from the current front to the second initial front in the propagation of the second front.

17. A method as claimed in claim 1, wherein the step of determining the flight path comprises the step of first determining a surface at the maximum, constant distance from the walls of the structure and then determining the shortest path along the surface which will constitute the flight path.

18. A method as claimed in claim 17, wherein the step of determining the surface at the maximum, constant distance from the walls of the structure comprises the steps of determining edges of a central region which are equidistant from edges of the structure.