[发明专利]一种遥感地表温度产品角度效应校正方法

权利要求书

1.一种遥感地表温度产品角度效应校正方法，其特征在于，包括：

步骤1：计算像元MODIS 31,32波段的亮度温度T31、T32；

步骤2：获得所述像元大气下界温度T_air、大气总水蒸汽含量CWV、传感器观测天顶角θv；

步骤3：使用像元大气下界温度T_air、大气总水蒸汽含量CWV、传感器观测天顶角θv从劈窗算法系数查找表中获得所述像元传感器观测方向通用劈窗算法系数C，A₁，A₂，A₃，B₁，B₂，B₃；

步骤4：确定所述像元的天顶发射率ε₀；

步骤5：计算所述像元的坡度α和坡向

步骤6：计算本地反射角θ₀；

步骤7：获得地表不同覆盖类型纯像元在传感器观测角度﹣65°到﹢65°之间的131个观测角度每个角度的平均发射率数据集；

步骤8：将步骤7中结果都减去获得到不同观测天顶角相对天顶方向的发射率差数据集；

步骤9：将所述像元在观测方向相对天顶方向发射率变异的表达式假设为：

Vε(θ)＝a₁θ+a₂θ²+a₃θ³+...+a_nθⁿ

使用数据集确定多项式中的系数a₁，a₂，a₃，...a_n以及n的值；

步骤10：确定所述像元的在所述观测方向的发射率表达式；

步骤11：确定所述观测方向的发射率，结合亮度温度，劈窗算法系数计算地表温度T_s。

2.根据权利要求1所述的方法，其特征在于，在步骤1中使用MOD021KM数据和普朗克函数近似公式计算所述像元MODIS 31,32波段的亮度温度T31、T32，

T31＝1304.413871/ln(1+729.541636/b31)

T32＝1196.978785/ln(1+474.684780/b32)。

3.根据权利要求1所述的方法，其特征在于，在步骤2中，使用MOD07_L2数据获得所述像元大气下界温度T_air、大气总水蒸汽含量CWV，以及从MOD021KM数据中获得传感器观测天顶角θv。

4.根据权利要求1所述的方法，其特征在于，在步骤4中，使用MOD21A22001年至2010年共十年数据计算所述像元十年历史月平均发射率，将结果作为地表天顶发射率数据集，根据所述像元获取月份确定所述像元的天顶发射率ε₀。

5.根据权利要求1所述的方法，其特征在于，在步骤5中，使用NOAA GLOBE 1KM DEM数据和ArcMap软件中aspect和slope功能计算所述像元的坡度α和坡向

6.根据权利要求1所述的方法，其特征在于，在步骤6中，使用如下公式计算本地反射角

其中，式中φ_s，θv分别是传感器观测方位角和观测天顶角。

7.根据权利要求1所述的方法，其特征在于，在步骤7中，使用2006年至2010年的MODIS地表覆盖类型数据MOD12Q1和MODIS地表发射率数据MOD11B1，地表覆盖数据经过聚合升尺度为和发射率数据相同的分辨率，获得地表不同覆盖类型纯像元在传感器观测角度﹣65°到﹢65°之间的131个观测角度每个角度的平均发射率数据集。