Apparatus and method for generation of attenuation map

Patent metadata

Field	Value
Publication number	US-9619905-B2
Application number	US-86424608-A
Country	US
Kind code	B2
Filing date	Oct 29, 2008
Priority date	Nov 9, 2007
Publication date	Apr 11, 2017
Grant date	Apr 11, 2017

Abstract

A method of image reconstruction corrected for attenuation is provided for use with radioemission-based imaging, such as SPECT and PET. This method includes collecting measured emission projection data. The emission projection data, a reconstruction of the emission projection data, and a priori organ information are collectively analyzed to generate a body region estimate of the imaged subject. Each voxel of the body region estimate is then homogenously assigned an attenuation coefficient to generate an initial attenuation map estimate. An initial emission assumption is also generated based on a reconstruction of the emission projection data. The initial emission assumption and initial attenuation map estimate are then processed and refined to produce an image reconstruction.

First claim

Having thus described the preferred embodiments, the invention is now claimed to be:

1. A method of generating an attenuation map estimate for use in attenuation correction with an image reconstruction of emission projection data of an imaged subject using of a medical diagnostic imaging apparatus, the method comprising:
generating, by an imaging processor, an outer body contour estimate of the imaged subject based on the emission projection data;
generating, by the imaging processor, an organ boundary estimate of the imaged subject based on the emission projection data;
generating, by the imaging processor, a three-dimensional body region estimate of the imaged subject based upon the outer body contour estimate and the organ boundary estimate, wherein the body region estimate is made up of a plurality of voxels, each voxel being assigned one of a plurality of body region indicators;
assigning, by the imaging processor, each voxel of the body region estimate an attenuation coefficient based upon the body region indicator of that voxel, wherein the same attenuation coefficient is assigned to each voxel with the same body region indicator or wherein the attenuation coefficient assigned to each voxel is selected from a range of predetermined attenuation coefficients, wherein the range of predetermined attenuation coefficients corresponds to the body region indicator of the voxel;
generating, by the imaging processor, an attenuation map estimate based on the assigned attenuation coefficients;
iteratively updating, by the imaging processor, the attenuation map estimate during a reconstruction of the emission projection data;
generating, by the imaging processor, a reconstructed image of the imaged subject based on the updated attenuation map; and
displaying, by a display device, the reconstructed image of the imaged subject.

2. The method of claim 1 , wherein the step of generating an organ boundary estimate comprises the step of analyzing a reconstructed image of the emission projection data for areas of activity.

3. The method of claim 2 , wherein the step of generating a three-dimensional body region estimate further comprises the step of refining the organ boundary estimate by analyzing a priori organ knowledge.

4. The method of claim 1 , wherein the same attenuation coefficient is assigned to each voxel with the same body region indicator.

5. The method of claim 1 , wherein the attenuation coefficient assigned to each voxel is selected from a range of predetermined attenuation coefficients, wherein the range of predetermined attenuation coefficients corresponds to the body region indicator of the voxel.

6. The method of claim 1 , wherein the attenuation coefficient assigned to each voxel is selected based upon at least one physical characteristic of the imaged subject.

7. The method of claim 1 , wherein the attenuation map is generated for use in attenuation correction with a SPECT image reconstruction or a PET image reconstruction.

8. The method of claim 1 , wherein the step of generating a three-dimensional body region estimate further comprises the step of refining the organ boundary estimate by analyzing a priori organ knowledge.

9. The method of claim 8 , wherein the step of generating a three-dimensional body region estimate further comprises the step of refining the organ boundary estimate by spatially relating the organ boundary estimate with the outer body contour estimate.

10. The method of claim 1 , further comprising:
collecting the emission projection data;
generating, by the imaging processor, the emission map estimate from a reconstruction of the collected emission projection data; and
iteratively reconstructing, by the imaging processor, the emission map estimate.

11. The method of claim 10 , further comprising simultaneous iterative reconstruction of the emission map estimate and iterative updating of the attenuation map estimate.

12. The method of claim 11 , wherein each iteration of the step of iteratively reconstructing the emission map estimate and iteratively updating the attenuation map estimate comprises:
forward projecting, by the imaging processor, the emission map estimate and the attenuation map estimate to generate an estimated emission projection;
comparing, by the imaging processor, the estimated emission projection with the emission projection data to generate estimated emission projection corrections;
back projecting, by the imaging processor, the estimated emission projection corrections and the attenuation map estimate to generate back-projected emission projection corrections; and
generating, by the imaging processor, an updated attenuation map estimate by analyzing the emission projection data, the attenuation map estimate and effective source information to calculate corrections for the attenuation map estimate.

13. The method of claim 10 , wherein the step of generating an organ boundary estimate comprises the step of analyzing a reconstructed image of the emission projection data for areas of activity.

14. The method of claim 13 , wherein the step of generating a three-dimensional body region estimate further comprises the step of refining the organ boundary estimate by analyzing a priori organ knowledge.

15. The method of claim 10 , wherein each body region indicator corresponds to a single attenuation coefficient.

16. The method of claim 10 , wherein the attenuation coefficient assigned to each voxel is selected from a range of predetermined attenuation coefficients corresponding to the body region indicator of the voxel.

17. The method of claim 10 , wherein the attenuation coefficient assigned to each voxel is selected based upon at least one physical characteristic of the imaged subject.

18. The method of claim 10 , wherein the medical diagnostic imaging apparatus comprises a SPECT imaging apparatus or a PET imaging apparatus.

19. The method of claim 10 , wherein the step of generating a three-dimensional body region estimate further comprises the step of refining the organ boundary estimate by analyzing a priori organ knowledge.

20. The method of claim 19 , wherein the step of generating a three-dimensional body region estimate further comprises the step of refining the organ boundary estimate by spatially relating the organ boundary estimate with the outer body contour estimate.

21. A medical diagnostic imaging apparatus, comprising:
an emission memory configured to store emission projection data;
an image memory configured to store an emission map estimate generated from a reconstruction of the emission projection data;
an attenuation memory configured to store an attenuation map estimate;
a reconstruction data processor configured to iteratively reconstruct the emission map estimate and iteratively update the attenuation map estimate to obtain a final emission map, comprising logic that instructs the reconstruction data processor to:
generate an outer body contour estimate of an imaged subject based on the emission projection data;
generate an organ boundary estimate of the imaged subject based on the emission projection data;
generate a three-dimensional body region estimate of the imaged subject based upon the outer body contour estimate and the organ boundary estimate, wherein the body region estimate is made up of a plurality of voxels, each voxel being assigned one of a plurality of body region indicators;
assign each voxel of the body region estimate an attenuation coefficient based upon the body region indicator of that vox

Assignees

• Salomon Andre
• Goedicke Andreas
• Koninklijke Philips Nv

Inventors

• Salomon Andre
• Goedicke Andreas

Classifications

• G06T12/10 Primary

  Image preprocessing, e.g. calibration, positioning of sources or scatter correction · CPC title

• G06T7/149

  involving deformable models, e.g. active contour models · CPC title

• G06T7/11

  Region-based segmentation · CPC title

• G06T2211/424

  Iterative · CPC title

• G06T2207/30004

  Biomedical image processing · CPC title