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

Abstract

We consider the problem of estimating a spatially varying density function, motivated by problems that arise in large-scale radiological survey and anomaly detection. In this context, the density functions to be estimated are the background gamma-ray energy spectra at sites spread across a large geographical area, such as nuclear production and waste-storage sites, military bases, medical facilities, university campuses, or the downtown of a city. Several challenges combine to make this a difficult problem. First, the spectral density at any given spatial location may have both smooth and non-smooth features. Second, the spatial correlation in these density functions is neither stationary nor locally isotropic. Finally, at some spatial locations, there is very little data. We present a method called multiscale spatial density smoothing that successfully addresses these challenges. The method is based on recursive dyadic partition of the sample space, and therefore shares much in common with other multiscale methods, such as wavelets and Pólya-tree priors. We describe an efficient algorithm for finding a maximum a posteriori (MAP) estimate that leverages recent advances in convex optimization for non-smooth functions. We apply multiscale spatial density smoothing to real data collected on the background gamma-ray spectra at locations across a large university campus. The method exhibits state-of-the-art performance for spatial smoothing in density estimation, and it leads to substantial improvements in power when used in conjunction with existing methods for detecting the kinds of radiological anomalies that may have important consequences for public health and safety.

36 pages, 10 figures

References in corpus (8)

• A multivariate semiparametric Bayesian spatial modeling framework for hurricane surface wind fields
• Fast and Flexible ADMM Algorithms for Trend Filtering
• Total Variation Classes Beyond 1d: Minimax Rates, and the Limitations of Linear Smoothers
• Trend Filtering on Graphs
• Detecting changes in maps of gamma spectra with Kolmogorov-Smirnov tests
• Spatially-Aware Temporal Anomaly Mapping of Gamma Spectra
• Bayesian Trend Filtering
• Locally adaptive smoothing with Markov random fields and shrinkage priors

Cited by in corpus (7)

• Reconstructing the Position and Intensity of Multiple Gamma-Ray Point Sources with a Sparse Parametric Algorithm
• Detecting changes in maps of gamma spectra with Kolmogorov-Smirnov tests
• Interpretable Low-Dimensional Regression via Data-Adaptive Smoothing
• Modelling Heterogeneity Using Bayesian Structured Sparsity
• How Likely are Ride-share Drivers to Earn a Living Wage? Large-scale Spatio-temporal Density Smoothing with the Graph-fused Elastic Net
• Deep Nonparametric Estimation of Discrete Conditional Distributions via Smoothed Dyadic Partitioning
• GapTV: Accurate and Interpretable Low-Dimensional Regression and Classification