Determining cardiac wall thickness and motion by imaging and three-dimensional modeling

• US 5,601,084 A
• Filed: 05/26/1995
• Issued: 02/11/1997
• Est. Priority Date: 06/23/1993
• Status: Expired due to Term

• Current Assignee: University of Washington, Washington University In St Louis
• Original Assignee: University of Washington
• Inventors: Bolson, Edward L., Jin, Huaichuan, Sheehan, Florence H.
• Primary Examiner(s): Lateef, Marvin M.
• Assistant Examiner(s): FIELDS, DERRICK LAMAR
• Application Number:
  US08/451,757
• Time in Patent Office: 627 Days
• Field of Search: 128/653.1, 128/661.04, 128/663.02, 128/660.07
• US Class Current: 600/450
• CPC Class Codes: B82Y 15/00 Nanotechnology for interacting, sensing or actuating, e.g. quantum dots as markers in protein assays or molecular motors
• CPC Class Codes: G06T 17/20 Finite element generation, e.g. wire-frame surface description, tesselation
• CPC Class Codes: G06T 2207/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/10132 Ultrasound image
• CPC Class Codes: G06T 2207/30048 Heart; Cardiac
• CPC Class Codes: G06T 7/12 Edge-based segmentation
• CPC Class Codes: G06T 7/62 of area, perimeter, diameter or volume
• 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

Abstract

A method for imaging and three-dimensional modeling portions of the heart, particularly the left ventricular endocardial and epicardial surfaces, using imaging data. Preferably, a transesophageal ultrasound probe is inserted into an esophagus of a patient to provide multiple plane imaging data at end systole and end diastole during a cardiac cycle. The image planes are then traced along the boundaries of the epicardial and endocardial surfaces to produce sets of data points, which are further processed and expanded through interpolation. These data points are used for modeling the endocardial and epicardial surface at the end systole and end diastole extremes of the cardiac cycle. A center surface is constructed between an inner and outer surface of the modeled surfaces of the left ventricle and either an average template of tiled sections is mapped onto this center surface, or a mesh of triangular tiled sections representing the inner surface is projected on the center surface to define corresponding numbered tiled sections on the center surface. The tiled sections of the center surface are then projected onto the inner and outer surfaces, defining triangular prisms. By determining the volume of the triangular prisms and an average area for their ends, the range of movement, which is equal to the volume divided by the average area, is determined. Similarly, for changes in wall thickness between the endocardial and epicardial surfaces at end diastole and end systole, the thickening of the cardiac wall is determined.

Claims

1. A method for analyzing characteristic cardiac parameters of a patient'"'"'s heart, comprising the steps of:

(a) imaging a heart to produce imaging data;

(b) modeling the heart using the imaging data, producing modeling data corresponding to an inner surface and an outer surface of the heart, at a plurality of times during a cardiac cycle, including at an end systole and at an end diastole of the cardiac cycle;

(c) using the inner surface of the heart produced with the modeling data, creating a mesh of connected polygons that represent the inner surface of the heart, said polygons including a plurality of connected edges that meet at vertices of the polygons;

(d) determining a normal for a face of each of the polygons;

(e) determining a vertex normal for each of the vertices by averaging the normals for the faces of the polygons meeting at the vertices;

(f) extending the vertex normals as a plurality of line segments that intersect the outer surface at a plurality of points, said line segments extending between the vertices of the polygons on the inner surface and the points on the outer surface;

(g) connecting midpoints of the line segments to form a mesh having a plurality of connected polygons with vertices disposed at the midpoints of the line segments, said mesh representing the center surface; and

(h) using the center surface represented by the mesh formed by connecting the midpoints of the line segments, determining characteristic cardiac parameters for the heart.

2. The method of claim 1, further comprising the step of reiteratively adjusting the center surface before it is used in step (h) so that the center surface meets predefined criteria, said predefined criteria indicating that the center surface is substantially centered between the inner and the outer surfaces.

3. The method of claim 2, wherein the step of reiteratively adjusting comprises the steps of:

(a) for each vertex of the connected polygons comprising the mesh representing the center surface;

(i) determining an average normal by averaging normals to the polygons surrounding the vertex with a vertex normal for the vertex, said average normal being then extended as a test line segment that intersects the inner surface and the outer surface;

(ii) comparing a first distance measured between the center surface and one of the inner surface and the outer surface with a sum of the first distance and a second distance, said second distance being measured between the center surface and the outer surface, where both the first distance and the second distance are measured along the test line segment; and

(iii) moving the vertex of the center surface to a midpoint of the test line segment; and

(b) accepting the center surface as thus adjusted, for use in step (h) if a last adjustment of the center surface made has enabled the center surface to meet the predefined criteria, and if not, then repeating steps (a)(i) through (a)(iii).

4. The method of claim 3, wherein the predefined criteria are met if:

(a) the last adjustment of the center surface has not resulted in a connectivity change of the connected polygons extending over the center surface; and

(b) a mean of a ratio of the first distance to the sum of the first distance and the second distance for all of the vertices on the center surface as last adjusted is substantially equal to 0.5.

5. The method of claim 1, wherein the characteristic cardiac parameters for the heart include at least one of:

(a) a range of motion for a wall of the heart; and

(b) a thickness of the wall of the heart.

6. The method of claim 1, wherein the step of imaging the heart includes the steps of:

(a) digitizing an analog output signal of an imaging device to produce a digital signal;

(b) graphically displaying images of the inner and the outer surfaces of the heart, produced using the digital signal; and

(c) tracing borders of the inner and the outer surfaces of the heart on the graphically displayed images to define the imaging data at an end systole and at an end diastole of the cardiac cycle, said center surface being determined using the modeled data for both the end systole and the end diastole.

7. The method of claim 1, wherein the connected polygons comprising both the mesh representing the inner surface and the mesh representing the center surface comprise triangles.

8. The method of claim 1, further comprising the step of determining if a fold exists on the center surface.

9. The method of claim 8, wherein the step of determining if a fold exists comprises the step of determining if a dihedral angle between adjacent polygons on the center surface that are joined along a common edge is less than a predefined value.

10. The method of claim 8, wherein if a fold is found to exist, said method further comprises the steps of:

(a) for each vertex of the polygons comprising the mesh representing the inner surface;

(i) determine a smoothing normal by averaging the vertex normal of said vertex with the vertex normals of surrounding vertices; and

(ii) replacing the vertex normal of said vertex with the smoothing normal; and

(b) using said smoothing normals as the vertex normals for the mesh representing the inner surface in steps (f) through (h).

11. The method of claim 2, further comprising the step of determining if a fold exists in the center surface after the center surface is reiteratively adjusted.

12. The method of claim 11, wherein the step of determining if a fold exists comprises the step of determining if a dihedral angle between adjacent polygons on the center surface that are joined along a common edge is less than a predefined value.

13. The method of claim 11, wherein if a fold is found to exist in the center surface, said method further comprises the steps of:

(a) undoing a last adjustment to the center surface;

(b) for each vertex of the polygons comprising the mesh representing the center surface;

(i) determining a smoothing normal by averaging the vertex normal of the vertex with the vertex normals of adjacent vertices; and

(ii) replacing the vertex normal of the vertex with the smoothing normal; and

(c) using said smoothing normals as the vertex normals for the mesh representing the center surface in determining a subsequent adjustment to the center surface.

14. The method of claim 1, wherein the step of using the center surface to determine characteristic cardiac parameters for the heart comprises the steps of:

(a) for each vertex on the center surface;

(i) determining a normal to each connected polygon surrounding said vertex;

(ii) determining an average vertex normal for said vertex by averaging the normals for the polygons surrounding said vertex;

(iii) extending the average vertex normals to intersect the inner surface and the outer surface at a plurality of points;

(iv) connecting the plurality of points on the inner surface to form a plurality of tiled sections on the inner surface; and

(v) connecting the plurality of points on the outer surface to form a plurality of tiled sections on the outer surface, said tiled sections on the inner surface forming first ends of volumetric elements that extend through the center surface, between the inner surface and the outer surface, and said tiled sections on the outer surface forming second ends of said volumetric elements;

(b) determining a mean area for the first and the second ends of each of the volumetric elements;

(c) determining a volume for each of the volumetric elements; and

(d) as a function of the mean area of the first and second ends and the volume of said volumetric elements, determining a thickness of the cardiac wall at each of the volumetric elements.

15. The method of claim 14, further comprising the step of determining a thickening of the cardiac wall at a volumetric element by analyzing a difference between the thickness of the cardiac wall at the volumetric element at an end diastole and an end systole of the cardiac cycle, where the center surface at the end systole and at the end diastole is used to determine the volumetric elements at the end systole and at the end diastole.

16. A method for determining a center surface that is centered between a first surface and a second surface represents a cardiac wall of a patient'"'"'s heart, for use in analyzing characteristic cardiac parameters, said method comprising the steps of:

(a) defining a first mesh of connected polygons representing the first surface, said connected polygons having vertices and edges along which they are connected;

(b) for each vertex, determining a vertex normal by averaging normals to all of the connected polygons surrounding the vertex;

(c) extending the vertex normals as line segments that intercept the second surface at a plurality of points;

(d) creating a second mesh of connected polygons to represent the center surface, by connecting midpoints of the lines segments, said midpoints comprising vertices of said plurality of connected polygons representing the center surface, and

(e) using said center surface for analyzing characteristic cardiac parameters.

17. The method of claim 16, further comprising the step of reiteratively adjusting the center surface until criteria that determine a quality of the center surface converge to within predefined limits.

18. The method of claim 17, wherein the step of reiteratively adjusting comprises the steps of:

(a) for each vertex of the connected polygons comprising the second mesh representing the center surface;

(i) determining an average normal by averaging normals to faces of the connected polygons surrounding said vertex, said average normal being then extended as a test line segment, to intersect the first surface and the second surface;

(ii) comparing a first distance between the center surface and the first surface with a second distance between the first surface and the second surface, where both the first distance and the second distance are measured along the test line segment; and

(iii) moving said vertex of the second mesh to a midpoint of the test line segment; and

(b) accepting the center surface as thus adjusted for use in analyzing the characteristic cardiac parameters, if a last adjustment of the center surface made has enabled the center surface to meet the criteria, and if not, repeating steps (a)(i) through (a)(iii).

19. The method of claim 18, wherein the criteria for convergence are met if:

(a) the last adjustment of the center surface has not resulted in a connectivity change of the connected polygons comprising the second mesh; and

(b) a mean of a ratio of the first distance to the second distance for all of the vertices on the center surface as last adjusted is substantially equal to 0.5.

20. The method of claim 16, wherein the connected polygons on the first surface and on the center surface comprise triangles.

21. The method of claim 16, further comprising the step of determining if a fold exists in the second mesh representing the center surface.

22. The method of claim 21, wherein the step of determining if a fold exists comprises the step of determining if a dihedral angle between faces of adjacent polygons comprising the second mesh that are joined along a common edge is less than a predefined value.

23. The method of claim 22, wherein if a fold is found to exist, said method further comprises the steps of:

(a) for each vertex of the polygons comprising the first mesh representing the first surface;

(i) determine a smoothing normal by averaging the vertex normals of said vertex with the vertex normal of surrounding vertices; and

(ii) replacing the vertex normal of said vertex with the smoothing normal; and

(b) using said smoothing normals as the vertex normals for the first mesh representing the first surface in steps (c) through (d).

24. The method of claim 17, further comprising the step of determining if a fold exists in the center surface after a last adjustment to the center surface was made.

25. The method of claim 24, wherein the step of determining if a fold exists comprises the step of determining if a dihedral angle between faces of adjacent polygons on the center surface that are joined along a common edge is less than a predefined value.

26. The method of claim 24, wherein if a fold exists in the center surface, said method further comprises the steps of:

(a) undoing a last adjustment to the center surface;

(b) for each vertex of the polygons comprising the second mesh representing the center surface;

(i) determining a smoothing normal by averaging the vertex normal of said vertex with the vertex normals of surrounding vertices; and

(ii) replacing the vertex normal of said vertex with the smoothing normal; and

(c) using said smoothing normals as the vertex normals for the second mesh that represents the center surface in a subsequent adjustment of the center surface.