System and method for highly attenuating material artifact reduction in x-ray computed tomography

• US 8,280,135 B2
• Filed: 01/20/2010
• Issued: 10/02/2012
• Est. Priority Date: 01/20/2009
• Status: Active Grant

• Current Assignee: Mayo Foundation For Medical Education And Research (Mayo Clinic)
• Original Assignee: Mayo Foundation For Medical Education And Research (Mayo Clinic)
• Inventors: Yu, Lifeng, McCollough, Cynthia H
• Primary Examiner(s): Seth, Manav
• Application Number:
  US12/690,765
• Publication Number:
  US 20100183214A1
• Time in Patent Office: 986 Days
• Field of Search: None
• US Class Current: 382/128
• CPC Class Codes: A61B 6/027 characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/12 Arrangements for detecting or locating foreign bodies
• CPC Class Codes: A61B 6/5258 involving detection or reduction of artifacts or noise
• CPC Class Codes: G06T 11/005 Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
• CPC Class Codes: G06T 2207/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/20036 Morphological image processing
• CPC Class Codes: G06T 2207/30004 Biomedical image processing
• CPC Class Codes: G06T 5/77 Retouching; Inpainting; Scratch removal
• CPC Class Codes: G06T 7/13 Edge detection

Abstract

The present invention is a method for reducing artifacts caused by highly attenuating materials in x-ray computed tomography (“CT”) images. The method includes combining projection views acquired at equivalent view angles to generate a projection plane data set, from which a reformatted projection is produced. The reformatted projection is then processed to detect and segment regions corresponding to objects composed of metals, metal alloys, or other highly attenuating materials. These segmented regions are then removed from the reformatted projection and the removed portions replaced by attenuation information interpolated from portions of the reformatted projection adjacent the removed portions. The interpolated reformatted projection is then mapped back to a projection plane data set, and an image of the subject is reconstructed from the projection views contained in that data set. The reconstructed image, therefore, is one in which artifacts caused by highly attenuating materials are substantially suppressed.

Claims

1. A method for producing, with an x-ray computed tomography (“

CT”

) system, an image of a subject in which artifacts associated with a highly attenuating material are substantially suppressed, the steps of the method comprising;

a) acquiring, with the CT system and using a helical-scan with a selected pitch factor, projection data that includes a plurality of projection views of a subject at a plurality of view angles over a selected longitudinal scan range;

b) producing, from the acquired projection data, a projection plane data set by combining a plurality of projection views from the projection data that were acquired at the same view angle;

c) producing a reformatted projection by mapping attenuation information in the projection plane data set to an image matrix;

d) correcting the reformatted projection for errors resulting from the selected pitch factor;

e) segmenting a region of the corrected reformatted projection that corresponds to an object composed of a highly attenuating material;

f) removing attenuation information in the segmented region from the corrected reformatted projection;

g) interpolating attenuation information from portions of the corrected reformatted projection into the segmented region to produce an interpolated reformatted projection; and

h) producing an image of the subject, in which artifacts caused by highly attenuating materials are substantially suppressed, by mapping the interpolated reformatted projection to a corrected projection plane data set and performing image reconstruction on the projection plane data set.

2. The method as recited in claim 1 in which step e) includes:

e)i) producing an edge-detected image by performing edge detection on the corrected reformatted projection;

e)ii) producing a dilated image by performing morphological dilation on the edge-detected image; and

e)iii) producing a hole-filled image by performing hole-filling on the dilated image.

3. The method as recited in claim 2 in which the hole-filled image is segmented by comparing regions in the hole-filled image with corresponding regions in the corrected reformatted projection.

4. The method as recited in claim 3 in which the regions in the hole-filled image that correspond to regions in the corrected reformatted projection having attenuation values above a selected threshold are selected as the detected region.

5. A method for producing, with an x-ray computed tomography (“

CT”

) system, an image of a subject in which artifacts associated with a highly attenuating material are substantially suppressed, the steps of the method comprising;

a) acquiring, with the CT system, projection data that includes a plurality of projection views of a subject at a plurality of view angles over a selected longitudinal scan range;

b) producing, from the acquired projection data, a reformatted projection by combining a plurality of projection views from the projection data that were acquired at a common view angle;

c) detecting a region in the produced reformatted projection that corresponds to an object composed of a highly attenuating material;

d) producing a corrected reformatted projection from the produced reformatted projection using the detected region; and

e) reconstructing, from the corrected reformatted projection, an image of the subject in which artifacts caused by highly attenuating materials are substantially suppressed.

6. The method as recited in claim 5 in which step c) includes producing a pitch-corrected reformatted projection by correcting for variations in the produced reformatted projection resulting from a pitch of the CT system utilized when acquiring the projection data.

7. The method as recited in claim 5 in which step a) is performed by a helical-scan.

8. The method as recited in claim 5 in which step a) is performed by:

a)i) acquiring a plurality of projection views of the subject at the plurality of view angles at a longitudinal position; and

a)ii) repeating step a)i) a plurality of times in order to acquire the projection data, during each repetition acquiring a different plurality of projection views of the subject at the plurality of view angles at a different longitudinal position.

9. The method as recited in claim 5 wherein step c) includes performing at least one of edge detection, morphological dilation, hole-filling, and segmentation.

10. The method as recited in claim 9 in which step c) includes:

c)i) producing an edge-detected image by performing edge detection on the reformatted projection;

c)ii) producing a dilated image by performing morphological dilation on the edge-detected image; and

c)iii) producing a hole-filled image by performing hole-filling on the dilated image.

11. The method as recited in claim 10 in which edge detection is performed in step c)i) using a Sobel operator.

12. The method as recited in claim 10 in which the detected region is detected by segmenting the hole-filled image.

13. The method as recited in claim 12 in which the hole-filled image is segmented by comparing regions in the hole-filled image with corresponding regions in the reformatted projection.

14. The method as recited in claim 13 in which the regions in the hole-filled image that correspond to regions in the reformatted projection having attenuation values above a selected threshold are selected as the detected region.

15. The method as recited in claim 14 in which the selected threshold is 1000 Hounsfield units.

16. The method as recited in claim 5 in which step d) includes removing attenuation information associated with the detected region from the reformatted projection and replacing the removed attenuation information with attenuation information interpolated from portions of the reformatted projection adjacent the detected region.

17. The method as recited in claim 16 in which step d) includes performing interpolation using Delaunay triangulation.

18. The method as recited in claim 5 in which step b) includes:

b)i) forming a projection plane data set by combining the plurality of projection views from the projection data that were acquired at the same view angle; and

b)ii) mapping attenuation intensity information in the projection plane data set to an image matrix having a first image matrix index corresponding to a longitudinal position and a second image matrix index corresponding to a detector element.

19. The method as recited in claim 18 in which step e) includes:

e)i) mapping attenuation information in the corrected reformatted projection to a corrected projection plane data set that includes a plurality of corrected projection views; and

e)ii) reconstructing the image of the subject from the plurality of corrected projection views.

20. The method as recited in claim 5 in which step e) includes mapping the corrected reformatted projection into a plurality of projection views and reconstructing the image of the subject from the plurality of projection views.