Search Images Maps Play YouTube News Gmail Drive More »
Sign in
Screen reader users: click this link for accessible mode. Accessible mode has the same essential features but works better with your reader.

Patents

  1. Advanced Patent Search
Publication numberUSRE37899 E1
Publication typeGrant
Application numberUS 09/456,636
Publication dateNov 5, 2002
Filing dateDec 8, 1999
Priority dateMar 11, 1996
Fee statusPaid
Also published asUS5696806
Publication number09456636, 456636, US RE37899 E1, US RE37899E1, US-E1-RE37899, USRE37899 E1, USRE37899E1
InventorsLee Grodzins, Charles G. Parsons
Original AssigneeAmerican Science And Engineering, Inc.
Export CitationBiBTeX, EndNote, RefMan
External Links: USPTO, USPTO Assignment, Espacenet
Tomographic method of x-ray imaging
US RE37899 E1
Abstract
A method is disclosed for obtaining the density distributions of three-dimensional elements that compose objects or groups of objects, by examining the objects with beams of x-rays or gamma radiation that are transmitted through the object in a plurality of approximately parallel paths and measuring the intensity of the radiation, scattered approximately perpendicular to the parallel paths, in arrays of detectors around the object. The energy of the x-rays or gamma rays is such that dominant interaction in the object is Compton scattering. The density of each element is determined from the totality of measurements by standard mathematical tomographic or relaxation techniques of data manipulation.
Images(5)
Previous page
Next page
Claims(36)
We claim:
1. A method for determining a three-dimensional density, distribution among volume elements in a volume containing at least one material, the method comprising:
a. producing a beam of energetic photons for penetrating the volume;
b. scanning the beam sequentially across incremental positions of the volume in a plurality, of paths having substantially parallel directions separated by incremental steps;
c. detecting scattered photons of substantially all energies scattered by the material in said volume with a detector having a spatial resolution in a direction substantially parallel to the paths of the beam;
d. measuring the intensity of scattered photons scattered approximately perpendicular to each substantially parallel direction of the beam at each incremental position of the beam to derive an independent measurement of intensity, of scattered photons scattered from the volume elements in the volume;
e. identifying a volume element as the approximate origin of scattering along the path of the beam giving rise to the measurement of intensity of scattered photons; and
f. calculating an independent density for each volume element of the material in the volume.
2. A method according to claim 1, in which the size of each volume element in the volume is determined by the diameter of the beam that traverses said volume element, the step size of the incremental steps of said beam, and the spatial resolution of the detector substantially parallel to the beam path.
3. A method according to claim 1, in which the number of independent measurements of the intensity of scattered photons scattered from the volume elements in said volume is at least approximately equal to the number of volume elements.
4. A method according to claim 1, in which the intensity, of each beam transmitted through the volume is measured together with the intensity of said scattered photons scattered by the material.
5. A method according to claim 1, in which the energies of the energetic photons in each beam are in an energy range where the interaction of the energetic photons with the material is dominated by the probability for Compton scattering.
6. A method according to claim 1, in which the step of calculating the density of each volume element of the material in the volume includes a mathematical relaxation procedure in which a first trial function in the mathematical relaxation procedure uses the independent measurements of the intensity, of scattered photons scattered from the volume elements in the volume.
7. A method according to claim 1, in which the density distribution among the volume elements in the volume is determined by a computerized tomographic reconstruction procedure having a number of independent measurements at least equal to the number of independent densities calculated in the step of calculating.
8. A method according to claim 1, in which the step of scanning the beam sequentially across incremental position of the volume includes scanning a first beam and a second beam, the first beam having a first mean energy high enough so that the Compton interaction dominates the interaction of the first beam with the material, and the second beam having a mean energy low enough so that the photoelectric interaction makes a significant contribution to the interaction of the second beam with at least some material in the volume, such that the intensity of scattered radiation from the first and second beams are used to determine the effective atomic number of the at least one material in the volume by a mathematical reconstruction technique.
9. A method according to claim 8, wherein the mathematical reconstruction technique includes relaxation methods.
10. A method according to claim 8, wherein the mathematical reconstruction technique includes methods of computerized tomographic reconstruction.
11. A device for determining densities in volume elements in a material present in an assembly of objects, the device comprising:
a. a source for producing a beam of energetic photons having a direction for penetrating the material;
b. an arrangement for scanning the beam of energetic photons in a sequence of sequential beams across the assembly in a manner such that successive directions of the sequential beams are substantially parallel to each other, the sequential beams passing through every volume of the assembly;
c. a detector disposed substantially parallel to the direction of the beam of energetic photons for providing measurements of the intensity of scattered photons of substantially all energies scattered approximately perpendicular to the direction of the beam at each incremental position of the beam; and
d. a computer for determining the densities in the volume elements of the material in the assembly from a totality of measurements of the intensity of scattered photons scattered approximately perpendicular to the direction of the beam at each incremental position of the beam by a mathematical reconstruction technique.
12. A device according to claim 11, further comprising a transmission detector for detecting energetic photons transmitted through the assembly of objects.
13. A device according to claim 11, wherein the beam of energetic photons comprises energetic photons in the range where the interaction of the energetic photons with the material in the assembly is dominated by the probability for Compton scattering.
14. A device according to claim 11, wherein the detector is segmented in a direction having a component substantially parallel to the beam of energetic photons.
15. A device according to claim 11, wherein the detector is collimated with respect to directions substantially perpendicular to the beam of energetic photons.
16. A device according to claim 11, wherein the mathematical reconstruction technique includes relaxation methods.
17. A device according to claim 11, wherein the mathematical reconstruction technique includes methods of computerized tomographic reconstruction.
18. A tomography system for analyzing a material concealed within an enveloping surface, the system comprising:
A. at least one source of penetrating radiation for emitting a beam along a beam axis having an orientation disposed with respect to the enveloping surface;
B. a scanner for varying the orientation of the beam axis with respect to the material in incremental steps; and
C. at least one array of segmented detectors disposed along a detector axis disposed substantially parallel to the beam axis for detecting scattered radiation of substantially all energies of scattered radiation and producing signals corresponding at least to the scattered radiation.
19. The tomography system in accordance with claim 18, additionally comprising at least one transmission detector disposed along the beam axis for measuring penetrating radiation transmitted through the material.
20. The tomography system in accordance with claim 18, additionally comprising a conveyor for transporting the material.
21. The tomography system in accordance with claim 18, wherein the scanner arrangement includes a scanner for raster-scanning the beam axis in a plane transverse to the beam axis.
22. The tomography system in accordance with claim 18, further comprising a plurality of collimators disposed in directions substantially perpendicular to the beam axis for limiting the field of view of each segmented detector.
23. A method for analyzing material concealed within an enveloping surface, the method comprising:
a. illuminating the enveloping surface with penetrating radiation propagating substantially along a beam axis, the penetrating radiation characterized by a first incident energy;
b. measuring a profile of penetrating radiation characterized by a first incident energy that is scattered by the concealed material;
c. illuminating the enveloping surface with penetrating radiation propagating substantially along the beam axis, the penetrating radiation characterized by a second incident energy;
d. measuring a profile of penetrating radiation characterized by a first incident energy that is scattered by the concealed material; and
e. determining at least one of the density and atomic number associated with each of a plurality of voxels based on the profiles of penetrating radiation characterized by the first and second incident energies and scattered by the concealed material.
24. A method for analyzing a material concealed within an enveloping surface, the method comprising:
A. producing a beam of energetic photons for penetrating the volume;
B. scanning the beam across the material incrementally in a plurality of beam directions;
C. detecting scattered photons of substantially all energies scattered by the material with a detector array having a spatial resolution in a direction substantially parallel to the beam directions;
D. measuring the intensity of scattered photons scattered approximately perpendicular to each direction of the beam; and
E. characterizing the material based on the intensity of scattered photons scattered approximately perpendicular to each incremental direction of the beam.
25. The method according to claim 24, further including measuring penetrating radiation transmitted through the material.
26. The method according to claim 24, further including deriving an attenuation per unit length characteristic of the material.
27. The method according to claim 24, further comprising the step of raster-scanning the beam axis in a plane transverse to the beam axis.
28. The method according to claim 24, further comprising the step of conveying the enveloping surface in a direction perpendicular to the beam axis.
29. An apparatus for generating sequential beams of penetrating electromagnetic radiation comprising:
a. a source for producing a beam of charged particles;
b. a target having a surface which receives the beam of charged particles and emits electromagnetic waves in response thereto;
c. a beam director that directs the beam of charged particles to the plurality of specified locations on the target; and
d. a collimator, the collimator having an array of transmitting regions and being disposed proximal to the target such that electromagnetic waves emitted from the target pass through the collimator and emerge from the array of transmitting regions in a series of parallel beams as the beam of charged particles is directed at a plurality of specified locations on the target.
30. An apparatus according to claim 29, wherein the collimator includes an array of apertures.
31. An apparatus according to claim 29, wherein the collimator includes a plurality of apertures in an absorbing matrix.
32. An apparatus according to claim 31, wherein the absorbing matrix is lead.
33. An apparatus according to claim 29, wherein the penetrating electromagnetic radiation includes x-rays.
34. The apparatus according claim 29, wherein the source for producing a beam of charged particles includes an cathode ray tube.
35. A method for generating sequential beams of penetrating electromagnetic radiation comprising:
a. providing a source for producing charged particles;
b. directing the particles to a plurality of specified locations on a target having a surface to produce radiation; and
c. collimating the radiation.
36. A method according to claim 35, wherein the step of collimating the radiation includes directing the radiation through apertures in an absorbing matrix.
Description
BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to the examination of a body by means of x-radiation or gamma radiation.

2. Description of the Prior Art

The creation of images of x-ray attention coefficients through objects began with the discovery of x-rays by Roentgen in 1895 and has developed continuously. A major advance was the invention in the early 1970's, by GN Hounsfielch G. N. Hounsfield of computerized axial tomography, CAT. (U.S. Pat. Nos. 3,778,614, 4,035,647). Hounsfield showed that the linear attenuation coefficients of individual elements in an object could be reconstructed from the measurements of the intensities of x-ray beams that pass in a plurality of independent paths through the object. Since Hounsfield's invention, the technology of CAT scanning has developed in a number of ways, but all x-ray methods deduce the density distributions from measurements of the radiation transmitted through the object.

SUMMARY OF THE INVENTION

This invention is a method for determining the densities of objects or the distribution of densities interior to an object by measuring the intensities of x-rays that are Compton scattered from the volume elements (voxels) of the object. The term x-rays is used throughout the descriptions since it is anticipated that most applications will use an x-ray beam generated by energetic electrons inn an x-ray tube, but it should be evident that all types of energetic photons can be used, including monoenergetic gamma rays, that satisfy the criterion that the energy of the photons are such that Compton scattering dominates the interactions of the photons in the object. The use of scattered radiation to determine densities is fundamentally different in theoretical underpinning, in methodology and in implementing apparatus from the standard methods that use transmitted x-rays to determine densities of interior voxels.

The incident x-rays are rastered across at least one face of an object. The x-rays that are Compton scattered approximately perpendicular to the beam directions are detected in arrays of collimated detectors each of which is sensitive to radiation scattered from a specific portion of the incident x-ray beam.

The distinctive features of the preferred embodiment of this invention are: 1) The energy of the x-rays is high enough so that the interactions in the object are dominated by the Compton effect. In particular, the energy is high enough so that the photoelectric interaction makes a minor contribution in the analysis but it is not so high that pair production is significant. 2) The incident x-rays are collimated into a beam that is scanned through the object in a series of contiguous, approximately parallel paths; the scanning may be accomplished by moving the beam or the container or a combination of both so that the incident x-ray beam passes through every voxel that the detected scattered radiation passes through. 3) The scattered x-rays are detected by arrays of counters that distinguish x-rays that are scattered approximately perpendicular to the incident radiation. The detectors must be capable of sensing the direction of the scattered radiation. Methods for sensing the direction of an incident x-ray are well known; gamma cameras, for example, do so with collimators and position-sensitive detectors. 4) The volume element resolved by this invention is determined by the cross sectional area of the incident x-ray beam times the spatial resolution along the beam path of the origins of the scattered x-rays. The total number of independent measurements is at least equal to the total volume being examined divided by the volume element of spatial resolution 5. The densities of the voxels is rapidly and accurately determined from the totality of measurements by standard mathematical relaxation methods, without the need for transformation into frequency space or the use of back projection, though both of those techniques can be used.

This invention, which we will refer to as Compton Scatter Tomography or CST, is a new modality for tomography, quite distinct from the conventional method of computerized axial tomography, CAT, in which the linear attenuation coefficients in voxel elements in an object are determined from transmission measurements, or Single Photon Emission Computer Tomography, SPECT, in which the directions of gamma rays emitted from a radioactive source distributed in a body is used to measure the distribution of radioactivity. The invention is described in its application for the inspection of containers for contraband such as explosives. It should be appreciated, however, that this invention may be useful for a broad range applications in which a non-destructive method is needed to determine the density distributions of objects.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1. A schematic drawing of the main elements of a preferred embodiment of the invention. A beam of x or γ-rays 2 directed along the Y axis is stepped in a series of approximately parallel paths in the YZ plane so as to intersect every voxel in the plane 8 of the container 3. The examined container moves in the X direction. The x-rays 4 scattered through approximately 90° by elements in the container are counted by an array of collimated detectors 5 and 6, above and below, respectively, the container. A detector 7 of the transmitted x-rays is also shown.

FIG. 2. A cross section of the container showing one plane of the container, such as plane 8. The area of the plane is divided into 25 imaginary voxels. The incident x-ray beam 12 passes along the fourth row and is shown scattered in voxel 43 4,3into the upper detector 18 and the lower detector 20. The transmitted beam is shown stalking striking the detector 12.

FIG. 3. A cross section view of the main elements of the invention showing an x-ray machine 31 that generates a raster scanned beam of electrons that strike the anode 33. The electron beam 32 is shown striking the top of the anode 33; the electron beam 34 is shown striking near the mid-point of the anode 33. The electron beam 34 generates x-rays that are collimated into a beam 35 by the collimator 36. The x-ray beam 34 passes through the container 41. Some of the x-ray beam is scattered through approximately 90° into arrays of collimated detectors 42 and 43 that are above and below the container, respectively. The x-ray beam that transmits the container 41 is detected in a segmented detector 38.

FIG. 4. A plot of the ratio of the Compton mass attenuation factor to the total mass attenuation factor as a function of x-ray energy, for oxygen, silicon and iron.

EMBODIMENT EMPLOYING MONOENERGETIC GAMMA RAYS

The invention is first described using a monoenergic gamma ray source, in particular the 662 keV gamma ray from the decay of 137Cs. Other radioactive sources, such as the 356 keV gamma ray from the decay of 133Ba, or the 1117 keV and 1332 keV gamma rays from 60Co, or monochromatic x-ray sources might be appropriate for specific applications but 137Cs is an especially appropriate choice because of its long 30 year half-life, low cost, high specific activity and simplicity of its radiation spectrum.

FIG. 1 shows the essential features of the invention as it might be applied to determining the density distributions of materials in luggage. The 662 keV gamma rays from the radioactive source 1 of 137Cs, are collimated into a beam 2, which is aimed in the Y direction into the luggage 3. The luggage 3 is conveyed in the X direction, i.e., perpendicular to the direction of the beam 2. The intensity of gamma rays 4 that are Compton scattered through approximately 90° in the ±Z directions are measured in top and bottom detectors 5 and 6. The collimated beam 2 is incrementally moved in the Z direction to an adjacent row of voxels in the slice 8 and the measurements repeated until the x-ray beam has interacted with every voxel in the full YZ slice 8 of the luggage 3. The luggage is then moved incrementally in the X direction and the adjacent slice 9 is investigated. In this way, the entire container is examined. It will be appreciated that any combination of relative motions of the gamma ray beam with respect to the luggage is acceptable, including stationary luggage with all of the relative motion supplied by a raster scanned beam and the converse, a stationary beam with all of the relative motion supplied by a luggage conveyance. It will also be appreciated that the relative motions may be incremental or continuous depending on the application. The central requirement is that every voxel in the volume being interrogated must be traversed at least once by both the incident and scattered radiation. The detectors of scattered radiation 5 and 6 may each be single large volume detectors that have the energy and angle dispersive power to determine both the energy of the Compton scattered radiation and the direction from which it came. Such detectors are being developed with these capabilities, but at this time a more cost-effective solution is to use segment detectors 5 and 6 into arrays of collimated detectors each of which is sensitive to radiation from a particular voxel along the beam 2.

FIG. 2 shows one slice 8 of the luggage imagined to be divided into 5 rows and 5 columns to make a total of 25 voxels. The voxels are numbered sequentially, 11, 12, 13, 14, 15, 21, 22, 1,1; 1,2; 1,3; 1,4; 1,5; 2,1; 2,2, etc. The 662 keV beam of gamma rays 11 is shown passing along the fourth row, through voxels, 41, 42, 43, 44 and 45 4,1; 4,2; 4,3; 4,4; and 4,5to the transmission detector 21. The gamma rays that are scattered in the voxels of row 4 are counted in the top detector array 13 and the bottom detector array 14 which have appropriate collimators 15 and 16 to ensure that only scatterings through ˜90° are detected. The voxels are assumed to be cubic; in practice, the size of the voxels will depend on the beam cross section, beam divergence and the spatial resolution of the detector arrays. FIG. 2 shows an example of Compton scattered radiation 17 being scattered in the Z direction from pixel 43 4,3into the single detector 18 and another scattering 19 in the −Z-direction into detector 20.

The scattering of the 662 keV gamma ray through 90° results in 288 keV radiation. The intensity of the 288 keV signal, I17, in detector 18 is given by: I 17 ( 288 ) = I o ( 662 ) e - ( λ 41 s 41 + λ 42 s 42 ) 662 [ λ 43 Comp ( 662 ) Ω Ω ( π 2 ) r 43 ] e - ( λ 33 r 33 + λ 23 r 23 + λ 13 r 13 ) 288 1 )

where Io(662) is the intensity of the incident beam 11, the λ values in the first exponential are the total linear attenuation coefficients for 662 keV radiation, the λ, values in the second exponential term are the total linear attenuation coefficients for the scattered radiation of −288 keV, the t values are the linear dimensions of the pixels, and the square bracket term is the probability for Compton scattering in pixel 43. The linear attenuation coefficients are defined in terms of the cross section σ, the atomic weight A, Avogadro's number No, and the density ρ; λ = σ N o ρ A 2 )

An analogous equation to 1 can be written for the intensity scattered into detector 20; only the last exponential term and the solid angle factor changes. There will be a total of 10 equations describing the scattering from material along row 4 into the segmented detectors in the top 13 and bottom 14 arrays. Each row that the 662 beam traverses will produce 10 more independent equations. A total of 50 equations will be generated in a full scan of the slice 8. An additional 5 equations will be generated by the intensities in the transmission detector 12 21; the equation for the ray 12-22 shown in FIG. 2 is given by,

I22(662)=I0(662)e −(80 43 t 41 42 t 42 43 t 43 44 t 44 )I22( 662 )=I 0( 662 )e −(λ 41 t 41 42 t 42 43 t 43 44 t 44 )   3)

The transmission intensities exampled by Equations 3 are not necessary for solving Equations 1 for the densities in each of the volumes of the container 12, but they give important additional information that can speed up and make more secure the analytic procedures.

In the preferred embodiment using 137Cs, the interactions in the container are dominated by the Compton effect. FIG. 4 shows that the Compton effect in iron, the heaviest of the materials found in luggage in substantial quantities, accounts for 98% of the interactions for the incoming radiation of 662 keV and 88% of the interactions for the 288 keV scattered radiation. Equations 1, 2 and 3 then simplify enormously since,

λ=λComptonphotoeffectRayleighpair production≅λCompton   4

and the differential Compton attenuation λ 43 Comp ( 662 ) Ω

for 90° scattering of 662 keV gamma rays is simply related to the total linear attenuation, λ 43 Comp ( 662 ) Ω = 0.04 λ Comp ( 662 ) . 5 )

The consequence of these simplifications, all of which follow from the use of sufficiently high energy photons, is a set of 50 scattering equations and 5 transmission equations that contain just 25 unknowns, the 25 Compton linear attenuation coefficients. Equation 1, for example, becomes, I 17 ( 288 ) = I o ( 662 ) e - ( λ 41 e + λ 42 e ) 662 [ .04 λ e ( 662 ) Ω ( π 2 ) t 43 ] e - ( λ 33 e r 33 + λ 23 e r 23 + λ 13 e r 13 ) 6 )

where all of the unknown quantities in the inspection of a container are total linear Compton attenuation coefficients which depend primarily on the densities in the voxels.

The 50 scattering equations in this example can be rapidly solved by relaxation techniques; it is not necessary to use matrix inversions or convert to frequency space. To emphasize the simplicity we note that the equations describing the interaction of the beam with the top row of voxels of FIG. 2 (voxels 11 1,1through 15 1,5) result immediately in the densities of each of the voxels since the scattering from voxel 11 is described by an equation with only 1 unknown, the density of the voxel. The scattering from voxel 12 1,2is described by an equation with only 2 unknowns, one of which has been determined from the scattering from voxel 11 1,1, and so forth. Thus a series of simple iterations results in the linear Compton attenuation values for all of the voxels. In practice, one would use the values of the 50 intensifies intensities of scattered radiation and 5 intensities of transmitted radiation to determine a first-order map of the λ values in the 50 voxels and then use mathematical relaxation techniques to obtain a best set of density values. It should be noted that the method results in at least twice as many independent equations as the minimum required for a full tomographic analysis. These extra equations, plus the transmission equations 3, can be used to make the corrections to the scattering equations, such as Equation 6, to take into account small contributions of photo-electric absorption or coherent scattering that are present when the voxels contain higher Z material.

The linear Compton attenuation coefficients are directly proportional to the electron densities in the voxels; i.e., Equation 2 simplifies to Equation 8. 7,

where σe, the Compton scattering per electron, is a constant, and Z is the number of electrons per atom. The electron densities are, in turn, very closely related to the matter densities since, for most materials in luggage, Z/A≅0.5. (The avenge value of Z/A for plastics, explosives and other light materials is a few percent greater than 0.5; Z/A for heavier materials such as iron are a few percent less.)

EMBODIMENT EMPLOYING POLYENERGETIC X-RAYS

It is anticipated that the invention will usually be implemented using a spectrum of x-rays generated by an electron beam, The source of the x-rays could be a conventional x-ray tube with a fixed electron beam striking a fixed or rotating anode. Our preferred embodiment uses a raster scanned electron beam shown schematically in FIG. 3. The x-rays are generated in a x-ray tube 31, which produces a raster-scanned beam of x-rays by scanning the electron beam 32 and placing an appropriate collimator 36 in front of the anode 33. The x-ray tube 31 is similar to a conventional cathode ray tube with an appropriate heavy element anode 33 replacing the traditional phosphor screen. As the electron beam sweeps across the anode 33, x-rays are generated that pass through successive parallel holes in the collimator 36. The result is a rastering of approximately parallel beams of x-rays through the container 41. The anode potential determines the maximum energy of the x-ray beams. An absorber 37 eliminates the softer components of the x-ray beam and determines the effective lower energy of the x-rays that interact in the container.

It should be emphasized that the choice of anode voltage and x-ray strength depends on the application. Large containers might warrant anode voltages as high as 2 MeV (pair production is still negligible at this energy) in order to produce sufficiently penetrating x-rays, while small containers with primarily low Z components might be studied effectively with x-ray energies below 150 keV.

A practical choice of parameters for scanning airport luggage would be an anode potential of 450 keV and an electron current of 4 milliamps. The length of the anode 33 might be 20″, i.e., about the height of the anode of a 30″ TV tube; a power density of 2 kilowatts is easily handled by modest cooling of the large-area. The collimator 36 might be a set of parallel holes in a 4″ thick lead block (attenuation by the lead>1010). The holes should be appropriately designed to minimize internal scattering in the collimator. The absorber 37 might be 1 mm of tungsten that would reduce the 300 keV x-rays by a factor of ˜2 while killing 100 keV components by factors of 104.

Above and below the container are detectors or detector arrays 42 and 43, respectively that measures the scattered x-rays 46 as a function of position of scattering along the beam direction. To do this, we propose to use collimating slits 44 and 45 such as the Soller plates used extensively in x-ray diffraction. These slits restrict the direction of x-rays seen by the detector; their function is similar to the collimators used in Single Photon Emission Tomography (SPECT) in which the origins of the emission of gamma rays from radioactive sources is determined by the SPECT detector. Many options are available for the detectors including the hodoscopes of NaI(TI), BGO and CdZnTe now used for SPECT and Positron Emission Tomogaphy.

Equations 1 and 3 must now be written in terms of weighted integrals over the energy spectra. Exact expressions can be taken into account in the analysis, though we anticipate that in most practical cases it will be sufficient to use appropriate averages of the incident and scattered energies as well as the differential and integral linear attenuation coefficients in Equations 1 and 3, since the Compton cross sections vary slowly with energy, atomic number, and scattering angle around 90°. Specifically, for elements from carbon to iron, the total Compton cross section varies by only 25% from 150 keV to 450 keV; for a given x-ray energy, it varies by only 10%. Moreover, the differential Compton cross section is almost independent of angle from 80° to 110°. It should also be noted that beam hardening—the changing energy spectrum in the container due to absorption—will not be significant when the invention is applied to airline baggage since the high energies of the incident beam are not much attenuated traversing an airline suitcase.

We have carded out computer simulation studies, assuming a rastered x-ray beam with dimensions 5 mm×5 mm generated by a 2 kilowatt, 450 keV electron beam. We estimate that efficient side scattered detectors can determine the origin of the scattered x-rays to within 2 cm along the beam path. (Note that SPECT hodoscopes of 150 keV radiation have spatial resolutions of approximately 0.5 cm.) Each voxel thus has a volume of 0.5 cc so that 100 grams of explosives would occupy about 300 voxels. The simulation studies show that the interrogation of a piece of luggage, 1 meter×60 cm×20 cm can be carried out in 6 seconds, resulting in the determination of the linear attention coefficients of each voxel in the luggage to an accuracy of 30%. The mean values of the densities of any contiguous 300 voxels (˜100 g of explosives) would then be known to an accuracy of 2%. The simulation studies show that CST should have a minimum detection limit below 100 g of explosives.

It should also be noted that the CST method of tomographic analysis is very effective for finding sheet bombs, one of the most difficult of the explosive configurations to investigate by x-ray means.

A logical extension of the invention is to make scattering measurements at two incident energies, one at the preferred high energy where the Compton effect is dominant and the other at a lower energy where the photo-electric effect makes a substantial contribution to the interactions of the x-rays in those voxels with high Z component. This so-called dual-energy method is well know for transmission tomography where it is used to determine the effective atomic number of the elements in the voxels and we anticipate that the dual energy method could have applications in which the measurement of the effective atomic number of the voxels as well as the density is important. Referring to FIG. 4, one might choose a high x-ray energy greater than 300 keV where the Compton effect is more than 90% of the total interaction in iron, and a low energy x-ray in the 100 keV range where the photo-electric effect makes up about 80% of the cross section in iron. Dual energy, however, is not so easily applied to CST where the scattered radiation is substantially lower than the incident energy; e.g., the energy of Compton scattered 300 keV x-ray is 188 keV, and the energy of the Compton scattered 100 keV x-ray is 83 keV. The change in the energy in the former case makes only a minor, and easily accounted for, complication to the analysis. The change in the latter case is not so easily taken into account since small changes in the scattered energy make large changes in the attenuation coefficients. The effects of beam hardening—that is, the increase in the average energy of the beam as the lower energy components become absorbed in the container—are very difficult to take into account. Nevertheless, there may be applications for this extension, particularly when only a relatively few voxels contain high Z material.

The invention stresses that the incident beams should be rastered across the container in approximately parallel paths and that the detected radiation should be limited to those x-rays that are scattered approximately perpendicular to the incident beam direction. The allowable deviations from these conditions depend on the applications. For all applications we expect that the deviations can be at least ±20°, since the cos 20° deviates by only 6% from unity. For some applications, especially those in which the high density regions make up a small portion of the container, the deviations from ideal could be considerably larger. Simulation studies indicate that the CST method is robust with respect to deviations from parallel beams and 90° scattering but that the closer the rastered beams are to being parallel, and the closer the scatter angle of detected rays are to 90°, the simpler and more accurate will be the analytic tomographic procedures for determining the densities.

Patent Citations
Cited PatentFiling datePublication dateApplicantTitle
US3955089Oct 21, 1974May 4, 1976Varian AssociatesAutomatic steering of a high velocity beam of charged particles
US4002917Aug 27, 1975Jan 11, 1977Emi LimitedSources of X-radiation
US4144457Apr 5, 1976Mar 13, 1979Albert Richard DTomographic X-ray scanning system
US4149076Jun 24, 1977Apr 10, 1979Albert Richard DMethod and apparatus producing plural images of different contrast range by X-ray scanning
US4194123May 12, 1978Mar 18, 1980Rockwell International CorporationLithographic apparatus
US4196351Jun 3, 1977Apr 1, 1980Albert Richard DavidScanning radiographic apparatus
US4535243Mar 17, 1983Aug 13, 1985Imatron AssociatesX-ray detector for high speed X-ray scanning system
US4598415Aug 15, 1983Jul 1, 1986Imaging Sciences Associates Limited PartnershipMethod and apparatus for producing X-rays
US4672615Jul 19, 1983Jun 9, 1987Unisearch LimitedIon and electron beam steering and focussing system
US4694457Dec 26, 1985Sep 15, 1987Unisearch LimitedMethods of steering and focusing ion and electron beams
US4730350 *Apr 21, 1986Mar 8, 1988Albert Richard DMethod and apparatus for scanning X-ray tomography
US4864142Jan 11, 1988Sep 5, 1989Penetron, Inc.Method and apparatus for the noninvasive interrogation of objects
US5022062Sep 13, 1989Jun 4, 1991American Science And Engineering, Inc.Automatic threat detection based on illumination by penetrating radiant energy using histogram processing
US5097492Feb 12, 1990Mar 17, 1992Four Pi Systems CorporationSolder connection
US5153900Sep 5, 1990Oct 6, 1992Photoelectron CorporationMiniaturized low power x-ray source
US5179581 *Jun 4, 1991Jan 12, 1993American Science And Engineering, Inc.Automatic threat detection based on illumination by penetrating radiant energy
US5247561Jan 2, 1991Sep 21, 1993Kotowski Andreas FLuggage inspection device
US5420905 *Oct 21, 1993May 30, 1995Massachusetts Institute Of TechnologyDetection of explosives and other materials using resonance fluorescence, resonance absorption, and other electromagnetic processes with bremsstrahlung radiation
US5430787 *Dec 3, 1992Jul 4, 1995The United States Of America As Represented By The Secretary Of CommerceCompton scattering tomography
US5442678Jan 21, 1994Aug 15, 1995Photoelectron CorporationX-ray source with improved beam steering
US5504796Nov 30, 1994Apr 2, 1996Da Silveira; Enio F.Method and apparatus for producing x-rays
US5548630Aug 31, 1995Aug 20, 1996Siemens AktiengesellschaftX-ray radiator having an electron source for sending a beam of electrons along an elongated anode
US5682412Sep 20, 1996Oct 28, 1997Cardiac Mariners, IncorporatedX-ray source
US5712889Apr 11, 1997Jan 27, 1998Lanzara; GiovanniScanned volume CT scanner
US5805662 *Dec 7, 1995Sep 8, 1998Quanta Vision, Inc.Using deflected penetrating radiation to image an object's internal structure
US5841831Apr 28, 1997Nov 24, 1998Siemens AktiengesellschaftX-ray computed tomography apparatus
US5930326 *Jul 8, 1997Jul 27, 1999American Science And Engineering, Inc.For analyzing a material concealed within an enveloping surface
US6111974 *Feb 11, 1998Aug 29, 2000Analogic CorporationApparatus and method for detecting sheet objects in computed tomography data
DE19710222A1Mar 12, 1997Sep 17, 1998Siemens AgX=ray beam generator especially for fast computer tomography in medicine
Referenced by
Citing PatentFiling datePublication dateApplicantTitle
US6996265 *Nov 7, 2000Feb 7, 2006Teradyne, Inc.Inspection method utilizing vertical slice imaging
US7023956Nov 11, 2003Apr 4, 2006Lockheed Martin CorporaitonDetection methods and system using sequenced technologies
US7050533 *Oct 15, 2004May 23, 2006Siemens AktiengesellschaftMethod and device for determining the type of fluid in a fluid mass in an object
US7310407 *Sep 2, 2005Dec 18, 2007Juni Jack ENuclear medical imaging device
US7322745Aug 9, 2004Jan 29, 2008Rapiscan Security Products, Inc.Single boom cargo scanning system
US7327870Mar 14, 2006Feb 5, 2008Teradyne, Inc.Method for inspecting a region of interest
US7369643Jan 12, 2007May 6, 2008Rapiscan Security Products, Inc.Single boom cargo scanning system
US7461032Jan 31, 2006Dec 2, 2008Lockheed Martin CorporationDetection methods and systems using sequenced technologies
US7486768Sep 13, 2004Feb 3, 2009Rapiscan Security Products, Inc.Self-contained mobile inspection system and method
US7517149Nov 30, 2007Apr 14, 2009Rapiscan Security Products, Inc.Cargo scanning system
US7519148Mar 20, 2008Apr 14, 2009Rapiscan Security Products, Inc.Single boom cargo scanning system
US7526064May 4, 2007Apr 28, 2009Rapiscan Security Products, Inc.Multiple pass cargo inspection system
US7551718Aug 22, 2007Jun 23, 2009American Science And Engineering, Inc.Scatter attenuation tomography
US7702069Feb 27, 2006Apr 20, 2010Rapiscan Security Products, Inc.X-ray security inspection machine
US7711089Apr 11, 2007May 4, 2010The Invention Science Fund I, LlcScintillator aspects of compton scattered X-ray visualization, imaging, or information providing
US7720195Jan 7, 2009May 18, 2010Rapiscan Security Products, Inc.Self-contained mobile inspection system and method
US7724871Apr 11, 2007May 25, 2010The Invention Science Fund I, LlcCompton scattered X-ray visualization, imaging, or information provider in soft matter such as tissue, organs, or blood, and/or in hard matter such as bones or teeth
US7734012Apr 11, 2007Jun 8, 2010The Invention Science Fund I, LlcVolumetric type compton scattered X-ray visualization, imaging, or information provider
US7742567Apr 11, 2007Jun 22, 2010Searete LlcCompton scattered X-ray visualization, imaging, or information provider with time of flight computation
US7769133Dec 19, 2008Aug 3, 2010Rapiscan Systems, Inc.Relocatable X-ray imaging system and method for inspecting commercial vehicles and cargo containers
US7783004Oct 31, 2008Aug 24, 2010Rapiscan Systems, Inc.Cargo scanning system
US7817776Mar 3, 2009Oct 19, 2010Rapiscan Systems, Inc.Cargo scanning system
US7856081Mar 13, 2008Dec 21, 2010Rapiscan Systems, Inc.Methods and systems for rapid detection of concealed objects using fluorescence
US7860213Mar 16, 2009Dec 28, 2010Rapiscan Systems, Inc.Multiple pass cargo inspection system
US7876880Mar 2, 2009Jan 25, 2011Rapiscan Systems, Inc.Single boom cargo scanning system
US7924979Sep 1, 2009Apr 12, 2011American Science And Engineering, Inc.Scatter attenuation tomography
US7963695Dec 19, 2008Jun 21, 2011Rapiscan Systems, Inc.Rotatable boom cargo scanning system
US7991113Jun 24, 2010Aug 2, 2011Rapiscan Security Products, Inc.Relocatable x-ray imaging system and method for inspecting commercial vehicles and cargo containers
US7995705Apr 5, 2010Aug 9, 2011Rapiscan Security Products, Inc.Self-contained mobile inspection system and method
US8041006Aug 1, 2007Oct 18, 2011The Invention Science Fund I LlcAspects of compton scattered X-ray visualization, imaging, or information providing
US8059781Aug 2, 2010Nov 15, 2011Rapiscan Systems, Inc.Cargo scanning system
US8090150Oct 10, 2008Jan 3, 2012Morpho Detection, Inc.Method and system for identifying a containment vessel
US8170177Nov 22, 2010May 1, 2012Rapiscan Systems, Inc.Multiple pass cargo inspection system
US8213570Mar 2, 2010Jul 3, 2012Rapiscan Systems, Inc.X-ray security inspection machine
US8275091May 21, 2010Sep 25, 2012Rapiscan Systems, Inc.Compact mobile cargo scanning system
US8356937May 16, 2011Jan 22, 2013Rapiscan Systems, Inc.Rotatable boom cargo scanning system
US8385501Jul 1, 2011Feb 26, 2013Rapiscan Systems, Inc.Self contained mobile inspection system and method
US8389941Dec 22, 2010Mar 5, 2013Rapiscan Systems, Inc.Composite gamma-neutron detection system
US8389942Jun 11, 2009Mar 5, 2013Rapiscan Systems, Inc.Photomultiplier and detection systems
US8428217Nov 15, 2010Apr 23, 2013Rapiscan Systems, Inc.Methods and systems for rapid detection of concealed objects
US8433036Feb 25, 2009Apr 30, 2013Rapiscan Systems, Inc.Scanning systems
US8457275Mar 28, 2012Jun 4, 2013Rapiscan Systems, Inc.Multiple pass cargo inspection system
US8491189Oct 26, 2011Jul 23, 2013Rapiscan Systems, Inc.Radiation source apparatus
US8503605Jul 12, 2010Aug 6, 2013Rapiscan Systems, Inc.Four sided imaging system and method for detection of contraband
US8579506May 20, 2009Nov 12, 2013Rapiscan Systems, Inc.Gantry scanner systems
US8644453Feb 25, 2009Feb 4, 2014Rapiscan Systems, Inc.Scanning systems
US8668386Aug 30, 2012Mar 11, 2014Rapiscan Systems, Inc.Compact mobile cargo scanning system
US8687765Dec 18, 2012Apr 1, 2014Rapiscan Systems, Inc.Cargo scanning system with boom structure
US8712138Oct 5, 2010Apr 29, 2014Koninklijke Philips N.V.Device and method for generating soft tissue contrast images
US8774357Apr 4, 2013Jul 8, 2014Rapiscan Systems, Inc.Scanning systems
US8831176May 20, 2009Sep 9, 2014Rapiscan Systems, Inc.High energy X-ray inspection system using a fan-shaped beam and collimated backscatter detectors
US8837677Aug 1, 2007Sep 16, 2014The Invention Science Fund I LlcMethod and system for compton scattered X-ray depth visualization, imaging, or information provider
US8840303May 20, 2009Sep 23, 2014Rapiscan Systems, Inc.Scanner systems
US8842808Jun 20, 2011Sep 23, 2014American Science And Engineering, Inc.Scatter attenuation tomography using a monochromatic radiation source
US20130032715 *Aug 2, 2012Feb 7, 2013Georgia Tech Research CorporationX-ray compton scatter imaging on volumetric ct systems
WO2011045705A1 *Oct 5, 2010Apr 21, 2011Koninklijke Philips Electronics N.V.Device and method for generating soft tissue contrast images
Classifications
U.S. Classification378/86, 378/87, 378/88
International ClassificationG01N23/20, G01N23/04
Cooperative ClassificationG01N23/20083, G01N23/046
European ClassificationG01N23/20K, G01N23/04D
Legal Events
DateCodeEventDescription
Nov 23, 2009ASAssignment
Owner name: AMERICAN SCIENCE AND ENGINEERING, INC., MASSACHUSE
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SILICON VALLEY BANK;REEL/FRAME:023556/0062
Effective date: 20091118
May 12, 2009FPAYFee payment
Year of fee payment: 12
Mar 11, 2005FPAYFee payment
Year of fee payment: 8
Sep 25, 2003ASAssignment
Owner name: SILICON VALLEY BANK DBA SILICON VALLEY EAST, CALIF
Free format text: SEURITY AGREEMENT;ASSIGNOR:AMERICAN SCIENCE AND ENGINEERING, INC.;REEL/FRAME:014007/0604
Effective date: 20030811
Owner name: SILICON VALLEY BANK DBA SILICON VALLEY EAST 3003 T
Free format text: SEURITY AGREEMENT;ASSIGNOR:AMERICAN SCIENCE AND ENGINEERING, INC. /AR;REEL/FRAME:014007/0604