# Patent

## Title
Identification and localization of explosives and other material

## Abstract

A neutron source illuminates suspect material leading to emission of gamma rays characteristic of the isotopes present. The system measures Compton scattering of emitted gamma rays using detectors with three dimension event localization capability. Detection does not require full energy deposition. A spatial correlation of projection vectors is computed by a reconstruction that searches for solutions that generate spatial correlation. Identification and location for contraband material is determined from solutions that generate spatial correlation.

## Background
RELATED APPLICATIONS
This application is a continuation-in-part of application Ser. No. 12/436,987, filed May 7, 2009 and titled “IDENTIFICATION AND LOCALIZATION OF RADIOISOTOPES” which is a continuation-in-part of application Ser. No. 11/413,480, filed Apr. 28, 2006, and titled “NUCLEAR MATERIAL IDENTIFICATION AND LOCALIZATION”, which claims priority to U.S. Patent Application No. 60/676,082 filed on Apr. 28, 2005, all of which are incorporated herein by specific reference.


TECHNICAL FIELD
The present invention relates methods and devices to identify and localize compounds such as explosives.
BACKGROUND
Neutron activitation is a valuable methodology to determine the identity of materials having a characteristic elemental composition. For example explosives can be characterized and by their unique ratios of Carbon, Ozygen and Nitrogen. Many neutron activitation applications, including identification for explosives relies on the emission of high energy gamma rays. In real situations with shielding material surrounding material the high energy gamma rays are more likely to emerge for possible detection than low energy gamma rays. Unfortunately these same high energy gamma rays are hard to detect with full energy deposition making collection of the energy spectra inefficient. A neutron activation and detection system is desired that can both identify and locate a shielded source of unknown type and quantity using only a few detected gamma rays.
The identification of the illicit movement of suspect materials requires a system capable of identifying the material and locating the source. Innocuous monitoring of movement in open staging or storage areas is highly desirable but difficult to achieve. Conventional neutron activation monitoring techniques involve illumination of the suspect area with neutrons and subsequent detection of emitted characteristic radiation, the identification of the energy of those emitted radiations and correlation of the measured energy or energies to the known emissions spectra of isotopes. Gamma ray spectra are useful for identifying an isotope, but not for locating it. Various imaging techniques can be used in conjunction with the energy spectroscopy including physical collimation of the incident neutrons or directional information of the emitted gamma-ray radiation by collimation, coded aperture arrays or conventional Compton cameras. Current methods for locating material by neutron activation typically rely on collimators and require long data acquisition times when the location of the material is not know.
SUMMARY OF THE INVENTION
Disclosed are embodiments of methods for identifying the type and location of materials using the neutron activation technique. In one embodiment, multiple radiation detectors surround the area of inspection and record interactions involving Compton scattered events. A spatial correlation method is employed to identify the location and type of material based on a geometric convergence of possible backprojected paths. Searching is accomplished using a list of known gamma ray decay energies for isotopes of interest.
The system utilizes a spatial correlation of projection vectors arising from Compton scattering of any, but especially of, the high energy emissions. Source identification and location is supplied by a reconstruction that searches for solutions with isotopes of unknown type. Detection does not require full energy deposition. Identification and location for contraband material somewhere in a large active area of interrogation is possible with many fewer detected events than conventional methods.
A neutron source is used to activate nuclei in the suspect area. Gamma ray emissions from the decay of neutron activated nuclei interact in the detector system. The detector system includes an array of detectors allowing an initial interaction via Compton scattering and a second interaction. The second interaction can be Compton scatter, photoelectric absorption, pair production or any other interaction that deposits energy in the detector. The deposited energies and the 3D location of the interactions is collected by a data acquisition and analysis system. A search mechanism looks for possible solutions for multiple events that will results in a geometric correlation. The resulting geometric correlation identifies the type of material and location.



DESCRIPTION OF THE FIGURES
FIG. 1 shows diagram of conical surfaces from three detector events backprojected from measurement or scattered photon energy at a single-point-of-observation detector system indicating the most probable location of the source at the intersection of the circles.
FIG. 2 shows diagram of a conical surface backprojected from measurement of scattered photon energy at a single-point-of-observation detector system.
FIG. 3 shows the detector system dispersed around control area. Backprojected Compton cones intersect at the radioisotope source location showing how the intersection of backprojected cones determines the location of a radioactive material.
FIG. 4 is a diagram of Compton scatter kinematics.
FIG. 5 is a flow chart showing the data analysis and processing for creating an image space corresponding to true and random backprojection cones.
FIG. 6 is a flow chart showing the data acquisition, analysis and decision making paths for identification of a radioisotope.
FIG. 7 is a diagram of the overview of an embodiment of how the system would be configured.
FIG. 8 is a diagram of the architecture of an embodiment of how the detector array would be configured.
FIG. 9 is a diagram of the architecture of an embodiment of how a detector node would be configured.
FIG. 10 is a diagram of the system configuration for illuminating the material in question with a neutron beam while observing the neutron activation gamma-ray emissions from a single detector position.
FIG. 11 is a diagram of the system configuration for illuminating the material in question with a neutron beam while observing the neutron activation gamma-ray emissions from multiple detector positions.



DETAILED DESCRIPTION OF THE INVENTION
A material detection, localization and identification system is disclosed. An array of radiation detectors is networked to act as a unified detection system. This array can have a wide area of view and high sensitivity to incident radiation. The system utilizes the properties of Compton scattering of gamma radiation. Rather than reliance on measuring full energy interactions for identification and possible localization, the system uses a geometric correlation to identify and localize radiation sources.
Neutron activation analysis is a valuable technique for the quantitative multi-element analysis of major, minor, trace and rare elements. Neutron activation analysis begins with neutron bombardment 130 of a sample 120 to convert stable isotopes to radioactive isotopes. The procedure involves irradiating the samples 120 to be analyzed in the neutron field produced by a neutron source 110. Radioisotopes created during the irradiation will decay with time. A portion of the energy released during decay is often in the form of gamma radiation 140, which is capable of traveling out of the sample. The gamma rays possess unique energies that are characteristic of the isotope undergoing decay. Gamma rays detected at a particular energy are indicative of a specific radionuclide's presence.
Data reduction of gamma ray spectra yields the concentrations of various elements in samples being studied. Neutron activation analysis sensitivities and accuracy are dependent on the concentration of a particular element and radionuclide parameters (i.e., parent isotope abundance, neutron cross-section, half-life, and gamma ray abundance).
Neutrons may be produced 110 using a number of techniques including radioactive isotopic sources, electrophysical neutron generators and accelerators.
Isotopic neutron sources produce continuous fluxes of neutrons. Typical isotopic sources are Californium-252 (252Cf), with a half-life of about 2.6 years (one mg 252Cf produces about 2.3·106 n/s), or Americium-Beryllium (AmBe), which produces neutrons via the 9Be(a,n)12C reaction (241Am has a half-life of 458 years). Isotopic neutron sources have the advantage of having a long useful life and producing a relatively constant flux of neutrons. They are relatively inexpensive for low flux (<108 neutrons per second) sources. However, isotopic sources have several disadvantages. The neutron output cannot be turned off; requiring that they be contained within bulky shielding at all times. Isotopic neutron sources cannot be pulsed and the energy spectrum of the emitted neutrons is broad and peaks at energies below the threshold for some important reactions.
Small neutron generators using the deuterium (2H)-tritium (3H) reaction are the most common accelerator based (as opposed to isotopic) neutron sources. Creating deuterium ions and accelerating these ions into a tritium or deuterium target produces neutrons. Deuterium atoms in the beam fuse with deuterium and tritium atoms in the target to produce neutrons.

d+t→n+4He En=14.2 MeV

d+d→n+3He En=2.5 MeV

d+t reaction has the largest maximum cross-section of 5.0 Barn (10−24 cm2) of all fusion reactions. Maximum cross-section of this reaction for energies of incoming particle below 1 MeV is reached at the energy of 130 keV.
Neutrons produced from the d-t reaction are emitted isotropically from the target. Neutron emission from the d-d reaction is slightly peaked in the forward (along the axis of the ion beam) direction
Most small d-t accelerators are sealed tube neutron generators. The ion source, ion optics and the accelerator target are enclosed in within a vacuum tight enclosure. Either glass or ceramic insulators provide high voltage insulation between the ion optical elements of the tube. The neutron tube is, in turn, enclosed in a metal housing, the accelerator head, which is filled with a dielectric media to insulate the high voltage elements of the tube from the surroundings. The accelerator voltage is typically between 80 and 180 kilovolts.
The accelerated ions strike the target. The target is usually made of titanium, scandium, or zirconium that form stable chemical compounds called metal hydrides when combined with hydrogen or its isotopes. These metal hydrides are made up of two hydrogen (deuterium or tritium) atoms per metal atom and allow the target to have extremely high densities of hydrogen. This is important for maximizing the neutron yield of the neutron tube. The gas reservoir element also uses metal hydrides as the active material.
The neutron generator does not create any radiation when it is switched-off. They may be operated either as continuous or pulsed neutron sources. The neutrons produced are mono-energetic (2.5 MeV or 14 MeV).
An array of radiation detectors is networked to act as a unified detection system. This network can have a wide area of view and high sensitivity to radiation sources. The system utilizes the properties of Compton scattering, well suited for high energy gamma radiation. Rather than the traditional reliance on measuring full energy interactions for identification and possible localization, the system uses a geometric correlation to identify and localize sources.
The disclosed system works for all gamma ray emitting radioisotope sources but especially well for sources with high energy emissions emerging, such as those commonly produced in Neutron Activation Analysis. Current systems using total energy identification work best for lower energy gamma rays, where full energy absorption is more probable.
In one embodiment the deployment of the system is accomplished by stationing a neutron source and a plurality of radiation detectors with three-dimensional and multiple-interaction-point functionality around the area to be monitored. These detectors 101 are relatively inexpensive to produce compared with large detector systems typically used for full energy measurements. System sensitivity and area of coverage can be scaled by adjusting the number of individual detector stations 101. Data from the multiple locations are collected in an analysis station. The system will continue to operate (at slightly reduced sensitivity) even if individual detector modules 101 become inoperable.
In another embodiment the area to be interrogated is irradiated by a neutron source and viewed from one location by a detector system with 3-D event localization capability. In another embodiment a neutron source and a mobile detector system or systems can move around the area of interest and record events. In yet another embodiment a combination of mobile and stationary neutron sources and detectors can be used to acquire the data.
The system operates is a very different mode than conventional radiation detection and localization systems. Each detector 101 has three dimensional event localization capability and the capability to measure individually multiple interaction points occurring simultaneously in the detector module. Each 3-D radiation detector 101 detects Compton scattered events, identifying the primary 405 and secondary 406 interaction vertices. While the detector module 101 identifies the energy and three dimensional positions for the two interaction vertices, the detector module 101 does not need to identify which vertex is the primary and which is the secondary. The detected coordinates are used by the analysis engine to search for combinations of primary and secondary that are consistent with known physical scattering properties of gamma ray emissions from suspect neutron activated material. The disclosed detection and localization system does not require detection of the full energy of each event to identify the neutron activated source, localize the source and identify the isotope.
The system benefits from the fact that Compton events are easier to detect than full energy events for high energy emissions. We obtain information on the direction of the original gamma ray if the primary 405 and secondary 406 interaction points for Compton events are measured. The interaction point 405 is a first interaction of the incoming gamma ray and interacts with the detector 402 by a Compton scattering interaction. A scattered gamma ray 404 (the Compton scatter secondary) emerges from the detector 402 and interacts in a second detector 401. The interaction in detector 401 may be any interaction that deposits energy in the detector, i.e. may be a Compton interaction of the scattered gamma ray, may be a photoelectric absorption with or without x-ray release from the detector or may interact by pair production. In interaction 401 the event may have interacted by leaving either full or partial energy of the Compton scatter secondary. While measuring full energy of the gamma ray is difficult, we start with the exact knowledge of the potential energies' of gamma rays we are interested in detecting, e.g., emissions from the activation of Carbon, Nitrogen or Oxygen. The use of coincidence in a single detector 101 or among detectors 401, 402 greatly reduces the random background. Looking for additional geometric correlation further reduces the false positive rate and provides a very reliable detection system.
Detector systems 101 to measure energy, timing and three dimensional position from gamma ray interactions are well know to those skilled in the art. Detector systems 101 vary in their energy and position resolution, timing resolution, in their stopping power characteristics and in their cost. Choice of detectors 101 will vary as one selects the target resolution and sensitivity desired while balancing cost considerations. The system described here will work for any of the detector choices made while performance results will vary accordingly.
Possible detector systems 101 are: CsI(T1) segmented detector using Photo Multiplier Tubes (PMTs), CsI(T1) segmented detector using pin diodes, CsI(T1) segmented detector using avalanche diodes, NaI(T1) segmented detector using PMTs, NaI(T1) single crystal planes using PMTs (similar to a conventional gamma camera), CsI(T1) single crystal planes using PMT's (similar to a conventional gamma camera), BGO segmented detector using PMTs, LSO segmented detector using PMTs, Silicon pixilated detector, Intrinsic Ge pixilated detector, CZT pixilated detector planes, CZT volumetric detector, gas scintillation detectors or gas proportional detectors. Other similar gamma ray detectors can be used within the scope of this invention as would be understood by one of ordinary skill in the art of gamma ray detection.
The detectors 101 provide measurement 601 of energy and position and event timing information. Timing information is used to measure relative time between two or more distinct interactions in the detector. If the two or more distinct interactions are from a single Compton scatter event, the relative timing of the events will be within the time for photons to traverse the distance between the event positions, usually on the order of a nanosecond. The timing resolution of the detectors is usually greater than this. Possible detector timing effects, such as charge transport or electronics response can make the coincidence window larger. One sets a coincidence time window to allow simultaneous events to be registered as possible multiple vertices from a single Compton scattering event. Multiple vertices can arise from two separate events that happen to interact in coincidence. The probability that accidental coincidences occur is a function of detector count rate. For most applications the count rate is very small, making accidental coincidences very rare.
The measured energy, timing or coincidence, three dimensional positions are recorded for further processing in a data storage device 602, 501. Most typically this data storage device will be the storage memory of a computer, such as a hard disk drive or into other digital memory of the computer. Alternative means to store data may be used if it provides the capability to further process the data as described in this specification. Timing information is used to determine if there is coincidence between two or more interaction vertices. Timing information for each vertex can be stored or simply an indication that the vertices are connected in coincidence may be used. The timing information establishes that the vertices are part of the same Compton scatter event. Energy data is typically produced as the channel output from an ADC and the digital channel is stored or can be converted to another scale such as a calibrated energy value. Energy calibration may occur here but can also occur later in the processing in connection with or prior to looking up hypothetical neutron activation emission energies 502. When multiple detectors are used as the gamma ray detector then a detector identifier or equivalent will be recorded with the event information. If a mobile detector is used, either alone or with other detectors (mobile or stationary) then the detector identifier (if more than one detector) and position of mobile detectors will be recorded. Position measurement of the mobile detectors can be accomplished by several techniques including GPS signals, encoded position readouts, RF tagging, optical systems and other means known to those skilled in the art or yet to be developed.
Events are selected from the date pool 602 following a variety of mechanisms. In one embodiment the events may be analyzed in the time sequence that they are acquired. In another embodiment the events may be analyzed following a predetermined algorithm that allocates certain energies or energy ranges more significant than others. The more significant events are selected and analyzed first in anticipation of more quickly ascertaining the identity of the radioisotope or radioisotopes. Alternatively the events may be selected in a totally random mode, or in a mix of methods. Other means of selecting event order that can accomplish the same function are within the scope of this invention.
This system uses the knowledge of the known neutron activation emissions from suspected isotopes 603. The neutron activation emission gamma rays are cataloged and available to those skilled in the art or nuclear spectroscopy or radioisotope identification. In addition to the energy of the gamma rays the relative emission strengths of the emission gamma rays can be included if available. One such depository of emission gamma rays is the National Nuclear Data Center of the National Institute of Standards and Technology (NIST) of the United States Department of Commerce. This neutron activation emission gamma ray data is loaded into the system in a readable form 603. The most common means to store this information is into the storage memory of a computer such as a hard disk drive or into other digital memory of the computer. Other means that can accomplish the same function are within the scope of this invention.
Within the calculation of possible backprojection cones arising from a set of vertices, the energy of the original gamma ray is typically unknown. If energy information is available for both vertices then the energy of the original gamma ray associated with a set of Compton vertices is equal to (in consideration of detector energy resolution limitations) or greater than the sum of the energies corresponding to the set of vertices. From the set of possible neutron activation emission gamma-ray energies under consideration for the particular isotope identification and localization task engaged a possible neutron activation emission gamma-ray energy is chosen 605, 502 whose energy is equal to or greater than the measured energy for the event. From the subset of neutron activation emission energies with acceptable energy, multiple methods of selection are possible. In one method, all acceptable neutron activation emission energies are chosen for backprojection calculation. Then each cone is calculated and backprojected with equal weights to form a total image. In another method using all the acceptable neutron activation emission energies the cones are backprojected using weights proportional to emission branching ratios and energy dependant detector efficiency. Another method will use emission energies from a class of isotopes (e.g. those found in explosives) and backproject them either with or without a weighting factor.
The equation for Compton scattering relates the angle the photon is deflected from its path (φ) to the energies of the scattered (hv) and incident (hv0) photons:




  
    
      
        φ
        =
        
          
            
              cos
              
                -
                1
              
            
            ⁡
            
              [
              
                1
                -
                
                  
                    m
                    0
                  
                  ⁢
                  
                    
                      c
                      2
                    
                    ⁡
                    
                      (
                      
                        
                          1
                          hv
                        
                        -
                        
                          1
                          
                            hv
                            0
                          
                        
                      
                      )
                    
                  
                
              
              ]
            
          
          .
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        1
      
    
  




This angle provides the basis for all reconstruction algorithms involving Compton scattering. If the scattered photon energy is known, the angle of deflection can be obtained from knowledge of the incident photon energy. Since any azimuthal deflection is possible the direction of the initial photon in three-dimensional space can be limited to a cone 403. Note the direction of the final photon also is to be determined, which in this example is accomplished with the use of a detector array 401. In general, three or more measurements are acquired by the detector to pinpoint the source 404, 103 of the initial gamma photon. This is accomplished with the use of reconstruction methods that essentially determine the intersection of the backprojected cones. This intersection represents the most probable location of the radioactive material.
For each coincidence event, the spatial coordinates for two or more vertices are identified along with the energy deposited in the detector at the associated vertices. In some cases, the first point of interaction is distinguishable from the second point of interaction due to the energy partitioning according to the physics of Compton scattering. In other cases it will not be possible to distinguish which is the first interaction vertex and which is the second interaction vertex. In cases with more than two interaction vertices in coincidence it is unknown what vertex is the first, second, third, etc. For those cases where it is possible to assign an unambiguous first interaction point, that unambiguous first interaction point is used as the origin of the backprojection cone. The axis of the cone is established as the line from the first vertex to the second interaction vertex. The backprojection cone 403 is then the cone originating from the origin with an axis to cone surface angle equal to φ as calculated from equation 1. The energy of the incident photon (hv0) is chosen 605 from the table of possible emission energies as the hypothetical incident energy and the energy of the scattered photon (hv) is the energy of the incident photon minus the measured interaction energy at the first vertex.
For those cases where it is not possible to assign an unambiguous first interaction point, two projections are calculated with each vertex used as the origin of the backprojection cone. The axis of the cone is established as the line from the assumed first vertex to the remaining second interaction vertex. The backprojection cone is then the cone originating from the origin with an axis to cone surface angle equal to φ as calculated from equation 1. The energy of the incident photon (hv0) is chosen from the table of possible neutron activation emission energies as the hypothetical incident energy and the energy of the scattered photon (hv) is the energy of the incident photon minus the measured interaction energy at the first vertex. For events with more than two coincidence vertices, each vertex will be used as the origin with the line to each other vertex as the axis. In configurations where not all the detectors produce an energy measurement, only the detector events that have energy measurement are used as origins for backprojection cones. One means of measuring overlap of calculated backprojection cones is to fill the image space with backprojection cones and look for the region of maximum overlap 401, 405. From the measured data set the set of hypothetical neutron activation emission energies 502 from a single radioisotope or a set of known related radioisotopes (for example as might be found in a nuclear weapon with a mix of radioisotopes and daughter radioisotopes) compatible with the data is identified. The backprojection cones 503 are calculated for each of the hypothetical neutron activation emission energies for each of the event data. The backprojection cones 503 are added 504 to an image space resulting in an image with relatively uniform background. The relatively uniform background is a result of the backprojection cones that do not correlate with the actual energies and directions of the gamma rays that produced the event data. This background is spatially uncorrelated because on average the incorrect cones direction and angles will be random. If that radioisotope is present then the set of correct backprojection cones will intersect at a point corresponding to the location of the radioisotope. Examination of the image space for such a region of maximum intensity 505 corresponds with a positive identification of the presence and position of the radioisotope 506. If no region of enhanced image intensity is present another possible neutron activation emission isotope can be used to process the data. The system may examine all possible neutron activation isotopes in its library if necessary to find the source of radiation.
One embodiment of an evaluation process follows the flow chart of FIG. 6. In this process backprojection cones are projected into an image space. Regions of overlap are identified for the set of backprojected cones. Additional events are backprojected into the image space and areas of multiple overlap identified. This process looks for regions of overlap for multiple possible events 607. When regions are identified they are stored and event cones that do not produce overlap are recalculated 611 using an alternative neutron activation emission energy hypothesis 605. In this manner the data set is searched through for possible neutron activation emission energies that result in a region of multiple overlap. If no region of multiple overlap is identified, the process begins again with another isotope chosen for hypothetical neutron activation emission energies.
Additional means common in imaging techniques can be used to create or explore the overlap space of the backprojection cones. A reconstruction algorithm for the linear superposition method is as follows. First, the algorithm constructs conical shells arising from the backprojections. This procedure is executed using spherical harmonic expansions and coordinate transformations for different coordinate systems that are specific to each detector. Other geometric approaches for constructing the shells are known. The algorithm then assigns unit amplitudes to the shells. Linear superposition of these shells in the image space will sum the amplitudes where multiple shells intersect. The amplitudes will therefore be proportional to the number of events (cones) that coincide (intersect) in space. Identification of the 3D location of the source will therefore be produced.
Backprojection of cones is one representation of backprojection data. Backprojection data is data that represents the possible trajectory paths that the incoming gamma ray may have followed that are consistent with the measured event data. Geometric cones in Cartesian coordinates may be used as well as, for example, other representations from cylindrical or spherical coordinates, special mathematical functions or representations in less than three dimensions. If the source location is known to lie in a plane (e.g. know to be on the surface) then backprojection data may be lines rather than cones.
Matrix methods for Compton search array reconstruction may use multiple spherical coordinate systems with coordinate transformations between systems. Iterative solution methods are combined with a generalized voxel in the reconstruction algorithm.
Matrix methods are a straightforward approach for reconstructing images from tomographic or similar data sets. The reconstruction problem is formulated as a transformation between the space/source distribution to be imaged and the response of the detectors. This transformation is a function of the geometry of the photon paths and the detector positions. By expressing the transformation as a set of linear equations (i.e., the projection matrix), the reconstruction problem can be solved using matrix algebra (e.g., matrix inversion). Matrix reconstruction is sometimes called the algebraic reconstruction technique (ART).
Matrix inversion methods have been presented for Compton cameras. Symmetries in the reconstruction problem were used to simplify the algorithm. The sampling geometry consisted of nested hemispherical shells sectioned into a grid of voxels of uniform angular size (Δθ and Δφ). The voxels were approximated with polyhedrons to compute the relationship between the image space and the projection matrix. Finally, the 3D spatial computations were simplified by transforming the matrices into the 2D Fourier domain and solving the matrix equation with an iterative reconstruction algorithm. Iterative solution of matrix equations offers several advantages in reconstruction problems, particularly for very large or very sparse matrices where matrix inversion is difficult or inaccurate.
Direct linear superposition of backprojected images can be an effective reconstruction method, similar to a triangulation method. For Compton scattering, the backprojected images would be comprised of the conical shells representing regions with a high probability of containing a radioactive source. Superposition of these shells defines regions of higher probability (i.e., where the shells intersect) for location of the source distribution. Backprojection algorithms that use linear superposition (summation over all detected events) are known in the art for Compton scattering data.
One embodiment of the system operates in the following manner. A neutron source 110 irradiates the area of interest. A network of small radiation detectors 101 with 3-D position readout are dispersed preferably around the perimeter of the area of interest. These detectors 101 are networked together and feed their data to a central processing station. Emissions from the area of interest 140 caused by neutron 130 activation are measured. Compton scatter events (at least the primary and secondary interaction vertices and the energy of the primary interaction) are recorded from the network. In each detector 101 we do not know which interaction vertex is the primary and which interaction vertex is the secondary. Coincidence within each detector module 101 tells us that the events are part of the same Compton interaction event. The event rate can be very slow and events in multiple detectors are accumulated for analysis over any desired data acquisition time used to control the area.
When multiple events are recorded the data are analyzed to compute possible Compton projection cones from the vertices. An assumption is made that the events are from a particular class of isotopes (e.g. those found in explosives) 502 with emissions equal or greater energy than measured by the sum of the energies for the two coincident Compton vertices. For example, if a coincidence produces a sum of energy for an event of 700 keV, we compute possible trajectories for known neutron activation emission gamma-rays, for the class of isotope hypothesis (e.g. those found in explosives), greater than or equal to 700 keV. Using the proposed energies of the suspected gamma rays, the projection cones are analyzed to determine if a single point of intersection exists (to within error calculated from known energy and spatial resolutions of the detector). If a true single source of radiation is emitting the radiation at that energy, then a single point of intersection will be found by the cone projection analysis. If the events are random then the probability of producing a single interaction point will be very small. Additional events in any case can bring false positives to any desired lower limit. A threshold level can be set to trigger identification of the source.
Multiple energies from the analysis data base are used to see if the projection cones produce a single focal point of interaction. If multiple energies are known to be emitted from a suspected source then multiple energies can be tested (consistent with emission energy larger or equal to the detected energy). When a single focal point is detected, it indicates that a particular set of isotopes (that associated with the energies used to analyze the data) is present and the location of the source is identified as that detected point. The proportions of elements within the set of isotopes so identified can now be used to identify if the material is of a suspected class, such as explosives.
If the system is deployed around the perimeter of a control area, the system has the highest sensitivity for sources located within the perimeter but is still sensitive to sources outside the perimeter. The system has superior localization because the detectors are situated at various locations around the control area and angular dependent resolution is combined from those multiple positions. Conventional detectors situated at a single location will have difficulty localizing a source that is more distant from the detector. Because the detectors are small modules, the system can be expanded by adding more detectors modules at more locations to improve sensitivity or to monitor a larger area. Specific applications for this technology include; surveys of warehouses, bunkers, and other storage facilities for locating materials; and the monitoring of transportation facilities and other areas for contraband materials such as explosives.
In a similar manner, the Compton search array would use three or more widely spaced detectors to triangulate on the location of material. Note the intersection of two cones would form a line curving in three dimensions. The addition of a third cone produces additional lines, with the source location defined by the intersection of either the three conical surfaces or the corresponding set of two-cone lines. Backprojected cones from additional detectors would further refine the location determination.
Several reconstruction methods are available for analyzing data from the disclosed Compton search array. Two example approaches that have been used for Compton reconstruction are matrix methods and linear superposition. Any computational method for reconstructing the backprojected cones in three-dimensional space and determining the intersection points of the cones may be used. The reconstruction method accommodates arbitrarily positioned detectors, and therefore a variety of imaging volumes and multiple cone geometries. Spreading of the cones due to energy and position uncertainties would also require consideration, and a “probability of location” analysis may be used since multiple intersection points may result from events detected from gamma rays of different energies.
The initial energy of the neutron activation emission gamma-ray 140 is treated as unknown. Search algorithms based on the known physical distributions, attenuation and energy dependent detector response provide the most optimum searching methods for the actual gamma ray energies. These search algorithms involve searching for solutions causing multiple events backprojections to intersect.
Because of their similarity, reconstruction methods developed for Compton cameras and SPECT imaging systems may be adaptable to Compton search array data. The majority of reconstruction methods developed for Compton imagers use the backprojection method. The backprojection method is based on the physics of Compton scattering, and is often combined with other tomographic methods such as line integrals, Radon inversion, Fourier transform inversion, and iterative techniques. Methods have also been developed to facilitate reconstruction of backprojections using spherical harmonic expansions.
In backprojection, Compton scattering is analyzed in a reverse fashion. Measuring the angle and energy of the Compton-scattered photon defines a conical surface in space from which the incident photon originated. Uncertainties in both the angular and energy measurements, however, convert the idealized conical surface into a more realistic conical shell. This shell is the region in space defined by the probe to contain the photon source. The backprojection method reconstructs these conical surfaces or shells from multiple detection events (photons). The intersections of the surfaces or shell volumes pinpoint the most probable locations of the emission sources.
A search array would be configured to surround the area of inspection as much as practical. This may include the placement of detectors at different elevations (stair platforms, roof supports, and ceilings) as well as floor (or ground) positions. The reconstruction is fully three-dimensional with variable reconstruction geometry. In contrast to the Compton camera, the Compton search array will sample gamma rays from a much wider range of angles. In tomography, greater coverage typically results in more accurate and robust image reconstruction.
In another embodiment the neutron source 110 irradiates the area of interest 120 and the detector system 101 operates at a single-point-of-observation detector location with 3 dimensional event localization capability and event energy information. The detector can be any such system capable of delivering such event data. The detector or collection of detectors will be able to measure two interaction vertices in coincidence for a single event. A monolithic detector structure with this operational characteristic can be constructed using a variety of techniques known to those skilled in the art such as but not limited to; light segmentation, multiple charge collection electrodes or time of flight measurements. Similarly a segmented detector structure can perform the data collection using segmented scintillators, solid state detectors or gaseous media. In general the detectors can consist of scintillators, solids state detectors, liquid scintillators, liquid charge drift systems, gaseous scintillators or gaseous charge drift systems. It will be obvious to those having skill in the art that many changes may be made to the details of the above-described detector without departing from the underlying principles of the invention.
In general when using a detector system with a single point of observation the final localization information will be limited to direction to the source location. If the single-point-of-observation detector is mobile then additional localization will be possible if the event stream is sampled from multiple locations. The event stream can be sampled along a continuous path as the detector is moved or can be sampled from observation point to observation point. When observing with a mobile or transportable system, a Global Positioning System (GPS) for detector position and an inertial measurement unit (IMU) for detector orientation can be used to provide better absolute position information for events and for final source localization. An IMU works by detecting the current rate of acceleration, as well as changes in rotational attributes, including pitch, roll and yaw. This data may then be fed into a computer, which calculates the current speed and position, given a known initial speed and position. IMUs may be only one component of a navigation system. Other systems may be included and their information may be used in place of IMU data or used to correct the inaccuracies that IMUs inevitably suffer. Such systems are GPS, gravity sensors (for local vertical), external speed sensors (to compensate for velocity drift), a barometer for altitude correction and a magnetic compass.
The method of search can follow several possible strategies. From a single point of search it might be advantageous to project cones in polar or cylindrical coordinates. Any coordinate system can be used, but some will be more conducive for computation of different particular situations. Cartesian coordinates might be preferred for a system with multiple detector locations or a moving detector system.
Search methods can include the actual energy measured for an event as a possible solution along with hypothetical energies associated with a particular isotope used for other events. To enhance the localization and overlap, the angular error can be calculated for each event according to the hypothesis or the actual energy measured. In general if the actual energy measured is used the angular error will be larger due to statistical error in the total energy. For hypothetical energies it is assumed that the energy is either completely correct or completely wrong. Energies that are completely wrong will not lead to geometric correlation. In some cases it is possible that the photons could Compton scatter within a selected angular range. The search algorithm could use this to compute the possible cone projections using a range of potential scatter angles.
The image quality and resultant location accuracy from the Compton search array will be dependent on the number of events detected. Actual angular resolution is not equal to the angular resolution of the imaging systems for a single event. System resolution will be better than from a single event and will be statistically improving with number of true events detected. Since geometric correlation is used, the nominal direction determination of the system will be that of maximum correlation while the resolution will be determined by statistical variation of that data set.
In one embodiment of the search, the mix of unknown energies to use in a hypothetical search allows each data set to stand on its own in conjunction with hypothetical energies for other events. When doing this the overlap region may be blurred to account for different errors among the different possible solutions cones. The blurring algorithm may account for the general resolution inaccuracies and may be a probability based distribution of locations. Those skilled in the art would be able to utilize a variety of techniques to account for the limited resolution blurring.
Additional parameters may be incorporated into the search process based on the knowledge about their relations to the physics of the measurement process. For example additional information that may be added to the search analysis include: singles rates, singles rate detector dependent, single rate detector side dependent, energy values from singles, energy from coincidences, probability of gamma ray energy emission, and mix of energies in coincidence. In addition there could be physical constraints placed upon the search analysis such as requiring the source location to be limited to a certain area (e.g. limited to within a cargo container or vehicle). If more flux is on one side of the detector, some indication of the origin is possible. One searches on solutions that project into this direction rather than all directions. If the source will be known to be located on a plane, then solution searches will be restricted to solution sets that overlap that plane and projection of lines rather than cones will be sufficient. This will be the case of a detector in the open looking for materials when there are no means of elevation (hills, buildings, flying vehicles, etc.) In a more sophisticated system the ground contour could be mapped and a viable location space defined to limit searching that could be more sophisticated than a plane, could include curved surfaces like hills or partial 3D areas such as buildings. The examples and embodiments disclosed herein are to be construed as merely illustrative and not a limitation of the scope of the present invention in any way.
The equation for Compton scattering relates the angle the photon is deflected from its path (φ) to the energies of the scattered (hv) and incident (hv0) photons 140. The angle and event vertices provide the basis for all reconstruction algorithms involving Compton scattering. The data E1 and E2 measured energy are the energies measured by the detector system.
In one embodiment the search method operates by searching the possible source type and location solution space to find one that maximizes the likelihood that a source type and source location could account for the measured data. The search algorithm proposes a hypothesis for source material, and thus a collection of possible neutron activation gamma-ray energies. We then test the hypothesis on the data to see if there is any geometric correlation. We test other hypotheses and compare to determine what the most likely source material is and identify the location from overlap of the back-projected cones. We add additional data as it becomes available to either validate the leading hypothesis or it may give credence to an alternative hypothesis. When multiple events are correlated to a particular direction a solution certainty can be computed. Different thresholds of certainty can provide different alarms to users.
The search method is adaptable to incorporate additional peripheral information to increase accuracy of the detection system. For example, the detector could experience higher flux on one face of the detector than the others. This directional evidence may lead the algorithm to focus on a more confined solution space and thus a more efficient final solution. The known distribution of singles energy and spatial distribution may be used to modify the prioritization of the search algorithm. Also, if there is a specific direction of interest the system can probe that direction of search space through algorithmic methods while the detector stays active to possible sources in all directions. Sources from more than one location or more than one source material can also be identified.
Spatial resolution for this method is superior to conventional Compton camera methods. The errors in scattering angle associated with conventional Compton cameras are primarily determined by the uncertainty of the measurement of the initial gamma ray energy and the uncertainty in the scattered photon. As noted in Equation 1, the angle is dependent on knowledge of both the initial total energy of the incoming gamma ray and the energy of the scattered gamma ray, thus the angular error is influenced by both the error of the initial total energy of the incoming gamma ray and the error of the energy of the scattered gamma ray.
Conventional Compton cameras rely on the measured total energy, E1+E2, and measured scattered energy, E2. Error in the computed angle is dependent on the errors in both of these measured energies. In absolute terms (keV) the uncertainty in the full energy, E1+E2 is greater because the initial gamma ray is always higher energy than the scattered gamma ray, E2. The analysis does not require full energy events, and can use more detected events, for example, Compton then Compton rather than Compton then photoelectric, without the concerns of lost energy from traditional Compton imagers.
When desired, the angular field of view can be controlled electronically. When a restricted field of view is desired the search for solutions can be constrained to only within that field of view. For example if there is a suspected source location we can more quickly determine if radiation is being emitted from that location by requiring only proposed solutions that overlap (within resolution constraints) with that location. One further feature is that, should the suspect location not be the correct one, further analysis of the captured data without the imposed restriction could yield the true source location. That would not be possible for a system that had a physically restricted field of view.
If no hypothesis generated by the search algorithm will provide geometric correlation, the system can still operate as a conventional Compton camera to see if the lower resolution imaging will provide an intersection of overlap, indicating a possible source of either some not listed material or perhaps shielding, such that only secondary scattered gamma rays are emerging from the shielding. In either case, systems that rely on full energy capture for identification will fail to identify the material. If the source material is not in standard libraries, it will not be identified even if full energies are detected.
False alarm rates are very low. Not only is geometric correlation hard to achieve by random cone projections, but since the detector is 4 pi sensitive, even if a random cone projection passes threshold, it would still have to intersect with a viable object to be a false alarm, i.e. if the random cone projection was in the direction over the detector where no objects were located.
In one embodiment of the system, the detector system may contain a primary Operator Control Unit (OCU) 701 that may be based on a single board computer (SBC). The OCU 701 will allow the user to operate the system and receive radiation detection status and location information.
The system controller 703 may be based upon a motherboard such as the Systel IPC 4472-4U Fault Resilient Rack Mount Military Computer with a 7 Slot ATX Active Motherboard. The IPC 4472 features;

    
    
        Shock mounted drive bays and motherboard
        Positive pressure front cooling with one 102 CFM ball bearing fan and 2 rear I/O fans
        Firmly held plug-in cards
        Internal chassis bracing
        Wire bridges to eliminate the need for stick-on tyrap holders
    
    


The system controller interfaces 704 to the Detector Array 702 using, for example, a MOXA CP-134U V2 Universal PCI Multiport Serial Board that has 4 RS-422/485 serial ports, each capable of data rates up to 921.6 Kbps. For this system the ports are configured in multidrop, 4-wire, RS-485 modes capable of interfacing to multiple nodes per port.
The system controller may use a combined Global Positioning System (GPS) 705 and Inertial Measurement Unit (IMU) 706, for example the one provided by NovAtel, for location of the system. Using NovAtel's Synchronized Position Attitude Navigation (SPAN) Technology featuring tight integration of a GPS and IMU would allow the system to provide uninterrupted position and attitude information during short periods of GPS outage, or when less than four satellites are received. From this position and attitude information the system is able to determine the relative location of a radiation source.
Communication capability of the system could be provided by a MultiTech Systems MTCBA-G-U-F2 GPRS wireless modem 707. This standalone modem provides wireless data communication based on industry-standard open interfaces. Some of the key features of the modem are:

    
    
        GPRS Class 10
        Dual-band 850/1900 or 900/1800 MHz GSM/GPRS
        Packet data up to 85.6K bps
        Circuit-switched data up to 14.4K bps transparent and non-transparent
        Short Message Services (SMS)
        USB interfaces
        SMA antenna connector and SIM socket
        Serial interface supports DTE speeds to 115.2K bps
        AT command compatible
        MNP2 error correction, V.42bis compression
        PTCRB certified
        Carrier approved
    
    


In one embodiment the detector system has four detector arrays 702 each of which have 31 detector nodes. This gives the system a total of 992 individual scintillation detectors in this embodiment. The system may be built using Detector Nodes 801. In one embodiment the detector 902 is a Thallium-activated sodium iodide scintillation detector, i.e. an NaI(T1) crystal coupled to a photomultiplier tube (PMT) 903.
The node controller 901 may be based upon the Xilinx Spartan-3 FPGA. This allows for flexibility in the node design and allows for a distributed architecture. This distributed architecture is well suited for determining correlation of events in a single node or between multiple nodes, providing a 10 us time resolution for coincidence identification.
The above description discloses the invention including preferred embodiments thereof. The examples and embodiments disclosed herein are to be construed as merely illustrative and not a limitation of the scope of the present invention in any way. It will be obvious to those having skill in the art that many changes may be made to the details of the above-described embodiments without departing from the underlying principles of the invention.

## Claims

1. A material identification device comprising:
a neutron source;
said neutron source illuminated said material;
said illuminated material emitting characteristic neutron activation gamma-rays;
a gamma-ray detector having three dimensional vertex localization;
a gamma-ray event comprising Compton interaction within said gamma-ray detector;
said gamma-ray detector producing a first three dimensional vertex localization and a second three dimensional vertex localization corresponding to two interaction locations resulting from said gamma-ray event;
said gamma-ray detector having a first deposited energy output for said gamma-ray event;
wherein said gamma-ray detector is configured to produce an event datum containing said first deposited energy output, said first three dimensional vertex localization and said second three dimensional vertex localization from said gamma-ray event;
a computer-readable storage media, wherein a data set of neutron activation emission gamma-ray energies characteristic of suspect materials are stored;
a computational processor communicatively coupled to said computer-readable storage media;
wherein said computational processor is configured to select a first event datum and a second event datum;
wherein said computational processor is configured to select a first hypothetical emission gamma-ray energy for said first event datum and to select a second hypothetical emission gamma-ray energy for said second event datum;
wherein said computational processor is configured to calculate first backprojection data using value of said first hypothetical emission gamma-ray energy and value of said first deposited energy output, said first three dimensional vertex localization and said second three dimensional vertex localization for said first event datum and to calculate second backprojection data using value of said second hypothetical emission gamma-ray energy and value of said first deposited energy output, said first three dimensional vertex localization and said second three dimensional vertex localization for said second event datum;
wherein said computational processor is configured to calculate directional correlation between said first backprojection data for said first event datum and said second backprojection data for said second event datum;
wherein said computational processor is configured to associate said directional correlation with identity of said material; and
an output from said computational processor indicating presence of specific type of material.

2. The material identification device of claim 1 wherein:
said gamma-ray detector is mobile.

3. The material identification device of claim 2 further comprising:
a GPS system physically coupled to said gamma-ray detector and communicatively coupled to said computational processor.

4. The material identification device of claim 2 further comprising:
an IMU system physically coupled to said gamma-ray detector and communicatively coupled to said computational processor.

5. The material identification device of claim 1 further comprising:
a system controller comprising a processor wherein said system controller is communicatively coupled to said gamma-ray detector.

6. The material identification device of claim 5 wherein:
said gamma-ray detector is mobile.

7. The material identification device of claim 6 further comprising:
a GPS system physically coupled to said gamma-ray detector and communicatively coupled to said system controller.

8. The material identification device of claim 6 further comprising:
an IMU system physically coupled to said gamma-ray detector and communicatively coupled to said system controller.

9. The material identification device of claim 3 further comprising:
a second gamma-ray detector having three dimensional vertex localization and communicatively coupled to said system controller.

10. The material identification device of claim 1 wherein:
said gamma-ray detectors is a NaI type scintillator.

11. The material identification device of claim 1 wherein:
said gamma-ray detectors is a CZT type solid state detector.

12. The material identification device of claim 1 wherein:
said computational processor is communicatively coupled to said gamma ray detector through said system controller.

13. The material identification device of claim 1 wherein:
said computational processor is also said system controller.

14. An explosives detection system comprising:
a neutron source;
said neutron source illuminated suspect material;
said illuminated material emitting characteristic neutron activation gamma-rays;
a plurality of gamma-ray detectors;
a first gamma-ray detector measuring interaction energy from gamma-ray interaction within said first gamma-ray detector wherein said first gamma-ray detector is configured to produce a first event datum containing interaction energy output for gamma-ray interaction within said first gamma-ray detector;
a second gamma-ray detector measuring interaction energy from gamma-ray interaction within said second gamma-ray detector wherein said second gamma-ray detector is configured to produce a second event datum containing interaction energy output for gamma-ray interaction within said second gamma-ray detector;
a computer-readable storage media wherein a data set of neutron activation emission gamma-ray energies characteristic of suspect explosives is stored;
a computational processor communicatively coupled to said computer-readable storage media;
wherein said computational processor is configured to receive said first event datum and said second event datum;
wherein said computational processor is configured to select a first hypothetical emission gamma-ray energy for said first event datum and to select a second hypothetical emission gamma-ray energy for said second event datum;
wherein said computational processor is configured to calculate first backprojection data using value of said first hypothetical emission gamma-ray energy and identification of first gamma ray detector for said first event datum and to calculate second backprojection data using value of said second hypothetical emission gamma-ray energy and identification of second gamma ray detector for said second event datum;
wherein said computational processor is configured to calculate directional correlation between said first backprojection data and said second backprojection data; and
an output from said computational processor indicating presence of explosives.

15. The explosives detection system of claim 14 further comprising:
a system controller comprising a processor wherein said system controller is communicatively coupled to said plurality of gamma-ray detectors.

16. The explosives detection system of claim 14 wherein:
at least one of said plurality of gamma-ray detectors is mobile.

17. The explosives detection system of claim 16 further comprising:
a GPS system physically coupled to said at least one of said plurality of gamma-ray detectors and communicatively coupled to said computational processor.

18. A method of material identification comprising:
Illuminating said material with neutrons from a neutron source;
said illuminated material emitting characteristic neutron activation gamma-rays;
creating event data by measuring with a detector three dimensional vertex localization data of gamma-ray events;
each of said event data having a first Compton interaction vertex and interaction energy;
each of said event data having a second interaction vertex;
transferring said event data to a computer;
a) using said computer to select a first hypothetical emission gamma-ray energy for a first event datum;
b) using said computer to select a second hypothetical emission gamma-ray energy for a second event datum;
c) calculating on said computer first backprojection data for said first event datum using said first event datum combined with said first hypothetical emission gamma-ray energy;
d) calculating on said computer second backprojection data for said second event datum using said second event datum combined with said second hypothetical emission gamma-ray energy;
e) measuring on said computer the spatial correlation of said calculated first backprojection data and second backprojection data; and
f) associating said spatial correlation with identity of said material.

19. The method of claim 18 further comprising:
storing a data set of neutron activation emission gamma-ray energies characteristic of suspect materials on a computer-readable storage media communicatively coupled to said computer.

20. The method of claim 19 wherein:
said data set of neutron activation emission gamma-ray energies are characteristic of explosives.

21. The method of claim 18 further comprising:
steps a), b), c), d), e) and f) are iterated for alternative said hypothetical emission gamma-ray energies; and
identifying most significant said spatial correlation.