FY-4A气溶胶光学厚度的反演方法、装置、设备和介质制造方法及图纸

本申请提供了一种FY

【技术实现步骤摘要】
FY
‑
4A气溶胶光学厚度的反演方法、装置、设备和介质


[0001]本申请涉及卫星
，尤其是涉及一种FY
‑
4A气溶胶光学厚度的反演方法、装置、设备和介质。

技术介绍

[0002]目前现有的遥感卫星反演气溶胶光学厚度AOD(Aerosol Optical Depth)的方法主要分为两类，深蓝算法和暗像元算法。其中，暗像元AOD反演中会通过建立蓝(0.47μm)、红波段(0.64μm)与受大气影响非常小的近红外波段(2.12μm)的线性关系来获取前两个波段的地表反射率，比如通过两个受大气影响非常低的近红外波段计算NDWI
SWIR
，然后利用红、蓝波段与近红外波段的线性关系、加上与NDWI
SWIR
的关系，来确定两个波段在某个像元的地表反射率，再利用大气传输模型计算该像元的AOD。也有研究直接复用了MODIS现有的线性关系函数，利用其他卫星传感器进行反演。
[0003]目前针对风云四号A星(FY
‑
4A)开发的AOD反演研究较少，官方也无正式公布的AOD产品。对于现有的暗像元算法，估算地表反射率的方法主要依赖MODIS都算法，即寻找红、蓝波段与近红波段的线性关系，更有研究直接使用了MODIS地表反射率产品如MOD09GQ、MOD34A3等。然而不同的传感器在波段设置上不尽相同，相同的反演方法无法在不同传感器之间简单地复用。即使FY
‑
4A搭载的多通道扫描成像辐射计AGRI(Advanced Geosynchronous Radiation Imager)在红、蓝波段的中心波长设置上与MODIS极为相近，他们的光谱响应函数也有巨大的区别，因此无法直接复用MODIS上的线性关系套用在FY
‑
4A的AOD反演上。

技术实现思路

[0004]本申请的目的在于提供一种FY
‑
4A气溶胶光学厚度的反演方法、装置、设备和介质，减少了因不同卫星传感器而造成的地表反射率估算错误，提升了FY
‑
4A气溶胶光学厚度的反演的准确性，提升了实时气溶胶光学厚度的生产效率和速度。
[0005]第一方面，本专利技术提供一种FY
‑
4A气溶胶光学厚度的反演方法，方法包括：获取FY
‑
4A预设时段的全圆盘数据，并读取指定波段下的实测数据；基于预先配置的气溶胶光学厚度查找表对指定波段下的实测数据进行大气校正，得到指定波段对应的目标地表反射率；其中，预先配置的气溶胶光学厚度查找表为配置有FY
‑
4A光谱响应函数的6S大气校正查找表；根据指定波段对应的目标地表反射率进行气溶胶光厚度的多进程反演处理，得到指定波段对应的目标气溶胶光学厚度。
[0006]在可选的实施方式中，预先配置的气溶胶光学厚度查找表包括预设个数的太阳高度角、预设个数的太阳天顶角、预设个数的相对方位角与多种气溶胶光学厚度进行组合后的组合结果；其中，每种组合都对应有相应的大气程辐射、总透射率和半球反照率。
[0007]在可选的实施方式中，指定波段包括第2波段和第3波段，实测数据包括卫星图像；基于预先配置的气溶胶光学厚度查找表对指定波段下的实测数据进行大气校正，得到指定
波段对应的目标地表反射率，包括：对卫星图像中进行暗像元判断；基于预先配置的气溶胶光学厚度查找表和第一大气传输模型确定第3波段对应的多种气溶胶光学厚度分别对应的第一目标地表反射率；第一大气传输模型为：；其中，为地表反射率；为太阳天顶角；为卫星天顶角；为方位角；为表观反射率；为大气程辐射、为总透射率，为半球反照率；基于预先配置的气溶胶光学厚度查找表和第一目标地表反射率确定第2波段的暗像元对应的目标表观反射率；基于指定的背景气溶胶光学厚度和表观反射率对第2波段和第3波段进行大气校正，得到第2波段对应的第二目标地表反射率。
[0008]在可选的实施方式中，对卫星图像中进行暗像元判断，包括：对卫星图像中的预设晴空陆地像元的第6波段进行阈值判断；若小于预设阈值，则为暗像元。
[0009]在可选的实施方式中，基于预先配置的气溶胶光学厚度查找表和第一目标地表反射率确定第2波段的暗像元对应的目标表观反射率，包括：基于多种气溶胶光学厚度分别对应的第一地表反射率和波段比值，模拟第2波段对应的多种气溶胶光学厚度分别对应的第二地表反射率；其中，波段比值为第2波段和第3波段的比值；基于第二地表反射率、第二大气传输模型和预先配置的气溶胶光学厚度查找表中的参数，反演得到第2波段对应的多种气溶胶光学厚度分别对应的目标表观反射率；第二大气传输模型为：。
[0010]在可选的实施方式中，在得到目标表观反射率后，方法还包括：获取第2波段对应的多种气溶胶光学厚度分别对应的实测表观反射率，通过插值得到暗像元在第2波段下的气溶胶光学厚度。
[0011]在可选的实施方式中，方法还包括：基于全圆盘数据中第2波段和第5波段的实测数据计算归一化水指数；将第2波段中的反射率数据进行灰度处理得到云像元分类阈值；基于归一化水指数和云像元分类阈值，提取非云水像元；计算非云水像元对应的气溶胶光学厚度。
[0012]第二方面，本专利技术提供一种FY
‑
4A气溶胶光学厚度的反演装置，装置包括：数据获取模块，用于获取FY
‑
4A预设时段的全圆盘数据，并读取指定波段下的实测数据；大气校正模块，用于基于预先配置的气溶胶光学厚度查找表对指定波段下的实测数据进行大气校正，得到指定波段对应的目标地表反射率；其中，预先配置的气溶胶光学厚度查找表为配置有FY
‑
4A光谱响应函数的6S大气校正查找表；反演模块，用于根据指定波段对应的目标地表反射率进行气溶胶光厚度的多进程反演处理，得到指定波段对应的目标气溶胶光学厚度。
[0013]第三方面，本专利技术提供一种电子设备，包括处理器和存储器，存储器存储有能够被处理器执行的计算机可执行指令，处理器执行计算机可执行指令以实现前述实施方式任一项的FY
‑
4A气溶胶光学厚度的反演方法。
[0014]第四方面，本专利技术提供一种计算机可读存储介质，计算机可读存储介质存储有计算机可执行指令，计算机可执行指令在被处理器调用和执行时，计算机可执行指令促使处理器实现前述实施方式任一项的FY
‑
4A气溶胶光学厚度的反演方法。
[0015]本申请提供的FY
‑
4A气溶胶光学厚度的反演方法、装置、设备和介质，该方法首先
获取FY
‑
4A预设时段的全圆盘数据，并读取指定波段下的实测数据，然后基于预先配置的气溶胶光学厚度查找表（也即配置有FY
‑
4A光谱响应函数的6S大气校正查找表）对指定波段下的实测数据进行大气校正，得到指定波段对应的目标地表反射率，最后根据指定波段对应的目标地表反射率进行气溶胶光厚度的多进程反演处理，得到指定波段对应的目标气溶胶光学厚度。该方法本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种FY
‑
4A气溶胶光学厚度的反演方法，其特征在于，所述方法包括：获取FY
‑
4A预设时段的全圆盘数据，并读取指定波段下的实测数据；基于预先配置的气溶胶光学厚度查找表对所述指定波段下的实测数据进行大气校正，得到所述指定波段对应的目标地表反射率；其中，所述预先配置的气溶胶光学厚度查找表为配置有FY
‑
4A光谱响应函数的6S大气校正查找表；根据所述指定波段对应的目标地表反射率进行气溶胶光厚度的多进程反演处理，得到所述指定波段对应的目标气溶胶光学厚度。2.根据权利要求1所述的FY
‑
4A气溶胶光学厚度的反演方法，其特征在于，所述预先配置的气溶胶光学厚度查找表包括预设个数的太阳高度角、预设个数的太阳天顶角、预设个数的相对方位角与多种气溶胶光学厚度进行组合后的组合结果；每种组合都对应有相应的大气程辐射、总透射率和半球反照率。3.根据权利要求1所述的FY
‑
4A气溶胶光学厚度的反演方法，其特征在于，所述指定波段包括第2波段和第3波段，所述实测数据包括卫星图像；基于预先配置的气溶胶光学厚度查找表对所述指定波段下的实测数据进行大气校正，得到所述指定波段对应的目标地表反射率，包括：对所述卫星图像中进行暗像元判断；基于预先配置的气溶胶光学厚度查找表和第一大气传输模型确定所述第3波段对应的多种气溶胶光学厚度分别对应的第一目标地表反射率；所述第一大气传输模型为：；其中，为地表反射率；为太阳天顶角；为卫星天顶角；为方位角；为表观反射率；为大气程辐射、为总透射率，为半球反照率；基于预先配置的气溶胶光学厚度查找表和所述第一目标地表反射率确定所述第2波段的所述暗像元对应的目标表观反射率；基于指定的背景气溶胶光学厚度和所述表观反射率对所述第2波段和所述第3波段进行大气校正，得到所述第2波段对应的第二目标地表反射率。4.根据权利要求3所述的FY
‑
4A气溶胶光学厚度的反演方法，其特征在于，对所述卫星图像中进行暗像元判断，包括：对卫星图像中的预设晴空陆地像元的第6波段进行阈值判断；若小于预设阈值，则为暗像元。5.根据权利要求3所述的FY
‑
4A气溶胶光学厚度的反演方法，其特征在于，基于预先配置的气溶胶光学厚度查找表和所述第一目标地表反射率确定所述第2波段的所述暗像元对应的目标表...

【专利技术属性】
技术研发人员：吴达周，黄葵，王宇翔，徐欢，宋蕾，刘雨生，汤琼，王振刚，夏丽华，
申请(专利权)人：航天宏图信息技术股份有限公司，
类型：发明
国别省市：