Quantifying the Uncertainty of Contour Maps

ABSTRACT Contour maps are widely used to display estimates of spatial fields. Instead of showing the estimated field, a contour map only shows a fixed number of contour lines for different levels. However, despite the ubiquitous use of these maps, the uncertainty associated with them has been given a surprisingly small amount of attention. We derive measures of the statistical uncertainty, or quality, of contour maps, and use these to decide an appropriate number of contour lines, which relates to the uncertainty in the estimated spatial field. For practical use in geostatistics and medical imaging, computational methods are constructed, that can be applied to Gaussian Markov random fields, and in particular be used in combination with integrated nested Laplace approximations for latent Gaussian models. The methods are demonstrated on simulated data and an application to temperature estimation is presented.

[1]  Finn Lindgren,et al.  Bayesian Spatial Modelling with R-INLA , 2015 .

[2]  J. French.  Confidence regions for the level curves of spatial data , 2014 .

[3]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[4]  David Bolin,et al.  A comparison between Markov approximations and other methods for large spatial data sets , 2013, Comput. Stat. Data Anal..

[5]  Finn Lindgren,et al.  Bayesian computing with INLA: New features , 2012, Comput. Stat. Data Anal..

[6]  F. Lindgren,et al.  Excursion and contour uncertainty regions for latent Gaussian models , 2012, 1211.3946.

[7]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach , 2011 .

[8]  F. Lindgren,et al.  Spatial models generated by nested stochastic partial differential equations, with an application to global ozone mapping , 2011, 1104.3436.

[9]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[10]  C. Amante,et al.  ETOPO1 arc-minute global relief model : procedures, data sources and analysis , 2009 .

[11]  S. Martino.  Approximate Bayesian Inference for Latent Gaussian Models , 2007 .

[12]  A. Wameling.  Vertical and Horizontal Spatial Variation of Geostatistical Prediction , 2006 .

[13]  Thomas Polfeldt.  On the quality of contour maps , 1999 .

[14]  N. Thrower,et al.  Maps and Civilization: Cartography in Culture and Society , 1996 .

[15]  Igor Rychlik,et al.  How reliable are contour curves? Confidence sets for level contours , 1995 .