A multi-resolution approximation for massive spatial datasets

Abstract

Automated sensing instruments on satellites and aircraft have enabled the collection of massive amounts of high-resolution observations of spatial fields over large spatial regions. If these datasets can be efficiently exploited, they can provide new insights on a wide variety of issues. However, traditional spatial-statistical techniques such as kriging are not computationally feasible for big datasets. We propose a multi-resolution approximation (M-RA) of Gaussian processes observed at irregular locations in space. The M-RA process is specified as a linear combination of basis functions at multiple levels of spatial resolution, which can capture spatial structure from very fine to very large scales. The basis functions are automatically chosen to approximate a given covariance function, which can be nonstationary. All computations involving the M-RA, including parameter inference and prediction, are highly scalable for massive datasets. Crucially, the inference algorithms can also be parallelized to take full advantage of large distributed-memory computing environments. In comparisons using simulated data and a large satellite dataset, the M-RA outperforms a related state-of-the-art method.

23 pages; to be published in Journal of the American Statistical Association

References in corpus (2)

• Covariance approximation for large multivariate spatial data sets with an application to multiple climate model errors
• When does the screening effect hold?

Cited by in corpus (21)

• A general framework for Vecchia approximations of Gaussian processes
• A class of multi-resolution approximations for large spatial datasets
• Vecchia approximations of Gaussian-process predictions
• The rational SPDE approach for Gaussian random fields with general smoothness
• Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains
• Practical Bayesian Modeling and Inference for Massive Spatial Datasets On Modest Computing Environments
• Vecchia-Laplace approximations of generalized Gaussian processes for big non-Gaussian spatial data
• On statistical approaches to generate Level 3 products from satellite remote sensing retrievals
• Online Joint State Inference and Learning of Partially Unknown State-Space Models
• Scalable Bayesian modeling for smoothing disease risks in large spatial data sets
• Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data
• Adaptive LASSO estimation for functional hidden dynamic geostatistical model
• Modeling Massive Spatial Datasets Using a Conjugate Bayesian Linear Regression Framework
• A Large Scale Spatio-temporal Binomial Regression Model for Estimating Seroprevalence Trends
• Big problems in spatio-temporal disease mapping: methods and software
• High-dimensional order-free multivariate spatial disease mapping
• Correlation-based sparse inverse Cholesky factorization for fast Gaussian-process inference
• Identification of Dominant Features in Spatial Data
• The SPDE approach for Gaussian and non-Gaussian fields: 10 years and still running
• Inside-out cross-covariance for spatial multivariate data
• Estimating Basis Functions in Massive Fields under the Spatial Mixed Effects Model