[发明专利]基于TSEVI的森林覆盖变化检测方法

说明书

本发明提出一种基于TSEVI的森林覆盖变化检测方法，包括以下步骤：TSEVI优化计算、TSEVI标准化、ROC_TSEVI（TSEVI变化率）计算、ROC_TSEVI可视化、森林增长与损失提取。本发明利用两期以上光学遥感影像，无需其它外部异源数据辅助，就可以有效检测出复杂地形山区森林覆盖细微变化（包括森林增长和森林损失）。

技术领域

本发明涉及植被和环境监测领域，尤其涉及一种基于TSEVI的森林覆盖变化检测方法。

背景技术

现有森林覆盖变化检测方法主要有差值比较法、分类后比较法、空间模型法等。由于森林主要分布在复杂地形山区，地形阴影效应会影响森林覆盖信息的准确提取并最终影响变化检测结果。

发明内容

有鉴于此，本发明的目的在于提供一种基于TSEVI的森林覆盖变化检测方法，消除地形阴影效应对森林覆盖信息的干扰，并可以灵敏检测森林覆盖微小变化情况。

本发明具体采用以下技术方案：

一种基于TSEVI的森林覆盖变化检测方法，其特征在于，包括以下步骤：

步骤S1：选择至少两期光学遥感影像，计算TSEVI：

其中：TSEVI为地形阴影消除植被指数，B_nir为遥感影像近红外波段表观反射率数据，B_r为遥感影像红光波段表观反射率数据，f(Δ)为调节因子；

步骤S2：以光学遥感影像上完全无植被区域的TSEVI值作为极小值，以去除噪声的影像TSEVI最大值为极大值，进行TSEVI标准化计算：

TSEVI′＝(TSEVI-TSEVI_rock)/(TSEVI_max′-TSEVI_rock)其中：TSEVI’为标准化后的TSEVI值，TSEVI_rock为完全无植被区域的TSEVI值，TSEVI_max’为去除噪声的影像TSEVI最大值；

步骤S3：计算ROC_TSEVI，即TSEVI变化率：

ROC_y2-y1＝(TSEVI′_y2-TSEVI′_y1)/TSEVI′_y1

其中：ROC_y2-y1为TSEVI从时期y₁到时期y₂的变化率,TSEVI'_y1为时期y₁标准化后的TSEVI值，TSEVI'_y2为时期y₂标准化后的TSEVI值；

步骤S4：根据ROC_TSEVI的值获得森林覆盖变化情况。

优选地，在步骤S1中，调节因子f(Δ)采用匹配寻优法(中国专利200910111688X提供的现有技术方案)或极大值优化法(中国专利201010180895.3提供的现有技术方案)或相关系数法(中国专利2015108077580提供的现有技术方案)或变异系数法(中国专利2018112055412提供的技术方案)计算获得。

优选地，在步骤S2中，所述完全无植被区域包括岩石、沙滩区域；前0.001％的TSEVI最大值作为噪声剔除。

优选地，在步骤S4中，ROC_TSEVI的值为正代表森林覆盖增长；ROC_TSEVI的值为负代表森林覆盖损失。

优选地，步骤S4具体包括以下步骤：

步骤S41：将步骤S3计算获得的ROC_TSEVI与地图进行匹配出图；