Gamma-ray imaging attempts to reconstruct the spatial and intensity distribution of gamma-emitting radionuclides from a set of measurements. Generally, this problem is solved by discretizing the spatial dimensions and employing the maximum likelihood expectation maximization (ML-EM) algorithm, with or without some form of regularization. While the generality of this formulation enables use in a wide variety of scenarios, it is susceptible to overfitting, limited by the discretization of spatial coordinates, and can be computationally expensive. We present a novel approach to 3D gamma-ray image reconstruction for scenarios where sparsity may be assumed, for example, radiological source search. In this paper, we first formulate a point-source localization (PSL) approach as an optimization problem, where both position and source intensity are continuous variables. We then extend and generalize this formulation to an iterative algorithm, called additive PSL (APSL), for sparse parametric image reconstruction. A set of simulated source search scenarios using a single non-directional detector are considered, finding improved image accuracy and computational efficiency with APSL over traditional grid-based approaches.

Snowpack is a vital component of Earth's hydrological cycle and one of the most sensitive to global warming. In order to develop a predictive understanding of the hydrological and biogeochemical dynamics in snow-dominated watersheds, scientists and water resources managers need spatially and tempo- rally dense data sets of snowpack parameters, which could be obtained by networks of low-cost, low-power, distributed sensors. We investigate how 60 GHz frequency modulated continuous wave (FMCW) radar systems-on-chips (SoCs) can be used for measuring a variety of snowpack parameters. For snow depth measurements, we present a dedicated radar detection algorithm that detects the top of the observed snowpack, and we compare its performance to the established cell averaging constant false alarm rate (CA-CFAR) technique. We also evaluate the impact of non-coherent integration of radar frames over receive channels and over time. For measuring snow density and snow water equivalent (SWE), we rely on an adjustment for the signal's propagation speed in the snowpack based on prior knowledge about the radar's true distance to the ground. Our lab and field experiments show that snow layers can be detected, bulk snow density can be calculated, and the 90th percentile of snow depth measurement errors is 25 mm.

Static and mobile radiation detectors can be deployed in urban environments for a range of nuclear security applications, including radiological source search-and-tracking scenarios. Modeling detector performance for such applications is challenging, as it does not depend solely on the detector capabilities themselves. Many factors must be taken into consideration, including specific source and background signatures, the topology and constraints of the deployment environment, the presence of nuisance sources, and whether detectors are mobile or static. When considering the simultaneous deployment of multiple, heterogeneous detectors, assessment of the system-wide performance requires the simulation of the individual detectors, and a system-level analysis of the detection performance. In radiological source search-and-tracking scenarios, performance is mostly dominated by the probability of encounter, which depends on the specifics of a given deployment, e.g., static vs. mobile detectors or a combination of both modalities, the number of detectors deployed, the dynamic vs. static setting of false alarm rates, and individual vs. networked operation. The Urban Deployment Model (UDM) toolset was specifically developed to cover the gap in the available generic frameworks for the simulation of radiation detector deployments at city scales. UDM provides a unified and modular framework to support the simulation and performance characterization of heterogeneous detector deployments in urban environments. This paper presents the key components along the UDM workflow.

We propose a complete framework for Bayesian image reconstruction and uncertainty quantification based on a Gaussian process prior (GPP) to overcome limitations of maximum likelihood expectation maximization (ML-EM) image reconstruction algorithm. The prior distribution is constructed with a zero-mean Gaussian process (GP) with a choice of a covariance function, and a link function is used to map the Gaussian process to an image. Unlike many other maximum a posteriori approaches, our method offers highly interpretable hyperparamters that are selected automatically with the empirical Bayes method. Furthermore, the GP covariance function can be modified to incorporate a priori structural priors, enabling multi-modality imaging or contextual data fusion. Lastly, we illustrate that our approach lends itself to Bayesian uncertainty quantification techniques, such as the preconditioned Crank-Nicolson method and the Laplace approximation. The proposed framework is general and can be employed in most radiation image reconstruction problems, and we demonstrate it with simulated free-moving single detector radiation source imaging scenarios. We compare the reconstruction results from GPP and ML-EM, and show that the proposed method can significantly improve the image quality over ML-EM, all the while providing greater understanding of the source distribution via the uncertainty quantification capability. Furthermore, significant improvement of the image quality by incorporating a structural prior is illustrated.

We present an experimental demonstration of additive point source localization (APSL), a sparse parametric imaging algorithm that reconstructs the 3-D positions and activities of multiple gamma-ray point sources. Using a handheld gamma-ray detector array and up to four 8 mu Ci 137Cs gamma-ray sources, we performed both source-search and source-separation experiments in an indoor laboratory environment. In the majority of the source-search measurements, APSL reconstructed the correct number of sources with position accuracies of 20 cm and activity accuracies (unsigned) of 20%, given measurement times of 2 to 3 min and distances of closest approach (to any source) of 20 cm. In source-separation measurements where the detector could be moved freely about the environment, APSL was able to resolve two sources separated by 75 cm or more given only 60 s of measurement time. In these source-separation measurements, APSL produced larger total activity errors of 40%, but obtained source-separation distances accurate to within 15 cm. We also compare our APSL results against traditional maximum likelihood-expectation maximization (ML-EM) reconstructions and demonstrate improved image accuracy and interpretability using APSL over ML-EM. These results indicate that APSL is capable of accurately reconstructing gamma-ray source positions and activities using measurements from existing detector hardware.

The ability to map and estimate the activity of radiological source distributions in unknown three-dimensional environments has applications in the prevention and response to radiological accidents or threats as well as the enforcement and verification of international nuclear non-proliferation agreements. Such a capability requires well-characterized detector response functions, accurate time-dependent detector position and orientation data, a digitized representation of the surrounding 3D environment, and appropriate image reconstruction and uncertainty quantification methods. We have previously demonstrated 3D mapping of gamma-ray emitters with free-moving detector systems on a relative intensity scale using a technique called Scene Data Fusion (SDF). Here we characterize the detector response of a multi-element gamma-ray imaging system using experimentally benchmarked Monte Carlo simulations and perform 3D mapping on an absolute intensity scale. We present experimental reconstruction results from hand-carried and airborne measurements with point-like and distributed sources in known configurations, demonstrating quantitative SDF in complex 3D environments.

By combining radiation detection technologies with robotics sensing, the ability to continuously conduct gamma-ray imaging using freely-moving systems was demonstrated in 2015.1 This new method, which was named free-moving 3D Scene Data Fusion (SDF), was then applied to mapping radioactive contamination and to contextualizing the extent of contamination and the efficacy of radiological clean-up efforts.2,3 Since then, further studies into the types of radiation detection systems to which SDF could be applied resulted in the discovery and demonstration that neutron activity could be mapped using neutron-sensitive CLLBC scintillators, arrays of pix-elated CZT detectors could be used to create multi-modal imagers, and more rudimentary detector systems such as arrays of four CsI modules could still achieve good-quality mapping by inferring source positioning through the encoded modulation of source-to-detector distance. This paper provides an overview of the SDF technology, highlights recent measurements leveraging SDF-equipped systems, discusses the continued development of quantitative algorithms4,5 and their ramifications for developing autonomous SDF-capabilities, and summarizes future directions of research and application development for free moving radiation detection systems.