Normalizing Basis Functions: Approximate Stationary Models for Large Spatial Data

Antony Sikorski, Daniel McKenzie, Douglas Nychka
{"title":"Normalizing Basis Functions: Approximate Stationary Models for Large Spatial Data","authors":"Antony Sikorski, Daniel McKenzie, Douglas Nychka","doi":"arxiv-2405.13821","DOIUrl":null,"url":null,"abstract":"In geostatistics, traditional spatial models often rely on the Gaussian\nProcess (GP) to fit stationary covariances to data. It is well known that this\napproach becomes computationally infeasible when dealing with large data\nvolumes, necessitating the use of approximate methods. A powerful class of\nmethods approximate the GP as a sum of basis functions with random\ncoefficients. Although this technique offers computational efficiency, it does\nnot inherently guarantee a stationary covariance. To mitigate this issue, the\nbasis functions can be \"normalized\" to maintain a constant marginal variance,\navoiding unwanted artifacts and edge effects. This allows for the fitting of\nnearly stationary models to large, potentially non-stationary datasets,\nproviding a rigorous base to extend to more complex problems. Unfortunately,\nthe process of normalizing these basis functions is computationally demanding.\nTo address this, we introduce two fast and accurate algorithms to the\nnormalization step, allowing for efficient prediction on fine grids. The\npractical value of these algorithms is showcased in the context of a spatial\nanalysis on a large dataset, where significant computational speedups are\nachieved. While implementation and testing are done specifically within the\nLatticeKrig framework, these algorithms can be adapted to other basis function\nmethods operating on regular grids.","PeriodicalId":501215,"journal":{"name":"arXiv - STAT - Computation","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Computation","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2405.13821","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

Abstract

In geostatistics, traditional spatial models often rely on the Gaussian Process (GP) to fit stationary covariances to data. It is well known that this approach becomes computationally infeasible when dealing with large data volumes, necessitating the use of approximate methods. A powerful class of methods approximate the GP as a sum of basis functions with random coefficients. Although this technique offers computational efficiency, it does not inherently guarantee a stationary covariance. To mitigate this issue, the basis functions can be "normalized" to maintain a constant marginal variance, avoiding unwanted artifacts and edge effects. This allows for the fitting of nearly stationary models to large, potentially non-stationary datasets, providing a rigorous base to extend to more complex problems. Unfortunately, the process of normalizing these basis functions is computationally demanding. To address this, we introduce two fast and accurate algorithms to the normalization step, allowing for efficient prediction on fine grids. The practical value of these algorithms is showcased in the context of a spatial analysis on a large dataset, where significant computational speedups are achieved. While implementation and testing are done specifically within the LatticeKrig framework, these algorithms can be adapted to other basis function methods operating on regular grids.
归一化基函数:大型空间数据的近似静态模型
在地理统计中,传统的空间模型通常依靠高斯过程(GP)来拟合数据的静态协方差。众所周知,当处理大量数据时,这种方法在计算上变得不可行,因此必须使用近似方法。有一类功能强大的方法将 GP 近似为具有随机系数的基函数之和。虽然这种技术具有计算效率高的特点,但本质上并不能保证协方差的稳定。为了缓解这一问题,可以对基值函数进行 "归一化 "处理,以保持恒定的边际方差,避免不必要的假象和边缘效应。这样就可以将接近静态的模型拟合到大型的、可能是非静态的数据集上,为扩展到更复杂的问题提供了一个严谨的基础。为了解决这个问题,我们为归一化步骤引入了两种快速而精确的算法,从而可以在精细网格上进行高效预测。在对大型数据集进行空间分析时,我们展示了这些算法的实用价值,计算速度明显加快。虽然这些算法是专门在 LatticeKrig 框架内实施和测试的,但它们也可适用于在常规网格上运行的其他基函数方法。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
自引率
0.00%
发文量
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信