Reconstruction of images from three-dimensional cone beam data

• US 6,219,441 B1
• Filed: 06/22/1993
• Issued: 04/17/2001
• Est. Priority Date: 06/22/1993
• Status: Expired due to Term

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Hu, Hui
• Primary Examiner(s): CHANG, JON CARLTON
• Application Number:
  US08/081,040
• Time in Patent Office: 2,856 Days
• Field of Search: 382/6, 382/131, 382/254, 382/275, 364/413.14, 364/413.15, 364/413.16, 364/413.17, 364/413.18, 364/413.19, 364/413.21, 378/4, 378/14, 378/15, 378/901
• US Class Current: 382/131
• CPC Class Codes: G06T 11/008 Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
• CPC Class Codes: G06T 2211/421 Filtered back projection [FBP]
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

A computed tomography x-ray imaging system acquires a three-dimensional array of x-ray attenuation values using a cone beam x-ray source and a curved two-dimensional array of detector elements. Two-dimensional image slices are reconstructed using a filtered back projection method, and corrections are made to the images to account for incomplete data due to the cone beam geometry.

Claims

1. A computed tomography imaging system which comprises:

a two-dimensional array of detector elements for receiving photons emanating in a cone beam from a source;

a digital acquisition system for acquiring two-dimensional arrays of cone beam data from the array of detector elements at a series of views in which the array of detector elements revolves around a central axis;

a filter for receiving the cone beam data and filtering the same;

means for back projecting the filtered cone beam data to produce image data f_d({right arrow over (r)});

means for receiving the cone beam data and producing from the received cone beam data estimated cone beam values not provided by the received cone beam data;

means for receiving the estimated cone beam values and calculating from the estimated cone beam values corrected image data f_C({right arrow over (r)}); and

summing means for combining the image data f_D({right arrow over (r)}) with the corrected image data f_C({right arrow over (r)}) to produce an image slice.

2. The system as recited in claim 1 in which the means for estimating cone beam values includes means for determining the boundaries in Radon space between regions therein for which complete cone beam data is received and regions therein for which cone beam data is incomplete, and means for interpolating between complete cone beam data at the Radon space boundaries to provide said estimated values.

3. The system as recited in claim 1 which includes means for weighting the cone beam data, based on the distance from the source to the detector elements, prior to its application to the filter.

4. The system as recited in claim 1 in which the source produces x-rays and is located on the side opposite the array of detector elements from the central axis, and both the array of detector elements and x-ray source are revolved around the central axis during the acquisition of said series of views.

5. The system as recited in claim 4 in which the acquired cone beam data is preprocessed to provide x-ray attenuation data for a three-dimensional region about the central axis and between the x-ray source and the array of detector elements, and the image data f_D({right arrow over (r)}) and correction image data f_C({right arrow over (r)}) are summed to produce a plurality of two-dimensional image slices through said three-dimensional region.

6. A method for producing an image with a computed tomography imaging system having a two-dimensional array of detector elements, the steps comprising:

a) acquiring two-dimensional arrays of cone beam data from the array of detector elements at a series of views in which the array of detector elements revolves around a central axis;

b) filtering the acquired cone beam data;

c) back projecting the filtered cone beam data to produce image data f_D({right arrow over (r)});

d) producing estimated cone beam values not provided by the acquired cone beam data, the estimated cone beam values being calculated from the acquired cone beam data;

e) calculating corrected image data f_C({right arrow over (r)}) from the estimated cone beam values; and

f) producing an image slice by combining the image data f_D({right arrow over (r)}) with the corrected image data f_C({right arrow over (r)}).

7. The method as recited in claim 6 in which step d) includes:

i) determining boundaries in Radon space between regions therein for which complete cone beam data was acquired and regions therein for which cone beam data was not acquired; and

ii) interpolating between complete cone beam data at the Radon space boundaries to produce said estimated cone beam values.

8. The method as recited in claim 6 which includes:

weighting the acquired cone beam data prior to filtering in step b).

9. The method as recited in claim 6 which includes:

preprocessing the acquired cone beam data to produce x-ray attenuation data for a three-dimensional region about the central axis; and

producing a plurality of two-dimensional image slices through said three-dimensional region in step f).