[发明专利]基于环境一号卫星的陆地气溶胶光学厚度反演方法

权利要求书

1.一种基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，该方法包括步骤：

S1.根据环境一号HJ-1A/B星的CCD相机的特征，模拟各像元的观测几何、大气状况，进行辐射传输计算，获得相应的大气参数，存储于数据文件中形成查找表；

S2.根据环境一号HJ-1A/B星的CCD数据的辅助数据，提取各像元的观测天顶角以及观测方位角，并计算太阳天顶角以及太阳方位角，获得各像元的观测几何；

S3.对环境一号HJ-1A/B星的CCD数据进行如下预处理：辐射定标获得相机各波段的表观反射率、去除云覆盖像元、以及海陆分离提取陆地像元；

S4.利用暗目标算法获得各波段地表反射率信息，将所述查找表插值到各像元的观测几何，根据经过预处理得到的表观反射率计算不同陆地气溶胶光学厚度下的地表反射率，拟合得到陆地气溶胶光学厚度。

2.如权利要求1所述的基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，步骤S1进一步包括：

S1.1设定不同的观测情况：9个太阳天顶角，12个观测天顶角，16个相对方位角，气溶胶模式为大陆型气溶胶，相对于0.55μm波长处的气溶胶光学厚度设为6个等级，查找表计算的波段为环境一号HJ-1A/B星的CCD的蓝波段以及红波段，海拔设置为0m，所述相对方位角为太阳方位角与观测方位角之差；

S1.2调用辐射传输模式进行模拟运算，得到相应的大气程辐射、大气整层双程透过率、大气向下的半球反射率的三个参数，并存储在数据文件中形成查找表。

3.如权利要求2所述的基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，步骤S2进一步包括：

S2.1根据环境一号HJ-1A/B星的CCD数据中的几何定位信息，计算每个像元的经纬度，得到太阳天顶角和太阳方位角，太阳天顶角θ_s为：

cosθ_s＝sin(lat)sinδ+cos(lat)cosδcost

其中，lat为地理纬度，δ为太阳光与地球赤道平面的夹角，t为太阳的时角，定义纬地方时12点的时间为0，6点时为-π/2，18点时为π/2；

太阳方位角φ为：

cos(φ)=tan(π/2-θ)tan(lat)-sinδcos(π/2-θ)*cos(lat);]]>

S2.2根据环境一号HJ-1A/B星的CCD数据的辅助数据中的观测方位角以及所述太阳方位角，计算所述相对方位角。