Method of and system for adaptive scatter correction in multi-energy computed tomography

• US 7,136,450 B2
• Filed: 05/26/2004
• Issued: 11/14/2006
• Est. Priority Date: 05/26/2004
• Status: Expired due to Fees

• Current Assignee: Analogic Corporation
• Original Assignee: Analogic Corporation
• Inventors: Crawford, Carl R., Simanovsky, Sergey, Ying, Zhengrong, Naidu, Ram
• Primary Examiner(s): Glick, Edward J.
• Assistant Examiner(s): TANINGCO, ALEXANDER H
• Application Number:
  US10/853,942
• Publication Number:
  US 20050276373A1
• Time in Patent Office: 902 Days
• Field of Search: 378/4, 378/7, 382/131, 382/261
• US Class Current: 378/7
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/5282 due to scatter
• CPC Class Codes: G01N 2223/419 computed tomograph
• CPC Class Codes: G01N 23/046 using tomography, e.g. computed tomography [CT]
• 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 2211/408 Dual energy

Abstract

Method of and system for adaptive scatter correction in the absence of scatter detectors in multi-energy computed tomography are provided, wherein input projection data acquired using at least two x-ray spectra for scanned objects may include a set of low energy projections and a set of high energy projections; wherein a low-pass filter of variable size is provided; the method comprises estimating the size of the low-pass filter; computing amounts of scatter; and correcting both sets of projections for scatter. The estimation of low-pass filter size comprises thresholding high energy projections into binary projections; filtering the binary projections; finding the maximum of the filtered binary projections; calculating the low-pass filter size from the found maximum. The computation of amounts of scatter comprises exponentiating input projections; low-pass filtering the exponentiated projections with the estimated filter size; computing the amounts of scatter from the filtered projections.

Claims

1. A method of adaptive scatter correction in the absence of scatter detectors for the input projection data in multi-energy X-ray computed tomography, wherein the input projection data include a set of low energy projections and a set of high energy projections acquired by scanning a set of objects using at least two x-ray spectra, wherein a low-pass filter of variable size is provided, comprising:

A. Estimating the filter size of the low pass filter as a function of the input projection data, including;

A1. Thresholding at least one set of projections into binary projections;

A2. Filtering the binary projections;

A3. Finding maximum of the filtered binary projections;

A4. Calculating the filter size from the maximum found in A3;

B. Computing amounts of scatter in the input projection data using the low-pass filter with the calculated filter size from A4;

C. Correcting the input projection data for the scatter using the amounts of scatter computed from B.

2. The method of claim 1, wherein C includes subtracting the computed amounts of scatter from both sets of projections to yield scatter corrected projections.

3. The method of claim 1, wherein the low-pass filter may be a one-dimensional odd-number-point long rectangular filter so that the estimating step A includes estimating the filter size of such a filter.

4. The method of claim 1, wherein A includes using high energy projections for estimating the filter size.

5. The method of claim 1, wherein A includes setting the calculated filter size to a pre-defined maximal value if the calculated filter size exceeds the predefined maximal value.

6. A method of adaptive scatter correction in the absence of scatter detectors for the input projection data in multi-energy X-ray commuted tomography, wherein the input projection data include a set of low energy projections and a set of high energy projections acquired by scanning a set of objects using at least two x-ray spectra, wherein a low-pass filter of variable size is provided, comprising:

A. Estimating the filter size of the low pass filter as a function of the input projection data;

B. Computing amounts of scatter in the input projection data using the low-pass filter with the estimated filter size from A, including;

B1. Exponentiating both sets of projections so as to provide exponentiated projections;

B2. Filtering the exponentiated projections with the estimated filter size;

B3. Computing amounts of scatter from the filtered projections; and

C. Correcting the input projection data for the scatter using the amounts of scatter computed from B3.

7. The method of claim 6, wherein B3 includes multiplying the filtered projections with a larger scale factor for low energy projections and a smaller scale factor for high energy projections.

8. A system for adaptive scatter correction in the absence of scatter detectors for the input projection data in multi-energy X-ray computed tomography, wherein the input projection data may include a set of low energy projections and a set of high energy projections acquired by scanning a set of objects using at least two x-ray spectra, wherein a low-pass filter of variable size is provided, comprising:

A. A subsystem constructed and arranged so as to estimate the filter size of the low-pass filter as a function of the input projection data, including;

A1. A module constructed and arranged so as to threshold one set of projections into binary projections;

A2. A module constructed and arranged so as to filter the binary projections;

A3. A module constructed and arranged so as to find maximum of the filtered binary projections;

A4. A module constructed and arranged so as to calculate the filter size from the maximum found by the module of A3.

B. A subsystem constructed and arranged so as to compute amounts of scatter in the input projection data using the low-pass filter with the calculated filter size from the module of A4;

C. A subsystem constructed and arranged so as to correct the input projection data for scatter using the amounts of scatter computed from the subsystem of B.

9. The system of claim 8, wherein subsystem C includes a subtractor constructed and arranged so as to subtract the computed amounts of scatter from both sets of projection data so as to yield scatter corrected projections.

10. The system of claim 8, wherein the low-pass filter may be a one-dimensional odd-number-point long rectangular filter.

11. The system of claim 8, wherein said subsystem is constructed and arranged so as to estimate the filter size of the low-pass filter as a function of the high energy input projection data.

12. The system of claim 8, wherein said module A4 includes a feature for setting the calculated filter size to a pre-defined maximal value if the calculated filter size exceeds the pre-defined maximal value.

13. A system for adaptive scatter correction in the absence of scatter dotectors for the input projection data in multi-energy X-ray computed tomography, wherein the input projection data may include a set of low energy projections and a set of high energy projections acquired by scamming a set of objects using at least two x-ray spectra, wherein a low-pass filter of variable size is provided, comprising:

A.A subsystem constructed and arranged so as to estimate the filter size of the low-pass filter as a function of the input projection data;

B. Computing amounts of scatter in the input projection data using the low-pass filter with the estimated filter size from the subsystem of A, including;

B1. A module constructed and arranged so as to exponentiate both sets of input projection data;

B2. A module constructed and arranged so as to filter the exponentiated projections with the estimated filter size;

B3. A module constructed and arranged so as to compute amounts of scatter from the filtered projection data; and

C. A subsystem constructed and arranged so as to correct the input projection data for scatter using the amounts of scatter computed by the module of B3.

14. The system of claim 13, wherein said module B3 includes a multiplier constructed and arranged so as to multiply the filtered projections with a larger scale factor for low energy projections and a smaller scale factor for high energy projections.