基岩露头采样密度对黄土古地形重建的影响

doi:10.3724/SP.J.1047.2016.00461

摘要/Abstract

摘要：

下伏古地形对黄土地貌的形成、演化与发育具有重要的控制作用。基于下伏地层露头的采样数据,通过空间插值方法进行古地形数字模型(DEM)重建是研究黄土地貌的重要方法。其中,采样点密度是下伏地形DEM重建精度的主要影响因子。本文以1:20万绥德幅地质图所涉及区域为实验样区,研究采样点密度及样条函数方法对下伏古地形重建精度的影响。实验结果表明,在高密度样本条件下,采用规则样条插值方法进行古地形重建时,插值结果呈现显著的“龙格现象”,即多项式插值结果出现剧烈波动现象,且随着样本密度增加,样条插值结果的误差呈逐步趋缓的降低趋势,并逐步趋于稳定。同时,基于古地形DEM提取的特征点,其局部最高点和局部最低点的数目呈直线上升趋势,说明在一定程度上基于有限数据采用样条函数进行地下三维建模并不一定能获得平滑曲面。研究结果对如何选择合理的空间采样密度进行古地形DEM重建具有重要的指导意义。

关键词: DEM, 黄土地貌, 下伏古地形, 采样点密度, 样条空间插值

Abstract:

The Pre-Quaternary underlying terrain profoundly controls the evolution and formation of loess landform. Obvious relationships, i.e. the geomorphological inheritance, could be found between the underlying terrain and the modern terrain. As a consequence, the Pre-Quaternary underlying terrain in the Loess Plateau should be regarded as the key factor for the understanding of the loess landform evolution. Among numerous numerical calculation methods, spatial interpolation has been regarded as an important method to reconstruct the DEM of underlying terrain by using the sampled bedrock outcrop points selected from a geological map. However, the sampling density has a great impact on the accuracy of the reconstructed underlying terrain. In this paper, the Suide geological map area (1:200 000) was selected as the study area, and then the influence of sampling density on the accuracy of the reconstructed underlying terrain was investigated using spline method. By adopting cross-validation method to evaluate reconstructed underlying terrain, the result shows that, different interpolation methods cause uncertainties to different degrees during the reconstruction of underlying terrain, particularly the spline method. On a basis of high density outcrop points and spline function interpolation process, the morphology of underlying terrain exhibits a typical “Runge phenomenon”. This phenomenon was always resulted from a polynomial interpolation process. With an increased sampling density, the error in underlying terrain appears a slowly decrease tendency firstly, and then it keeps stable. Meanwhile, the number of the extracted features has a linear upward trend. The result also shows that the sampling density of 1.7-2.0 points per square kilometer could achieve a good balance between the accuracy and underlying terrain feature reservation. The aforementioned results adjust our previous understandings that spline function could smooth the interpolated surface to some extent. And the result also provides guidance for the selection of a reasonable spatial sampling density.

Key words: DEM, Loess underlying terrain, sampling density, spline interpolation

段家朕, 熊礼阳, 汤国安. 基岩露头采样密度对黄土古地形重建的影响[J]. 地球信息科学学报, 2016, 18(4): 461-468.DOI:10.3724/SP.J.1047.2016.00461

DUAN Jiazhen,XIONG Liyang,TANG Guoan. Effect of Outcrop Sampling Density on the Underlying Terrain Reconstruction[J]. Journal of Geo-information Science, 2016, 18(4): 461-468.DOI:10.3724/SP.J.1047.2016.00461

图/表 14

图1

图2

图3

表1

图4

图5

图6

表2

图7

图8

表3

图9

图10

图11

参考文献 24

[1]	刘东生. 黄土与环境[M].北京:科学出版社,1985.
	[ Liu D S, Loess and the environment[M]. Beijing Science Press,1985. ]
[2]	夏正楷. 黄土高原第四纪期间水土流失的地质记录和基本规律[J].水土保持研究,1999,6(4):49-53. doi: 10.3969/j.issn.1005-3409.1999.04.010
	[ Xia Z K.The records of quaternary soil erosion in the Loess Plateau[J]. Research of Soil and Water Conservation, 1999,6(4):49-53. ] doi: 10.3969/j.issn.1005-3409.1999.04.010
[3]	桑广书,陈雄,陈小宁,等.黄土丘陵地貌形成模式与地貌演变[J].干旱区地理,2007,30(3):375-380. doi: 10.3321/j.issn:1000-6060.2007.03.009
	[ Sang G S, Chen X, Chen X N, et al.Formation model and geomorphic evolution of loess hilly landforms[J]. Arid Land Geography, 2007,30(3):375-380. ] doi: 10.3321/j.issn:1000-6060.2007.03.009
[4]	程彦培,石建省,杨振京,等.古地形对黄土区岩土侵蚀趋势的控制作用[J].干旱区地理,2010,33(3):334-339.
	[ Cheng Y P, Shi J S, Yang Z J, et al.Control of ancient landform on rock-soil erosion in loess area[J]. Arid Land Geography, 2010,33(3):334-339. ]
[5]	袁宝印,郭正堂,郝青振,等.天水-秦安一带中新世黄土堆积区沉积-地貌演化[J].第四纪研究,2007,27(2):161-171.
	[ Yuan B Y, Guo Z T, Hao Q Z, et al.Cenozoic evolution of geomorphic and sedimentary environments in the tianshui-qin’an regions[J]. Quaternary Sciences, 2007,27(2):161-171. ]
[6]	李新,程国栋,卢玲.空间内插方法比较[J].地球科学进展,2000,15(3):260-265.
	[ Li X, Cheng G D, Lu L.Comparison of spatial interpolation methods[J]. Advance in Earthences, 2000,15(3):260-265. ]
[7]	李志林,朱庆.数字高程模型[M].武汉:武汉大学出版社,2000:125-139.
	[ Li Z L, Zhu Q.Digital elevation model[M]. Wuhan: Wuhan University Press, 2000:125-139. ]
[8]	Lam N S.Spatial interpolation methods: A review[J]. The American Cartographer, 1983,10(2):129-149. doi: 10.1559/152304083783914958
[9]	Robeson S M.Spherical methods for spatial interpolation: Review and evaluation[J]. Cartography and Geographic Information Systems, 1997,24(1):3-20. doi: 10.1559/152304097782438746
[10]	Burrough P A, McDonnell R A. Principles of geographical information systems[M]. New York: Oxford University Press, 1998:333.
[11]	Franke R.Scattered data interpolation: Tests of some methods[J]. Mathematics of Computation, 1982,38(157):181-200. doi: 10.1090/S0025-5718-1982-0637296-4
[12]	Weber D, Englund E.Evaluation and comparison of spatial interpolators II[J]. Mathematical Geology, 1994,26(5):589-603. doi: 10.1007/BF02089243
[13]	Declercq F A N. Interpolation methods for scattered sample data: Accuracy, spatial patterns, processing time[J]. Cartography and Geographic Information Systems, 1996,23(3):128-144. doi: 10.1559/152304096782438882
[14]	Yang X, Hodler T.Visual and statistical comparisons of surface modeling techniques for point-based environmental data[J]. Cartography and Geographic Information Science, 2000,27(2):165-175. doi: 10.1559/152304000783547911
[15]	张锦明,游雄,万刚.DEM插值参数优选的试验研究[J].测绘学报,2014,43(2):178-185. doi: 10.13485/j.cnki.11-2089.2014.0026
	[ Zhang J M, You X, Wan G.Experimental research on optimization of DEM interpolation parameters[J]. Acta Geodaetica et Cartographica Sinica, 2014,43(2):178-185. ] doi: 10.13485/j.cnki.11-2089.2014.0026
[16]	熊礼阳,汤国安,袁宝印,等.基于DEM的黄土高原(重点流失区)地貌演化的继承性研究[J].中国科学:地球科学,2014,44(2):313-321.
	[ Xiong L Y, Tang G A, Yuan B Y, et al.Geomorphological inheritance for loess landform evolution in a severe soil erosion region of Loess Plateau of China based on digital elevation models[J]. Science China: Earth Sciences, 2014,44(2):313-321. ]
[17]	王春,江岭,徐静,等.DEM地面形态重构方法的精度差异特征研究[J].地理与地理信息科学,2014,30(4):18-21,26. doi: 10.3969/j.issn.1672-0504.2014.04.004
	[ Wang C, Jiang L, Jing X U, et al.Accuracy differences among approaches of DEM surface morphology reconstruction[J]. Geography and Geo-Information Science, 2014,30(4):18-21,26. ] doi: 10.3969/j.issn.1672-0504.2014.04.004
[18]	Aguilar F J, Agüera F, Aguilar M A, et al.Effects of terrain morphology, sampling density, and interpolation methods on grid DEM accuracy[J]. Photogrammetric Engineering & Remote Sensing, 2005,71(7):805-816. doi: 10.1002/ppp.529
[19]	王春,刘学军,汤国安,等.格网DEM地形模拟的形态保真度研究[J].武汉大学学报(信息科学版),2009,34(2):146-149.
	[ Wang C, Liu X J, Tang G A, et al.Morphologic fidelity of grid digital elevation model[J]. Geomatics & Information Science of Wuhan University, 2009,34(2):146-149. ]
[20]	黄秉维. 编制黄河中游流域土壤侵蚀分区图的经验教训[J].科学通报,1955,12:15-21,14.
	[ Huang B W.The Experience and lessons of the compilation of soil erosion zoning map in the middle reaches of the Yellow River[J]. Chinese Science Bulletin, 1955,12:15-21,14. ]
[21]	田剑,汤国安,周毅,等.黄土高原沟谷密度空间分异特征研究[J].地理科学,2013,33(5):622-628.
	[ Tian J, Tang G A, Zhou Y, et al.Spatial variation of gully density in the Loess Plateau[J]. Scientia Geographica Sinica, 2013,33(5):622-628. ]
[22]	邬伦,吴小娟,肖晨超,等.五种常用降水量插值方法误差时空分布特征研究——以深圳市为例[J].地理与地理信息科学,2010,26(3):19-24.
	[ Wu L, Wu X J, Xiao C C, et al.On temporal and spatial error distributions of five precipitation interpolation models: a case of Shenzhen[J]. Geography and Geo-Information Science, 2010,26(3):19-24. ]
[23]	王春,顾留碗,陶旸,等.DEM地形描述误差(Et)计算模型研究[J].地球信息科学学报,2014,16(5):699-706. doi: 10.3724/SP.J.1047.2014.00699
	[ Wang C, Gu L W, Tao Y, et al.A new calculation model of DEM terrain description error[J]. Journal of Geo-information Science, 2014,16(5):699-706. ] doi: 10.3724/SP.J.1047.2014.00699
[24]	仲腾,汤国安,周毅,等.基于反地形DEM的山顶点自动提取[J].测绘通报,2009(4):35-37.
	[ Zhong T, Tang G A, Zhou Y, et al.Method of extracting surface peaks based on reverse DEMs[J]. Bulletin of Surveying & Mapping,2009,4:35-37. ]

样本集	n1	n2	n3	n4	n5	n6	n7	n8	n9	n10
点数	2531	5061	7592	10 122	12 653	15 183	17 714	20 244	22 775	25 305
等效密度/（个/km²）	0.346	0.693	1.040	1.386	1.733	2.080	2.427	2.773	3.120	3.466

样本集	n1	n2	n3	n4	n5	n6	n7	n8	n9	n10
MAE	16.132	13.325	11.295	10.912	9.032	9.032	8.004	7.343	6.718	6.587
RMSE	31.647	26.701	21.030	25.904	17.708	18.607	16.358	13.654	12.612	13.060
MRE	0.018	0.015	0.013	0.012	0.010	0.010	0.009	0.008	0.007	0.007

样本集	n1	n2	n3	n4	n5	n6	n7	n8	n9	n10
LHP	300	589	895	1234	1515	1853	2158	2457	2762	3097
LLP	370	761	1105	1493	1895	2308	2650	2968	3279	3641