X-ray scatter image reconstruction by balancing of discrepancies between detector responses, and apparatus therefor

• US 7,203,276 B2
• Filed: 07/27/2005
• Issued: 04/10/2007
• Est. Priority Date: 08/27/2004
• Status: Active Grant

• Current Assignee: University of New Brunswick
• Original Assignee: University of New Brunswick
• Inventors: Arsenault, Paul Jacob, Hussein, Esam
• Primary Examiner(s): THOMAS, COURTNEY D
• Application Number:
  US11/189,950
• Publication Number:
  US 20060043310A1
• Time in Patent Office: 622 Days
• Field of Search: 378/9, 378/50, 378/51, 378/54, 378/57, 378/62, 378/87, 250/393
• US Class Current: 378/87
• CPC Class Codes: G01N 23/20008 Constructional details of analysers, e.g. characterised by X-ray source, detector or optical system; Accessories therefor; Preparing specimens therefor monochromators for X- rays using crystals G21K1/06

Abstract

The present invention pertains to a method and an apparatus to generate a density image of an object using fan or cone beams of radiation and collimated detectors on one side of the object. The method consists of irradiating an object with a plurality of pairs of non-parallel radiation beams wherein the beams in each pair intersect a same segment along the axis of the detector. Compton-scattering radiation from the beams are then measured, and corrected attenuation coefficients along each beam are obtained. This latter step is effected by taking a first ratio of detector measurements for the beams in each pair; comparing the first ratio with a second ratio of corresponding calculated detector measurements and balancing discrepancies between the first and second ratios in a forward-inverse numerical analysis algorithm. Taking ratios of attenuation coefficients along related incident beams eliminates non-linearity problems whereby the aforesaid algorithm can be solved.

Claims

1. A method for imaging an object using x-ray, gamma rays or fast neutrons;

and a forward-inverse numerical analysis algorithm, comprising the steps of;

irradiating said object with radiation beams selected from fan beams or cone beams, along a plurality of beam paths;

obtaining a plurality of measurements of radiation scattered out of said object from said beam paths, wherein each of said measurements has a linear influence, a first exponential influence and a second exponential influence related to a density of said object;

selecting pairs from said measurements having a similar origin within said object and similar first or second exponential influences;

taking a ratio of each of said pairs of said measurements;

iteratively solving said forward-inverse numerical analysis algorithm using said ratios instead of said measurements; and

reconstructing an image of said object using corrected attenuation coefficients corresponding to said first or second exponential influences.

2. The method as claimed in claim 1, wherein said measurements of radiation scattered out of the object are taken with a detector having a sight axis, and each of said beam paths forms an acute angle with said sight axis.

3. The method as claimed in claim 1, wherein each of said pair of measurements have similar linear influences.

4. The method as claimed in claim 1, wherein said second exponential influences are related to incident radiation penetrating said object, and said first exponential influences are related to scattering radiation exiting said object.

5. The method as claimed in claim 4, wherein said step of selecting is effected with pair of said measurements having similar said first exponential influences and said step of reconstructing is effected with attenuation coefficients corresponding to said second exponential influences.

6. The method as claimed in claim 4, wherein said step of obtaining is effected using a collimated detector, and said similar origin is a segment along a sight axis of said collimated detector.

7. A method for imaging object using radiation selected from x-rays, gamma rays and fast neutrons;

a source of said radiation, a radiation detector having a sight axis extending through said object and a forward-inverse numerical analysis algorithm;

comprising the steps of;

irradiating said object with a first beam of radiation extending along a first angle relative to said sight axis, and intersecting said sight axis along a first scattering segment;

using said radiation detector, obtaining a first detector measurement of scattering radiation caused by said first beam;

irradiating said object with a second beam of radiation extending along a second angle relative to said sight axis and intersecting said sight axis along said first scattering segment, said second angle being different from said first angle;

using said radiation detector, obtaining a second detector measurement of scattering radiation caused by said second beam;

formulating a first ratio with said first and second detector measurements;

calculating first and second calculated detector responses;

making a second ratio with said calculated detector responses;

iteratively balancing relative errors between said first and second ratios and calculating corrected radiation attenuation coefficients of said object alone said first and second beams and said segment; and

relating said corrected radiation attenuation coefficients to visual indicators and constructing an attenuation-coefficient image of said object along said first and second beams and said segment.

8. The method as claimed in claim 7, wherein said first and second angles are mirror angles relative to said sight axis.

9. The method as claimed in claim 7 wherein said step of balancing relative errors includes the step of adding a regularization constraint thereto causing each of said voxels in a same region of said object to have similar properties.

10. The method as claimed in claim 7 further including the steps of dividing said object into voxels and repeating said steps of irradiating said object along first and second angles;

obtaining first and second detector measurements; and

formulating a ratio with said first and second detector measurements;

with second and subsequent beam pairs and second and subsequent scattering segments until each of said voxels is intersected at least four times by said beams of radiation.

11. The method as claimed in claim 10, further including the steps of identifying all voxels in a field of view of said detector during each of said detector measurements, and constructing an image of said object using attenuation coefficient through each of said voxels.

12. A method for imaging an object using radiation selected from x-rays or gamma rays, comprising the steps of:

a) providing a collimated radiation detector pointing at said object;

said collimated detector having a sight axis extending through said object;

b) providing means for irradiating said object with one or more radiation beams having a shape selected from a pencil beam, a fan beam or a cone beam;

c) irradiating said object with a first radiation beam having a first beam path intersecting said sight axis of said detector along a first segment of said sight axis;

said first segment being selected so that said detector detects a non-zero measurement;

d) using said detector, obtaining a first actual detector measurement of radiation scattered by said object from said first radiation beam;

e) irradiating said object with a second radiation beam having a second beam path different from said first beam path;

said second beam path intersecting said sight axis of said detector along said first segment of said sight axis;

f) using said detector, obtaining a second actual detector measurement of radiation scattered by said object from said second radiation beam;

g) associating said first and second actual detector measurements with said first and second beam paths respectively;

h) moving said object relative to said detector, or moving said beam paths relative to said sight axis of said detector or relative to each other such that an intersection thereof with said sight axis is a second common segment different from said first segment;

i) repeating said steps of irradiating, obtaining, associating and moving, and obtaining first and second actual detector measurements for each of a second and subsequent pairs of actual detector measurements of radiation scattered by said object corresponding respectively to first and second beam paths in each of a second and subsequent pairs of beam paths, until all voxels in said object have been intersected at least twice by said beam paths;

while obtaining each of said pairs of actual detector measurements;

defining first and second guessed radiation attenuation coefficients through said object along a corresponding pair of beam paths;

using said first and second guessed radiation attenuation coefficients, calculating first and second estimated detector measurements corresponding to said guessed radiation attenuation coefficients along said corresponding pair of beam paths;

comparing a first ratio of said first and second estimated detector measurements with a second ratio of said pair of actual detector measurements for said corresponding pair of beam paths or, comparing a first ratio of said first estimated detector measurement and said first actual detector measurement for said first beam path in said corresponding pair of beam paths with a second ratio of said second estimated detector measurement and said second actual detector measurement for said second beam path in said corresponding pair of beam paths, and obtaining a difference between said first and second ratios;

using said difference, correcting said first and second guessed radiation attenuation coefficients,

repeating said steps of calculating, comparing and correcting, and obtaining first and second corrected radiation attenuation coefficients along said first and second beam paths in said corresponding pair of beam paths respectively; and

j) using matrix manipulation and image reconstruction algorithms, relating said corrected radiation attenuation coefficients to density values and constructing a density image of said object, or relating said corrected radiation attenuation coefficients to visual indicators and constructing an attenuation-coefficient image of said object.

13. The method as claimed in claim 12, wherein said first and second beam paths form symmetrical angles relative to said detector axis.

14. The method as claimed in claim 12, wherein said detector measurements having a zero value are eliminated from said step of repeating said steps of calculating, comparing and correcting, and obtaining.

15. An apparatus for measuring radiation attenuation coefficients of material inside an object, comprising:

a first collimated detector having a first sight axis;

a first radiation source having a first beam axis, said first beam axis and said first sight axis extending along a same plane;

a second radiation source having a second beam axis, said second beam axis also extending along said same plane and said first and second radiation sources and said first collimated detector being on a same side of said object;

said beam axes being non-parallel with each other and being aligned to intersect said first sight axis along a same segment of said first sight axis.

16. The apparatus as claimed in claim 15, further comprising means for varying orientations of said first and second radiation beams within said same plane.

17. The apparatus as claimed in claim 16, further comprising a plurality of collimated detectors having a second and subsequent sight axes parallel with said first sight axis and extending along said same plane.

18. The apparatus as claimed in claim 17, further including means for moving said detectors, said radiation sources and said plane in unison.

19. The apparatus as claimed in claim 15, further including means for moving said radiation sources along said plane.