Method and device for determining a contour in a space of image parameter values

• US 5,559,901 A
• Filed: 06/28/1994
• Issued: 09/24/1996
• Est. Priority Date: 06/29/1993
• Status: Expired due to Fees

• Current Assignee: US Philips Corporation (Koninklijke Philips N.V.)
• Original Assignee: US Philips Corporation (Koninklijke Philips N.V.)
• Inventors: Lobregt, Steven
• Primary Examiner(s): Couso, Jose L.
• Assistant Examiner(s): Bella, Matthew C.
• Application Number:
  US08/268,080
• Time in Patent Office: 819 Days
• Field of Search: 382/199, 382/256, 382/257, 382/258, 382/259, 382/266, 382/269, 382/197, 382/298
• US Class Current: 382/256
• CPC Class Codes: G06T 2207/10072 Tomographic images
• CPC Class Codes: G06T 2207/20192 Edge enhancement; Edge preservation
• CPC Class Codes: G06T 2207/30101 Blood vessel; Artery; Vein; Vascular
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/12 Edge-based segmentation
• CPC Class Codes: G06V 10/755 Deformable models or variational models, e.g. snakes or active contours
• CPC Class Codes: G06V 2201/03 Recognition of patterns in medical or anatomical images

Abstract

In a space having image parameter values, a contour (10) is defined consisting of vertices (11, 12, 13) connected by edges (14, 15, 16). By variation of the position and number of the vertices as a function of the variation of image parameter values in the space, the contour (10) is determined to follow features in the density variation. The movement of vertices is restricted to a direction (P_i) perpendicular to the local direction (t_i) of the contour (10). For smoothing the contour (10), an internal force or energy is defined at each of the vertices (11, 12, 13) which force or energy depends on the angle (c_i) between the edges at the vertex, or of the variation of the angles along a sequence of adjacent vertices. A direct interaction between vertices is thereby avoided.

Claims

1. Method for determining a contour in a space of image parameter values, comprising the following steps implemented on a computer:

accessing a memory storing image parameter values in a space in the form of a matrix of pixel values;

defining a plurality of initial vertices being located on a seed contour in said space, said vertices being connected by edges to neighbouring vertices;

evaluating an energy function comprising an internal portion which is a function of a curvature at each of the vertices and an external portion determined by the variation of the image parameter values at the location of each of said vertices;

determining a final contour by continuously varying the position of the vertices till the energy function reaches a minimum, characterised in that,

for each of the vertices the internal portion of the energy function is dependent on the angle between the edges at the respective vertex and in that the external portion of the energy function is determined by the variation of the image parameter values in a single direction that bisects the angle between the edges at the respective vertex.

2. Method according to claim 1, characterised in that the internal energy portion is derived from the difference in angles of the edges at a series of adjacent vertices.

3. Method according to claim 2, characterised in that said internal energy portion is obtained by convolving said difference in angles with a symmetric discrete filter having a zero frequency component equal to zero.

4. Method according to claim 1, characterised in that the seed contour is open.

5. Method according to claim 1, characterised in that if during minimization of the energy function an edge length exceeds a first predetermined threshold, the edge is replaced by an additional vertex and connecting edges.

6. Method according to claim 1, characterised in that if during minimization of the energy function an edge length falls below a second predetermined threshold, the edge and the connected vertices are replaced by a single vertex.

7. Method according to claim 1, characterised in that a user defined energy distribution is added to the external portion of the energy function.

8. Method according to claim 1, characterised in that firstly a contour is determined in a first space and that, secondly, a contour is determined in a second space having a variation in image parameter values having similar features to the variation of the image parameter values in the first space whereby the seed contour in the second space is defined by transfer of the final contour obtained in the first space.

10. Method according to claim 2, characterized in that the seed contour is open.

11. Method according to claim 3, characterized in that the seed contour is open.

12. Method according to claim 2, characterized in that if during minimization of the energy function an edge length exceeds a first predetermined threshold, the edge is replaced by an additional vertex and connecting edges.

13. Method according to claim 3, characterized in that if during minimization of the energy function an edge length exceeds a first predetermined threshold, the edge is replaced by an additional vertex and connecting edges.

14. Method according to claim 4, characterized in that if during minimization of the energy function an edge length exceeds a first predetermined threshold, the edge is replaced by an additional vertex and connecting edges.

15. Method according to claim 11, characterized in that if during minimization of the energy function an edge length exceeds a first predetermined threshold, the edge is replaced by an additional vertex and connecting edges.

16. Method according to claim 2, characterized in that a user defined energy distribution is added to the external portion of the energy function.

17. Method according to claim 3, characterized in that a user defined energy distribution is added to the external portion of the energy function.

18. Method according to claim 15, characterized in that a user defined energy distribution is added to the external portion of the energy function.

19. The method of claim 1 wherein said image parameter values are obtained from a medical imaging modality.

9. Arrangement for determining a contour in a space of image parameter values, the arrangement comprising

a memory for storing image parameter values in said space in the form of a matrix of pixel values;

a display for displaying a two-dimensional image determined from said pixel values;

a pointing device for indicating initial vertex positions in the displayed image for defining a seed contour;

a computer, having access to said memory and said vertex positions, for evaluating an energy function comprising an internal portion which is a function of a curvature at each of the vertices and an external portion determined by a variation of image parameter values at the location of each of said vertices;

the computer being arranged for determining a final contour by continuously varying the position of the vertices till the energy function reaches a minimum, while

for each of the vertices the internal portion of the energy function is dependent on the angle between the edges at the respective vertex and in that the external portion of the energy function is determined by the variation of image parameter values in a single direction that bisects the angle between the edges at the respective vertex.

20. The arrangement of claim 9 wherein said image parameter values are obtained from a medical imaging modality.