[发明专利]地表反射率计算方法和装置

权利要求书

1.一种用于HJ-1A/B卫星的地表反射率计算方法，其特征在于，包括步骤：

1)基于暗目标法或不变目标方法得到大气气溶胶光学厚度；

2)基于HJ-1A/B XML文件和影像数据计算每日周期的太阳天顶角、太阳方位角、观测天顶角和观测方位角；

3)基于上述步骤1)和2)中得到的参数，利用MODTRAN模拟星上幅亮度L_m；

4)通过步骤3)建立地表反射率查找表；

5)根据大气参数和待校正图像，利用查找表实现大气校正。

2.根据权利要求1所述的地表反射率计算方法，其特征在于，在步骤3)之前，还包括将气溶胶光学厚度转换为水平气象视距。

3.根据权利要求1所述的地表反射率计算方法，其特征在于，在步骤3)之前，还包括通过表观辐亮度及反射率的反演，得到大气上界的表观反射率。

4.根据权利要求1所述的地表反射率计算方法，其特征在于，在步骤4)中，查找表的自变量是大气模式、水平气象视距、太阳天顶角、卫星天顶角，因变量是求地表反射率的三个参数。

5.根据权利要求4所述的地表反射率计算方法，其特征在于，还包括步骤6)：参考MODIS大气校正算法中水平气象视距、太阳天顶角这两个变量的变化间距，然后根据MODTRAN运行结果来判定水平气象视距、太阳天顶角的变化间距上是否符合精度，若不符合精度，则进行插值。

6.根据权利要求1所述的地表反射率计算方法，其特征在于，通过输入HJ-1A/B CCD蓝光、红光和近红外地表反射率数据和高光谱反射率数据，实现植被指数的反演。 

7.根据权利要求1所述的地表反射率计算方法，其特征在于，所述生态参数反演模块通过植被指数和叶面积指数LAI，实现植被覆盖度的反演。

8.一种用于HJ-1A/B卫星的地表反射率计算装置，其特征在于，包括：

大气参数反演模块，其基于暗目标法或不变目标方法得到大气气溶胶光学厚度；

太阳及观测几何参数反演模块，其基于HJ-1A/B XML文件和影像数据计算每日周期的太阳天顶角、太阳方位角、观测天顶角和观测方位角；

地表反射率反演模块，其利用大气参数反演模块和太阳及观测几何参数反演模块计算得到的参数，利用MODTRAN模拟星上幅亮度L_m；并建立地表反射率查找表；然后根据大气参数和待校正图像，利用查找表实现大气校正。

9.根据权利要求8所述的地表反射率计算装置，其特征在于，还包括生态参数反演模块，其通过输入HJ-1A/B CCD蓝光、红光和近红外地表反射率数据和高光谱反射率数据，实现植被指数的反演。

10.根据权利要求9所述的地表反射率计算装置，其特征在于，所述生态参数反演模块通过土地利用/土地覆盖数据和HJ-1A/HJ-BCCD第3、第4地表反射率数据，实现叶面积指数LAI的反演。

11.根据权利要求10所述的地表反射率计算装置，其特征在于，所述生态参数反演模块通过植被指数和叶面积指数LAI，实现植被覆盖度的反演。

12.根据权利要求8所述的地表反射率计算装置，其特征在于，所述地表反射率反演模块还将气溶胶光学厚度转换为水平气象视距。 

13.根据权利要求12所述的地表反射率计算装置，其特征在于，所述查找表的自变量是大气模式、水平气象视距、太阳天顶角、卫星天顶角，因变量是求地表反射率的三个参数。

14.根据权利要求13所述的地表反射率计算装置，其特征在于，所述地表反射率反演模块参考MODIS大气校正算法中水平气象视距、太阳天顶角这两个变量的变化间距，然后根据MODTRAN运行结果来判定水平气象视距、太阳天顶角的变化间距上是否符合精度，若不符合精度，则进行插值。