Scatter correction methods and systems in single photon emission computed tomography

• US 5,903,008 A
• Filed: 07/02/1997
• Issued: 05/11/1999
• Est. Priority Date: 07/02/1997
• Status: Expired due to Term

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Li, Jianying
• Primary Examiner(s): Hannaher, Constantine
• Assistant Examiner(s): Jiron, Darren M.
• Application Number:
  US08/888,341
• Time in Patent Office: 678 Days
• Field of Search: 250/363.02, 250/363.04, 250/363.07, 250/369, 378/7, 378/86
• US Class Current: 250/363.04
• CPC Class Codes: G01T 1/1647 Processing of scintigraphic data not related to a particular imaging system G01T1/2992

Abstract

An emission tomographic system for imaging an object of interest is described. The system, in one form, includes a gantry and a patient table. A detector including a collimator is secured to the gantry, and a computer is coupled to the gantry and to the detector to detect and control the position of the detector relative to the table. The system is configured to determine a transmission measurement and generate a scatter fraction utilizing the transmission measurement. A dual energy window data acquisition algorithm then determines non-scatter photons in a primary energy window utilizing the scatter fraction.

Claims

1. An emission tomographic system for imaging an object of interest, said system comprising a gantry having at least one detector secured thereto, a patient table, at least one transmitting source, and a computer coupled to the gantry and to the detector to detect and control the position of the detector relative to the table, said transmitting source configured to transmit gamma rays toward said detector, said system configured to:

determine a transmission measurement; and

generate a three-dimensional scatter fraction utilizing the transmission measurement.

2. An emission tomographic system in accordance with claim 1 wherein said transmission measurement is determined in accordance with:

space="preserve" listing-type="equation">T/T.sub.0 =exp(-μ

.sub.TCT D)

where;

T is a projection of a transmission scan;

T₀ is a projection of a blank scan;

D is a water-equivalent total body width that the transmission source passes through; and

μ

_TCT is a transmission source water attenuation coefficient.

3. An emission tomographic system in accordance with claim 2 wherein said scatter fraction is generated in accordance with:

space="preserve" listing-type="equation">k(i,j,l)=A-B* T/T.sub.0 !.sup.C/2

where;

k(i,j,l) is said scatter fraction;

A is a first constant;

B is a second constant; and

C is a third constant.

4. An emission tomographic system in accordance with claim 3 wherein if T/T₀ equals zero, then said scatter fraction is generated in accordance with:

space="preserve" listing-type="equation">T.sub.0 *exp (-μ

.sub.TCT *2*R),

where R is a radius of rotation of the detector.

5. An emission tomographic system in accordance with claim 3 wherein at least one of said first constant A, said second constant B, and said third constant C is stored in a memory of said computer.

6. An emission tomographic system in accordance with claim 1 further configured to reduce photon scatter utilizing a dual energy window data acquisition algorithm.

7. A method for performing a tomographic scan of an object with an emission tomographic system including a gantry having at least one detector and at least one transmitting source secured thereto, a patient table, and a processor coupled to the gantry and to the detector to detect and control the position of the detector relative to the table, said method comprising the steps of:

determining a transmission measurement; and

generating a three-dimensional scatter fraction utilizing the transmission measurement.

8. A method in accordance with claim 7 wherein said transmission measurement is determined in accordance with:

space="preserve" listing-type="equation">T/T.sub.0 =exp(-μ

.sub.TCT D)

where;

T is a projection of a transmission scan;

T₀ is a projection of a blank scan;

D is a water-equivalent total body width that the transmission source passes through; and

μ

_TCT is a transmission source water attenuation coefficient.

9. A method in accordance with claim 8 wherein said scatter fraction is generated in accordance with:

space="preserve" listing-type="equation">k(i,j,l)=A-B* T/T.sub.0 !.sup.C/2

where;

k(i,j,l) is said scatter fraction;

A is a first constant;

B is a second constant; and

C is a third constant.

10. A method in accordance with claim 9 wherein if T/T₀ equals zero, then said scatter fraction is generated in accordance with:

space="preserve" listing-type="equation">T.sub.0 *exp(-μ

.sub.TCT *2*R),

where R is a radius of rotation of the detector.

11. A method in accordance with claim 9 wherein at least one of said first constant A, said second constant B, and said third constant C is stored in a memory of said computer.

12. A method in accordance with claim 7 wherein the emission tomographic system is configured to reduce photon scatter utilizing a dual energy window data acquisition algorithm.

13. A method for performing scatter correction in radiographic imaging, said method comprising the steps of:

acquiring projection data for a primary window;

generating a three-dimensional scatter fraction; and

generating estimated scattered projections in the primary window using the three-dimensional scatter fraction.

14. A method in accordance with claim 13 further comprising the step of determining a transmission measurement.

15. A method in accordance with claim 14 wherein said transmission measurement is determined in accordance with:

space="preserve" listing-type="equation">T/T.sub.0 =exp(-μ

.sub.TCT D)

where;

T is a projection of a transmission scan;

T₀ is a projection of a blank scan;

D is a water-equivalent total body width that the transmission source passes through; and

μ

_TCT is a transmission source water attenuation coefficient.

16. A method in accordance with claim 15 wherein said scatter fraction is generated in accordance with:

space="preserve" listing-type="equation">k(i,j,l)=A-B* T/T.sub.0 !.sup.C/2

where;

k(i,j,l) is said scatter fraction;

A is a first constant;

B is a second constant; and

C is a third constant.

17. A method in accordance with claim 16 wherein if T/T₀ equals zero, then said scatter fraction is generated in accordance with:

space="preserve" listing-type="equation">T.sub.0 *exp(-μ

.sub.TCT *2*R),

where R is a radius of rotation of the detector.

18. A method in accordance with claim 16 wherein at least one of said first constant A, said second constant B, and said third constant C is stored in a memory of said computer.

19. A method in accordance with claim 13 further comprising the step of subtracting the estimated scattered projections from the projection data acquired for the primary window.

20. A method in accordance with claim 13 further comprising the step of performing iterative reconstructions utilizing the estimated scattered projections in the primary window.