[发明专利]一种关于新TAVI模型中s因子的算法

摘要

本发明涉及一种关于新TAVI模型中s因子的算法。该方法包括以下步骤影像分析、计算新TAVI、计算太阳入射角余弦值cosi、计算TAVI与cosi相关系数r、计算有效区间集、计算s因子有效区间集、计算s因子结果。本发明的新TAVI模型中s因子的算法，弥补了目前s因子没有计算方法的不足，为基于太阳高度角的新TAVI模型实际应用提供了有效的计算方法，对新TAVI从经验模型向物理模型演进具有重要的科学意义与实用价值。

主权项

一种关于新TAVI模型中s因子的算法，其特征在于：包括如下步骤，步骤S1、影像分析：选择研究区，收集研究区时间序列遥感影像和坡度、坡向数据；对遥感影像进行包括辐射定标、几何校正的预处理，生成研究区影像表观反射率数据集；对影像数据集进行质量分析，统计遥感影像红光波段和近红外波段表观反射率均值、方差、极大值、极小值；分析山区植被在红光波段和近红外波段的反射率数据是否合理，决定影像是否正常可用；步骤S2、采用以下公式计算新TAVI：TAVI=BnirBr+f(Δ)·1Br]]>f(Δ)＝s‑sin(α)式中，Bnir为遥感影像近红外波段表观反射率数据，Br为遥感影像红光波段表观反射率数据，f(Δ)为调节因子，s为传感器因子，α为卫星过境时太阳高度角；步骤S3、采用以下公式计算太阳入射角余弦值cosi：cosi＝cosδ·cosθ+sinδ·sinθ·cos(β‑ω)式中，i为太阳入射角，δ是坡度角，β是坡向角，θ、ω是卫星过境时的太阳天顶角与太阳方位角；δ和β可以从坡度、坡向图中获取，θ和ω可从遥感影像头文件中获取；步骤S4、采用以下公式计算研究区时序影像TAVI与cosi相关系数r：r=nΣxy-Σx·ΣynΣx2-(Σx)2·nΣy2-(Σy)2]]>式中，r为TAVI与cosi的相关系数，x、y分别为所述样区对应时相的TAVI、cosi数据，n为影像数据的像元数；步骤S5、计算f(Δ)有效区间集：假设r位于(‑0.05,+0.05)时，TAVI计算结果满足精度要求；令f(Δ)从0开始，以0.001为步长，依次递增，计算时序影像的TAVI，并计算TAVI与cosi的相关系数r，当r＝0.05时，停止计算，得到相应f(Δ)值f1及其集合F1；再令f(Δ)从1.0开始，以0.001为步长，依次递减，计算TAVI，并计算TAVI与cosi的相关系数r，当r＝‑0.05时，停止计算，得到相应f(Δ)值f2及其集合F2；步骤S6、采用以下公式计算s因子有效区间集S1和S2：s＝f(Δ)+sin(α)式中，s为传感器因子；α为卫星过境时太阳高度角，从遥感影像头文件读取；f(Δ)为调节因子，根据步骤S5计算结果取相应数据集F1和F2；计算结果得到对应f(Δ)数据集F1和F2的s因子s1及其集合S1和s2及其集合S2；步骤S7、计算s因子结果：计算S1的极大值MAXS1和S2的极小值MINS2，当MAXS1<＝MINS2时得到结果s∈(MAXS1,MINS2)；步骤S8、如果MAXS1>MINS2,返回步骤S5调整r可接受精度范围，重新计算步骤S5到步骤S7；r原则上不超出(‑0.10,+0.10)；否则，剔除异常遥感影像重新计算。