Methods and apparatus for three-dimensional and maximum intensity projection image reconstruction in a computed tomography system

• US 5,818,896 A
• Filed: 11/18/1996
• Issued: 10/06/1998
• Est. Priority Date: 11/18/1996
• Status: Expired due to Fees

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Hsieh, Jiang
• Primary Examiner(s): Porta, David P.
• Assistant Examiner(s): Bruce, David Vernon
• Application Number:
  US08/751,648
• Time in Patent Office: 687 Days
• Field of Search: 378/15, 378/901
• US Class Current: 378/15
• 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: G01N 2223/419 computed tomograph
• CPC Class Codes: G01N 23/046 using tomography, e.g. computed tomography [CT]
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

The present invention, in one form, is a system for reducing noise artifacts in three-dimensional image reconstruction using data acquired in a helical scan. More specifically, a standard deviation ratio of reconstructed images with and without helical weighting is identified. Such standard deviation ratio is then used to generate filter coefficients, which are then applied to the data in an adaptive smoothing algorithm.

Claims

1. A system for generating an image of an object using data acquired in a helical scan, said system comprising an x-ray source and a detector array, said system configured to:

determine a helical weighting function;

determine a standard deviation ratio of reconstructed images with and without helical weighting;

identify a filter coefficient based at least in part on the determined standard deviation ratio; and

apply the identified filter coefficient to the data.

2. A system in accordance with claim 1 wherein to identify the filter coefficient, said system is further configured to segment the data.

3. A system in accordance with claim 2 wherein to segment the data, said system is configured to perform grey scale thresholding.

4. A system in accordance with claim 1 further configured to generate a Maximum Intensity Projection image.

5. A system in accordance with claim 1 wherein said standard deviation ratio of reconstructed images with and without helical weighting is:

##EQU7## where;

ω

(γ

,β

) is the helical weighting function;

σ

(γ

,β

) is a standard deviation distribution function of the projection;

β

is a gantry angle;

γ

is a detector angle;

L'"'"' is a projection distance between an x-ray source and a point of interest, (r,φ

).

6. A system in accordance with claim 5 further comprising a 3×

3 kernel, and wherein the filter coefficient is identified according to;

##EQU8##

7. A system in accordance with claim 1 further configured to generate a three-dimensional image.

8. A method for generating an image of an object using data acquired in a helical scan, said method comprising the steps of:

determining a helical weighting function;

determining a standard deviation ratio of reconstructed images with and without helical weighting;

identifying a filter coefficient based at least in part on the determined standard deviation ratio; and

applying the identified filter coefficient to the data.

9. A method in accordance with claim 8 further comprising the step of generating a Maximum Intensity Projection image.

10. A method in accordance with claim 8 wherein identifying the filter coefficient comprises the step of segmenting the data.

11. A method in accordance with claim 10 wherein segmenting the data comprises the step of performing grey scale thresholding.

12. A method in accordance with claim 8 further comprising the step of generating a three-dimensional image.

13. A method in accordance with claim 8 wherein the standard deviation ratio of reconstructed images with and without helical weighting is:

##EQU9## where;

ω

(γ

,β

) is the helical weighting function;

σ

(γ

,β

) is a standard deviation distribution function of the projection;

β

is a gantry angle;

γ

is a detector angle;

L'"'"' is a projection distance between an x-ray source and a point of interest, (r,φ

).

14. A method in accordance with claim 13 wherein the filter coefficient is identified according to:

##EQU10##