An improved coarse-grained parallel algorithm for computational acceleration of ordinary Kriging interpolation

Heavy computation limits the use of Kriging interpolation methods in many real-time applications, especially with the ever-increasing problem size. Many researchers have realized that parallel processing techniques are critical to fully exploit computational resources and feasibly solve computation-intensive problems like Kriging. Much research has addressed the parallelization of traditional approach to Kriging, but this computation-intensive procedure may not be suitable for high-resolution interpolation of spatial data. On the basis of a more effective serial approach, we propose an improved coarse-grained parallel algorithm to accelerate ordinary Kriging interpolation. In particular, the interpolation task of each unobserved point is considered as a basic parallel unit. To reduce time complexity and memory consumption, the large right hand side matrix in the Kriging linear system is transformed and fixed at only two columns and therefore no longer directly relevant to the number of unobserved points. The MPI (Message Passing Interface) model is employed to implement our parallel programs in a homogeneous distributed memory system. Experimentally, the improved parallel algorithm performs better than the traditional one in spatial interpolation of annual average precipitation in Victoria, Australia. For example, when the number of processors is 24, the improved algorithm keeps speed-up at 20.8 while the speed-up of the traditional algorithm only reaches 9.3. Likewise, the weak scaling efficiency of the improved algorithm is nearly 90% while that of the traditional algorithm almost drops to 40% with 16 processors. Experimental results also demonstrate that the performance of the improved algorithm is enhanced by increasing the problem size. We propose an improved coarse-grain parallel algorithm of ordinary Kriging.We compare two parallel algorithms derived from improved and traditional approaches.Parallel programs based on MPI are implemented in a distributed memory system.Improved algorithm performances better in both parallel efficiency and scalability.Increasing the problem size enhances the performance of improved algorithm.

[1]  Alexander Litvinenko,et al.  Kriging and Spatial Design Accelerated by Orders of Magnitude: Combining Low-Rank Covariance Approximations with FFT-Techniques , 2013, Mathematical Geosciences.

[2]  Albrecht Gebhardt,et al.  PVM kriging with R , 2003 .

[3]  G. Heuvelink,et al.  A generic framework for spatial prediction of soil variables based on regression-kriging , 2004 .

[4]  D. Nychka,et al.  Covariance Tapering for Interpolation of Large Spatial Datasets , 2006 .

[5]  Jian Wang,et al.  Explorations of the implementation of a parallel IDW interpolation algorithm in a Linux cluster-based parallel GIS , 2011, Comput. Geosci..

[6]  P. Goovaerts.  Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall , 2000 .

[7]  N. Memarsadeghi,et al.  Efficient Kriging via Fast Matrix-Vector Products , 2008, 2008 IEEE Aerospace Conference.

[8]  Kenneth A. Hawick,et al.  Kriging Interpolation on High-Performance Computers , 1998, HPCN Europe.

[9]  J. Jones,et al.  A parallel implementation of kriging with a trend , 1997 .

[10]  Ola Hössjer,et al.  Fast kriging of large data sets with Gaussian Markov random fields , 2008, Comput. Stat. Data Anal..

[11]  Barry Wilkinson,et al.  Parallel programming , 1998 .

[12]  Ana Cortés,et al.  Parallel ordinary kriging interpolation incorporating automatic variogram fitting , 2011, Comput. Geosci..

[13]  Peter S. Pacheco.  An Introduction to Parallel Programming , 2011 .

[14]  Xuan Shi,et al.  Kriging interpolation over heterogeneous computer architectures and systems , 2013 .

[15]  Mark Richards,et al.  A regionalized national universal kriging model using Partial Least Squares regression for estimating annual PM2.5 concentrations in epidemiology. , 2013, Atmospheric environment.

[16]  Marc P. Armstrong,et al.  Massively parallel strategies for local spatial interpolation , 1997 .

[17]  Huayi Wu,et al.  Leveraging the power of multi-core platforms for large-scale geospatial data processing: Exemplified by generating DEM from massive LiDAR point clouds , 2010, Comput. Geosci..

[18]  Francisco J. Jiménez-Hornero,et al.  Using general-purpose computing on graphics processing units (GPGPU) to accelerate the ordinary kriging algorithm , 2014, Comput. Geosci..

[19]  Michael W.D. Davis,et al.  Kriging in a global neighborhood , 1984 .

[20]  Marc P. Armstrong,et al.  Local Interpolation Using a Distributed Parallel Supercomputer , 1996, Int. J. Geogr. Inf. Sci..

[21]  Douglas W. Nychka,et al.  Covariance Tapering for Likelihood-Based Estimation in Large Spatial Data Sets , 2008 .

[22]  Zbigniew J. Czech,et al.  Introduction to Parallel Computing , 2017 .

[23]  N. Cressie,et al.  Fixed rank kriging for very large spatial data sets , 2008 .

[24]  Michael F. Goodchild,et al.  A parallel computing approach to fast geostatistical areal interpolation , 2011, Int. J. Geogr. Inf. Sci..

[25]  W. Nowak,et al.  Application of FFT-based Algorithms for Large-Scale Universal Kriging Problems , 2009 .

[26]  Holdaway.  Spatial modeling and interpolation of monthly temperature using kriging , 1996 .

[27]  Akula Venkatram,et al.  On the use of kriging in the spatial analysis of acid precipitation data , 1988 .

[28]  Michael J. Quinn,et al.  Parallel programming in C with MPI and OpenMP , 2003 .

[29]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[30]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[31]  Tangpei Cheng,et al.  Accelerating universal Kriging interpolation algorithm using CUDA-enabled GPU , 2013, Comput. Geosci..

[32]  Qun Wang,et al.  On Parallelizing Universal Kriging Interpolation Based on OpenMP , 2010, 2010 Ninth International Symposium on Distributed Computing and Applications to Business, Engineering and Science.

[33]  Shaowen Wang,et al.  A quadtree approach to domain decomposition for spatial interpolation in Grid computing environments , 2003, Parallel Comput..