[发明专利]一种基于遥感影像的草地退化程度自动提取方法

说明书

技术领域

本发明涉及草原退化监测领域，尤其涉及一种基于遥感影像的草地退化程度自动提取方法。

背景技术

近年来，由于气候变化及人为干扰等原因，导致草地出现不同程度的退化，草地生产力明显下降，而对退化情况进行及时的诊断与提取是草地退化管理的重要前提，有助于针对不同区域退化情况的分布特征实施相应的恢复措施及管理手段。

目前的研究方法主要有目视综合判读法和结合实测数据反演草地退化指标法。目视综合判读法以遥感图像计算机屏幕人机交互判读为核心，保证了高精度，但对于解译人员和遥感影像数据的选取、处理等要求严格，工作量最大。而退化指标法中指标主要有盖度、生物量、净第一生产力(NPP)等生物学指标和土壤含水量、水土保持能力等地表特征指标，可以通过实测数据与退化指标拟合直接反演整个区域的退化情况，省时省力。但该方法中需要较多的实测数据来与相应的指标进行拟合，同时存在经验参数较多，退化指标不统一，参照系统模糊等问题，故反演精度和推广性有限。

因此，在认识到现有的提取方法客观性、合理性上方面所具有的局限性的基础上，鉴于草地退化所呈现的变化特征，建立了一种基于遥感影像的草地退化程度自动提取方法，以快速、准确的提取不同退化程度的草地。

发明内容

本发明的目的在于提供一种基于遥感影像的草地退化程度自动提取方法。

本发明的基于遥感影像的草地退化程度自动提取方法，包括如下步骤：

步骤(1)数据获取：获取待监测区域的卫星遥感影像数据，包括两景前后时相的原始影像数据；

步骤(2)数据预处理：对各个波段的光谱影像进行目视检查，确定由于大气或其他因素对光谱产生较大噪声的异常波段，对异常波段直接剔除，之后对影像数据依次进行辐射定标、镶嵌、裁剪和大气校正；

步骤(3)最小噪声分离变换：将预处理后的前后时相影像数据进行叠加处理，对叠加后的影像进行最小噪声分离变换(MNF变换)；

步骤(4)多尺度分割：对步骤(3)处理后得到的数据进行多尺度分割，根据设定的波段数、光谱与形状标准参数，选择20，25，30，35，40，45，50，55，60，65，70，75，80，85，90共15个分割尺度对影像进行分割，最后目视判定观察对象内光谱均一性最好的则为最优分割尺度；

步骤(5)在最优分割尺度下，对于每一个对象，如果光谱特征阈值满足逻辑或条件，则属于非草地对象集NG，则剩余对象集G均为草地，这样可以剔除非草地的影响，从而提高分类精度。

步骤(6)在影像对象集G上，根据先验知识，在每个退化类别区域内均匀选取一定比例的训练样本，通过最近邻分类器，结合光谱、形状和纹理特征提取得到不同退化程度的草地，其中包含未退化、轻度退化、中度退化和严重退化等4种退化类型。

步骤(3)所述的MNF变换具体实现步骤如下：

G.把混有噪声的信号进行高通滤波处理，得到噪声协方差矩阵∑_n。

H.计算噪声协方差矩阵∑_n的特征值λ和特征向量P，由λ组成的对角矩阵记为Δ，由P组成的矩阵记为E。

I.求∑_n的白化矩阵F：

F＝EΔ^-1/2

J.计算噪声白化之后的观测数据F^TZ协方差矩阵:

∑_w＝cov(F^TZ)＝F^T∑_zF

其中，Z为影像数据，∑_z为影像的协方差矩阵。

K.由∑_w求得特征值和特征向量，由所有特征值构成的对角矩阵记为Λ，相应特征向量构成的矩阵记为G。

L.将G作用到F^TZ的每个象元上，G^TF^TZ即为MNF变换后得到的图像。

步骤(4)所述的波段数、光谱与形状标准参数的设定，主要是通过专家知识结合目视观察获得的。

步骤(5)所述的逻辑或条件为：设定NDVI阈值为0.2-0.3，MNF处理后第四波段均值阈值为8-9，当原始影像波段NDVI值均小于NDVI阈值，MNF处理后第四波段均值均大于上述均值阈值，则将该对象归入非草地，通过NDVI或MNF变换第四波段均值的阈值调节可以剔除多数非草地对象，提高分类精度。