US20080302967A1 - Tracking of moving radioactive sources - Google Patents
Tracking of moving radioactive sources Download PDFInfo
- Publication number
- US20080302967A1 US20080302967A1 US11/413,698 US41369806A US2008302967A1 US 20080302967 A1 US20080302967 A1 US 20080302967A1 US 41369806 A US41369806 A US 41369806A US 2008302967 A1 US2008302967 A1 US 2008302967A1
- Authority
- US
- United States
- Prior art keywords
- radiation
- source
- detection
- detector
- detectors
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000002285 radioactive effect Effects 0.000 title claims abstract description 15
- 230000005855 radiation Effects 0.000 claims abstract description 74
- 238000001514 detection method Methods 0.000 claims abstract description 57
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims description 38
- 238000012360 testing method Methods 0.000 claims description 20
- 238000012549 training Methods 0.000 claims description 18
- 238000007476 Maximum Likelihood Methods 0.000 claims description 12
- 238000013480 data collection Methods 0.000 claims description 11
- 238000012882 sequential analysis Methods 0.000 claims description 7
- 238000004891 communication Methods 0.000 claims description 5
- 238000013500 data storage Methods 0.000 claims description 2
- 238000004590 computer program Methods 0.000 claims 2
- 230000003287 optical effect Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 19
- 230000033001 locomotion Effects 0.000 abstract description 8
- 239000000463 material Substances 0.000 abstract description 5
- 230000006855 networking Effects 0.000 abstract description 2
- 239000011824 nuclear material Substances 0.000 abstract 2
- 230000006870 function Effects 0.000 description 34
- 238000005259 measurement Methods 0.000 description 19
- 238000009826 distribution Methods 0.000 description 11
- 230000035945 sensitivity Effects 0.000 description 11
- 230000008569 process Effects 0.000 description 10
- 238000012544 monitoring process Methods 0.000 description 9
- 238000003745 diagnosis Methods 0.000 description 8
- 230000004044 response Effects 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 8
- 238000013507 mapping Methods 0.000 description 7
- 230000014509 gene expression Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000004807 localization Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000005251 gamma ray Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000037361 pathway Effects 0.000 description 2
- 230000005258 radioactive decay Effects 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000003657 Likelihood-ratio test Methods 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000013101 initial test Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
-
- G01V5/26—
-
- G01V5/271—
Definitions
- This invention relates to a portable monitor for tracking a source of radiation. More specifically this invention relates to a portable highly-sensitive radioactive source tracking system.
- the first scenario involves source interception along primary transportation corridors leading into major metropolitan centers. Choke points (such as rest stops, toll collection plazas, truck stops, weigh stations, and bridges) are logical hosting points where one could take advantage of existing infrastructure and visual clutter and where reduced traffic speed improves counting statistics. Since, presumably, these sites would not be among the intended site of ultimate dispersal and since the purveyor would be unaware that he has been detected, the element of surprise would exist.
- the second scenario involves source interception at border crossings and at maritime ports of entry in the plaza where traffic first congregates before being funneled into parallel lanes already equipped with deterrent-based portal monitors.
- An early-warning signal in advance of a signal from the existing system would provide increased response time.
- the third scenario involves personnel and small package monitoring i.e. mass transit stations, convention centers, stadiums, and airports.
- the tracking system could also be integrated with mobile systems with appropriate technology.
- a source tracking system which 1) detects the introduction of a source into the monitored area, 2) tracks its position in the monitored area, and 3) identifies the type of radioisotope and the source strength.
- the present invention relates to systems and methods for 1) detecting the introduction of a source into the monitored area, 2) tracking its position in the monitored area, and 3) identifying the type of radioisotope and the source strength.
- the present invention provides these capabilities while minimizing the probability of false positives and the probability of a missed positive (false negative) even in environments where the signal-to-noise ratio may be lower than previously achievable by prior systems.
- the present invention relates to a portable highly-sensitive radioactive source tracking system.
- the system can include one or more of gamma radiation detectors, networking, power, signal processing; algorithms, and central processing and control could be temporarily installed in an area to monitor for radioactive sources, to localize the source, and to track the movement of the source.
- the source may be shielded to reduce its dose rate for personnel protection or shielded in an attempt to hide its detection.
- the detectable gamma radiation will be distributed lower in energy due to Compton scattering. This results in a lower full energy photopeak and a larger continuum of lower energy photons.
- the likelihood of a signal of comparable intensity to background and overlapping in spectrum requires a sensitive test and one that incorporates the characteristics of the local background radiation environment. This suggests the use of statistical methods that take into account the variability of both the source and the background radiation fields and uncertainties in the electronic detection hardware. The need for timely near real-time predictions argue similarly. These uncertainties ultimately limit the precision with which the source location can be identified.
- the solution must also provide a measure of the uncertainty in the position of the source.
- the invention includes a portable, highly-sensitive radioactive source tracking system.
- the system is intended for temporary installation to monitor an area for radioactive sources, to localize the source, and to track movement of the source.
- the system of the present invention could also be used at fixed installations and in fixed configurations for the routine monitoring of source movements.
- the system exploits the mathematical relationships that exist among multiple detectors deployed in a known spatial configuration, the inherent sensitivity characteristics of the detectors, the properties of radioactive sources, and the correlation among signals from multiple detectors.
- the portable system of the present invention provides three distinct functions during monitoring: alarming when a source is detected by any single detector node, tracking of the location of the source in real-time, and collection of a gamma-ray energy spectrum for performing source identification.
- the present invention allows for: 1) alarm on source presence (by detector), 2) source location, 3) source tracking, 4) source identification, 5) source strength determination.
- the present invention offers the advantage of locating and tracking a source without having a person in close proximity to the source.
- the present invention is applicable to moving detectors when the detector and position data can be transmitted in near real-time.
- FIG. 1 is a flowchart overview of the one embodiment of the present invention
- FIG. 2 is a flow chart showing the steps in one embodiment for calibrating radiation detectors with respect to energy functions
- FIG. 3 is a flow chart depicting one embodiment of a train detection algorithm function
- FIG. 4 is a flow chart depicting one embodiment of an on-line detection and tracking algorithm function
- FIG. 5 is a source-tracking system architecture in accordance with the principles of the present invention.
- FIG. 6 is a comparison of Gaussian and Poisson distributions
- FIG. 7 is a general layout for a main panel for one embodiment of the present invention.
- FIG. 8 is a screen shot of a high voltage input GUI for one embodiment of the present invention.
- FIG. 9 is a screen shot of a detector calibration GUI for one embodiment of the present invention.
- FIG. 10 is a screen shot of a detector layout specification GUI for one embodiment of the present invention.
- FIG. 11 is a screen shot of a data acquisition checkout GUI for one embodiment of the present invention.
- FIG. 12 is a screen shot of a data viewing and analysis GUI for one embodiment of the present invention.
- FIG. 13 is a screen shot of a learning venue and detector specific characteristics GUI for one embodiment of the present invention.
- FIG. 14 is a screen shot of a source detection and tracking GUI for one embodiment of the present invention.
- FIG. 15 is graph depicting a location of source in accordance with the principles of the present invention.
- the present invention relates to a portable system and methods for detection of radioactive sources and subsequent tracking.
- the principal components of the present invention are detector electronics, data collection/archiving, operator console, functional capabilities, self diagnostics, supporting signal processing, detection, and backend electronics.
- the function capabilities of the present invention include detection of source, identification of source location, identification of source type, signal strength, and site-specific configuration and learning.
- the self-diagnostics include a degraded detector diagnosis, a failed detector diagnosis, and a shielded detector diagnosis.
- the supporting signal processing includes SAT, DDFM, and cluster analysis.
- One embodiment of the present invention utilizes NaI detectors with 14 pin connectors.
- FIG. 1 depicts a flow chart for one general method of operation for a device 101 ( FIG. 5 ) of the present invention.
- the process begins with a startup step 110 .
- the radiation detectors are calibrated in the calibration step 112 .
- a check step 114 determines if the device is appropriately calibrated. If the device isn't calibrated properly, it re-enters step 112 for calibration again. If the device is calibrated properly, then it enters a training step 116 . Following training 116 , a training check 118 is performed to determine if training was accomplished. If not, then the training step 116 is repeated. However, if training was correct, then the device progresses to a online detection step 120 .
- the results of the online detection step 120 are shown in a display step 122 .
- the device alternates between detection 120 and display 122 , in another exemplary embodiment the device performs real-time or near simultaneous detection 120 and display 122 .
- FIG. 2 illustrates one embodiment of a calibration step 112 .
- the device 101 enters at entry step 124 and proceeds to an input step for selection of a detector type and spatial location.
- the voltage threshold and gain for the detectors is set at step 128 .
- a calibration collection step 130 collects count rates from a weak radiation source.
- a display step 132 the spectra for the weak radiation counted in step 130 is displayed.
- a gain and threshold check step 134 is performed following collection 130 . If the threshold and gain assignments were not completed or correct, the device 101 proceeds again through steps 126 , 128 , 130 , and 132 . If the assignments were completed, then the calibration step 112 is exited at step 136 .
- FIG. 3 illustrates one embodiment of a training step 116 .
- the device 101 enters the training step 116 at an entry step 128 .
- a background collection step is performed to collect data regarding background radiation.
- An analysis step 142 is performed to determine if the background data exhibits Gaussian distribution. If not, the background collection step 140 is performed again. If the background data does exhibit a Gaussian distribution, then an estimation step 144 is performed to estimate the decay parameter, in one embodiment using the maximum likelihood method.
- a directional response step 146 is performed. In one embodiment, a weak radiation source is introduced to obtain the detector's directional response. In another embodiment, a Monte Carlo simulation is performed to obtain the directional response. Following acquisition of the directional response in step 146 , the non-linear mapping is determined in a mapping step 148 for the direction response of step 146 . Once the non-linear mapping is learned, the device 101 exits the training step 116 at step 150 .
- FIG. 4 illustrates one embodiment of a detection step 120 .
- the device 101 enters the detection step at 152 .
- a sampling step 154 takes a sample of the environment at a predetermined interval, in an exemplary embodiment a one second interval.
- the sample is tested at step 156 for the presence of a radiation source, in an exemplary embodiment the testing utilizes sequential analysis.
- a presence determination is made at step 158 , if no radiation source is detected, then steps 154 and 156 are performed again, if a source is detected, then the device 101 proceeds to step 160 for alarm display on a GUI.
- an autoassociative mapping step 162 is performed to generate surrogate signals for shield detectors.
- a location step 164 locates the source and determines its strength, in one embodiment the maximizing likelihood function is used to make this determination.
- a display step 166 provides a graphical display of the probability or density of the environment to indicate the radiation source to a user.
- the display step 120 is exited at step 168 .
- FIG. 5 depicts one embodiment of the present invention.
- the portable detection device 101 includes data collection 210 , decision making 212 , and user interface 214 subsystems.
- the data collection subsystem 210 includes at least one radiation detector 216 , in communication with a data archiving/collection device 218 .
- the data collection subsystem 210 further includes a video device 220 .
- the decision making subsystem 212 includes an information bus 224 for transporting data from the data collection subsystem 212 to be processed.
- the processing includes in an exemplary embodiment detection of a source 226 , location of the source 228 , strength of the source 230 , failed or degraded detector diagnosis 232 , shield detector diagnosis 234 .
- the portable detection device 101 further includes an operator interface subsystem 214 .
- the operator interface subsystem 214 provides a user with information regarding the outputs of the other subsystems 210 , 212 and provides the user with the ability to control the device 101 .
- An operator console 238 is provided as a means for input by a user for controlling the device 101 .
- an auto diagnosis scheduler 242 is provided which controls the degraded/failed detector diagnosis 232 and the shielded detector diagnosis 234 for determination of the status of the at least one detector 216 .
- a function schedule 240 provides the corresponding controls for the detection of source 226 , location of source 228 , and strength of source 230 processing.
- a model for the emission of gamma rays from a radioactive source in accordance with the practice of the present invention is based upon a few basic assumptions. One of these is that the probability of a single emission per unit time is ⁇ It has been shown previously that the number of gamma rays, n, emitted in the time interval ⁇ t by a source is a random variable with Poisson distribution
- ⁇ is treated as a free parameter whose value depends on the type of source and will be estimated from measurements.
- the counts registered by a detector near the source are a subset of all the emissions.
- the detector has a view which is a fraction of the total solid angle surrounding the source and only a fraction of the gamma rays that enter the detector actually register as events. That is, the detector has efficiency less than unity.
- the probability of the detector registering an emission per unit time is then ⁇ Then the probability of an observed count rate n is from Eq. (1)
- FIG. 6 compares Eqs. (2) and (3) for different values of the average number of counts. Essentially the figure shows that if the average is greater than 20, then Eq. (3) is a very good approximation to Eq. (1).
- the count window interval ⁇ t for the present invention is preferably large enough that on average at least 20 counts register.
- the window size of the present invention can be decreased and the Gaussian distribution still be preserved.
- the introduction of a source can be detected by comparing the current count rate with the normal background count rate without a source present. A statistically significant deviation in total count rate from the historical norm for the background would suggest the presence of a new radiation source.
- the Sequential Analysis Test (SAT) algorithm provides a highly sensitive test for detecting the presence of a source and for detecting when a detector is temporarily shielded or failing. It accomplishes this by forming the residual between an actual measurement and a value based on a probabilistic model of the process. Successive deviations in the residual from those implied by the model indicate that the model is no longer valid. That is, a source has entered the venue, a detector has become shielded, or a detector is failing.
- the Data-Driven Functional Mapping (“DDFM”) algorithm provides a capability to recall an association among signals that is implicit rather than explicit. That is, there is no requirement for declaration of independent variables as in the definition of a single valued function. This capability is used to create an auto-associative mapping that can be used to sense which out of a number of variables has deviated as a result of a new process unique to that variable that was not present at the time the mapping was learned. Important instances are a temporarily shielded or failing detector.
- the Maximum Likelihood (ML) algorithm provides a capability to determine from sensor data values for unmeasured process model parameters such that resulting model is the one that most likely generated the data. This capability is used to determine the most likely position of a source given 1) detector readings and 2) a probabilistic model for the processes that give rise to detector counts from source emissions.
- the State Estimation (SE) algorithm (with the Kalman Filter as a special case) provides a capability to predict the value of an unmeasured variable through a combination of measurements and a dynamic model for the process underlying the variable.
- the algorithm includes a representation for noise in the measurements and uncertainties in the model.
- the dynamic model provides the ability to include past measurements in addition to the current measurements (i.e., temporal data) to estimate a current value for the unmeasured variable.
- source position is an unmeasured variable.
- the dynamic model can be the equations of motion for the source or constraints on its spatially rate of change. Including past count data is a way to improve estimated source position over that obtained by only using the count from the current observation window.
- SAT is based on the likelihood ratio test of Wald, A. 1947. Sequential Analysis. John Wiley & Sons, New York, herein incorporated by reference.
- SAT supposes that a random process r uncorrelated in time is normally distributed with zero mean and depends on the parameter ⁇ .
- the likelihood of observing r i given the parameter value ⁇ 1 is denoted p( ⁇ /r i ).
- the task is to decide as each new observation becomes available what is the value of ⁇ . If the choice is between two values, 61 and 62 then form the following test statistic:
- H o that the parameter value is ⁇ 1
- H 1 that the parameter value is ⁇ 2 .
- ⁇ is the probability that the test will result in the selection of H 1 when hypothesis 0 is correct and ⁇ is the probability that the test will result in the selection of H 0 when hypothesis 1 is correct.
- this ratio test is used to decide with the arrival of each new measurement sample when the model is no longer a statistically valid representation of the background radiation.
- the present invention looks specifically for a drift in measurement count away from the model, characterizing the drift in terms of the parameter ⁇ defined to be the difference between the mean of the model output and the measurement counts. Immediately following training and during subsequent normal operation its value will be zero. If, however, there is a change in the measurement counts, then ⁇ will no longer be zero. To monitor its value, a residual vector is used:
- the likelihood function for the residual vector is defined by:
- the decision as to whether or not a failure has occurred is then given by Eq. (5) where the likelihood function is given by Eq. (5).
- the probability that a decision is incorrect is given by Eq. (6).
- the SAT algorithm is used to detect a change in detector counts from normal background thereby signaling the introduction of a source.
- the values of the two parameters, ⁇ and K, which are specific to the site, are all that is needed to complete the description of background given by Eq. (8).
- these values can be estimated in one of two ways. In the maximum likelihood method, the values are those that maximize the probability that Eq. (8) generated an observed data set of background counts. Alternatively, if the number of these data is large, then the usual statistical estimators for mean and variance will yield good values in a computationally more straightforward manner than the maximum likelihood method.
- Testing refers to the actual on-line real-time monitoring of total counts to detect statistically significant deviation of the count rate from normal background to estimate the position of the source. Normal background characteristics are established during training. There is a mathematical link as described above between the sensitivity of the test and the probability that a wrong decision as to the presence of a source is made.
- a user has control of the sensitivity through setting the values of the parameters A and B in Eq. (9), but as described above, the values then fix the probability of a wrong SAT decision.
- radioactive source that is concentrated at a point can be located in space through a combination of triangulation and the knowledge that the count rate falls off as one over the distance squared.
- the count rate is a probability density function, the exact location will be subject to uncertainty.
- both the distance law and emission probability are used to derive the likelihood that the source is located at a particular position.
- One embodiment of the present invention involves a Maximum Likelihood Algorithm.
- the probability distribution function of the counts that register at a detector is a function of the detection efficiency (intrinsic detector efficiency and solid angle) and the distance from the source.
- the source strength is a factor in the overall count rate but the probability density function is normalized per unit source emission.
- Let the subscript 0 designate a reference detector that is located a distance r 0 from a source with a total detection efficiency ⁇ 0 . As before the probability of observing n counts while sampling for ⁇ t is
- the parameter ⁇ is exclusive to the radioactive decay process and has no dependence on detector location or detection efficiency.
- the number of counts is the sum of the contribution from the source and from background. It is assumed the background can be characterized by a Gaussian distribution.
- the sum of the background counts and source counts, both normally distributed random variables, is a normal random variable with probability density function
- the above expression reflects the fact that the probability of a single count being observed in unit time is proportional to the detector efficiency and inversely proportional to the square of the distance to the source.
- the probability that two detectors will see a particular set of counts can be derived from this.
- the subscripts i and j indicate the two detectors. Since a count seen at detector i is independent of a count at detector j (each count is the consequence of two separate disintegrations and the disintegrations are independent of each other) the joint probability is the product p i ( ⁇ 0 ⁇ t s , ⁇ 0 ⁇ t b , r 0 /r i , ⁇ i / ⁇ 0
- the probability density function for seeing counts n 0 , n 1 , . . . , n N-1 at detectors 0 , 1 , . . . , N ⁇ 1, respectively, is
- the location of the source is taken to be the point where the likelihood of finding the source is greatest. This corresponds to the point where the joint probability density function given by Eq. (13) has its maximum value.
- a measure of the uncertainty in the estimate for the source location is provided by the behavior of the likelihood function in the region that surrounds the most likely location.
- p(n 0 , n 1 , . . . , n N-1 , r 0 , r 1 , . . . , r N-1 ) denote the probability that the source is at r 0 , r 1 , . . . , r N-1 and let the point where p has maximum value, the most likely location of the source, be designated by r 0 ′, r 1 ′, . . . , r N-1 ′.
- This function plotted as contours of constant value over two-dimensional space provides a graphical image of how well the source can be located in space. Closely spaced contour lies would indicate the likelihood that the source being at a particular point falls off rapidly as one moves away from the most likely point.
- ⁇ is the parameter with over whose space the optimization is to be performed and ⁇ is the mean and ⁇ 2 is the variance.
- ⁇ is the parameter with over whose space the optimization is to be performed and ⁇ is the mean and ⁇ 2 is the variance.
- the value of ⁇ at the maximum of a product of such terms is the same as the value at the minimum of the negative of the natural logarithm. That is
- the likelihood that a source at a particular position produced the observed measurements is computed.
- the likelihood is given by Eq. (13).
- There are several parameters that appear in Eq. (13), that must first be estimated: ⁇ 0 ⁇ b ⁇ t, ⁇ 0 ⁇ s , ⁇ t, and ⁇ i / ⁇ 0 , i 1, N ⁇ 1
- Equation (13) assumes a single count observation window from each detector. This extends to M contiguous observation windows for each detector. Let the counts obtained for detector 0 in the successive windows 0 , 1 , 2 , . . . , M ⁇ 1 be denoted by n 00 , n 01 , . . . , n 0M-1 , respectively. Since the total count in a window is independent of that in any other window, the probability of getting a particular sequence of window counts is the product of the probabilities for the count in each of the windows. Thus, the likelihood that our model produced the data is
- ⁇ i / ⁇ 0 is estimated by maximizing with respect to ⁇ i / ⁇ 0 the likelihood our model produced the detector data.
- the assumption is that the background count rate registered by each detector differs from the others only as a result of differing total detection efficiency. These differences can be due to a variety of factors including the individual detector characteristics, the detector sizes and shapes, the lower level discriminator settings, and the detector calibrations. The likelihood is then
- n ik ( 18 )
- n ik 19
- the value for ⁇ 0 ⁇ b ⁇ t is the value obtained in Step 1 and the value for ⁇ i / ⁇ 0 is the value obtained in Step 2 .
- the likelihood must be maximized with respect to not only ⁇ 0 ⁇ s ⁇ t but also r 0 /r i .
- Testing refers to the actual on-line real-time monitoring of counts to estimate the position of the source.
- initial data acquisition, collection, and storage relies on the MCA, such as that sold by ORTEC under the trade name digiBase digital MCA, coupled to NaI detectors.
- Data acquisition software for the command, control, data collection, and storage may be performed using the LabView application. This includes functions for displaying the spectra, changing detector parameters (like gain and HV), and even energy calibrating the detectors (See Figures and discussion below regarding GUI).
- data acquisition is designed for a system of up to 6 detectors.
- Post-processing tasks and data transmission times associated with the list mode collection can cause a practical limit on the count rate that is unacceptable. So, in one embodiment, the data collection allows the system to collect a spectrum for each detector for a preselected time interval. In one exemplary embodiment, one second intervals are used during all of the initial tests, but it can be reduced to a practical limit of about 100 milliseconds or increased.
- a Sequential Analysis Test is used to detect a source.
- SAT is a statistical test for consistency between a time-sampled signal and a reference Gaussian distribution of known mean and variance. With each new sample, the SAT algorithm issues good, bad, or need more data decision.
- a Training Phase is used in one embodiment to learn the mean and variance of reference distribution. Behind the main window are tabs for all of the system functions like detector and system set-up, detector calibration, detector spectral display, and parameter and display controls. Some of these windows were shown in the previous figures.
- the present invention includes a method and apparatus for improving performance for “low” signal strength cases.
- the present invention overcomes poor signal-to-noise ratio, and the degeneracy it leads to, by introducing physically-based constraints on the solution.
- the laws of motion are used to constrain predicted position given detector data across successive time frames.
- source strength parameter is estimated using data across multiple time frames rather than a single time frame as is presently done.
- the present invention provides for the consideration of the shielding effect from people.
- the tracking method of the present invention is extended to account for distributed shielding materials.
- the shielding of the source by individuals passing between the source and detector can be thought as falling into three classes, each amenable to solution: point shield, uniformly distributed shielding, and spatially skewed shielding.
- point shield a single person shielding a single detector can be detected as a loss of correlation among detector signals.
- a surrogate can be generated using the DDFM algorithm for the signal from the shielded detector.
- uniformly distributed shielding the presence of uniformly distributed shielding material introduces an additional parameter whose value can be estimated if there are four or more detectors.
- a uniform distribution would approximate the case where the monitored area has a constant flux of individuals such that the collective shielding effect can be considered relatively constant.
- the population density in the monitored space is not reasonably uniform.
- other techniques can be employed (such as video image diagnostics) to determine spatially-dependent shielding coefficients.
- the present invention accounts for the efficiency of numerical algorithms. Avenues which increase data processing speeds of the signal processing algorithms are utilized in one embodiment of the present invention. Real-time implementation requires rapid execution of the signal processing algorithms.
- the present invention utilizes a detector count probability density function which is a function of efficiency and distance:
- localization is derived using a Multiple Detector Probability Density Function (MDPDF)
- MDPDF Multiple Detector Probability Density Function
- n k
- FIG. 15 is an graphical depiction of the likelihood contour of a source located within a boundary having detectors at the four corners.
- Detector response in general is non-uniform over a terrain map.
- fixed obstructions block view and the detector response function has a spatial/directional dependence.
- DDFM is used to learn these response relations for a specific site.
- the GUI provides a display of the present invention's results.
- the GUI is being comprises a laptop computer using LabViewTM.
- the user interface is set up to contain one high level window which contains the detector readings and alarm indicators, the PDF map of the source location, video display of the monitored area, and control buttons for system command and configuration.
- the present invention relates to a method for testing an area for a radiation source.
- the detector is setup and calibrated.
- the user inputs the high voltage threshold in the FIG. 8 GUI 403 .
- the user sets the energy to channel number correspondence via the GUI 405 in FIG. 9 using a known source.
- the detector layout is specified.
- the detector configuration is provided to the system as a collection of points on a two-dimensional grid.
- the GUI 407 for inputting this data is shown in FIG. 10 .
- data acquisition occurs.
- the basic operability of the data acquisition task is assessed by viewing data from individual detectors.
- the GUI 409 in FIG. 11 provides a capability for selecting individual detectors, controlling the time window over which they are sampled, and then viewing the spectral content of that detector. When used in conjunction with a known source, the operability can be confirmed.
- the strength-distance detection threshold is determined.
- the SAT alarming capability is used to establish the size of the Detection Perimeter (DP) that surrounds the detectors for which a source of Specified Minimum Strength (SMS) will be detected.
- the detection sensitivity will be a function of parameters that control the probability of obtaining a false positive and missing a true positive. These parameters will in general not be operator settable.
- the SMS will be a function of the types of perceived threats and in general will not be operator settable.
- the perimeter may be established by walking a source with strength equal to the SMS about the venue, entering the locations where the SAT alarms into the venue grid and building up a map of the perimeter for display in this manner.
- One specific embodiment of the general panel layout of FIG. 7 is shown as the GUI 401 of FIG. 14 .
- the present invention is provided with Connectivity to computer through USB port(s).
- a computer e.g., a laptop, hard drive used for data collection either as permanent or temporary storage on hard drive or for download to central server or CD-ROM.
- One embodiment of the present invention includes a graphical user interface. While various functionality of the such a GUI would be understood by one skilled in the art, such functionality can include the following. For system set-up, venue specific information, detector information, calibration, surveillance and video. For monitor system status and display, functionality includes alarm states, detector states, and a video feed. Reporting functionality may also be provided on the GUI.
- the method of learning the background statistics comprises, in one embodiment, dividing time into contiguous windows, accumulating the counts in each window, representing statistics of the background by the PDF formed from the count rate in each window over all windows.
- the background is assumed to be Gaussian with Constant ⁇ and ⁇ .
- SAT is used to detect introduction of a new source, wherein the present invention looks for a sequence of time windows exhibiting count rates inconsistent with PDF of background.
- the time window may be designed in accordance with the principles of the present invention to be of various lengths. The choice of the time windowing is based on many factors including: the number of detectors and sizes, the minimum detectable level, the uncertainty on location (speed of target), and the population density
- gamma rays are the most readily detectable radioactive decay products. They are energetic, have large interaction cross sections with detector materials, and require relatively large shielding masses to hinder chance of detection.
- the energy of gamma rays of interest from an unshielded source is typically above 100 keV up to 3 MeV.
- Background gamma ray energies are concentrated predominantly at low energy, typically below about 200 keV, with a few weak peaks at higher energies due to known isotopes, such as 40 K. Background radiation tends to be at different energies than full energy photopeaks from source emissions—although this is a generic statement that does not apply to all sources.
- the different spectral characteristics of source emissions compared to background emissions can be used to increase detector sensitivity.
- shielding of a source introduces Compton scattering that can downshift the spectrum emanating from the source.
- the signal-to-noise ratio for a detector signal can be significantly improved by only considering a small energy band around the full energy photopeak energy.
- the real challenge for detection is to perform source detection and location tracking in situations where the source strength as seen at the detector is not much above background. It follows that one should take the statistics of the source and uncertainties in the detectors into account to get best detection sensitivity. The longer the sample, the more counts obtained and the better characterized of the radiation environment. This applies to both background radiation signal and source radiation signal. In the case of background, it is possible to sample for an extended time period prior to the equipment going into service to obtain a well characterized background. But in the case of in-service operation, the source may be moving, thus requiring a constant update on its estimated position, and so this is not given to accumulating good statistics if the signal is just above background. Also, because the physics of point sources is dominated by the 1/r 2 effect, the ability of the present invention to localize the source based on smaller samples allows for a larger detection range, lower minimum detectable sources, and improved tracking for faster moving sources.
- the present invention allows for pinpointing and tracking of the location of the source in real-time. This differs from existing systems which perform a crude localization function based on which detectors are alarming.
- an operator interface console runs on laptop for ease of deployment.
- the operator interface supports initial setup including system calibration, data acquisition for training of signal processing algorithms, and viewing results of source detection and tracking tasks.
- detectors can be self-powered and wireless. Rechargeable battery provides for 24 hours continuous operation. Detectors can also be operated using Cabling (USB and Cat-5) and powered from standard AC power through USB hubs.
- the present invention is a low profile system for unobtrusive monitoring.
- the number of detectors is user selectable, subject to sensitivity requirements.
- confidence in predicted source location is reflected in likelihood contours.
- the present invention includes the learning of site specific characteristics including background radiation level, detector spatial arrangement, and detector response function including effect of shielding from permanent obstructions for increased overall system sensitivity and reliability.
- the predicted location of source is rendered as ghost-image that is superimposed on video camera image of venue. Operator uses composite image as an aid for visual identification of purveyor of source.
- the present invention offers the advantage of locating and tracking a source without having a person in close proximity to the source.
- the present invention is applicable to moving detectors when the detector and position data can be transmitted in near real-time.
- fixed and mobile detectors are integrated into same display.
- the system provides three distinct functions during monitoring: alarming when a source is detected by any single detector node, tracking of the location of the source in real-time, and collection of a gamma-ray energy spectrum for performing source identification.
Abstract
Description
- The United States Government has rights in this invention pursuant to Contract No. W-31-109-ENG-38 between the United States Government and Argonne National Laboratory.
- This invention relates to a portable monitor for tracking a source of radiation. More specifically this invention relates to a portable highly-sensitive radioactive source tracking system.
- Current devices for radioactive source detection exist but their capabilities do not include a precise identification of the position of the source nor its movement. Data collection is independent, not correlated, and not stored. There are currently no systems which integrate signals from multiple sensors or which provide directional information or discrete source location information. In addition, prior systems allow a target to know that they are being tracked.
- These prior art devices tend to be limited with respect to the range of geometrical configurations in which they can be deployed. Generally, they must be deployed in a single fixed geometrical arrangement of detectors and thus are limited to applications where the flow of traffic can be directed to meet the spatial requirements imposed by the detection device. This generally involves defining a surveillance pathway and having to alter the free flow of traffic so it is directed into the pathway. However, there is a need for a device capable of monitoring applications where the venue is defined by pre-existing boundaries such as walls or other barriers.
- In particular, there are three deployment scenarios for which state-of-the-art methods are insufficient and can be significantly improved. The first scenario involves source interception along primary transportation corridors leading into major metropolitan centers. Choke points (such as rest stops, toll collection plazas, truck stops, weigh stations, and bridges) are logical hosting points where one could take advantage of existing infrastructure and visual clutter and where reduced traffic speed improves counting statistics. Since, presumably, these sites would not be among the intended site of ultimate dispersal and since the purveyor would be unaware that he has been detected, the element of surprise would exist. The second scenario involves source interception at border crossings and at maritime ports of entry in the plaza where traffic first congregates before being funneled into parallel lanes already equipped with deterrent-based portal monitors. An early-warning signal in advance of a signal from the existing system would provide increased response time. The third scenario involves personnel and small package monitoring i.e. mass transit stations, convention centers, stadiums, and airports. The tracking system could also be integrated with mobile systems with appropriate technology.
- Thus, there is a need for a source tracking system which 1) detects the introduction of a source into the monitored area, 2) tracks its position in the monitored area, and 3) identifies the type of radioisotope and the source strength. These capabilities need to be provided in a way that minimizes the probability of false positives and the probability of a missed positive (false negative). In addition, for reasons related to equipment portability and cost, this must be achieved in environments where the signal-to-noise ratio may be lower than previously achievable.
- The present invention relates to systems and methods for 1) detecting the introduction of a source into the monitored area, 2) tracking its position in the monitored area, and 3) identifying the type of radioisotope and the source strength. The present invention provides these capabilities while minimizing the probability of false positives and the probability of a missed positive (false negative) even in environments where the signal-to-noise ratio may be lower than previously achievable by prior systems.
- In one embodiment, the present invention relates to a portable highly-sensitive radioactive source tracking system. In the scenario envisioned, the system can include one or more of gamma radiation detectors, networking, power, signal processing; algorithms, and central processing and control could be temporarily installed in an area to monitor for radioactive sources, to localize the source, and to track the movement of the source.
- A solution to this problem must take account of a number of factors. The source may be shielded to reduce its dose rate for personnel protection or shielded in an attempt to hide its detection. With shielding materials present, the detectable gamma radiation will be distributed lower in energy due to Compton scattering. This results in a lower full energy photopeak and a larger continuum of lower energy photons. The likelihood of a signal of comparable intensity to background and overlapping in spectrum requires a sensitive test and one that incorporates the characteristics of the local background radiation environment. This suggests the use of statistical methods that take into account the variability of both the source and the background radiation fields and uncertainties in the electronic detection hardware. The need for timely near real-time predictions argue similarly. These uncertainties ultimately limit the precision with which the source location can be identified. The solution must also provide a measure of the uncertainty in the position of the source.
- In one exemplary embodiment, the invention includes a portable, highly-sensitive radioactive source tracking system. In one embodiment, the system is intended for temporary installation to monitor an area for radioactive sources, to localize the source, and to track movement of the source. In another embodiment, the system of the present invention could also be used at fixed installations and in fixed configurations for the routine monitoring of source movements. The system exploits the mathematical relationships that exist among multiple detectors deployed in a known spatial configuration, the inherent sensitivity characteristics of the detectors, the properties of radioactive sources, and the correlation among signals from multiple detectors.
- In one exemplary embodiment, the portable system of the present invention provides three distinct functions during monitoring: alarming when a source is detected by any single detector node, tracking of the location of the source in real-time, and collection of a gamma-ray energy spectrum for performing source identification.
- The present invention allows for: 1) alarm on source presence (by detector), 2) source location, 3) source tracking, 4) source identification, 5) source strength determination. In one embodiment, the present invention offers the advantage of locating and tracking a source without having a person in close proximity to the source. In addition, the present invention is applicable to moving detectors when the detector and position data can be transmitted in near real-time.
- These and other objects, advantages, and features of the invention, together with the organization and manner of operation thereof, will become apparent from the following detailed description when taken in conjunction with the accompanying drawings, wherein like elements have like numerals throughout the several drawings described below.
-
FIG. 1 is a flowchart overview of the one embodiment of the present invention; -
FIG. 2 is a flow chart showing the steps in one embodiment for calibrating radiation detectors with respect to energy functions; -
FIG. 3 is a flow chart depicting one embodiment of a train detection algorithm function; -
FIG. 4 is a flow chart depicting one embodiment of an on-line detection and tracking algorithm function; -
FIG. 5 is a source-tracking system architecture in accordance with the principles of the present invention; -
FIG. 6 is a comparison of Gaussian and Poisson distributions; -
FIG. 7 is a general layout for a main panel for one embodiment of the present invention; -
FIG. 8 is a screen shot of a high voltage input GUI for one embodiment of the present invention; -
FIG. 9 is a screen shot of a detector calibration GUI for one embodiment of the present invention; -
FIG. 10 is a screen shot of a detector layout specification GUI for one embodiment of the present invention; -
FIG. 11 is a screen shot of a data acquisition checkout GUI for one embodiment of the present invention; -
FIG. 12 is a screen shot of a data viewing and analysis GUI for one embodiment of the present invention; -
FIG. 13 is a screen shot of a learning venue and detector specific characteristics GUI for one embodiment of the present invention; -
FIG. 14 is a screen shot of a source detection and tracking GUI for one embodiment of the present invention; -
FIG. 15 is graph depicting a location of source in accordance with the principles of the present invention. - The present invention relates to a portable system and methods for detection of radioactive sources and subsequent tracking. In general, the principal components of the present invention are detector electronics, data collection/archiving, operator console, functional capabilities, self diagnostics, supporting signal processing, detection, and backend electronics.
- In one embodiment, the function capabilities of the present invention include detection of source, identification of source location, identification of source type, signal strength, and site-specific configuration and learning. In one embodiment, the self-diagnostics include a degraded detector diagnosis, a failed detector diagnosis, and a shielded detector diagnosis. In addition, in on embodiment, the supporting signal processing includes SAT, DDFM, and cluster analysis. One embodiment of the present invention utilizes NaI detectors with 14 pin connectors. One of ordinary skill in the art will appreciate that a variety of different sized detectors and different types could be used as known in the art.
-
FIG. 1 depicts a flow chart for one general method of operation for a device 101 (FIG. 5 ) of the present invention. The process begins with astartup step 110. The radiation detectors are calibrated in thecalibration step 112. Acheck step 114 determines if the device is appropriately calibrated. If the device isn't calibrated properly, it re-enters step 112 for calibration again. If the device is calibrated properly, then it enters atraining step 116. Followingtraining 116, atraining check 118 is performed to determine if training was accomplished. If not, then thetraining step 116 is repeated. However, if training was correct, then the device progresses to aonline detection step 120. In one embodiment, the results of theonline detection step 120 are shown in adisplay step 122. In an exemplary embodiment, the device alternates betweendetection 120 anddisplay 122, in another exemplary embodiment the device performs real-time or nearsimultaneous detection 120 anddisplay 122. -
FIG. 2 illustrates one embodiment of acalibration step 112. Thedevice 101 enters atentry step 124 and proceeds to an input step for selection of a detector type and spatial location. Next the voltage threshold and gain for the detectors is set atstep 128. Acalibration collection step 130 collects count rates from a weak radiation source. Next adisplay step 132, the spectra for the weak radiation counted instep 130 is displayed. A gain andthreshold check step 134 is performed followingcollection 130. If the threshold and gain assignments were not completed or correct, thedevice 101 proceeds again throughsteps calibration step 112 is exited atstep 136. -
FIG. 3 illustrates one embodiment of atraining step 116. Thedevice 101 enters thetraining step 116 at anentry step 128. A background collection step is performed to collect data regarding background radiation. Ananalysis step 142 is performed to determine if the background data exhibits Gaussian distribution. If not, thebackground collection step 140 is performed again. If the background data does exhibit a Gaussian distribution, then anestimation step 144 is performed to estimate the decay parameter, in one embodiment using the maximum likelihood method. Next, adirectional response step 146 is performed. In one embodiment, a weak radiation source is introduced to obtain the detector's directional response. In another embodiment, a Monte Carlo simulation is performed to obtain the directional response. Following acquisition of the directional response instep 146, the non-linear mapping is determined in amapping step 148 for the direction response ofstep 146. Once the non-linear mapping is learned, thedevice 101 exits thetraining step 116 atstep 150. -
FIG. 4 illustrates one embodiment of adetection step 120. Thedevice 101 enters the detection step at 152. Asampling step 154 takes a sample of the environment at a predetermined interval, in an exemplary embodiment a one second interval. Next the sample is tested atstep 156 for the presence of a radiation source, in an exemplary embodiment the testing utilizes sequential analysis. A presence determination is made atstep 158, if no radiation source is detected, then steps 154 and 156 are performed again, if a source is detected, then thedevice 101 proceeds to step 160 for alarm display on a GUI. In one embodiment, followingalarm display 160, anautoassociative mapping step 162 is performed to generate surrogate signals for shield detectors. Alocation step 164 locates the source and determines its strength, in one embodiment the maximizing likelihood function is used to make this determination. Following the determination in thelocation step 164, adisplay step 166 provides a graphical display of the probability or density of the environment to indicate the radiation source to a user. Thedisplay step 120 is exited atstep 168. -
FIG. 5 depicts one embodiment of the present invention. Theportable detection device 101 includesdata collection 210, decision making 212, anduser interface 214 subsystems. Thedata collection subsystem 210 includes at least oneradiation detector 216, in communication with a data archiving/collection device 218. In an exemplary embodiment, thedata collection subsystem 210 further includes avideo device 220. Thedecision making subsystem 212 includes aninformation bus 224 for transporting data from thedata collection subsystem 212 to be processed. In one embodiment, the processing includes in an exemplary embodiment detection of asource 226, location of thesource 228, strength of thesource 230, failed or degradeddetector diagnosis 232,shield detector diagnosis 234. In one exemplary embodiment, the data from thedata collection subsystem 210 may be processed to account for sitespecific learning 222. Theportable detection device 101 further includes anoperator interface subsystem 214. Theoperator interface subsystem 214 provides a user with information regarding the outputs of theother subsystems device 101. Anoperator console 238 is provided as a means for input by a user for controlling thedevice 101. In one embodiment, anauto diagnosis scheduler 242 is provided which controls the degraded/faileddetector diagnosis 232 and the shieldeddetector diagnosis 234 for determination of the status of the at least onedetector 216. In addition, afunction schedule 240 provides the corresponding controls for the detection ofsource 226, location ofsource 228, and strength ofsource 230 processing. - A model for the emission of gamma rays from a radioactive source in accordance with the practice of the present invention is based upon a few basic assumptions. One of these is that the probability of a single emission per unit time is λ It has been shown previously that the number of gamma rays, n, emitted in the time interval Δt by a source is a random variable with Poisson distribution
-
- In the present invention, λ is treated as a free parameter whose value depends on the type of source and will be estimated from measurements.
- Thus, the counts registered by a detector near the source are a subset of all the emissions. The detector has a view which is a fraction of the total solid angle surrounding the source and only a fraction of the gamma rays that enter the detector actually register as events. That is, the detector has efficiency less than unity. Thus, for a given detector in a given position relative to a source only a fraction η of the total emissions from the source register as events in the detector. The probability of the detector registering an emission per unit time is then ηλ Then the probability of an observed count rate n is from Eq. (1)
-
- When the average number of counts, ηλΔt, is large the Poisson distribution is well approximated by the Gaussian distribution
-
-
FIG. 6 compares Eqs. (2) and (3) for different values of the average number of counts. Essentially the figure shows that if the average is greater than 20, then Eq. (3) is a very good approximation to Eq. (1). Thus, the count window interval Δt for the present invention is preferably large enough that on average at least 20 counts register. Similarly, as the count rate increases due to the introduction of a source, the window size of the present invention can be decreased and the Gaussian distribution still be preserved. - In one embodiment of the present invention, the introduction of a source can be detected by comparing the current count rate with the normal background count rate without a source present. A statistically significant deviation in total count rate from the historical norm for the background would suggest the presence of a new radiation source.
- The Sequential Analysis Test (SAT) algorithm provides a highly sensitive test for detecting the presence of a source and for detecting when a detector is temporarily shielded or failing. It accomplishes this by forming the residual between an actual measurement and a value based on a probabilistic model of the process. Successive deviations in the residual from those implied by the model indicate that the model is no longer valid. That is, a source has entered the venue, a detector has become shielded, or a detector is failing.
- The Data-Driven Functional Mapping (“DDFM”) algorithm provides a capability to recall an association among signals that is implicit rather than explicit. That is, there is no requirement for declaration of independent variables as in the definition of a single valued function. This capability is used to create an auto-associative mapping that can be used to sense which out of a number of variables has deviated as a result of a new process unique to that variable that was not present at the time the mapping was learned. Important instances are a temporarily shielded or failing detector.
- The Maximum Likelihood (ML) algorithm provides a capability to determine from sensor data values for unmeasured process model parameters such that resulting model is the one that most likely generated the data. This capability is used to determine the most likely position of a source given 1) detector readings and 2) a probabilistic model for the processes that give rise to detector counts from source emissions.
- The State Estimation (SE) algorithm (with the Kalman Filter as a special case) provides a capability to predict the value of an unmeasured variable through a combination of measurements and a dynamic model for the process underlying the variable. The algorithm includes a representation for noise in the measurements and uncertainties in the model. The dynamic model provides the ability to include past measurements in addition to the current measurements (i.e., temporal data) to estimate a current value for the unmeasured variable. In our work, source position is an unmeasured variable. The dynamic model can be the equations of motion for the source or constraints on its spatially rate of change. Including past count data is a way to improve estimated source position over that obtained by only using the count from the current observation window.
- The SAT is based on the likelihood ratio test of Wald, A. 1947. Sequential Analysis. John Wiley & Sons, New York, herein incorporated by reference. SAT supposes that a random process r uncorrelated in time is normally distributed with zero mean and depends on the parameter δ. The likelihood of observing ri given the parameter value δ1 is denoted p(δ/ri). The task is to decide as each new observation becomes available what is the value of δ. If the choice is between two values, 61 and 62 then form the following test statistic:
-
- based on the sequence of observations r1, r2, . . . rm. One can form two hypotheses: Ho that the parameter value is δ1, and H1 that the parameter value is δ2. The decision process is
-
If Sm#B terminate and accept H0, else (5) -
If B<Sm<A continue sampling, else -
If Sm∃A terminate and accept H1. - with the probability that the wrong decision will be made given by
-
- where α is the probability that the test will result in the selection of H1 when
hypothesis 0 is correct and β is the probability that the test will result in the selection of H0 whenhypothesis 1 is correct. - In accordance with the principles of the present invention, this ratio test is used to decide with the arrival of each new measurement sample when the model is no longer a statistically valid representation of the background radiation. The present invention looks specifically for a drift in measurement count away from the model, characterizing the drift in terms of the parameter δ defined to be the difference between the mean of the model output and the measurement counts. Immediately following training and during subsequent normal operation its value will be zero. If, however, there is a change in the measurement counts, then δ will no longer be zero. To monitor its value, a residual vector is used:
-
r i =y measi −y i (7) - defined as the difference between the measurement counts ymeas i and model output yi for input ui. Both the model and measurement counts have the same variance Ki and the mean value of the model output is μ and the mean value of the measurement counts is μ+δ. Then the likelihood function for the residual vector is defined by:
-
- In one embodiment, failure is to be declared when the mean of the measurement counts has drifted away from the model by an amount δ=δbad. Thus, for two hypotheses: H0, where δ=0; and H1 where δ=δbad, the decision as to whether or not a failure has occurred is then given by Eq. (5) where the likelihood function is given by Eq. (5). The probability that a decision is incorrect is given by Eq. (6).
- If the test for measurement count drift away from the model is to be symmetric, then a second test is needed for the case where δ=δbad. If the test is to be applied independently to each output and is to be symmetric, then a total of 2N tests is needed where N is the number of outputs. If Ki is a constant, then it is more convenient to work with the natural logarithm of Sm given by
-
- The SAT algorithm is used to detect a change in detector counts from normal background thereby signaling the introduction of a source. An assumption is that Eq. (8) with δ=0 represents the probability distribution of background counts. The values of the two parameters, μ and K, which are specific to the site, are all that is needed to complete the description of background given by Eq. (8). In accordance with the principles of the present invention, these values can be estimated in one of two ways. In the maximum likelihood method, the values are those that maximize the probability that Eq. (8) generated an observed data set of background counts. Alternatively, if the number of these data is large, then the usual statistical estimators for mean and variance will yield good values in a computationally more straightforward manner than the maximum likelihood method.
- Testing refers to the actual on-line real-time monitoring of total counts to detect statistically significant deviation of the count rate from normal background to estimate the position of the source. Normal background characteristics are established during training. There is a mathematical link as described above between the sensitivity of the test and the probability that a wrong decision as to the presence of a source is made. In one embodiment, a user has control of the sensitivity through setting the values of the parameters A and B in Eq. (9), but as described above, the values then fix the probability of a wrong SAT decision. The maximum likelihood method is used. From previously estimated values of η0λbΔt, ηi/η0, i=1, . . . , N−1, and η0λ5Δt, the most likely location of the source is located by maximizing with respect to location in two-dimensional space the likelihood Eq. (13) produced the detector data.
- In accordance with the principles of the present invention, radioactive source that is concentrated at a point can be located in space through a combination of triangulation and the knowledge that the count rate falls off as one over the distance squared. However, since the count rate is a probability density function, the exact location will be subject to uncertainty. In the method of the present invention described below both the distance law and emission probability are used to derive the likelihood that the source is located at a particular position.
- One embodiment of the present invention involves a Maximum Likelihood Algorithm. The probability distribution function of the counts that register at a detector is a function of the detection efficiency (intrinsic detector efficiency and solid angle) and the distance from the source. The source strength is a factor in the overall count rate but the probability density function is normalized per unit source emission. Let the
subscript 0 designate a reference detector that is located a distance r0 from a source with a total detection efficiency η0. As before the probability of observing n counts while sampling for Δt is -
- The parameter λ is exclusive to the radioactive decay process and has no dependence on detector location or detection efficiency.
- In the presence of background radiation the number of counts is the sum of the contribution from the source and from background. It is assumed the background can be characterized by a Gaussian distribution. The sum of the background counts and source counts, both normally distributed random variables, is a normal random variable with probability density function
-
- where subscripts s and b denote source and background, respectively.
- Now consider a detector i with efficiency ηi at a distance ri from the same source. The probability distribution function of the counts registered by the detector is
-
- The above expression reflects the fact that the probability of a single count being observed in unit time is proportional to the detector efficiency and inversely proportional to the square of the distance to the source.
- The probability that two detectors will see a particular set of counts can be derived from this. Let the subscripts i and j indicate the two detectors. Since a count seen at detector i is independent of a count at detector j (each count is the consequence of two separate disintegrations and the disintegrations are independent of each other) the joint probability is the product pi(η0λΔts, η0λΔtb, r0/ri, ηi/η0|ni) pj(η0λΔts, η0λΔtb, r0/rj, ηj|η0). Thus, for N detectors, the probability density function for seeing counts n0, n1, . . . , nN-1 at
detectors -
- where dropped from the left-hand side is the explicit reference to detector and source parameters. If the source is not in motion, then the above expression can be extended to take in multiple sample intervals.
- The location of the source is taken to be the point where the likelihood of finding the source is greatest. This corresponds to the point where the joint probability density function given by Eq. (13) has its maximum value.
- In one embodiment, a measure of the uncertainty in the estimate for the source location is provided by the behavior of the likelihood function in the region that surrounds the most likely location. Let p(n0, n1, . . . , nN-1, r0, r1, . . . , rN-1) denote the probability that the source is at r0, r1, . . . , rN-1 and let the point where p has maximum value, the most likely location of the source, be designated by r0′, r1′, . . . , rN-1′. Form the likelihood ratio
-
- This function plotted as contours of constant value over two-dimensional space provides a graphical image of how well the source can be located in space. Closely spaced contour lies would indicate the likelihood that the source being at a particular point falls off rapidly as one moves away from the most likely point.
- In the case where the likelihood function is to be maximized it is customary to operate on the negative of the natural logarithm. The resulting expression is simpler and the small numbers associated with a negative exponential are avoided. All the expressions for probability given above consist of products of the general form
-
- where γ is the parameter with over whose space the optimization is to be performed and μ is the mean and σ2 is the variance. The value of γ at the maximum of a product of such terms is the same as the value at the minimum of the negative of the natural logarithm. That is
-
- In one embodiment, the likelihood that a source at a particular position produced the observed measurements is computed. The likelihood is given by Eq. (13). There are several parameters that appear in Eq. (13), that must first be estimated: η0λbΔt, η0λs,Δt, and ηi/η0, i=1, N−1
- For the reference detector without the presence of source, estimate η0λbΔt by maximizing with respect to η0λbΔt the likelihood that Eq. (13) produced the detector data. For this special case, in Eq. (13), N=1 and since there is no source, η0λs=Δt=0. Equation (13) assumes a single count observation window from each detector. This extends to M contiguous observation windows for each detector. Let the counts obtained for
detector 0 in thesuccessive windows -
- where the right-hand side is evaluated according to Eq. (12)
- For the remaining detectors one-by-one without presence of a source, ηi/η0 is estimated by maximizing with respect to ηi/η0 the likelihood our model produced the detector data. The assumption is that the background count rate registered by each detector differs from the others only as a result of differing total detection efficiency. These differences can be due to a variety of factors including the individual detector characteristics, the detector sizes and shapes, the lower level discriminator settings, and the detector calibrations. The likelihood is then
-
- Note that in the above expression, the value for η0ζbΔt is the value obtained in
Step 1 - For all detectors simultaneously in presence of a source, estimate η0λsΔt by maximizing with respect to it and r0/ri (as they identify a hypothesized position of the source) the likelihood the present invention's model produced the detector data. This likelihood is
-
- Note that in the above expression, the value for η0λbΔt is the value obtained in
Step 1 and the value for ηi/η0 is the value obtained inStep 2. The likelihood must be maximized with respect to not only η0λsΔt but also r0/ri. - Testing refers to the actual on-line real-time monitoring of counts to estimate the position of the source. The maximum likelihood method is used. From previously estimated values of η0λbΔt, ηi/η0, i=1, . . . , N−1, and η0λsΔt the most likely location of the source is located by maximizing with respect to location in two-dimensional space the likelihood Eq. (13) produced the detector data.
- In one embodiment, initial data acquisition, collection, and storage relies on the MCA, such as that sold by ORTEC under the trade name digiBase digital MCA, coupled to NaI detectors. Data acquisition software for the command, control, data collection, and storage may be performed using the LabView application. This includes functions for displaying the spectra, changing detector parameters (like gain and HV), and even energy calibrating the detectors (See Figures and discussion below regarding GUI).
- In one embodiment, data acquisition is designed for a system of up to 6 detectors. Post-processing tasks and data transmission times associated with the list mode collection can cause a practical limit on the count rate that is unacceptable. So, in one embodiment, the data collection allows the system to collect a spectrum for each detector for a preselected time interval. In one exemplary embodiment, one second intervals are used during all of the initial tests, but it can be reduced to a practical limit of about 100 milliseconds or increased.
- In one embodiment, a Sequential Analysis Test (SAT) is used to detect a source. SAT is a statistical test for consistency between a time-sampled signal and a reference Gaussian distribution of known mean and variance. With each new sample, the SAT algorithm issues good, bad, or need more data decision. A Training Phase is used in one embodiment to learn the mean and variance of reference distribution. Behind the main window are tabs for all of the system functions like detector and system set-up, detector calibration, detector spectral display, and parameter and display controls. Some of these windows were shown in the previous figures.
- Rigorous Treatment of Uncertainties—A rigorous means for predicting from detector signals when the underlying physics model is becoming degenerate as a result of “low” signal strength on some detectors (i.e. signal strength is approaching background) will be developed. In essence the tracking method breaks down when noise or uncertainties are at a level where they obscure the underlying deterministic physics. This condition manifests itself as the most probable location occupying large regions rather than a single point. It can be detected through a stochastic analysis of the equations that are used to generate the likelihood function.
- In one embodiment, the present invention includes a method and apparatus for improving performance for “low” signal strength cases. The present invention overcomes poor signal-to-noise ratio, and the degeneracy it leads to, by introducing physically-based constraints on the solution. In one embodiment, the laws of motion are used to constrain predicted position given detector data across successive time frames. In another exemplary embodiment, the subdivision of energy spectrum into bands and limiting localization data analysis to only those bands where source energy is present. In yet another exemplary embodiment, source strength parameter is estimated using data across multiple time frames rather than a single time frame as is presently done.
- In one embodiment, the present invention provides for the consideration of the shielding effect from people. The tracking method of the present invention is extended to account for distributed shielding materials. The shielding of the source by individuals passing between the source and detector can be thought as falling into three classes, each amenable to solution: point shield, uniformly distributed shielding, and spatially skewed shielding. In point shield, a single person shielding a single detector can be detected as a loss of correlation among detector signals. A surrogate can be generated using the DDFM algorithm for the signal from the shielded detector. In uniformly distributed shielding the presence of uniformly distributed shielding material introduces an additional parameter whose value can be estimated if there are four or more detectors. A uniform distribution would approximate the case where the monitored area has a constant flux of individuals such that the collective shielding effect can be considered relatively constant. In spatially-skewed shielding the population density in the monitored space is not reasonably uniform. In addition, one of ordinary skill in the art would appreciate that other techniques can be employed (such as video image diagnostics) to determine spatially-dependent shielding coefficients.
- In one embodiment, the present invention accounts for the efficiency of numerical algorithms. Avenues which increase data processing speeds of the signal processing algorithms are utilized in one embodiment of the present invention. Real-time implementation requires rapid execution of the signal processing algorithms.
- In one embodiment, the present invention utilizes a detector count probability density function which is a function of efficiency and distance:
-
- In on embodiment, localization is derived using a Multiple Detector Probability Density Function (MDPDF)
-
- Likelihood Function from MDPDF provides the highest probability location and degree of uncertainty.
FIG. 15 is an graphical depiction of the likelihood contour of a source located within a boundary having detectors at the four corners. - In addition, to general location awareness, the present invention contemplates site-specific learning. Detector response in general is non-uniform over a terrain map. For example, fixed obstructions block view and the detector response function has a spatial/directional dependence. In one embodiment, DDFM is used to learn these response relations for a specific site.
- In one embodiment, the present invention provides for a shielded or failed detector. In one scenario the detector produces inconsistent or abnormal signals such as where a object or person temporarily shields a detector or where there is a detector failure. In one scenario, there is no source present. In this scenario, the present invention will poll detectors one by one for abnormal statistics compared to detector historical data. Detection with the SAT and sensor logic tree as discussed above is utilized. In addition, in one embodiment, the detector is dropped until statistics return to normal. In the scenario where a source is present, the detectors are polled one by one for a deterministic signal that is inconsistent with other detectors. In one embodiment, the inconsistency is detected among an over-determined set of detectors. In another embodiment, the detector is dropped until signal is again consistent with other detectors, wherein site learning data with time-correlated data may be used.
- The sensitivity of the tracking system is primarily a function of the number of detectors and their sizes. In one embodiment, the system comprise one of three different NaI detector combinations: four detectors each two inches by two inches, four detectors each three inches by three inch inches, and four detectors each three inches by twelve inches.
- The GUI provides a display of the present invention's results. In one embodiment, the GUI is being comprises a laptop computer using LabView™. The user interface is set up to contain one high level window which contains the detector readings and alarm indicators, the PDF map of the source location, video display of the monitored area, and control buttons for system command and configuration.
-
FIG. 7 depicts one embodiment of amain panel layout 401 on a GUI.Panel 421 displays counting level graphs for all thedetectors 216, wherein counts for each time interval are shown.Panel 422 is the alarming system for monitoring abnormal counting level where the logical output from the alarming function will be shown by LEDs.Panel 423 displays the probability density function (PDF) scale map for localizing the source. In one embodiment, room schematics and detector spots are also shown in this 2-D map.Panel 424 is a video image showing the real picture of the detection area. In one embodiment, the video focuses to where the highest probability is in the PDF. Additional panels 425-428 are provided as needed for system control buttons leading to lower level screens, including but not limited to: video camera control, hardware configuration, detector calibration, training mode, measurement control, and software display. - In one embodiment, approximately twenty minutes are needed for software in accordance with the principles of the present invention to characterize the background radiation and for it to estimate the detector efficiency ratios. On completion of this phase the system is ready to enter service and perform the source detecting and tracking functions. In service the system runs in essentially real time. Strictly, since the operator interface is updated every second (by default) with the results from new measurements, the time to collect and process new measurements is about a second. The update rate is user settable. The system can also be operated in a batch mode where data is collected over a user-defined window of time.
- In one embodiment the present invention relates to a method for testing an area for a radiation source. First, the detector is setup and calibrated. The user inputs the high voltage threshold in the
FIG. 8 GUI 403. The user sets the energy to channel number correspondence via theGUI 405 inFIG. 9 using a known source. Next the detector layout is specified. The detector configuration is provided to the system as a collection of points on a two-dimensional grid. TheGUI 407 for inputting this data is shown inFIG. 10 . Next, data acquisition occurs. The basic operability of the data acquisition task is assessed by viewing data from individual detectors. TheGUI 409 inFIG. 11 provides a capability for selecting individual detectors, controlling the time window over which they are sampled, and then viewing the spectral content of that detector. When used in conjunction with a known source, the operability can be confirmed. - Once data has been acquired, it can be viewed and analyzed. This step provides the operator with a simple capability to examine the interdependence between source strength and its distance from detector as they relate to the actual signal observed at the detector. This step is provided to allow the operator to perform some basic sensitivity studies. By moving a known source about the venue he can develop a qualitative appreciation for those factors specific to his hardware and venue that impact system sensitivity. Any factors that might adversely affect performance in a gross sense would become apparent in this step. The
GUI 411 for this step is shown inFIG. 12 . - In one embodiment, venue and detector specific characteristics are learned. There are venue and detector specific parameters that the system must learn the values of through a so-called training session. These are background strength and the ratio of detector efficiencies. Both these parameters are used in the SAT algorithm for alarming and the maximum likelihood algorithm for source localization. The learning is performed with the
GUI 413 shown inFIG. 13 . - In another exemplary embodiment, the strength-distance detection threshold is determined. In this step the SAT alarming capability is used to establish the size of the Detection Perimeter (DP) that surrounds the detectors for which a source of Specified Minimum Strength (SMS) will be detected. The detection sensitivity will be a function of parameters that control the probability of obtaining a false positive and missing a true positive. These parameters will in general not be operator settable. The SMS will be a function of the types of perceived threats and in general will not be operator settable. The perimeter may be established by walking a source with strength equal to the SMS about the venue, entering the locations where the SAT alarms into the venue grid and building up a map of the perimeter for display in this manner. One specific embodiment of the general panel layout of
FIG. 7 is shown as theGUI 401 ofFIG. 14 . - In one embodiment, the present invention is provided with Connectivity to computer through USB port(s). In one embodiment a computer, e.g., a laptop, hard drive used for data collection either as permanent or temporary storage on hard drive or for download to central server or CD-ROM.
- One embodiment of the present invention includes a graphical user interface. While various functionality of the such a GUI would be understood by one skilled in the art, such functionality can include the following. For system set-up, venue specific information, detector information, calibration, surveillance and video. For monitor system status and display, functionality includes alarm states, detector states, and a video feed. Reporting functionality may also be provided on the GUI.
- For detection of a source, the background statistics must first be learned. The method of learning the background statistics comprises, in one embodiment, dividing time into contiguous windows, accumulating the counts in each window, representing statistics of the background by the PDF formed from the count rate in each window over all windows.
- The background is assumed to be Gaussian with Constant μ and σ. Then SAT is used to detect introduction of a new source, wherein the present invention looks for a sequence of time windows exhibiting count rates inconsistent with PDF of background. The time window may be designed in accordance with the principles of the present invention to be of various lengths. The choice of the time windowing is based on many factors including: the number of detectors and sizes, the minimum detectable level, the uncertainty on location (speed of target), and the population density
- Of the possible emission types one can monitor, gamma rays are the most readily detectable radioactive decay products. They are energetic, have large interaction cross sections with detector materials, and require relatively large shielding masses to hinder chance of detection. The energy of gamma rays of interest from an unshielded source is typically above 100 keV up to 3 MeV. Background gamma ray energies are concentrated predominantly at low energy, typically below about 200 keV, with a few weak peaks at higher energies due to known isotopes, such as 40K. Background radiation tends to be at different energies than full energy photopeaks from source emissions—although this is a generic statement that does not apply to all sources. However, in accordance with the principles of the present invention the different spectral characteristics of source emissions compared to background emissions, can be used to increase detector sensitivity. As a further complication to detection, shielding of a source introduces Compton scattering that can downshift the spectrum emanating from the source. In accordance with the principles of the present invention, the signal-to-noise ratio for a detector signal can be significantly improved by only considering a small energy band around the full energy photopeak energy.
- The real challenge for detection is to perform source detection and location tracking in situations where the source strength as seen at the detector is not much above background. It follows that one should take the statistics of the source and uncertainties in the detectors into account to get best detection sensitivity. The longer the sample, the more counts obtained and the better characterized of the radiation environment. This applies to both background radiation signal and source radiation signal. In the case of background, it is possible to sample for an extended time period prior to the equipment going into service to obtain a well characterized background. But in the case of in-service operation, the source may be moving, thus requiring a constant update on its estimated position, and so this is not given to accumulating good statistics if the signal is just above background. Also, because the physics of point sources is dominated by the 1/r2 effect, the ability of the present invention to localize the source based on smaller samples allows for a larger detection range, lower minimum detectable sources, and improved tracking for faster moving sources.
- In one exemplary embodiment, the present invention allows for pinpointing and tracking of the location of the source in real-time. This differs from existing systems which perform a crude localization function based on which detectors are alarming.
- In one exemplary embodiment, an operator interface console runs on laptop for ease of deployment. The operator interface supports initial setup including system calibration, data acquisition for training of signal processing algorithms, and viewing results of source detection and tracking tasks.
- In one exemplary embodiment, detectors can be self-powered and wireless. Rechargeable battery provides for 24 hours continuous operation. Detectors can also be operated using Cabling (USB and Cat-5) and powered from standard AC power through USB hubs.
- In one exemplary embodiment, the present invention is a low profile system for unobtrusive monitoring. In one exemplary embodiment, the number of detectors is user selectable, subject to sensitivity requirements.
- In one exemplary embodiment, confidence in predicted source location is reflected in likelihood contours.
- In one exemplary embodiment, the present invention includes the learning of site specific characteristics including background radiation level, detector spatial arrangement, and detector response function including effect of shielding from permanent obstructions for increased overall system sensitivity and reliability.
- In one exemplary embodiment, the predicted location of source is rendered as ghost-image that is superimposed on video camera image of venue. Operator uses composite image as an aid for visual identification of purveyor of source.
- In one embodiment, the present invention offers the advantage of locating and tracking a source without having a person in close proximity to the source. In addition, the present invention is applicable to moving detectors when the detector and position data can be transmitted in near real-time. In one exemplary embodiment, fixed and mobile detectors are integrated into same display.
- In one embodiment, the system provides three distinct functions during monitoring: alarming when a source is detected by any single detector node, tracking of the location of the source in real-time, and collection of a gamma-ray energy spectrum for performing source identification.
- The foregoing description of embodiments of the present invention have been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the present invention to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the present invention. The embodiments were chosen and described in order to explain the principles of the present invention and its practical application to enable one skilled in the art to utilize the present invention in various embodiments, and with various modifications, as are suited to the particular use contemplated.
Claims (25)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/413,698 US7465924B1 (en) | 2006-04-28 | 2006-04-28 | Tracking of moving radioactive sources |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/413,698 US7465924B1 (en) | 2006-04-28 | 2006-04-28 | Tracking of moving radioactive sources |
Publications (2)
Publication Number | Publication Date |
---|---|
US20080302967A1 true US20080302967A1 (en) | 2008-12-11 |
US7465924B1 US7465924B1 (en) | 2008-12-16 |
Family
ID=40094993
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/413,698 Expired - Fee Related US7465924B1 (en) | 2006-04-28 | 2006-04-28 | Tracking of moving radioactive sources |
Country Status (1)
Country | Link |
---|---|
US (1) | US7465924B1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110024611A1 (en) * | 2009-07-29 | 2011-02-03 | Ut-Battelle, Llc | Calibration method for video and radiation imagers |
EP2733663A1 (en) * | 2011-07-11 | 2014-05-21 | Nec Corporation | Work assistance system, terminal, method, and program |
EP2796898A4 (en) * | 2011-12-23 | 2015-09-02 | Nuctech Co Ltd | Method and device for inspecting and positioning radioactive substance |
US20170031034A1 (en) * | 2015-07-31 | 2017-02-02 | Nucare Medical Systems, Inc. | Radiation imaging method and system |
CN113468985A (en) * | 2021-06-16 | 2021-10-01 | 北京科技大学 | Method for locking suspicious radiation source carrying personnel |
CN113589758A (en) * | 2021-07-19 | 2021-11-02 | 华中科技大学 | Numerical control machine tool working space point clustering method based on modal mass distribution |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100224788A1 (en) * | 2001-10-26 | 2010-09-09 | Innovative American Technology, Inc. | Various arrangements of radiation and fissile materials detection systems using sensor arrays in spreader bars, gantry cranes, self-propelled frame structures, and transport vehicles |
US20120153162A1 (en) * | 2001-10-26 | 2012-06-21 | Innovative American Technology Inc. | High performance straddle carrier cbrne radiation verification system |
US20100226580A1 (en) * | 2005-12-01 | 2010-09-09 | Innovative American Technology Inc. | System and method for increased gamma/neutron detection |
US8247781B2 (en) | 2005-12-01 | 2012-08-21 | Innovative American Technology, Inc. | Fabrication of a high performance neutron detector with near zero gamma cross talk |
US8330115B2 (en) * | 2005-12-01 | 2012-12-11 | Innovative American Technology, Inc. | High performance neutron detector with near zero gamma cross talk |
WO2009143131A2 (en) | 2008-05-19 | 2009-11-26 | Innovative American Technology Inc. | Mobile frame structure with passive/active sensor arrays for non-invasive identification of hazardous materials |
US8930165B2 (en) * | 2010-03-31 | 2015-01-06 | Uchicago Argonne, Llc | Source localization using multiple units of a tight-pitched detector array |
US8494808B2 (en) | 2010-05-17 | 2013-07-23 | The Johns Hopkins University | Method for optimizing parameters for detection systems |
US9212905B2 (en) * | 2010-11-10 | 2015-12-15 | Uchicago Argonne, Llc | Method and system for determining radiation shielding thickness and gamma-ray energy |
US8842881B2 (en) * | 2012-04-26 | 2014-09-23 | General Electric Company | Real-time video tracking system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5324948A (en) * | 1992-10-27 | 1994-06-28 | The United States Of America As Represented By The United States Department Of Energy | Autonomous mobile robot for radiologic surveys |
US6080989A (en) * | 1995-06-24 | 2000-06-27 | British Nuclear Fuels Plc | Apparatus and methods for detecting and/or imaging gamma radiation |
-
2006
- 2006-04-28 US US11/413,698 patent/US7465924B1/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5324948A (en) * | 1992-10-27 | 1994-06-28 | The United States Of America As Represented By The United States Department Of Energy | Autonomous mobile robot for radiologic surveys |
US6080989A (en) * | 1995-06-24 | 2000-06-27 | British Nuclear Fuels Plc | Apparatus and methods for detecting and/or imaging gamma radiation |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110024611A1 (en) * | 2009-07-29 | 2011-02-03 | Ut-Battelle, Llc | Calibration method for video and radiation imagers |
US7973276B2 (en) * | 2009-07-29 | 2011-07-05 | Ut-Battelle, Llc | Calibration method for video and radiation imagers |
EP2733663A1 (en) * | 2011-07-11 | 2014-05-21 | Nec Corporation | Work assistance system, terminal, method, and program |
EP2733663A4 (en) * | 2011-07-11 | 2015-03-11 | Nec Corp | Work assistance system, terminal, method, and program |
US9024257B2 (en) | 2011-07-11 | 2015-05-05 | Nec Corporation | Work assistance system, terminal, method and program |
EP2796898A4 (en) * | 2011-12-23 | 2015-09-02 | Nuctech Co Ltd | Method and device for inspecting and positioning radioactive substance |
US9829601B2 (en) | 2011-12-23 | 2017-11-28 | Nuctech Company Limited | Method and device for patrol inspecting and locating a radioactive substance |
GB2512547B (en) * | 2011-12-23 | 2019-04-24 | Nuctech Co Ltd | Method and device for patrol inspecting and locating a radioactive substance |
US20170031034A1 (en) * | 2015-07-31 | 2017-02-02 | Nucare Medical Systems, Inc. | Radiation imaging method and system |
US9953463B2 (en) * | 2015-07-31 | 2018-04-24 | Nucare Medical Systems, Inc. | Radiation imaging method and system |
CN113468985A (en) * | 2021-06-16 | 2021-10-01 | 北京科技大学 | Method for locking suspicious radiation source carrying personnel |
CN113589758A (en) * | 2021-07-19 | 2021-11-02 | 华中科技大学 | Numerical control machine tool working space point clustering method based on modal mass distribution |
Also Published As
Publication number | Publication date |
---|---|
US7465924B1 (en) | 2008-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7465924B1 (en) | Tracking of moving radioactive sources | |
JP5465383B2 (en) | High energy particle detection method | |
US7391028B1 (en) | Apparatus and method for detection of radiation | |
CN101833116B (en) | Special nuclear material examination by neutrons and delayed neutrons | |
US8930165B2 (en) | Source localization using multiple units of a tight-pitched detector array | |
US9212905B2 (en) | Method and system for determining radiation shielding thickness and gamma-ray energy | |
US9239404B2 (en) | Device and method for continuous monitoring of persons, vehicles, containers or packets | |
US20140224964A1 (en) | Apparatus and method for detection of radiation | |
CN101019042B (en) | Object examination by dual energy radiation scanning and delayed neutron detection | |
Bukartas et al. | A Bayesian method to localize lost gamma sources | |
Mendes et al. | Radioactive hot-spot localisation and identification using deep learning | |
US10132943B2 (en) | Spectral segmentation for optimized sensitivity and computation in advanced radiation detectors | |
CN110673191A (en) | Channel type radioactive source positioning system and method | |
US10444385B1 (en) | Shieldless detector with one-dimensional directionality | |
Wagner et al. | Tests of a Gamma Spectrometer-Neutron Counter Relationship as a Neutron Alarm Metric in Mobile Radiation Search Systems | |
US20170302759A1 (en) | Intelligent server in a system of networked sensors | |
ROBERTO et al. | SYSTEM AND METHOD FOR THE | |
WO2023017045A1 (en) | Method for identifying a moving radiation source | |
Zhao | Radiation source detection from mobile sensor networks using principal component analysis | |
MX2007010503A (en) | Apparatus and method for detection of radioactive materials |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: UNIVERSITY OF CHICAGO, ILLINOIS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KLANN, RAYMOND;PARK, YOUNG SOO;VILIM, RICHARD B.;REEL/FRAME:017825/0217 Effective date: 20060428 |
|
AS | Assignment |
Owner name: U CHICAGO ARGONNE LLC, ILLINOIS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:UNIVERSITY OF CHICAGO, THE;REEL/FRAME:018385/0618 Effective date: 20060925 Owner name: U CHICAGO ARGONNE LLC,ILLINOIS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:UNIVERSITY OF CHICAGO, THE;REEL/FRAME:018385/0618 Effective date: 20060925 |
|
AS | Assignment |
Owner name: ENERGY, UNITED STATES DEPARTMENT OF, DISTRICT OF C Free format text: CONFIRMATORY LICENSE;ASSIGNOR:UNIVERSITY OF CHICAGO, THE;REEL/FRAME:018482/0123 Effective date: 20060915 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20201216 |