Fixed rank kriging for very large spatial data sets

Citation data:

Journal of the Royal Statistical Society. Series B: Statistical Methodology, ISSN: 1369-7412, Vol: 70, Issue: 1, Page: 209-226

Publication Year:
Usage 553
Abstract Views 432
Full Text Views 114
Link-outs 7
Captures 265
Readers 239
Exports-Saves 26
Citations 342
Citation Indexes 342
Repository URL:;
Cressie, Noel A; Johannesson, Gardar
Mathematics; Decision Sciences; spatial; fixed; rank; kriging; very; data; large; sets; Physical Sciences and Mathematics
article description
Spatial statistics for very large spatial data sets is challenging. The size of the data set, n, causes problems in computing optimal spatial predictors such as kriging, since its computational cost is of order. In addition, a large data set is often defined on a large spatial domain, so the spatial process of interest typically exhibits non-stationary behaviour over that domain. A flexible family of non-stationary covariance functions is defined by using a set of basis functions that is fixed in number, which leads to a spatial prediction method that we call fixed rank kriging. Specifically, fixed rank kriging is kriging within this class of non-stationary covariance functions. It relies on computational simplifications when n is very large, for obtaining the spatial best linear unbiased predictor and its mean-squared prediction error for a hidden spatial process. A method based on minimizing a weighted Frobenius norm yields best estimators of the covariance function parameters, which are then substituted into the fixed rank kriging equations. The new methodology is applied to a very large data set of total column ozone data, observed over the entire globe, where n is of the order of hundreds of thousands. © 2008 Royal Statistical Society.