基于Lansat―8LI影像的南伊内里切克冰川运动速度提取与分析.docVIP

基于Landsat — 8 OLI影像的南伊内里切克冰川运动速度提取与分析 当前全球气候变化问题被日益关注，冰川作为气候变化的敏 感指示器之一，对分析气候变化过程有重要科学意义。冰川表面运动 是冰川的一个重要基本体征，它能够及时快速地反映冰川动态变化信 息，为研究冰川对气候变化的响应提供了有效途径。本次研究以天山 南伊内里切克冰川为研究对象，利用Landsat-8 OLI卫星影像，采用 归一化互相关算法对2013-2016年的四期遥感影像进行分析处理，得 到冰川表面的速度时空分布。通过分析得出，该冰川具有以下运动特 征：冰川的轴部是冰川运动的主流线，在冰川轴线上，速度呈现先增 大后减小的趋势，最大速度出现在平衡线附近；流速在轴部最大，在 两侧呈现减小的趋势；2013-2016年间，该冰川主体流速约为34cm/d， 呈现出较稳定的态势。 ，供学习和研宂 使用， 和著作权归原 所有， 关键词：光学遥感；相关性分析；COSI-Corr；冰川流速 P343.6 A 2095-2945(2017) 20-0008-03 1概述 全球气候变化是当今科学研究的热点问题。冰川作为冰冻圈系统 的主要组成部分之一，也是天然的气候变化指示器，为了解水资源现 状提供了研究途径［1］。迄今为止，前人的研究区域大多集中在南极和 格陵兰两大冰盖等地区［2-4］，但只右少数研究是关于内陆山岳冰川 的。谢自楚等在珠穆朗玛峰地区科学考察报告中指出，我国大多数的 冰川数据资料都是于1960-1980年间获取的［5］;王欣等通过进一步研 究，发现目前我国对于冰川的连续监测，仅局限于天山乌鲁木齐河源 1号冰川等少数几条冰川［6］。而对于其他的冰川覆盖地区，则缺乏对 它们运动的深入研究。 目前，冰川表面运动速度的提取主要有两种方式。一种是实地花 杆数据测量［7］，这种方法的精度高，但是受时间和空间的局限性比较 大;另一种方法则基于遥感影像进行提取［8］，相比传统野外实测方法， 遥感获取冰川表面流速具有快速和成本低的优势，对于大面积的以及 人类难以踏足的冰川地区，利用遥感影像的方式获取冰川表面速度的 优势更加明显。 2研究区介绍 托木尔峰地区南伊内里切克冰川（South Inilchek Glacier，中 国冰川编目编号：5Y673K1； 42° 12. 5Z N, 80° 12.2 K）是我国典型 的树枝状山谷冰川，冰川长约60. 5千米，面积约567. 2平方千米，是 中低纬度区域长度超过50千米的八大冰川之一，图1是研究区域示意 图。 3研究方法 3. 1基本原理 对不同时期天山冰川的遥感影像进行相关性分析，进而提取同名 像点的位移。经过图像预处理之后，可以去除太阳辐射、轨道偏差、 地形影响等误差。再利用归一化互相关算法得到的同名像点的东西向 和南北向位移，经过合成，就可以得到冰川表面的运动总位移，可以 认为就是冰川运动的结果［9］，结合时间可以得到该阶段内冰川的曰平 均速度。 3.2相关性分析及算法介绍 归一化互相关算法（NCC，Normalized Cross Correlation） ［10］， 也被称为泊松系数法或者归一化相关系数法。对于归一化互相关算法 的表达式定义如下： 其中，f （X，y）搜索影像的搜索窗n，g （x，y）是原始影像的 模板窗口，u，V是坐标的偏移景，U，V）是两个窗口各自的灰度平均 值。如图2所示，m, n为模板窗口的长和宽，M, N是搜索窗口的长和 宽，搜索窗口一般大于模板窗口。归一化处理能够增强该算法的鲁棒 性，减少其他因素对相关系数带来的影响，同时降低相关系数范围对 像特征的依赖［11］ 在矩阵NCC （u，v）中，值为1的点表示搜索窗口中该点及其周 围像素与模板窗U完全相同；值为-1的点表示搜索窗U中该点及周围 像素与模板窗口完全相反。所以，通过选取值最大的点即可初步获得 特征点的同名点，之后可以根据图像的坐标信息，计算出关键点的位 移［12］。 4数据选择与处理 与其他遥感影像数据相比，Landsat系列数据价格较低，获取途 径多，分辨率较高时问序列较长的Landsat影像适合大区域长周期的 冰川变化研究［13］。本次研究选用了 4组Landsat-8 OLI影像，用来 反演冰川的表面速度。 孔繁司［14］等在对四款常用的光学影像运动软件进行比较后发 现，⑶Sl-Corr?件作为一款基于图像配准和相关系数计算的软件，该 软件的配准精度可达到0. 1个像元，因而获取到的冰川运动速度精度 较高，而且花费的时间也较短。在计算最大相关系数时，采用频率域 算法，参考窗口设为128,搜索窗口设为32。将位移数据分成了 3个 图层：东西向位移、南北向位移以及信噪比(Signal Noise Ratio, SNR)。为