[发明专利]一种气溶胶光学厚度反演方法

说明书

本申请提供了一种气溶胶光学厚度反演方法，涉及遥感图像处理技术领域，包括：对获取的静止卫星的遥感影像进行预处理，去除云像元、水像元以及冰雪像元，得到影像AOD反演区域；针对反演区域的各暗像元，通过对短波红外波段、红光波段和蓝光波段进行地气解耦，得到暗像元的红光地表反射率和蓝光地表反射率；读取可见红光波段查找表和蓝光波段查找表，通过求解暗像元的第一代价函数最小值，获取暗像元的最优气溶胶光学厚度；针对反演区域的各亮像元，基于地表反射率库计算亮像元的蓝光地表反射率；读取蓝光波段查找表，通过求解亮像元的第二代价函数最小值，获取亮像元的最优气溶胶光学厚度。从而能够提高静止卫星遥感影像的气溶胶光学厚度反演精度。

技术领域

本申请涉及遥感图像处理技术领域，尤其是涉及一种气溶胶光学厚度反演方法。

背景技术

气溶胶对全球大气环境、气候变化以及人类健康具有重要影响，气溶胶光学厚度（Aerosol Optical Depth ，AOD）是气溶胶最重要的参数之一，也是评价大气污染程度的关键因子。AOD的监测有两种途径：地基观测和卫星反演。地基监测仪器价格昂贵，并且仪器的操作和维护过程较为繁琐，因此地基观测能力十分有限，无法满足长时序、高频次、大范围动态监测AOD的要求。而卫星反演的AOD可以有效满足以上要求。

目前应用最广的极轨卫星MODIS_AOD产品主要关注区域的某一时刻的瞬时AOD，不具备观测区域全天的AOD分布和变化情况的能力。利用静止气象卫星数据反演的AOD可以大幅度提高观测频率，产品虽然能达到一定的精度水平，但受限于静止卫星光谱分辨率、空间分辨率以及定标精度较低等因素，很难对其反演的误差进行控制。

发明内容

有鉴于此，本申请提供了气溶胶光学厚度反演方法，能够解决现有卫星反演气溶胶光学厚度存在精度不高的技术问题。

本申请实施例提供了一种气溶胶光学厚度反演方法，包括：

对获取的静止卫星的遥感影像进行预处理，去除云像元、水像元以及冰雪像元，得到影像AOD反演区域；

对于影像AOD反演区域的各暗像元，通过对短波红外波段、红光波段和蓝光波段进行地气解耦，得到所述暗像元的红光地表反射率和蓝光地表反射率，由此计算模拟红光表观反射率和模拟蓝光表观反射率；读取预先建立的可见红光波段查找表和蓝光波段查找表，通过求解所述暗像元的第一代价函数最小值，获取所述暗像元的最优气溶胶光学厚度；所述第一代价函数值根据所述暗像元的实际红光表观反射率、实际蓝光表观反射率、模拟红光表观反射率和模拟蓝光表观反射率计算得到；

对于影像AOD反演区域的各亮像元，基于预先建立的地表反射率库计算所述亮像元的蓝光地表反射率，由此计算模拟蓝光表观反射率；读取所述蓝光波段查找表，通过求解所述亮像元的第二代价函数最小值，获取所述亮像元的最优气溶胶光学厚度，所述第二代价函数值根据所述亮像元的模拟蓝光表观反射率以及实际蓝光表观反射率计算得到。

进一步的，所述方法还包括：构建基于6S辐射传输模型的可见红光波段查找表和蓝光波段查找表，具体为：

在预设时间段内，调用6S辐射传输模型，在可见红光波段内，建立针对所有可能的太阳天顶角、卫星观测角、相对方位角和气溶胶光学厚度情况下得到的大气半球反射率、太阳-地面的大气路径透过率、地面-传感器的大气路径透过率和大气程辐射反射率，形成可见红光波段查找表；

在预设时间段内，调用6S辐射传输模型，在蓝光波段内，建立针对所有可能的太阳天顶角、卫星观测角、相对方位角、气溶胶光学厚度情况下得到的大气半球反射率、太阳-地面的大气路径透过率、地面-传感器的大气路径透过率和大气程辐射反射率，形成蓝光波段查找表。

进一步的，对遥感影像进行预处理，去除云像元、水像元以及冰雪像元，得到影像AOD反演区域；包括：

将遥感影像从RGB空间转换到HSI空间，构建遥感影像的显著性图像；