CN106483050A - 气溶胶光学厚度的反演方法及*** - Google Patents

气溶胶光学厚度的反演方法及*** Download PDF

Info

Publication number
CN106483050A
CN106483050A CN201510557877.5A CN201510557877A CN106483050A CN 106483050 A CN106483050 A CN 106483050A CN 201510557877 A CN201510557877 A CN 201510557877A CN 106483050 A CN106483050 A CN 106483050A
Authority
CN
China
Prior art keywords
aerosol
inverting
look
reflectivity
type
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510557877.5A
Other languages
English (en)
Other versions
CN106483050B (zh
Inventor
李莘莘
陈良富
范萌
司丹
司一丹
余超
苏林
陶金花
李小英
王子峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Remote Sensing and Digital Earth of CAS
Original Assignee
Institute of Remote Sensing and Digital Earth of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Remote Sensing and Digital Earth of CAS filed Critical Institute of Remote Sensing and Digital Earth of CAS
Priority to CN201510557877.5A priority Critical patent/CN106483050B/zh
Publication of CN106483050A publication Critical patent/CN106483050A/zh
Application granted granted Critical
Publication of CN106483050B publication Critical patent/CN106483050B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种气溶胶光学厚度的反演方法及***,其中的方法包括,建立气溶胶反演查找表以及获取地表反射率;基于高分一号卫星WFV数据,代入获取的地表反射率,并通过查找和插值建立的气溶胶反演查找表的大气参数,获得气溶胶光学厚度;其中,基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型;将前述的计算结果作为输入参数,利用6SV1辐射传输方程构建气溶胶反演查找表;以及,基于国外卫星MODIS的地表二向反射率产品,计算地表反射率。通过本发明能够快速准确的实现大范围的气溶胶监测,为颗粒物浓度的遥感估算提供必要前提。

Description

气溶胶光学厚度的反演方法及***
技术领域
本发明涉及大气遥感技术领域,更为具体地,涉及一种气溶胶光学厚度的反演方法及***。
背景技术
气溶胶是指悬浮在气体中的固体或液体微粒和气体载体共同组成的多相体系。气溶胶来源于自然界和人类活动,自然源包括火山灰、海洋碎末、风沙等。更大程度上,气溶胶是由工业生产、汽车尾气、燃煤取暖、秸秆焚烧等不完全燃烧产生并进入大气的颗粒物质。气溶胶含量虽然只占地球大气组分的很少一部分,但是由于其在大气中的广泛分布、巨大的空间变化、短暂的生命期及复杂的化学组成,对天气过程、全球气候和人类健康均有重要影响。因气溶胶问题一直是全球的大气科学研究工作者研究的热点问题。目前,传统的地面气溶胶监测仅能够提供“点”尺度的信息,因此,如何利用卫星遥感技术对大区域气溶胶进行快速监测具有重要的现实意义和科学价值。
高分一号卫星于2013年4月成功发射入轨,该星采用太阳同步轨道,平均轨道高度约为644km,属于光学成像遥感卫星,设计寿命为5-8年。高分一号卫星搭载有4个中分辨率相机,星下点能达到16m的空间分辨率,共有蓝、绿、红、近红外四个波段,每个相机幅宽能够达到200km,拼接4台相机(WFV1、WFV2、WFV3、WFV4)的观测数据可实现800km幅宽的观测,具有4天覆盖中国全境的能力。基于GF-1WFV数据,可以反演气溶胶光学厚度。
卫星遥感监测气溶胶最常用的参数是气溶胶光学厚度,它代表整层的气溶胶污染状况。气溶胶光学厚度的卫星遥感反演需要解决两大问题:气溶胶模式和地表噪声的去除。气溶胶模式即气溶胶类型,代表不同气溶胶粒子(如硫酸盐、铵盐、含碳类气溶胶、沙尘等)的组成情况和微物理特性。
目前,气溶胶模式主要基于地面观测数据进行模拟,此种方法的缺点是难以反映气溶胶模式的时空分布特性。而在地表噪声去除方面,目前应用最广泛的是美国国家航空航天局(National Aeronautics and Space Administration,简称NASA)的暗像元算法,暗像元算法利用在近红外波段受大气影响较小的特点来去除地表噪声。但是由于高分一号卫星WFV数据缺少近红外波段通道,无法构造地表反射率,从而无法实现气溶胶光学厚度的卫星遥感反演。
发明内容
鉴于上述问题,本发明的目的是提供一种气溶胶光学厚度的反演方法及***,以解决上述背景技术中所提出的问题。
本发明提供一种气溶胶光学厚度的反演方法,包括:
S1:建立气溶胶反演查找表以及获取地表反射率;
S2:基于高分一号卫星WFV数据,代入步骤S1中获取的地表反射率,并通过查找和插值步骤S12中建立的气溶胶反演查找表的大气参数,获得气溶胶光学厚度;
其中,
建立气溶胶反演查找表的过程包括:
S11:基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型;
S12:将步骤S11的计算结果作为输入参数,利用6SV1辐射传输方程构建气溶胶反演查找表;
获取地表反射率的过程包括:
基于国外卫星MODIS的地表二向反射率产品,计算地表反射率。
本发明还提供一种气溶胶光学厚度的反演***,包括:气溶胶反演查找表创建单元、地表反射率计算单元和气溶胶光学厚度反演单元;其中,
气溶胶反演查找表创建单元,用于建立气溶胶反演查找表;其中,气溶胶反演查找表创建单元包括:
气溶胶组分体积比计算模块,用于基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型;
气溶胶反演查找表构建模块,用于将气溶胶组分计算模块的计算结果作为输入参数,利用6SV1辐射传输方程构建气溶胶反演查找表;
地表反射率计算单元,用于基于国外卫星MODIS的地表二向反射率产品,计算地表反射率;
气溶胶光学厚度反演单元,用于基于高分一号卫星WFV数据,代入地表反射率计算单元计算出的地表反射率,并通过查找和插值气溶胶反演查找表构建模块构建的气溶胶反演查找表的大气参数,获得气溶胶光学厚度。
利用上述根据本发明提供的气溶胶光学厚度的反演方法及***,基于已有的地基颗粒物组分实验观测数据,利用相关的分类原则,通过实验观测结果和大气化学模式模拟结果的匹配验证,校正模式输出的沙尘型、水溶型、煤烟型和海盐型的气溶胶体积比,并作为辐射传输方程的输入参数构建反演查找表,再基于高分一号WFV数据精确的反演气溶胶光学厚度,为颗粒物卫星遥感定量估算提供必要的前提。
为了实现上述以及相关目的,本发明的一个或多个方面包括后面将详细说明并在权利要求中特别指出的特征。下面的说明以及附图详细说明了本发明的某些示例性方面。然而,这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外,本发明旨在包括所有这些方面以及它们的等同物。
附图说明
通过参考以下结合附图的说明及权利要求书的内容,并且随着对本发明的更全面理解,本发明的其它目的及结果将更加明白及易于理解。在附图中:
图1为根据本发明实施例的气溶胶光学厚度的反演方法的流程示意图;
图2为根据本发明实施例的气溶胶光学厚度的反演***的逻辑结构框图。
在所有附图中相同的标号指示相似或相应的特征或功能。
具体实施方式
在下面的描述中,出于说明的目的,为了提供对一个或多个实施例的全面理解,阐述了许多具体细节。然而,很明显,也可以在没有这些具体细节的情况下实现这些实施例。在其它例子中,为了便于描述一个或多个实施例,公知的结构和设备以方框图的形式示出。
本发明采用GEOS-Chem大气化学模式数据构建立格网尺度的气溶胶模式。GEOS-Chem大气化学模式数据通过NASA全球模式同化办公室的戈达得地球观测***(Goddard Earth Observing System,GEOS)同化气象数据驱动。GEOS-Chem大气化学模式数据提供的沙尘、硫酸盐、硝酸盐、铵盐、黑碳、有机碳、海洋型等颗粒物的体积比、浓度、光学特性、廓线等信息,可以有效支撑气溶胶模式的构建。在地表噪声去除方面,本发明通过利用其它卫星(如MODIS)构造辅助的地表反射率,来实现高分一号卫星数据的气溶胶光学厚度反演。
以下将结合附图对本发明的具体实施例进行详细描述。
图1示出了根据本发明实施例的气溶胶光学厚度的反演方法的流程。
如图1所示,本发明提供的气溶胶光学厚度的反演方法,包括:
S1:建立气溶胶反演查找表以及获取地表反射率。
建立气溶胶反演查找表以及获取地表反射率的过程不分先后顺序。
一、气溶胶反演查找表的建立
建立气溶胶反演查找表的过程包括:
S11:基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型。
基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比的过程,包括:
S111:使用GEOS-Chem大气化学模式数据来构建气溶胶模型,利用气溶胶模型估算硫酸盐气溶胶、硝酸盐气溶胶、铵盐气溶胶、沙尘气溶胶、黑碳气溶胶、有机碳气溶胶和海盐气溶胶的体积比。
需要说明的是,GEOS-Chem大气化学模式数据为GEOS-Chem全球大气化学模式的模拟结果。
S112:将硫酸盐气溶胶、硝酸盐气溶胶、铵盐气溶胶、沙尘气溶胶、黑碳气溶胶、有机碳气溶胶和海盐气溶胶按照组份划分为沙尘型、水溶型、烟尘型和海盐型,并通过如下公式计算每种类型气溶胶组份的体积比:
[Dust-like]=[DST1]+[DST2]+[DST3] (1)
[Water-Soluble]=[SO4]+[NH4]+[NO3] (2)
[Soot]=[BC]+[OC] (3)
[Oceanic]=[NaCl] (4)
其中,[Dust-like]表示沙尘型气溶胶的体积比;
[DST1]、[DST2]、[DST3]表示GEOS-Chem大气化学模式数据输出的三种不同粒径的沙尘气溶胶的体积比;
[Water-Soluble]表示水溶型气溶胶的体积比;
[SO4]、[NH4]、[NO3]分别表示GEOS-Chem大气化学模式数据输出的硫酸盐气溶胶、铵盐气溶胶和硝酸盐气溶胶的体积比;
[Soot]表示烟尘型气溶胶的体积比;
[BC]、[OC]分别表示GEOS-Chem大气化学模式数据输出的黑碳气溶胶和有机碳气溶胶的体积比;
[Oceanic]表示海洋型气溶胶的体积比;
[NaCl]表示GEOS-Chem大气化学模式数据输出的氯化钠气溶胶的体积比。
需要说明的是,沙尘型气溶胶包括沙尘气溶胶,沙尘型气溶胶的体积比为三种不同粒径的沙尘气溶胶的体积比之和;水溶型气溶胶包括硫酸盐气溶胶、铵盐气溶胶和硝酸盐气溶胶,水溶型气溶胶的体积比为硫酸盐气溶胶、铵盐气溶胶和硝酸盐气溶胶的体积比之和;烟尘型气溶胶包括黑碳气溶胶和有机碳气溶胶,烟尘型气溶胶的体积比为黑碳气溶胶和有机碳气溶胶的体积比之和;海洋型气溶胶包括氯化钠气溶胶,海洋型气溶胶体积比为氯化钠气溶胶的体积比。
S12:将步骤S11的计算结果作为输入参数,利用6SV1(Second Simulationof a Satellite Signal in the Solar Spectrum,Vector version 1)辐射传输方程构建气溶胶反演查找表。
6SV1辐射传输方程构建气溶胶反演查找表的过程,包括:
S121:设定6SV1辐射传输方程的输入参数,输入参数包括:9个太阳天顶角、12个观测天顶角、16个相对方位角、气溶胶模式,以及将550nm处的气溶胶光学厚度设定为10个等级、将气溶胶反演查找表计算的波段设置为高分一号卫星WFV传感器的波段响应函数;其中,气溶胶模式根据步骤S11中的四种类型气溶胶的组份确定。
S122:调用6SV1辐射传输方程,计算大气参数,并将计算出的大气参数存储在二进制文件中生成气溶胶反演查找表,大气参数包括:大气整层透过率、大气程辐射、大气向下的半球反射率。
二、地表反射率的获取
本发明实施例基于国外卫星MODIS的地表二向反射率产品,计算地表反射率。计算地表反射率的公式如下:
其中,r为地表反射率,kvol为体散射的核,kgeo为几何光学核,θS为光线入射角,θv为观测角,为观测与太阳间的相对方位角,三个角度可以通过高分数据获得,fiso、fgeo、fvol存储在国外卫星MODIS的地表二向反射率产品中,分别表示各向同性散射、几何光学散射、体散射在这三部分中所占的权重。
通过高分数据获得光线入射角,观测角,观测与太阳间的相对方位角,计算体散射的核和几何光学核;基于MODIS的地表二向反射率产品,计算各向同性散射、几何光学散射、体散射的权重。
S2:基于高分一号卫星WFV数据,代入步骤S1中获取的地表反射率,并通过查找和插值步骤S12中建立的气溶胶反演查找表的大气参数,获得气溶胶光学厚度。
气溶胶光学厚度通过如下公式反演获得:
其中,τ为气溶胶光学厚度,θ为天顶角,下标s表示太阳,v表示观测,φ为相对方位角,r为地表反射率,ρ0为大气程辐射的等效反射率,S为大气下界向下的半球反射率,T为大气整层透过率,其中,S、ρ0、T是表征大气状况的三个参数,存储在S12步骤建立的气溶胶反演查找表中,ρTOA为卫星观测到的表观反射率。
将S、ρ0、T三个大气参数和卫星观测到的表观反射率ρTOA、以及步骤S1步骤中获得的地表反射率,代入上述公式就可反演得到气溶胶光学厚度。
需要说明的是,对WFV数据匹配最为接近的卫星观测几何条件,输入卫星观测的近红外波段表观反射率,通过插值气溶胶光学厚度得到可见光波段的表观反射率,与卫星观测的可见光波段的表观反射率进行比较获得最小差值时,该插值得到的气溶胶光学厚度即为该点的气溶胶光学厚度。
与上述气溶胶光学厚度的反演方法相对应,本发明提供一种气溶胶光学厚度的反演***。图2示出了根据本发明实施例的气溶胶光学厚度的反演***的逻辑结构。
如图2所示,本发明实施例提供的气溶胶光学厚度的反演***包括气溶胶反演查找表创建单元210、地表反射率计算单元220和气溶胶光学厚度反演单元230。
气溶胶反演查找表创建单元210,用于建立气溶胶反演查找表;其中,气溶胶反演查找表创建单元210包括:
气溶胶组分体积比计算模块211,基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型。
气溶胶反演查找表构建模块212,将气溶胶组分体积比计算模块211的计算结果作为输入参数,利用6SV1辐射传输方程构建气溶胶反演查找表。
地表反射率计算单元220,基于国外卫星MODIS的地表二向反射率产品,计算地表反射率。
气溶胶光学厚度反演单元230,基于高分一号卫星WFV数据,代入地表反射率计算单元220计算出的地表反射率,并通过查找和插值气溶胶反演查找表构建模块212构建的气溶胶反演查找表的大气参数,获得气溶胶光学厚度。
其中,在气溶胶组分体积比计算模块211计算每种气溶胶组分体积比的过程中,
S111:根据GEOS-Chem大气化学模式数据估算硫酸盐气溶胶、硝酸盐气溶胶、铵盐气溶胶、沙尘气溶胶、黑碳气溶胶、有机碳气溶胶和海盐气溶胶的体积比;
S112:将硫酸盐气溶胶、硝酸盐气溶胶、铵盐气溶胶、沙尘气溶胶、黑碳气溶胶、有机碳气溶胶和海盐气溶胶按照组份划分为沙尘型、水溶型、烟尘型和海盐型,并通过如下公式计算每种类型气溶胶组份的体积比:
[Dust-like]=[DST1]+[DST2]+[DST3] (1)
[Water-Soluble]=[SO4]+[NH4]+[NO3] (2)
[Soot]=[BC]+[OC] (3)
[Oceanic]=[NaCl] (4)
其中,[Dust-like]表示沙尘型气溶胶的体积比;
[DST1]、[DST2]、[DST3]分别表示GEOS-Chem大气化学模式数据输出的三种不同粒径的沙尘气溶胶的体积比;
[Water-Soluble]表示水溶型气溶胶的体积比;
[SO4]、[NH4]、[NO3]分别表示GEOS-Chem大气化学模式数据输出的硫酸盐气溶胶、铵盐气溶胶和硝酸盐气溶胶的体积比;
[Soot]表示烟尘型气溶胶的体积比;
[BC]、[OC]分别表示GEOS-Chem大气化学模式数据输出的黑碳气溶胶和有机碳气溶胶的体积比;
[Oceanic]表示海洋型气溶胶的体积比;
[NaCl]表示GEOS-Chem大气化学模式数据输出的氯化钠气溶胶的体积比。
此外,在气溶胶反演查找表构建模块212构建气溶胶反演查找表的过程中,
S121:设定6SV1辐射传输方程的输入参数,所述输入参数包括:9个太阳天顶角、12个观测天顶角、16个相对方位角、气溶胶模式,以及将550nm处的气溶胶光学厚度设定为10个等级、将气溶胶反演查找表计算的波段设置为高分一号卫星WFV传感器的波段响应函数;其中,所述气溶胶模式根据步骤S11中的四种类型气溶胶的组份确定;
S122:调用6SV1辐射传输方程,计算大气参数并存储在二进制文件中生成气溶胶反演查找表,所述大气参数包括:大气整层透过率、大气程辐射、大气向下的半球反射率。
另外,地表反射率计算单元根据如下公式计算地表反射率:
其中,r为地表反射率,kvol为体散射的核,kgeo为几何光学核,θS为光线入射角,θv为观测角,为观测与太阳间的相对方位角,fiso、fgeo、fvol存储在国外卫星MODIS的地表二向反射率产品中,分别表示各向同性散射、几何光学散射、体散射的权重。
再者,气溶胶光学厚度反演单元230根据如下公式反演气溶胶光学厚度:
其中,τ为气溶胶光学厚度,θ为天顶角,下标s表示太阳,v表示观测,φ为相对方位角,r为地表反射率,ρ0为大气程辐射的等效反射率,S为大气向下的半球反射率,T为大气整层透过率,ρTOA为卫星观测到的表观反射率。
上述内容详细地说明了本发明实施例提供的气溶胶光学厚度的反演方法及***,基于已有的地基颗粒物组分实验观测数据,利用相关的分类原则,通过实验观测结果和大气化学模式模拟结果的匹配验证,校正模式输出的沙尘型、水溶型、煤烟型和海盐型等气溶胶体积比,并作为辐射传输方程的输入参数构建反演查找表,再基于高分一号WFV数据精确的反演气溶胶光学厚度,能够快速准确的实现大范围的气溶胶监测,为颗粒物浓度的遥感估算提供必要前提。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。

Claims (7)

1.一种气溶胶光学厚度的反演方法,包括:
S1:建立气溶胶反演查找表以及获取地表反射率;
S2:基于高分一号卫星WFV数据,代入步骤S1中获取的地表反射率,并通过查找和插值步骤S12中建立的所述气溶胶反演查找表的大气参数,获得气溶胶光学厚度;
其中,所述建立气溶胶反演查找表的过程包括:
S11:基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,所述四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型;
S12:将步骤S11的计算结果作为输入参数,利用6SV1辐射传输方程构建气溶胶反演查找表;
所述获取地表反射率的过程包括:
基于国外卫星MODIS的地表二向反射率产品,计算所述地表反射率。
2.如权利要求1所述的气溶胶光学厚度的反演方法,其中,步骤S11进一步包括:
S111:根据GEOS-Chem大气化学模式数据估算硫酸盐气溶胶、硝酸盐气溶胶、铵盐气溶胶、沙尘气溶胶、黑碳气溶胶、有机碳气溶胶和海盐气溶胶的体积比;
S112:将硫酸盐气溶胶、硝酸盐气溶胶、铵盐气溶胶、沙尘气溶胶、黑碳气溶胶、有机碳气溶胶和海盐气溶胶按照组份划分为沙尘型、水溶型、烟尘型和海盐型,并通过如下公式计算每种类型气溶胶组份的体积比:
[Dust-like]=[DST1]+[DST2]+[DST3] (1)
[Water-Soluble]=[SO4]+[NH4]+[NO3] (2)
[Soot]=[BC]+[OC] (3)
[Oceanic]=[NaCl] (4)
其中,[Dust-like]表示沙尘型气溶胶的体积比;
[DST1]、[DST2]、[DST3]分别表示GEOS-Chem大气化学模式数据输出的三种不同粒径的沙尘气溶胶的体积比;
[Water-Soluble]表示水溶型气溶胶的体积比;
[SO4]、[NH4]、[NO3]分别表示GEOS-Chem大气化学模式数据输出的硫酸盐气溶胶、铵盐气溶胶和硝酸盐气溶胶的体积比;
[Soot]表示烟尘型气溶胶的体积比;
[BC]、[OC]分别表示GEOS-Chem大气化学模式数据输出的黑碳气溶胶和有机碳气溶胶的体积比;
[Oceanic]表示海洋型气溶胶的体积比;
[NaCl]表示GEOS-Chem大气化学模式数据输出的氯化钠气溶胶的体积比。
3.如权利要求2所述的气溶胶光学厚度的反演方法,其中,步骤S12进一步包括:
S121:设定6SV1辐射传输方程的输入参数,所述输入参数包括:9个太阳天顶角、12个观测天顶角、16个相对方位角、气溶胶模式,以及将550nm处的气溶胶光学厚度设定为10个等级、将气溶胶反演查找表计算的波段设置为高分一号卫星WFV传感器的波段响应函数;其中,所述气溶胶模式根据步骤S11中的四种类型气溶胶的组份确定;
S122:调用所述6SV1辐射传输方程,计算大气参数并存储在二进制文件中生成气溶胶反演查找表,所述大气参数包括:大气整层透过率、大气程辐射、大气向下的半球反射率。
4.如权利要求1所述的气溶胶光学厚度的反演方法,其中,通过如下公式计算地表反射率:
其中,r为地表反射率,kvol为体散射的核,kgeo为几何光学核,θS为光线入射角,θv为观测角,为观测与太阳间的相对方位角,fiso、fgeo、fvol存储在国外卫星MODIS的地表二向反射率产品中,分别表示各向同性散射、几何光学散射、体散射的权重。
5.如权利要求1所述的气溶胶光学厚度的反演方法,在步骤S2中,
通过如下公式反演气溶胶光学厚度:
ρ T O A ( θ s , θ v , φ , τ ) = ρ o ( θ s , θ v , φ , τ ) + T ( θ s , θ v , τ ) r [ 1 - r S ( τ ) ]
其中,τ为气溶胶光学厚度,θ为天顶角,下标s表示太阳,v表示观测,φ为相对方位角,r为地表反射率,ρ0为大气程辐射的等效反射率,S为大气向下的半球反射率,T为大气整层透过率,ρTOA为卫星观测到的表观反射率。
6.如权利要求1-5中任一项所述的气溶胶光学厚度的反演方法,其中,GEOS-Chem大气化学模式数据为GEOS-Chem全球大气化学模式的模拟结果。
7.一种气溶胶光学厚度的反演***,包括:气溶胶反演查找表创建单元、地表反射率计算单元和气溶胶光学厚度反演单元;其中,
所述气溶胶反演查找表创建单元,用于建立气溶胶反演查找表;其中,所述气溶胶反演查找表创建单元包括:
气溶胶组分体积比计算模块,用于基于GEOS-Chem大气化学模式数据计算四种类型气溶胶组分的体积比,所述四种类型气溶胶分别为沙尘型、水溶型、烟尘型和海盐型;
气溶胶反演查找表构建模块,用于将所述气溶胶组分计算模块的计算结果作为输入参数,利用6SV1辐射传输方程构建气溶胶反演查找表;
所述地表反射率计算单元,用于基于国外卫星MODIS的地表二向反射率产品,计算地表反射率;
所述气溶胶光学厚度反演单元,用于基于高分一号卫星WFV数据,代入所述地表反射率计算单元计算出的地表反射率,并通过查找和插值所述气溶胶反演查找表构建模块构建的气溶胶反演查找表的大气参数,获得气溶胶光学厚度。
CN201510557877.5A 2015-09-02 2015-09-02 气溶胶光学厚度的反演方法及*** Active CN106483050B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510557877.5A CN106483050B (zh) 2015-09-02 2015-09-02 气溶胶光学厚度的反演方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510557877.5A CN106483050B (zh) 2015-09-02 2015-09-02 气溶胶光学厚度的反演方法及***

Publications (2)

Publication Number Publication Date
CN106483050A true CN106483050A (zh) 2017-03-08
CN106483050B CN106483050B (zh) 2019-08-13

Family

ID=58238658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510557877.5A Active CN106483050B (zh) 2015-09-02 2015-09-02 气溶胶光学厚度的反演方法及***

Country Status (1)

Country Link
CN (1) CN106483050B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107389617A (zh) * 2017-06-27 2017-11-24 环境保护部卫星环境应用中心 基于高分四号卫星的气溶胶光学厚度的反演方法及设备
CN108426815A (zh) * 2018-04-20 2018-08-21 中国科学院遥感与数字地球研究所 一种近地面细颗粒物组分浓度估算方法
CN108872093A (zh) * 2018-05-02 2018-11-23 复旦大学 基于被动遥感的o4吸收校正及气溶胶消光廓线反演方法
CN109272053A (zh) * 2018-10-12 2019-01-25 国网湖南省电力有限公司 极轨卫星监测气溶胶光学厚度的数据融合方法及***
CN109781593A (zh) * 2019-01-17 2019-05-21 南京泛在地理信息产业研究院有限公司 一种基于parasol多角度偏振数据的气溶胶二次反演方法
CN110186822A (zh) * 2019-05-13 2019-08-30 中国科学院遥感与数字地球研究所 一种气溶胶光学厚度遥感反演方法
CN110186823A (zh) * 2019-06-26 2019-08-30 中国科学院遥感与数字地球研究所 一种气溶胶光学厚度反演方法
CN110411927A (zh) * 2019-08-02 2019-11-05 中国科学院遥感与数字地球研究所 一种大气细粒子aod和地表偏振反射率协同反演方法
CN110455689A (zh) * 2019-08-01 2019-11-15 中国科学院遥感与数字地球研究所 一种模拟冰晶粒子的光散射特性的方法
CN111123382A (zh) * 2019-12-25 2020-05-08 中国科学院遥感与数字地球研究所 一种气溶胶和地表参数联合反演方法
CN111191380A (zh) * 2020-01-08 2020-05-22 北京大学 一种基于地基光谱仪测量数据的大气气溶胶光学厚度估算方法和装置
CN111709180A (zh) * 2020-06-24 2020-09-25 中国科学院空天信息创新研究院 耦合吸收性气溶胶模型的大气矢量辐射传输模拟方法
CN112964666A (zh) * 2021-02-03 2021-06-15 中国电子科技集团公司第二十八研究所 一种基于地表双向反射的大气二氧化碳含量计算方法
CN113466917A (zh) * 2020-03-31 2021-10-01 中国矿业大学(北京) 一种煤烟型气溶胶辐射强迫计算新方法
CN115524264A (zh) * 2022-11-04 2022-12-27 山东科技大学 一种基于激光雷达测量的气溶胶分类方法
CN116466368A (zh) * 2023-06-16 2023-07-21 成都远望科技有限责任公司 基于激光雷达和卫星资料的沙尘消光系数廓线估算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103927455A (zh) * 2014-04-24 2014-07-16 中国科学院遥感与数字地球研究所 基于高分一号卫星的陆地气溶胶光学性质反演方法
CN107273820A (zh) * 2017-05-26 2017-10-20 中国科学院遥感与数字地球研究所 一种土地覆盖遥感分类方法及***

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103927455A (zh) * 2014-04-24 2014-07-16 中国科学院遥感与数字地球研究所 基于高分一号卫星的陆地气溶胶光学性质反演方法
CN107273820A (zh) * 2017-05-26 2017-10-20 中国科学院遥感与数字地球研究所 一种土地覆盖遥感分类方法及***

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SHENSHEN LI等: ""Retrieval of the Haze Optical Thickness in North China Plain Using MODIS Data"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
王中挺等: ""利用MODIS 数据监测北京地区气溶胶"", 《遥感技术与应用》 *
茆佳佳: ""华东地区气溶胶特性的MODIS资料反演及其时空分布特征的研究"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107389617A (zh) * 2017-06-27 2017-11-24 环境保护部卫星环境应用中心 基于高分四号卫星的气溶胶光学厚度的反演方法及设备
CN108426815A (zh) * 2018-04-20 2018-08-21 中国科学院遥感与数字地球研究所 一种近地面细颗粒物组分浓度估算方法
CN108872093A (zh) * 2018-05-02 2018-11-23 复旦大学 基于被动遥感的o4吸收校正及气溶胶消光廓线反演方法
CN108872093B (zh) * 2018-05-02 2020-09-29 复旦大学 基于被动遥感的o4吸收校正及气溶胶消光廓线反演方法
CN109272053A (zh) * 2018-10-12 2019-01-25 国网湖南省电力有限公司 极轨卫星监测气溶胶光学厚度的数据融合方法及***
CN109272053B (zh) * 2018-10-12 2021-11-02 国网湖南省电力有限公司 极轨卫星监测气溶胶光学厚度的数据融合方法及***
CN109781593A (zh) * 2019-01-17 2019-05-21 南京泛在地理信息产业研究院有限公司 一种基于parasol多角度偏振数据的气溶胶二次反演方法
CN110186822B (zh) * 2019-05-13 2020-09-11 中国科学院遥感与数字地球研究所 一种气溶胶光学厚度遥感反演方法
CN110186822A (zh) * 2019-05-13 2019-08-30 中国科学院遥感与数字地球研究所 一种气溶胶光学厚度遥感反演方法
CN110186823A (zh) * 2019-06-26 2019-08-30 中国科学院遥感与数字地球研究所 一种气溶胶光学厚度反演方法
CN110455689A (zh) * 2019-08-01 2019-11-15 中国科学院遥感与数字地球研究所 一种模拟冰晶粒子的光散射特性的方法
CN110411927A (zh) * 2019-08-02 2019-11-05 中国科学院遥感与数字地球研究所 一种大气细粒子aod和地表偏振反射率协同反演方法
CN111123382A (zh) * 2019-12-25 2020-05-08 中国科学院遥感与数字地球研究所 一种气溶胶和地表参数联合反演方法
CN111123382B (zh) * 2019-12-25 2020-11-13 中国科学院遥感与数字地球研究所 一种气溶胶和地表参数联合反演方法
CN111191380B (zh) * 2020-01-08 2020-11-20 北京大学 一种基于地基光谱仪测量数据的大气气溶胶光学厚度估算方法和装置
CN111191380A (zh) * 2020-01-08 2020-05-22 北京大学 一种基于地基光谱仪测量数据的大气气溶胶光学厚度估算方法和装置
CN113466917A (zh) * 2020-03-31 2021-10-01 中国矿业大学(北京) 一种煤烟型气溶胶辐射强迫计算新方法
CN113466917B (zh) * 2020-03-31 2023-12-15 中国矿业大学(北京) 一种煤烟型气溶胶辐射强迫计算方法
CN111709180B (zh) * 2020-06-24 2023-10-13 中国科学院空天信息创新研究院 耦合吸收性气溶胶模型的大气矢量辐射传输模拟方法
CN111709180A (zh) * 2020-06-24 2020-09-25 中国科学院空天信息创新研究院 耦合吸收性气溶胶模型的大气矢量辐射传输模拟方法
CN112964666A (zh) * 2021-02-03 2021-06-15 中国电子科技集团公司第二十八研究所 一种基于地表双向反射的大气二氧化碳含量计算方法
CN115524264A (zh) * 2022-11-04 2022-12-27 山东科技大学 一种基于激光雷达测量的气溶胶分类方法
CN116466368B (zh) * 2023-06-16 2023-08-22 成都远望科技有限责任公司 基于激光雷达和卫星资料的沙尘消光系数廓线估算方法
CN116466368A (zh) * 2023-06-16 2023-07-21 成都远望科技有限责任公司 基于激光雷达和卫星资料的沙尘消光系数廓线估算方法

Also Published As

Publication number Publication date
CN106483050B (zh) 2019-08-13

Similar Documents

Publication Publication Date Title
CN106483050A (zh) 气溶胶光学厚度的反演方法及***
Banfield et al. The atmosphere of Mars as observed by InSight
Xu et al. Concentric gravity waves over northern China observed by an airglow imager network and satellites
CN110726653B (zh) 基于天地一体化信息的pm2.5浓度监测方法
Zhang et al. The impact of circulation patterns on regional transport pathways and air quality over Beijing and its surroundings
Kumar et al. Wintertime characteristics of aerosols at middle Indo-Gangetic Plain: Impacts of regional meteorology and long range transport
KR101307908B1 (ko) 지표면 도달 일사량 계산을 위한 단층 태양복사모델의 구성 방법 및 장치
CN107065036B (zh) 一种联合遥感和气象数据获取全天候蒸散发的方法
CN102435586A (zh) 地表反照率产品的生成方法及***
CN101936777A (zh) 一种基于热红外遥感反演近地层气温的方法
CN109406361B (zh) 一种基于遥感技术的干旱区灰霾污染预警方法
CN107389617A (zh) 基于高分四号卫星的气溶胶光学厚度的反演方法及设备
CN113324656A (zh) 无人机搭载红外遥感的地表热异常探测方法及***
Chaboureau et al. Fennec dust forecast intercomparison over the Sahara in June 2011
Zhang et al. Influence of ice particle model on satellite ice cloud retrieval: lessons learned from MODIS and POLDER cloud product comparison
Vasarevičius et al. Solar irradiance model for solar electric panels and solar thermal collectors in Lithuania
CN103927454B (zh) 一种基于环境卫星的灰霾污染监测方法
Dewanto et al. Application of web 3D GIS to display urban model and solar energy analysis using the unmanned aerial vehicle (UAV) data (Case study: National Cheng Kung university buildings)
CN103743488A (zh) 遥感卫星地球临边背景特性的红外成像仿真方法
CN111914396B (zh) 基于高分辨率dem数据的次网格地形三维地表太阳辐射强迫效应快速参数化方法
Shen et al. Identification of polluted clouds and composition analysis based on GF-5 DPC data
Wei et al. Black and odorous water monitoring by using gf series remote sensing data
Nesari et al. Aerosol optical depth retrieval using Landsat 8 satellite imagery-Case study: City of Tamanrasset, Algeria
Zoran et al. Impact of climate and anthropogenic changes on a periurban forest surface albedo derived from MODIS satellite data
Li et al. The impact of observation nudging on simulated meteorology and ozone concentrations during DISCOVER-AQ 2013 Texas campaign.

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant