Accéder directement au contenu Accéder directement à la navigation
Article dans une revue

Incomplete Cholesky Decomposition for the kriging of large datasets

Abstract : Kriging of very large spatial datasets is a challenging problem. The size nn of the dataset causes problems in computing the kriging estimate: solving the kriging equations directly involves inverting an n×nn×n covariance matrix. This operation requires O(n3) computations and a storage of O(n2). Under these circumstances, straightforward kriging of massive datasets is not possible. Several approaches have been proposed in the literature among which two main families exist: sparse approximation of the covariance function and low rank approaches. We propose here an approach that is built upon a low rank approximation of the covariance matrix obtained by incomplete Cholesky decomposition. This algorithm requires O(nk) storage and takes O(nk2) arithmetic operations, where kk is the rank of the approximation, whose accuracy is controlled by a parameter. We detail the main properties of this method and explore its links with existing methods. Its benefits are illustrated on simple examples and compared to those of existing approaches. Finally, we show that this low rank representation is also suited for inverse conditioning of Gaussian random fields.
Type de document :
Article dans une revue
Liste complète des métadonnées
Contributeur : Thomas Romary <>
Soumis le : vendredi 19 juillet 2013 - 18:15:38
Dernière modification le : jeudi 24 septembre 2020 - 16:34:08


  • HAL Id : hal-00846756, version 1


Thomas Romary. Incomplete Cholesky Decomposition for the kriging of large datasets. Spatial Statistics, Elsevier, 2013, 5 (85-99), ⟨hal-00846756⟩



Consultations de la notice