# [1]王金鑫,秦子龙,曹泽宁,等.基于八叉树的修正克里金空间插值算法[J].郑州大学学报(工学版),2021,42(06):22-28.[doi:10.13705/j.issn.1671-6833.2021.06.004] 　WANG Jinxin,QIN Zilong,CAO Zening,et al.Modified Kriging Spatial Interpolation Algorithm Based on Octree Mechanism[J].Journal of Zhengzhou University (Engineering Science),2021,42(06):22-28.[doi:10.13705/j.issn.1671-6833.2021.06.004] 点击复制 基于八叉树的修正克里金空间插值算法() 分享到： var jiathis_config = { data_track_clickback: true };

42卷

2021年06期

22-28

2021-11-10

## 文章信息/Info

Title:
Modified Kriging Spatial Interpolation Algorithm Based on Octree Mechanism

Author(s):
School of Earth Science and Technology of Zhengzhou University; School of Water Conservancy Science and Engineering, Zhengzhou University;

Keywords:
DOI:
10.13705/j.issn.1671-6833.2021.06.004

A

Abstract:
Neighborhood search is an important step in the spatial interpolation algorithm. Whether the neighborhood range is properly selected has a great impact on the interpolation efficiency and accuracy. Aiming at the problem that there were few studies on neighborhood search of interpolation algorithm, a neighborhood search strategy based on octree considering the spatial distribution of discrete points was proposed in this paper. Firstly, the minimal enclosing box of the sampling points was constructed and divided with octree, and the sampling points were grouped into each divided box. Then, the spatial distribution of the interpolating points was constrained by defining the point density. Finally, the above neighborhood search strategy is applied to the ordinary Kriging interpolation model. In order to verify the feasibility of the proposed method, in true 3D geological modeling, the proposed algorithm of this paper, the conventional Kriging interpolation based on fixed distance and fixed number strategy, and the inverse distance weighted interpolation were all used to calculate respectively, and the geological model was constructed from the data obtained from the interpolation. By comparing the method of this paper with the traditional spatial interpolation methods, it was concluded that the proposed method was superior to the traditional method in interpolation accuracy and efficiency when the same number of points were obtained under the same conditions, except for the fixed number method of 30 sample points. Although the fixed number method of 30 sample points had a slight advantage in accuracy, its calculation time is 6.6 times longer than that of the proposed method. In addition, the proposed method improved the time by 20% compared with the traditional method under the same conditions. Compared with the traditional method, the proposed method reduced the redundancy by nearly 1/3 when using the interpolated data to construct the geological model at the same level, thus improving the efficiency of calculation.

## 参考文献/References:

[1] TOBLER W R.A computer movie simulating urban growth in the Detroit region[J].Economic geography,1970,46(增刊1):234-240.

[2] 李小根,王安明.基于GIS的滑坡地质灾害预警预测系统研究[J].郑州大学学报(工学版),2015,36(1):114-118.
[3] ZHANG C S,MCGRATH D.Geostatistical and GIS analyses on soil organic carbon concentrations in grassland of southeastern Ireland from two different periods[J].Geoderma,2004,119(3/4):261-275.
[4] 王靖波,潘懋,张绪定.基于Kriging方法的空间散乱点插值[J].计算机辅助设计与图形学学报,1999,11(6):525-529.
[5] ADHIKARY S K,MUTTIL N,YILMAZ A G.Cokriging for enhanced spatial interpolation of rainfall in two Australian catchments[J].Hydrological processes,2017,31(12):2143-2161.
[6] 孙宗良.基于空间插值的三维近地表建模及可视化研究[D].成都:成都理工大学,2018.
[7] 黄蕾蕾.内蒙古乌努格吐山矿山高精度三维地质建模与评价[D].北京:中国地质大学(北京),2020.
[8] 王长鹏,梁勇,孙黎明,等.一种混合几何曲率和克里金插值的平滑地质曲面构建方法[J].测绘地理信息,2020,45(1):62-65.
[9] 冯波,陈明涛,岳冬冬,等.基于两种插值算法的三维地质建模对比[J].吉林大学学报(地球科学版),2019,49(4):1200-1208.
[10] 李慧晴,叶爱中.基于地形加权的降水空间插值方法研究[J].武汉大学学报(工学版),2021,54(1):28-37.
[11] 曹端广,张子民,王海,等.顾及风向和风速的气温空间插值方法[J].地理与地理信息科学,2021,37(1):47-52.
[12] 莫跃爽,索惠英,焦树林,等.喀斯特地区降水量空间插值方法对比:以贵州省为例[J].水土保持研究,2021,28(1):164-170.
[13] 孙立双,王恩德,王井利,等.基于空间分布权系数Kriging邻域选点算法[J].沈阳建筑大学学报(自然科学版),2007,23(3):423-426.
[14] KEBAILI BARGAOUI Z,CHEBBI A.Comparison of two Kriging interpolation methods applied to spatiotemporal rainfall[J].Journal of hydrology,2009,365(1/2):56-73.
[15] CRESSIE N,JOHANNESSON G.Fixed rank Kriging for very large spatial data sets[J].Journal of the royal statistical society:series B (statistical methodology),2008,70(1):209-226.
[16] MUHAMAD ALI M Z,OTHMAN F.Selection of variogram model for spatial rainfall mapping using analytical hierarchy procedure (AHP)[J].Scientia iranica,2017,24(1):28-39.
[17] 杜宇健,萧德云.Delaunay-固定距离滑动邻域Kriging算法[J].工程图学学报,2005,26(2):64-68.
[18] 徐建华.现代地理学中的数学方法[M].2版.北京:高等教育出版社,2002.
[19] 王政权.地统计学及在生态学中的应用[M].北京：科学出版社，1999.
[20] YARON A F. New methods for spatial statistics in geographic information systems [D]. Columbus: The Ohio State University, 2001.
[21] FRANKE R.Scattered data interpolation:tests of some methods[J].Mathematics of computation,1982,38(157):181.
[22] 唐泽圣.三维数据场可视化[M].北京:清华大学出版社, 1999.
[23] 王金鑫,郑亚圣,李耀辉,等.利用球体剖分瓦块构建真三维数字地球平台[J].测绘学报,2015,44(6):694-701.