[发明专利]一种计算高效的SEVI调节因子优化方法

权利要求书

1.一种计算高效的SEVI调节因子优化方法，其特征在于：包括以下步骤：

步骤S1：读取数据：读取一幅山地遥感图像，包括近红外波段反射率数据B_nir、红光波段反射率数据B_red，并记图像行列号为R×C；

步骤S2：控制窗口大小：设置初始窗口大小为A行B列，A∈(1,R)，B∈(1,C)；

步骤S3：预设中间变量：预设R行C列数值全是255的数组F_S与E_S0，同时预设S行S列的卷积核h_S，当A≥B时，S＝A，否则S＝B；

其中，F_S为预设包含所有窗口内f(Δ)_lopt的数组，f(Δ)_lopt为窗口内局部最优调节因子；E_S0为预设包含所有窗口内误差的数组；

步骤S4：植被指数计算：计算整幅图像的比值植被指数RVI、阴影植被指数SVI和阴影消除植被指数SEVI；

步骤S5：以h_S为卷积核对整幅图像计算出来的RVI、SVI、SEVI、RVI²、SVI²、SEVI²、(SEVI×RVI)、(SEVI×SVI)执行卷积操作，并计算当前f(Δ)对应下的E_S(f(Δ))；

其中，RVI²为RVI的平方、SVI²为SVI的平方、SEVI²为SEVI的平方、(SEVI×RVI)为SEVI与RVI的乘积、(SEVI×SVI)为SEVI与SVI的乘积，f(Δ)为调节因子，E_S(f(Δ))为当前f(Δ)包含所有窗口内误差的数组；

步骤S6：逼近优化，令f(Δ)从0到1，以L₁为间隔，依次递增迭代更新F_S与E_S0；

步骤S7：全局最优值确定，选用F_S中所有f(Δ)_lopt的平均值作为研究区遥感图像的全局最优调节因子f(Δ)_gopt。

2.根据权利要求1所述的一种计算高效的SEVI调节因子优化方法，其特征在于：步骤S3中所述F_S、E_S0、h_S为以下数组形式：

当A＝B时，窗口为正方形，h_S为元素全是1的数组：

当A＞B时，窗口为行数大于列数的矩形，h_S的第B+1列至S列的元素为0：

当A＜B时，窗口为行数小于列数的矩形，h_S的第A+1行至S行的元素为0：

3.根据权利要求1所述的一种计算高效的SEVI调节因子优化方法，其特征在于：步骤S4中，所述阴影消除植被指数SEVI、比值植被指数RVI和阴影植被指数SVI的计算公式为：

其中，f(Δ)为调节因子，B_red为遥感图像红光波段反射率数据，B_nir为遥感图像近红外波段反射率数据。