Automatically generating maps of a measured variable of interest can be problematic. In this work we focus on the monitoring network context where observations are collected and reported by a network of sensors, and are then transformed into interpolated maps for use in decision making. Using traditional geostatistical methods, estimating the covariance structure of data collected in an emergency situation can be difficult. Variogram determination, whether by method-of-moment estimators or by maximum likelihood, is very sensitive to extreme values. Even when a monitoring network is in a routine mode of operation, sensors can sporadically malfunction and report extreme values. If this extreme data destabilises the model, causing the covariance structure of the observed data to be incorrectly estimated, the generated maps will be of little value, and the uncertainty estimates in particular will be misleading. Marchant and Lark (2007) propose a REML estimator for the covariance, which is shown to work on small data sets with a manual selection of the damping parameter in the robust likelihood. We show how this can be extended to allow treatment of large data sets together with an automated approach to all parameter estimation. The projected process kriging framework of Ingram et al. (2008) is extended to allow the use of robust likelihood functions, including the two component Gaussian and the Huber function. We show how our algorithm is further refined to reduce the computational complexity while at the same time minimising any loss of information. To show the benefits of this method, we use data collected from radiation monitoring networks across Europe. We compare our results to those obtained from traditional kriging methodologies and include comparisons with Box–Cox transformations of the data. We discuss the issue of whether to treat or ignore extreme values, making the distinction between the robust methods which ignore outliers and transformation methods which treat them as part of the (transformed) process. Using a case study, based on an extreme radiological events over a large area, we show how radiation data collected from monitoring networks can be analysed automatically and then used to generate reliable maps to inform decision making. We show the limitations of the methods and discuss potential extensions to remedy these.
[1]
D. Cox,et al.
An Analysis of Transformations
,
1964
.
[2]
M. Opper.
Sparse Online Gaussian Processes
,
2008
.
[3]
J. Pilz,et al.
Bayesian Kriging with lognormal data and uncertain variogram parameters
,
2005
.
[4]
Marc G. Genton,et al.
Highly Robust Variogram Estimation
,
1998
.
[5]
P. Diggle,et al.
Model‐based geostatistics
,
2007
.
[6]
Lehel Csató,et al.
Sparse On-Line Gaussian Processes
,
2002,
Neural Computation.
[7]
Dan Cornford,et al.
Fast algorithms for automatic mapping with space-limited covariance functions
,
2008
.
[8]
Lehel Csat,et al.
Fast Spatial Interpolation using Sparse Gaussian Processes
,
2005
.
[9]
N. Cressie,et al.
Robust estimation of the variogram: I
,
1980
.
[10]
R. M. Lark,et al.
Robust estimation of the variogram by residual maximum likelihood
,
2007
.