System and method for fused PET-CT visualization for heart unfolding

• US 7,813,535 B2
• Filed: 04/11/2006
• Issued: 10/12/2010
• Est. Priority Date: 04/19/2005
• Status: Expired due to Fees

• Current Assignee: Siemens Medical Solutions USA Incorporated (Siemens AG)
• Original Assignee: Siemens Medical Solutions USA Incorporated (Siemens AG)
• Inventors: Moreau-Gobard, Romain, Florin, Charles
• Primary Examiner(s): Johns; Andrew W
• Assistant Examiner(s): ANSARI, TAHMINA N
• Application Number:
  US11/401,640
• Publication Number:
  US 20060235294A1
• Time in Patent Office: 1,645 Days
• Field of Search: 382/100, 382/128, 382/130, 382/131, 382/132, 382/168, 382/171, 382/173, 382/174, 382/181, 382/190, 382/195, 382/206, 382/276, 382/277, 382/285, 382/286, 382/288, 382/259, 382/293, 382/294, 382/295, 382/296, 378/4, 378/8, 378/18, 378/19, 378/21, 378/23, 378/24, 378/25, 378/26, 378/27, 378/162, 378/164, 378/901, 600/407, 600/425, 128/920, 434/188, 434/211, 434/214, 345/418, 345/419, 345/424, 345/581, 345/582, 345/587, 345/619, 345/649, 345/650, 345/653, 345/654, 33/301, 333/00, 333/18, 333/20
• US Class Current: 382/128
• CPC Class Codes: A61B 5/05 Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves impedance plethysmography A61B5/0295; measuring movement of the entire body or parts thereof A61B5/11; detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof A61B5/24

Abstract

A method for fusing two digitized image datasets for 2-dimensional visualization of the heart, includes providing a first digitized image and a second digitized image, each image acquired from a different imaging modality and each comprising a plurality of intensities corresponding to a domain of points on a 3-dimensional grid, and wherein the first and second images have been registered, segmenting the heart in one of said digitized images, forming a 2-dimensional maximum intensity projection of a surface of said heart from each of said digitized images, and fusing said images by overlaying one 2-dimensional projection of said surface over the other 2-dimensional projection of said surface.

Claims

1. A method of fusing two digitized image datasets for 2-dimensional visualization of the heart, comprising the steps of:

providing a first digitized image and a second digitized image, each image acquired from a different imaging modality and each comprising a plurality of intensities corresponding to a domain of points on a 3-dimensional grid, and wherein the first and second images have been registered;

segmenting the heart in one of said digitized images, wherein the heart segmentation identifies the surface of the heart;

forming a distance map that associates with each point in said heart a distance from said point to the surface of the heart;

using said distance map to model the heart surface by a 3-dimensional ellipsoid, wherein said 3-dimensional ellipsoid is a bi-axial ellipsoid wherein a center is computed as a center of gravity G of all points located on a distance map isosurface of distance 0, a major axis Z is a vector {right arrow over (G)}_M, where M is a point on said isosurface maximizing the length |{right arrow over (G)}_M|, a minor axis X is a vector {right arrow over (G)}_N, where N is a point on said isosurface minimizing the length |{right arrow over (G)}_N|, and a third axis Y orthogonal to said axes Z and X;

forming a 2-dimensional maximum intensity projection of a surface of said heart from each of said digitized images; and

fusing said images by overlaying one 2-dimensional projection of said surface over the other 2-dimensional projection of said surface.

2. The method of claim 1, wherein one of said images is acquired from a positron emission tomography (PET) scan, and the other of said two images is acquired from a computed tomography (CT) scan.

3. The method of claim 2, wherein the heart is segmented in the CT image.

4. The method of claim 1, wherein forming said 2-dimensional maximum intensity projection comprises casting rays from said ellipsoidal surface into said heart to form said 2-dimensional projection, wherein an intensity of a point I(X)=I(u,v) on said 2D projection is a normalized sum of each intensity, I_n along the ray

5. The method of claim 1, comprising using a transverse Mercator projection for forming said 2-dimensional maximum intensity projection of said heart surface from said ellipsoid.

6. The method of claim 3, wherein forming said 2-dimensional maximum intensity projection of said heart surface using said PET image comprises computing corresponding points on said heart surface in said PET image from said CT image using said registration, selecting points on the surface in said PET image, along with points slightly inside the heart and slightly outside the heart, and summing and averaging intensity values along rays cast through said selected points.

7. A method of fusing two digitized image datasets for 2-dimensional visualization of the heart, comprising the steps of:

providing a digitized computed tomography (CT) image and a digitized positron emission tomography (PET) image, each image comprising a plurality of intensities corresponding to a domain of points on a 3-dimensional grid;

segmenting the heart surface in said CT image;

modeling the heart surface with a 3-dimensional ellipsoid and forming a 2-dimensional maximum intensity projection of said heart surface from each of said digitized images, wherein forming a 2-dimensional projection of said heart surface from said CT image comprises casting rays from said ellipsoid into said heart to form a 2-dimensional projection of said heart surface, wherein an intensity of a point I(X)=I(u,v) on said 2D projection is a normalized sum of each intensity, I_n along the ray

8. The method of claim 7, further comprising registering said segmented heart surface from the CT image onto the PET image wherein correspondence formulae from the CT image to the PET image are obtained.

9. The method of claim 8, wherein coordinates of the rays are saved, wherein forming a 2-dimensional projection of said heart surface from said PET image comprises applying said correspondence formulae to said saved ray coordinates from said CT image to said PET image.

10. A program storage device readable by a computer, tangibly embodying a program of instructions executable by the computer to perform the method steps for fusing two digitized image datasets for 2-dimensional visualization of the heart, the method comprising the steps of:

providing a first digitized image and a second digitized image, each image acquired from a different imaging modality and each comprising a plurality of intensities corresponding to a domain of points on a 3-dimensional grid, and wherein the first and second images have been registered;

segmenting the heart in one of said digitized images, wherein the heart segmentation identifies the surface of the heart;

forming a distance map that associates with each point in said heart a distance from said point to the surface of the heart;

using said distance map to model the heart surface by a 3-dimensional ellipsoid;

using a transverse Mercator projection for forming said 2-dimensional maximum intensity projection of said heart surface from said ellipsoid for each of said digitized images; and

fusing said images by overlaying one 2-dimensional projection of said surface over the other 2-dimensional projection of said surface.

11. The computer readable program storage device of claim 10, wherein one of said images is acquired from a positron emission tomography (PET) scan, and the other of said two images is acquired from a computed tomography (CT) scan.

12. The computer readable program storage device of claim 11, wherein the heart is segmented in the CT image.

13. The computer readable program storage device of claim 10, wherein said 3-dimensional ellipsoid is a bi-axial ellipsoid wherein a center is computed as a center of gravity G of all points located on a distance map isosurface of distance 0, a major axis Z is a vector {right arrow over (G)}_M, where M is a point on said isosurface maximizing the length |{right arrow over (G)}_M|, a minor axis X is a vector {right arrow over (G)}_N, where N is a point on said isosurface minimizing the length |{right arrow over (G)}_N|, and a third axis Y orthogonal to said axes Z and X.

14. The computer readable program storage device of claim 12, wherein forming said 2-dimensional maximum intensity projection of said heart surface using said PET image comprises computing corresponding points on said heart surface in said PET image from said CT image using said registration, selecting points on the surface in said PET image, along with points slightly inside the heart and slightly outside the heart, and summing and averaging intensity values along rays cast through said selected points.