DEM分辨率与平均坡度的关系分析

doi:10.3724/SP.J.1047.2014.00524

摘要/Abstract

摘要：

由DEM提取坡度时,DEM分辨率的变化会造成局部栅格提取的坡度数值的变化,其变化规律对于深刻理解平均坡度,随DEM分辨率的变化而变化的规律具有现实意义。在黄土高原的研究中,选择了8个典型地貌类型的样区,获得了以1:10 000比例尺地形图为数据源的DEM(分辨率为5 m)。从DEM分辨率对栅格坡度影响规律的分析中,得到如下结论：(1)随着DEM分辨率粗略化,所提取的坡度平均值有可能升高;(2)在所有样区,随DEM分辨率粗略化,负误差的栅格均是正误差的栅格数量的1.13倍以上,保持负误差的栅格是保持正误差的栅格数量的1.81-4.04倍,同时坡度的负误差均值是正误差均值变化幅度的3.1倍以上;(3)结合结论(2)和本文提出的指标,可准确地解释8个样区中坡度平均值随DEM分辨率变化而变化的规律;(4)DEM的分辨率对正负误差的栅格数量没有明显影响,正误差(或负误差)栅格均集中于个别坡度的分级区间,而DEM分辨率对此集中的区间没有明显影响;(5)坡度总体误差可根据本研究所提出的经验公式由DEM分辨率和分辨率为5 m的 DEM所提取坡度的平均值计算出来。

关键词: 数字高程模型, 坡度, 误差, 分辨率

Abstract:

At present, slopes are usually derived from discrete DEMs. When deriving slope from DEM, the variation of DEM resolution will cause changes in the slope values which are extracted from the local grids. The variation will be significant in understating the general rule of the variation of average slope with respect to DEM resolution. 8 areas representing the typical geomorphology in Loess Plateau were chosen as the sample areas, and the DEMs with a 5m resolution were established according to the topographic map scaled at 1:10 000. In every sample area, 22 500 sampling points were selected, and the slopes which were derived according to the algorithm of third-order finite difference weighted by reciprocal of squared distance on the sampling points were taken as the research object. The regression analysis, single-factor variance analysis, contrast analysis were adopted to analyze the rule of the influence of DEM resolution on grid slopes. The following conclusions were obtained: (1) the average slope may increase as the resolution of DEM becomes coarser; (2) in all the sample areas, with the resolution becoming coarser, the number of the grids with negative slope error is at least 1.13 times of the grids with positive slope error. The amplitude of average slope of the grids with negative slope error is greater than 3.1 times of the latter. The number of grids with decreasing slopes is 1.81-4.04 times of the grids with increasing slopes; (3) the change of average slope can be explained in this research by combining the author proposed index with the previously mentioned conclusion (2); (4) the number of grids with varied slope values is not sensitive to the resolution of DEM. The slopes of the grids which have positive slope error or negative slope error are intensively distributed in certain intervals of slopes that derived from the DEMs with 5m resolution. The intervals are not sensitive to the resolutions of DEM; (5) the total slope error can be computed, according to the empirical formula proposed in this research and calculated from the DEM resolution and the average slope of DEM with a 5m resolution.

Key words: DEM, slope, error, resolution

陈楠. DEM分辨率与平均坡度的关系分析[J]. 地球信息科学学报, 2014, 16(4): 524-530.

CHEN Nan*. Relationship Between DEM Resolution and Average Slope Derived from DEM[J]. Journal of Geo-information Science, 2014, 16(4): 524-530.

图/表 10

图1

表1

图2

图3

图4

图5

图6

表2

表3

表4

参考文献 18

[1]	刘学军,张平.DEM坡度、坡向的有效尺度范围[J].武汉大学学报(信息科学版),2008,33(12):1254-1258.
[2]	汤国安,赵牡丹,李天文,等.DEM提取黄土高原地面坡度的不确定性[J].地理学报,2003,58(6):824-830.
[3]	邓仕虎,杨勤科.DEM采样间隔对地形描述精度的影响研究[J].地理与地理信息科学,2010,26(2):23-26.
[4]	Shi W Z, Tian Y.A Hybrid Interpolation method for the refinement of a regular grid digital elevation model[J]. International Journal of Geographical Information Science, 2006,20(1):53-67.
[5]	Jai V, Jin T, Georgina S.Impact of DEM accuracy and resolution on topographic indices[J]. Environmental Modeling & Software, 2010,25(10):1086-1098.
[6]	Erskine P H, Green T R, Ramirez J A, et al.Digital elevation accuracy and grid cell size: Effects on estimated terrain attributes[J]. Soil & Water Management & Conservation, 2007,71(4):1371-1380.
[7]	胡鹏. 新数字高程模型:理论、方法、标准和应用[M].北京:测绘出版社,2007:46-64.
[8]	周启明,刘学军.数字地形分析[M].北京:科学出版社,2006,266-290.
[9]	李志林,朱庆.数字高程模型[M].武汉:武汉测绘科技大学出版社,2000,60-72.
[10]	汤国安,刘学军,房亮,等.DEM及数字地形分析中尺度问题研究综述[J].武汉大学学报(信息科学版),2006,31(12):1059-1066.
[11]	Gao J, James E, Zhu A X.Neighborhood size and spatial scale in raster-based slope calculations[J]. International Journal of Geographical Information Science, 2012,26(10):1959-1978.
[12]	Zhang J X, Chang K T, Wu J Q.Effects of DEM resolution and source on soil erosion modeling: A case study using the WEPP model[J]. International Journal of Geographical Information Science, 2008,22(8):925-942.
[13]	Tang G A, Li F Y, Liu X J, et al.Research on the slope spectrum of the loess plateau[J].Science in China(E), 2008,51(1):175-185.
[14]	赵明伟,汤国安,张磊,等.基于DEM的5节点二阶差分坡度算法模型与实验分析[J].地球信息科学学报, 2012,14(3):286-291.
[15]	卢华兴,刘学军,王永君,等.插值条件下格网DEM 坡度计算模型的噪声误差分析[J].测绘学报,2012,41(6):926-932.
[16]	Livnea E,Svoraya T.Components of uncertainty in primary production model: The study of DEM, classification and location error[J]. International Journal of Geographical Information Science, 2011,25(3):473-488.
[17]	刘红艳,杨勤科,王春梅,等.坡度随水平分辨率变化及其空间格局研究[J].武汉大学学报(信息科学版),2012,37(1):105-109.
[18]	陈楠. DEM分辨率变化对坡度误差的影响[J].武汉大学学报(信息科学版),2013,38(5):594-598.

样区编号	样区名称	平均坡度(º)	海拔高度范围(m)	地貌类型
1	安塞	37.92	500~1846	黄土高原丘陵沟壑区
2	长武	13.32	161~1215	黄土塬
3	淳化	12.41	768~1188	黄土塬梁状丘陵区
4	富县	26.76	1147~1459	黄土梁峁状丘陵沟壑区
5	神木	9.45	1035~1323	风沙黄土过渡区域
6	绥德	26.39	820~1180	黄土梁峁状丘陵沟壑区
7	延川	31.30	923~1252	黄土梁峁状丘陵沟壑区
8	宜君	19.26	767~1158	黄土梁状丘陵区

比较项	样区
比较项	安塞	长武	淳化	富县	神木	绥德	延川	宜君
一直保持正误差的栅格个数(个)	5042	4137	4093	3548	3350	3644	2791	4237
一直保持正误差的栅格个数与10m DEM正误差栅格个数之比(%)	57.96	44.63	48.73	36.63	43.67	38.93	31.72	44.19
一直保持负误差的栅格个数(个)	9121	9256	10 235	10 078	11 311	10 551	11 288	9038
一直保持负误差的栅格个数与 10 m DEM负误差栅格个数之比(%)	92.27	72.40	72.68	78.66	76.31	85.33	82.42	70.27

分辨率(m)	坡度分级编号
分辨率(m)	1	2	3	4	5	6	7	8	9
10	0.00	0.00	0.00	0.00	66.04	32.35	1.61	0.00	0.00
15	0.00	0.00	0.00	0.00	64.88	33.45	1.67	0.00	0.00
20	0.00	0.00	0.00	0.00	62.94	35.30	1.76	0.00	0.00
25	0.00	0.00	0.00	0.00	60.10	38.00	1.90	0.00	0.00
30	0.00	0.00	0.00	0.00	56.59	41.35	2.06	0.00	0.00
35	0.00	0.00	0.00	0.00	52.76	45.00	2.24	0.00	0.00
40	0.00	0.00	0.00	0.00	48.44	49.12	2.45	0.00	0.00
45	0.00	0.00	0.00	0.00	44.14	53.21	2.65	0.00	0.00
50	0.00	0.00	0.00	0.00	39.62	57.51	2.87	0.00	0.00

分辨率(m)	坡度分级编号
分辨率(m)	1	2	3	4	5	6	7	8	9
10	82.06	0.10	0.39	6.76	10.69	0.00	0.00	0.00	0.00
15	71.72	0.12	0.56	9.13	18.47	0.00	0.00	0.00	0.00
20	63.54	0.20	0.90	11.46	23.90	0.00	0.00	0.00	0.00
25	55.74	0.37	1.36	13.02	29.51	0.00	0.00	0.00	0.00
30	48.78	0.50	2.05	15.34	33.32	0.00	0.00	0.00	0.00
35	43.14	0.69	2.83	17.97	35.38	0.00	0.00	0.00	0.00
40	38.27	1.02	3.41	20.51	36.80	0.00	0.00	0.00	0.00
45	34.49	1.36	3.98	23.38	36.79	0.00	0.00	0.00	0.00