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