基于多源高分辨率遥感影像的腾格里沙漠沙丘移动分析.docVIP

基于多源高分辨率遥感影像的腾格里沙漠沙丘移动分析 　　[摘 要]沙丘移动变化的研究是沙漠迁移过程研究的一个重要方面。地面观测手段难以满足大空间尺度范围的研究需要，而卫星遥感影像则为沙丘移动研究提供了长周期、大范围的数据保障。本文利用2010年ALOS、2013年国产高分一号和2014年WorldView-2共3期不同来源、不同分辨率的高分辨率遥感影像对腾格里沙漠沙丘移动进行研究分析，结果表明：（1）高分辨率遥感影像能有效反映1a尺度的沙丘变化研究；（2）研究区2010年～2013年间，沙丘平均位移4.36m；与研究区内风场风向基本一致；2013年至2014年间平均位移2.43m；平均移动方向指向方位角约36.5°SE；（3）个别沙丘局部表现出一定的西北―东南向的往复性运动，在2010年至2013年间样本沙丘最大逆向移动距离为2.55米。 　　[关键词]高分辨率遥感影像 几何配准 沙丘移动 　　中图分类号：TG333.7 文献标识码：A 文章编号：1009-914X（2016）02-03370-01 　　引言 　　沙漠扩张和沙漠化是人类面临的一个紧迫的环境难题，是社会经济发展中一个严重的制约因素。沙丘移动是判断沙漠扩张的一种重要依据。沙丘移动主要有野外观测法和遥感分析法。野外观测法，主要通过利用电子全站仪对沙丘形态测量及通过布设观测桩测定沙丘垂直断面的水平移动来确定移动距。该种方法最为直接，但是这种方法难以在大范围的空间上展开同步观测。遥感影像的出现克服了观测空间、时间的限制，给沙丘移动变化的观测带来了新的方式。遥感分析法，利用卫星、航空和雷达影像等分析对比沙丘在一定时期的变化状态。Bourke通过航空影像和雷达影像对维多利亚谷地1961年―2001年的沙丘变化进行了分析研究。陈芳利用Landsat TM多景影像对巴丹吉林沙漠的4个典型区域1989-2010年间的沙丘移动方向和距离以6a、8a、10a的时间尺度分别进行了实验分析，从而得到研究区沙丘移动的平均年移动距离。问题在于，中低分辨率遥感影像同样受空间分辨率的限制，无法做到更精确的量测。本文利用ALOS、中国高分一号、Worldview-2三种高分辨率遥感影像计算出腾格里沙漠东南缘区域沙丘在年际尺度上的移动速度和方向。 　　1.研究方法及数据准备 　　1.1.研究区概况 　　研究区位于腾格里沙漠东南边缘，属于内蒙古阿拉善左旗和宁夏中卫市的交界地带。地理坐标37°29′49.03″～37°33′15.40″N，105°01′33.14″～105°04′22.5″E。年平均空气湿度为36.51%～45.40%。年平均温度为11.21～12.65℃。年降水量71.44～116.60mm。年平均风速为2.48～2.79m/s，风向以西北风为主，其次为东北风。研究区南侧位于沙坡头自然保护区北部。2012年以来在北、东方向开始逐步进行光伏电场建设工程。研究区沙丘类型主要为仍然保持天然风貌的新月形沙丘链，至2015年尚未受到人类建设活动的干扰。 　　1.2 数据准备 　　1.2.1. 数据来源 　　遥感数据包括：（1）日本ALOS卫星10m分辨率的AVNIR-2多光谱数据和2.5m分辨率PRISM全色波段数据，1B2+RPC级产品；数据获取日期为2010年09月10日；（2）GF-1A级遥感数据，包括8m分辨率多光谱数据和2m分辨率全色波段数据；数据获取日期为2013年08月08日；（3）Worldview-2卫星2A级数据，包括2m分辨率多光谱影像数据和0.5m分辨率全色波段数据；数据获取日期为2014年10月04日。 　　1.2.2.影像配准 　　利用3种高分辨率影像数据所附带的RPC文件和GMTED2010数字高程模型数据进行正射影像校正。而后以Worldview-2影像为基准影像分别对ALOS和GF-1影像进行几何精纠正，使校正总精度保持在1个像元误差以内，单点的定位精度控制在0.5个像元误差以内，从而确保影像数据的配准精度。 　　1.2.3. 影像融合 　　影像融合是将低空间分辨率的多光谱图像或高光谱数据与高空间分辨率的单波段影像重采样，生成一幅高分辨率多光谱遥感影像的处理技术，使得处理后的影像既有较高的空间分辨率，又有多光谱特征。为便于更好的目视解译，本文应用ENVI5.1软件，对三期影像采用Gram-Schmidt方法进行融合。 　　1.2.4. 沙梁线的提取 　　沙丘由迎风坡和背风坡组成。沙丘在高分辨率遥感影像上表现出明显的纹理及灰度特性。两个坡面的交界线即是沙丘的沙梁线。分别在三期高分辨率遥感影像上勾绘沙梁线，实现对沙梁线的提取。 　　2.结果与讨论 　　沙梁的移动距离计算，采用了5点求平均法。即在每期影像的样本沙梁线上分别