Collimator-free photon tomography

• US 5,818,050 A
• Filed: 04/07/1997
• Issued: 10/06/1998
• Est. Priority Date: 04/07/1997
• Status: Expired due to Fees

• Current Assignee: Associated Universities, Inc., Brookhaven Science Associates, LLC (Government of the United States of America), The Research Foundation for The State University of New York (State University of New York)
• Original Assignee: Brookhaven Science Associates, LLC (Government of the United States of America)
• Inventors: Dilmanian, F. Avraham, Barbour, Randall L.
• Primary Examiner(s): Westin, Edward P.
• Assistant Examiner(s): HANIG, RICHARD E
• Application Number:
  US08/838,493
• Time in Patent Office: 547 Days
• Field of Search: 250/363.09, 250/363.01, 250/363.02, 250/363.03, 250/363.04, 250/252.1 R
• US Class Current: 250/363.09
• CPC Class Codes: A61B 6/037 Emission tomography
• CPC Class Codes: G01T 1/1642 using a scintillation crystal and position sensing photodetector arrays, e.g. ANGER cameras

Abstract

A method of uncollimated single photon emission computed tomography includes administering a radioisotope to a patient for producing gamma ray photons from a source inside the patient. Emissivity of the photons is measured externally of the patient with an uncollimated gamma camera at a plurality of measurement positions surrounding the patient for obtaining corresponding energy spectrums thereat. Photon emissivity at the plurality of measurement positions is predicted using an initial prediction of an image of the source. The predicted and measured photon emissivities are compared to obtain differences therebetween. Prediction and comparison is iterated by updating the image prediction until the differences are below a threshold for obtaining a final prediction of the source image.

Claims

1. A method of photon computed tomography comprising:

administering a radioisotope to a patient (14) for producing gamma ray photons (16) from a source (14a) inside said patient;

measuring emissivity of said photons (16) externally of said patient (14) at a plurality of measurement positions (24) surrounding said patient (14) for obtaining corresponding energy spectrums (E);

predicting photon emissivity (Q) at said plurality of measurement positions (24) using an initial prediction of an image (I) of said source (14a);

comparing said predicted and measured photon emissivities (Q,M) to obtain differences (d) therebetween; and

iterating said predicting and comparing steps by updating said image prediction (I) until said differences (d) are below a threshold for obtaining a final prediction of said image (I_n).

2. A method according to claim 1 further comprising measuring said photon emissivity uncollimated.

3. A method according to claim 2 wherein:

said energy spectrum (E) includes a plurality of incremental energy levels (1,2,3, . . . ,N); and

said measured emissivity (M) includes a corresponding intensity for each of said energy levels at each of said measurement positions.

4. A method according to claim 3 wherein said measurement positions each include a plurality of pixels (26) arranged in a two-dimensional spatial array, with each pixel (26) having a corresponding one of said energy spectrums (E).

5. A method according to claim 4 wherein said predicting step comprises applying a reference imaging operator (P) to said image prediction (I) to obtain said predicted photon emissivity (Q) at said plurality of pixels (26) for said plurality of measurement positions.

6. A method according to claim 5 wherein said reference imaging operator (P) is obtained by:

providing a reference source (30) of said radioisotope inside a phantom (28) representative of said patient (14) for producing gamma ray photons (16) therein;

indexing said reference source (30) incrementally inside said phantom (28) at a plurality of positions therein to define a respective plurality of volume element voxels (34) collectively defining said phantom (28); and

measuring emissivity of said photons (16) emitted from said phantom (28) externally thereof at each of said measurement positions for obtaining corresponding energy spectrums (E) thereat for each of said incrementally scanned voxels (34), with said imaging operator (P) being said measured voxel photon emissivity (E) due to each of said irradiated voxels (34) independently.

7. A method according to claim 6 wherein said phantom (28) comprises an anthropomorphic hollow vessel filled with liquid (28b).

8. A method according to claim 6 wherein said voxels (34) define said predicted image (I).

9. A method according to claim 8 further comprising:

perturbating said emissivity differences (d) in each of said iterating steps to obtain respective changes (D) in said predicted image (I) at each of said voxels (34); and

adjusting said predicted image (I) at each of said iterating steps using said image changes (D).

10. A method according to claim 9 wherein said perturbating step comprises Green'"'"'s function.

11. A photon computed tomography system (10) for imaging a source (14a) in a patient (14) producing gamma ray photons (16) comprising:

means (12) for measuring emissivity of said photons (16) externally of said patient (14) at a plurality of measurement positions (24) surrounding said patient (14) for obtaining corresponding energy spectrums (E);

means (20b) for predicting photon emissivity (Q) at said plurality of measurement positions (24) using an initial prediction of an image (I) of said source (14a);

means (20b) for comparing said predicted and measured photon emissivities (Q,M) to obtain differences (d) therebetween; and

means (20b) for iterating said predicting and comparing steps by updating said image prediction (I) until said differences (d) are below a threshold for obtaining a final prediction of said image (I_n).

12. A system according to claim 11 wherein said measuring means comprise an uncollimated gamma camera (12).

13. A system according to claim 12 wherein:

said energy spectrum (E) includes a plurality of incremental energy levels (1,2,3, . . . ,N); and

said measured emissivity (M) includes a corresponding intensity for each of said energy levels at each of said measurement positions.

14. A system according to claim 13 wherein said camera (12) includes a plurality of pixels (26) arranged in a two-dimensional spatial array, with each pixel (26) having a corresponding one of said energy spectrums (E).

15. A system according to claim 14 wherein said means for predicting, comparing, and iterating comprise a computer (20b) configured for performing said steps.

16. A system according to claim 15 wherein said computer (20b) is configured for applying a reference imaging operator (P) to said image prediction (I) to obtain said predicted photon emissivity (Q) at said plurality of pixels (26) for said plurality of measurement positions.

17. A system according to claim 16 further comprising:

a reference source (30) of said radioisotope inside a phantom (28) representative of said patient (14) for producing gamma ray photons (16) therein;

means (32) for indexing said reference source (30) incrementally inside said phantom (28) at a plurality of positions therein to define a respective plurality of volume element voxels (34) collectively defining said phantom (28); and

said camera (12) being effective for measuring emissivity of said photons (16) emitted from said phantom (28) externally thereof at each of said measurement positions for obtaining corresponding energy spectrums (E) thereat for each of said incrementally scanned voxels (34), with said imaging operator (P) being said measured voxel photon emissivity (E) due to each of said irradiated voxels (34) independently.

18. A system according to claim 17 wherein said phantom (28) comprises an anthropomorphic hollow vessel filled with liquid (28b).

19. A system according to claim 18 wherein said computer (20b) is further configured for:

perturbating said emissivity differences (d) in each of said iterating steps to obtain respective changes (D) in said predicted image (I) at each of voxels (34); and

adjusting said predicted image (I) at each of said iterating steps using said image changes (D).