基于时序多光谱影像的干旱草原区开采扰动信息提取方法

露天开采会彻底改变原有土地利用景观格局，直接破坏当地生态环境，甚至会影响附近居民的生产和生活，因此越来越多的学者开始关注开采扰动。先前有关利用时序多光谱影像提取开采扰动的研究区集中于扰动形式单一的森林区。而我国露天煤矿大多集中于草原区，且我国东北部的草原矿区因其脆弱的生态环境以及其他多种扰动形式的存在，使得开采扰动识别更加困难。为明确我国东北部生态脆弱区草原露采场的开采扰动，以胜利矿区为例，利用1986年-2017年27期Landsat多光谱遥感影像，基于归一化植被指数NDVI（normalized difference vegetation index）的长时间序列轨迹变化特征（为了去除物候、云和阴影等对时序多光谱影像的影响，利用BISE-WT滤波器对原始NDVI时间序列进行滤波处理, 有效地去除时序NDVI数据中的噪声并同时保留有效信息），经过样本点训练，获得CV阈值（变异系数coefficient of variation）和Max阈值（植被阈值），构建CV-Max扰动识别模型，提取研究区的扰动分布。并利用植被阈值，分析NDVI时序轨迹，获得扰动年际信息，重构扰动历史地图；进而通过分析研究区典型地物的光谱特征，构建裸煤提取规则，以此来提取研究区的裸煤分布；最后通过构建裸煤及扰动区两者间的拓扑关系，进行空间拓扑叠置分析，从而获得开采扰动信息。经过精度验证，开采扰动的提取精度达到93.17%（Kappa系数=0.85），扰动年际信息提取精度达到83.35%（Kappa系数=0.81）。结果表明：在研究期间，空间上，开采扰动面积占研究区总面积的8.90%；时间上，开采扰动的发生集中于2000年-2009年，期间开采扰动像元占开采扰动总像元的76.70%；1988年-1998年矿区属于土地损毁初始期，2000年-2005年矿区属于土地损毁加速期，2006年-2009年矿区属于土地损毁高峰期，2010年-2017年开采扰动像元占比趋势比较平缓且持续处于较低水平，矿区土地损毁范围基本稳定。所提出的针对我国东北部生态脆弱性草原矿区，基于时序多光谱影像，利用植被指数NDVI和裸煤光谱特征提取开采扰动信息的方法是可行的，该研究结果可为干旱、半干旱草原露天矿区的可持续发展提供数据和理论方法支撑。