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

• US 5,435,310 A
• Filed: 06/23/1993
• Issued: 07/25/1995
• Est. Priority Date: 06/23/1993
• Status: Expired due to Term

• Current Assignee: University of Washington
• Original Assignee: University of Washington
• Inventors: Bookstein, Fred L., Sampson, Paul D., Sheehan, Florence H., Bolson, Edward L.
• Primary Examiner(s): ZELE, KRISTA MICHELLE
• Application Number:
  US08/082,404
• Time in Patent Office: 762 Days
• Field of Search: 128/653.1, 128/653.2, 128/660.07, 382/6, 364/413.13, 364/413.22
• US Class Current: 600/416
• 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/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/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 an average template of tiled sections is mapped onto this center surface. The mapped tiled sections 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 at least one of an inner surface and an outer surface of the heart, at an end systole and at an end diastole of a cardiac cycle;

(c) defining a center surface that is generally midway between said inner surface and said outer surface, said center surface being defined by a plurality of tiled sections comprising polygons;

(d) projecting the tiled sections onto the inner surface and the outer surface to define volumetric elements having opposite ends on said surfaces, each of the volumetric elements containing a different one of the tiled sections;

(e) determining an area for each of the opposite ends of the volumetric elements;

(f) determining volumes for the volumetric elements;

(g) averaging the areas of the opposite ends of the volumetric elements; and

(h) determining characteristic cardiac parameters for the heart as a function of an average area of the ends and of the volume of the volumetric elements.

2. The method of claim 1, wherein the characteristic parameters are spatially related to the volumetric elements and 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.

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

(a) digitizing an analog output signal produced by 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.

4. The method of claim 3, further comprising the step of defining anatomical landmarks in the heart as borders in the imaging data to be used as spatial references.

5. The method of claim 4, wherein the step of defining the center surface comprises the step of mapping a tiled section template onto the center surface to define the tiled sections, the tiled sections on the center surface each having a shape and a size that depend upon individual characteristics of the heart, but being identical in number to the tiled section template, said tiled section template being derived by averaging the anatomical landmark and imaging data for a plurality of different hearts, said tiled section template being mapped onto the center surface of the patient'"'"'s heart, with respect to the spatial references in the heart.

6. The method of claim 1, wherein the tiled sections each comprise triangles and the volumetric elements are prism shaped.

7. The method of claim 1, wherein the step of projecting comprises the step of extending chords from corners of the tiled sections to the inner and outer surfaces, each of said chords extending in a direction normal to an average plane for all tiled sections having corners meeting at the chord, to ensure that the projections of the tiled sections on the inner and outer surfaces do not overlap.

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

(a) imaging at least a portion of a heart to produce a signal that includes first imaging data defining images of an inner surface and an outer surface of the portion of the heart at an end systole and second imaging data defining images of the inner surface and the outer surface of the portion of the heart at the end diastole;

(b) producing a first three-dimensional model of the portion of the heart that was imaged using the first and the second imaging data from the signal produced by the imaging, said first three-dimensional model comprising first model data that defines a first modeled inner surface at end systole and a first modeled inner surface at end diastole;

(c) producing a second three-dimensional model of the portion of the heart using the first and second imaging data, said second three-dimensional model comprising second model data that represent a second modeled inner surface of the portion of the heart at end diastole and a second modeled outer surface of the portion of the heart at end diastole;

(d) defining a first center surface that is disposed substantially midway between the first modeled inner surface and the first modeled outer surfaces, said first center surf;

ace being represented by a plurality of tiled sections;

(e) projecting the shape of the tiled sections onto the first modeled inner and outer surfaces to define a plurality of volumetric elements, by extending chords from vertices of the edges to said surfaces, each of said chords being transverse to a plane having an orientation that is an average of orientations for the tiled sections surrounding the chord, the projected shapes on said surfaces respectively corresponding to first and second ends of the volumetric element having sides aligned with the chords that define it, each tiled surface being contained in a different volumetric element;

(f) for each volumetric element of interest;

(i) determining a mean area of the first and second ends of the volumetric element;

(ii) determining a volume for the volumetric element; and

(iii) as a function of the mean area of the first and second ends and the volume of the volumetric element, determining a range of motion of the cardiac wall at that volumetric element;

(g) defining a second center surface that is disposed substantially midway between the second modeled inner surface and the second modeled outer surfaces, said second center surface being represented by a plurality of tiled sections;

(h) projecting the shape of the tiled sections of step (i) onto the second modeled inner and outer surfaces to define a plurality of volumetric elements, by extending chords from vertices of the edges to said surfaces, each of said chords being transverse to a plane having an orientation that is an average of orientations for said tiled sections surrounding the chord, the projected shapes on said surfaces respectively corresponding to first and second ends of the volumetric element having sides aligned with the chords that define it, each tiled surface being contained in a different volumetric element; and

(i) for each volumetric element of step (h) of interest;

(i) determining a mean area of the first and second ends of said volumetric element;

(ii) determining a volume for said volumetric element; and

(iii) as a function of the mean area of the first and second ends and the volume of said volumetric element, determining a thickness of the cardiac wall at that volumetric element.

9. The method of claim 8, wherein the step of imaging comprises the step of scanning the portion of the heart with one of an ultrasonic probe, a magnetic resonance imaging system, and a cine CT imaging system, to produce an imaging signal corresponding to the shape of a selected surface of the portion of the heart, at a selected time during the cardiac cycle.

10. The method of claim 9, wherein the step of imaging comprises the step of recording the imaging signal while scanning the portion of the heart to represent it as a plurality of cross-sectional slices.

11. The method of claim 9, wherein the step of imaging comprises the step of converting the imaging signal from an analog format into a digital format suitable for driving a display to present a visual image of the portion of the heart that is scanned.

12. The method of claim 11, wherein the step of imaging further comprises the step of tracing an outline of a selected surface of the portion of the heart presented as the visual image in multiple views to produce the first and the second imaging data, said step of tracing including the step of including borders representing anatomical landmarks in the heart in the imaging data for use as reference points.

13. The method of claim 12, wherein the anatomical landmarks include valves, papillary muscles, an apex of the heart, boundaries, and specific anatomical structures of the heart.

14. The method of claim 12, wherein the steps of defining the first and second center surfaces each comprise the step of mapping an average template of tiled sections onto a corresponding one of the first and second center surfaces, locating the average template relative to said corresponding one of the first and second center surfaces based upon the reference points.

15. The method of claim 14, wherein the average template is derived by averaging a plurality of center surface tiled sections for a plurality of other hearts to define a location for each of said tiled sections and a specified number of tiled sections for the average template, said step of mapping generally modifying the shape and size of the tiled sections from the average template to the center surface of the heart, but providing common spatial references for all such hearts.

16. The method of claim 8, wherein the step of defining a center surface comprises the steps of:

(a) determining a plurality of spheres that fill a space between an inner surface and an outer surface of one of the three-dimensional models, each sphere meeting predefined criteria relating to points on said surfaces that are in contact with the sphere; and

(b) connecting centers of each of the spheres to produce polygons that define the center surface for that three-dimensional model.

17. The method of claim 8, wherein the step of defining a center surface comprises the steps of:

(a) filling a space between an inner and an outer surfaces of one of the three-dimensional models with a plurality of tetrahedrons constructed to meet predefined criteria, said plurality of tetrahedrons extending between said two surfaces; and

(b) for each tetrahedron, determining a center polygon that is midway between said surfaces, said center polygon defining the center surface for that three-dimensional model.

18. The method of claim 8 wherein the step of defining a center surface includes the step of dividing the center surface into a predefined number of triangles, said triangles comprising the tiled sections.

19. The method of claim 18, wherein the volumetric elements comprise prisms having three sides extending between the first and second ends.

20. The method of claim 8, further comprising the step of displaying a three-dimensional graphic representation of the portion of the heart, said graphic image including color coding to represent features and the cardiac parameters of the heart, determined using the center surface.

21. The method of claim 20, wherein the color coding of the graphic image represents one of:

(a) a regional abnormality in a cardiac function; and

(b) a regional cardiac function.

22. The method of claim 8, further comprising the steps of:

(a) producing a third three-dimensional model of the portion of the heart using the first and second imaging data, said third three-dimensional model comprising third model data that represent a third modeled inner surface of the portion of the heart at end systole and a third modeled outer surface of the portion of the heart at end systole;

(b) defining a third center surface that is disposed substantially midway between the third modeled inner surface and the third modeled outer surfaces, said third center surface being represented by a plurality of tiled sections;

(c) projecting the shape of the tiled sections of step (b) in this claim onto the third modeled inner and outer surfaces to define a plurality of volumetric elements, by extending chords from vertices of the edges to said surfaces, each of said chords being transverse to a plane having an orientation that is an average of orientations for said tiled sections surrounding the chord, the projected shapes on said surfaces respectively corresponding to first and second ends of the volumetric element having sides aligned with the chords that define it, each tiled surface being contained in a different volumetric element; and

(d) for each volumetric element of step (c) in this claim, of interest;

(i) determining a mean area of the first and second ends of said volumetric element;

(ii) determining a volume for said volumetric element; and

(iii) as a function of the mean area of the first and second ends and the volume of said volumetric element, determining a thickness of the cardiac wall at that volumetric element.

23. The method of claim 22, further comprising the step of determining a thickening of the cardiac wall at a volumetric element by determining a difference between the thickness of the cardiac wall at the volumetric element at end diastole and end systole.

24. A method for determining at least one characteristic cardiac parameter of a patient'"'"'s heart using imaging data produced by imaging the heart, comprising the steps of:

(a) modeling the heart using the imaging data, producing modeling data corresponding to at least one of an inner surface and an outer surface of the heart, at an end systole and at an end diastole of a cardiac cycle;

(b) defining a center surface that is generally midway between the inner and the outer surfaces to which the modeling data correspond, said center surface being defined by a plurality of tiled sections comprising polygons;

(c) projecting the tiled sections onto the inner and the outer surfaces to define volumetric elements having opposite ends on the two surfaces, each volumetric element containing a different one of the tiled sections;

(d) determining an area for the opposite ends of at least one of the volumetric elements;

(e) determining a volume for said at least one volumetric element; and

(f) determining at least one characteristic cardiac parameter for the heart at said at least one volumetric element, as a function of the area of the opposite ends and of the volume of said at least one volumetric element.

25. A method for determining at least one characteristic cardiac parameter of a patient'"'"'s heart using modeling data corresponding to at least one of an inner surface and an outer surface of the heart, at an end systole and at an end diastole of a cardiac cycle, comprising the steps of:

(a) defining a center surface that is generally midway between the inner and the outer surfaces, said center surface being defined by a plurality of tiled sections comprising polygons;

(b) projecting the tiled sections onto the inner and the outer surfaces to define volumetric elements having opposite ends on the inner and the outer surfaces, each of said volumetric elements containing a different one of the tiled sections;

(c) determining an area for the opposite ends of at least one of the volumetric elements;

(d) determining a volume for said at least one volumetric element; and

(e) determining at least one characteristic cardiac parameter for the heart at said at least one volumetric element, as a function of the area of the opposite ends and of the volume of said at least one volumetric element.