Abstract:In this investigation, we propose several algorithms to recover the location and intensity of a radiation source located in a simulated 250 m x 180 m block in an urban center based on synthetic measurements. Radioactive decay and detection are Poisson random processes, so we employ likelihood functions based on this distribution. Due to the domain geometry and the proposed response model, the negative logarithm of the likelihood is only piecewise continuous differentiable, and it has multiple local minima. To address these difficulties, we investigate three hybrid algorithms comprised of mixed optimization techniques. For global optimization, we consider Simulated Annealing (SA), Particle Swarm (PS) and Genetic Algorithm (GA), which rely solely on objective function evaluations; i.e., they do not evaluate the gradient in the objective function. By employing early stopping criteria for the global optimization methods, a pseudo-optimum point is obtained. This is subsequently utilized as the initial value by the deterministic Implicit Filtering method (IF), which is able to find local extrema in non-smooth functions, to finish the search in a narrow domain. These new hybrid techniques combining global optimization and Implicit Filtering address difficulties associated with the non-smooth response, and their performances are shown to significantly decrease the computational time over the global optimization methods alone. To quantify uncertainties associated with the source location and intensity, we employ the Delayed Rejection Adaptive Metropolis (DRAM) and DiffeRential Evolution Adaptive Metropolis (DREAM) algorithms. Marginal densities of the source properties are obtained, and the means of the chains' compare accurately with the estimates produced by the hybrid algorithms.

Multiscale spatial density smoothing: an application to large-scale radiological survey and anomaly detection

Nonparametric Density Estimation for Data Scattered on Irregular Spatial Domains: A Likelihood-Based Approach Using Bivariate Penalized Spline Smoothing

Sensitivity-Aware Density Estimation in Multiple Dimensions

Flexible Basis Representations for Modeling Large Non-Gaussian Spatial Data

Hybrid optimization and Bayesian inference techniques for a non-smooth radiation detection problem

A Multilevel Stratified Spatial Sampling Approach for the Quality Assessment of Remote-Sensing-Derived Products

Spatial Deconvolution of Aerial Radiometric Survey and its application to the Fallout from a Radiological Dispersal Device

Multi-scale decomposition of astronomical maps -- a constrained diffusion method

A Scalable Partitioned Approach to Model Massive Nonstationary Non-Gaussian Spatial Datasets

A Weakly Supervised Approach for Estimating Spatial Density Functions from High-Resolution Satellite Imagery

Volumetric Data Exploration with Machine Learning-Aided Visualization in Neutron Science

Collective Spectral Density Estimation and Clustering for Spatially-Correlated Data

Maximizing Multi-scale Spatial Statistical Discrepancy.

Multiscale scattered data analysis in samplet coordinates

Recover Fine-Grained Spatial Data from Coarse Aggregation

Exact Bayesian Inference for Multivariate Spatial Data of Any Size with Application to Air Pollution Monitoring

Density-aware Stratified Sampling for Visualizing Large Volume Geo-Spatial Data

Thresholded Multiscale Gaussian Processes with Application to Bayesian Feature Selection for Massive Neuroimaging Data

Localizing Anomalies via Multiscale Score Matching Analysis

Multiscale Multi-Type Spatial Bayesian Analysis for High-Dimensional Data with Application to Wildfires and Migration

Kryging: Geostatistical analysis of large-scale datasets using Krylov subspace methods