[发明专利]一种遥感图像大气订正方法

权利要求书

1.一种遥感图像大气订正方法，具体为：1)查找表构建；2)表观反射率计算；3)角度信息获取；4)气溶胶光学厚度获取：基于纯像元方法利用HJ-1卫星CCD数据反演得到；5)由像元成像几何及光学厚度插值得到四个大气校正参数；6)逐像素进行大气校正，得到地表反射率。

2.如权利要求1所述的遥感图像大气订正方法，其特征在于，步骤1)中采用6S辐射传输模型构建关于大气校正参数R_r+a、T_g(θ_s，θ_v)、T(θ_s)·T(θ_v)的查找表。

3.如权利要求2所述的遥感图像大气订正方法，其特征在于，步骤2)根据定标系数将影像DN值转换为表观反射率。

4.如权利要求3所述的遥感图像大气订正方法，其特征在于，步骤3)中太阳天顶角、方位角通过查询数据描述文件XML文件得到；观测天顶角和方位角从Sat_Zenith_Azimuth.txt文件中得到。

5.如权利要求4所述的遥感图像大气订正方法，其特征在于，步骤4)中的反演方法具体为：a)根据辐射传输模型建立查找表；b)依据纯像元指数提取遥感图像中的纯像元，生成纯像元掩膜；c)确定纯像元地表反射率；d)根据查找表，将遥感观测的辐射亮度反演为气溶胶光学厚度。

6.如权利要求5所述的遥感图像大气订正方法，其特征在于，该反演方法还包括步骤e)对整幅影像进行内插，得到光学厚度图。

7.如权利要求6所述的遥感图像大气订正方法，其特征在于，步骤a)采用6S辐射传输程序进行辐射传输计算得出多组不同的大气参数R_r+a、T_g(θ_s，θ_v)、T(θ_s)·T(θ_v)和S组合而成的查找表，R_r+a是由分子散射加气溶胶散射所构成的路径程辐射反射率，T(θ_s)、T(θ_v)分别为太阳——目标、目标——传感器大气路径透射率，S为大气下届半球反照率。

8.如权利要求7所述的遥感图像大气订正方法，其特征在于，步骤c)具体为：选取一个时间段内的晴朗日影像和污染日影像，假设在这个时间段内纯像元地表反射率变化较小，对晴朗天气的影像进行大气校正，获取纯像元的地表反射率。

9.如权利要求8所述的遥感图像大气订正方法，其特征在于，步骤d)具体为：根据读取的太阳天顶角、太阳方位角、观测天顶角、观测方位角，在查找表选取相应的数据，进行线性插值，得到不同波段、不同气溶胶光学厚度的R_r+a、T_g(θ_s，θ_v)、T(θ_s)·T(θ_v)和S大气参数；结合步骤c)得到的纯像元地表反射率，按照公式

R*(θs,θv,φs-φv)=Tg(θs,θv)[Rr+a(θs,θv,φs-φv)+T(θs)T(θv)R1-RS]]]>计算每个预设气溶胶光学厚度对应的表观反射率，与卫星真实测量的表观反射率进行比较，相差最小的气溶胶光学厚度值即为所求，θ_s是太阳天顶角，φ_s传感器天顶角，θ_v是太阳方位角，φ_v是传感器方位角，T_g是气体吸收透过率，R_r+a是由分子散射加气溶胶散射所构成的路径程辐射反射率，T(θ_s)、T(θ_v)分别为太阳——目标、目标——传感器大气路径透射率，S为大气下届半球反照率，R是地物目标反射率。