基于泰森多边形算法的冰川长度自动提取方法研究 ——以三江源区为例.pdf

摘要 摘要 冰川长度是表征冰川几何形状的重要指标之一，然而冰川编目数据集大多都缺失 长度信息。传统的使用人工测量的方法获取冰川长度，费时费力。因此，迫切需要自 动化的提取方法来计算冰川长度。近几年，国外一些学者提出了几种冰川长度自动提 取方法。然而，国内尚无较成熟的自动提取冰川长度的方法。为此本文利用卫星资料 提出了一种自动提取冰川中心线的方法，并与其他冰川中心线提取方法进行比较分析， 同时将该方法应用于三江源地区冰川并对三江源地区冰川长度分布及变化特征进行 分析。 本文提出了一种基于泰森多边形理论的冰川中心线自动提取方法。该方法以冰川 轮廓数据和数字高程模型（DEM ）作为输入数据，主要分为三步，确定冰川的最低点 和局部最高点、构建泰森多边形、路径的选择与平滑处理。将该方法应用于阿拉斯加 南部的 1125 条冰川上，同现有研究中的两种冰川中心线自动提取方法进行了比较， 对比分析了三种提取方法的平均长度和不同方法之间的长度比。并利用本文提出的冰 川中心线自动提取方法对三江源地区2009 年和2016 年冰川长度信息进行提取，对三 江源地区冰川长度分布及变化特征进行分析。得到主要结论有： （1）三种方法提取的阿拉斯加南部 Stikine 冰原的冰川长度的平均值相差不大， 在 1.789-1.833 km2 之间。1125 条冰川的 R_(x/m) 的平均值为 1.05，1125 条冰川的 R_(x/k) 的平均值为1.01 ，与完全一致（R=1.00 ）相比，三种方法之间存在较小的正偏 差。92 % 的冰川长度比R_(x/k)在0.9-1.1 之间，84.2 %的冰川长度比R_(x/m) 在0.9- 1.1 之间。本文所提出的基于泰森多边形算法的冰川长度自动提取方法与其他自动提 取方法具有可比性，提取结果可信。 （2 ）三江源区冰川长度分布特征。2009 年和 2016 年，三江源地区的冰川平均 长度随冰川规模等级的递增而增加。2009-2016 年，随时间变化，不同规模等级的冰 川平均长度都在减小。2009 年平均长度为 1.48 km ，2016 年平均长度为 1.44 km ，这 也从侧面反映了研究区冰川萎缩的事实。将三江源地区按不同的山系进行空间划分， 研究区可分为阿尼玛卿山、曲阿加吉玛山峰地区、唐古拉山、巴颜喀拉山和昆仑山五 个山系，其中巴颜喀拉山没有分布冰川。阿尼玛卿山、曲阿加吉玛山峰地区、唐古拉 I 西北大学硕士学位论文 山和昆仑山的冰川平均长度分别为1.31 km 、0.91 km 、1.50 km 、1.80 km 。 （3 ）三江源区冰川长度变化特征。各规模等级冰川的面积和长度整体上都存在 不同程度的萎缩。2009-2016 年冰川面积年萎缩率为0.43 % ，冰川平均长度每年退缩 5 m 。对研究区阿尼玛卿山、曲阿加吉玛山峰地区、唐古拉山和昆仑山四个区域的冰 川长度进行统计分析。2009-2016 年，三江源区四个山系的冰川平均长度都减小了。 2009-2016 年，研究区按山系划分的四个区域冰川面积萎缩率和冰川长度变化速率一 一对应。其中，研究区曲阿加吉玛山峰地区的冰川长度变化速率最快，每年退缩8 m ， 其面积萎缩率为8.35 % ，在四个区域中冰川面积萎缩最快的。 （4 ）冰川面积和长度之间存在一定的相关性。由经验方程 = • , (k = [ ] [ ] 0.92, 1.70 , β = 0.467,0.606 ) ，得出2009 年和2016 年三江源地区冰川长度和面积 的经验公式分别为 = 1.65 0.53和 = 1.620.53 。 关键词：冰川， 长度， 中心线， 泰森多边形， 三江源区 II 目录 目录 摘要