本发明专利技术公开一种基于MODIS数据反演城市大气气溶胶光学厚度的方法。该方法包括：1)根据6S大气辐射传输方程，建立数学模型计算四种基本气溶胶沙尘性、水溶性、海洋性以及煤烟性的体积比浓度，得到自定义的气溶胶类型；2)根据MODIS数据中2.1μm通道的地表反射率与红、蓝通道的地表反射率的关系，得出红、蓝通道的地表反射率的估算值；3)根据6S大气辐射传输方程，结合步骤1中的气溶胶类型，分别建立红通道下气溶胶光学厚度查找表和蓝通道下气溶胶光学厚度查找表；4)对查找表中的数据进行三次样条插值，得出指定地表反射率下的气溶胶光学厚度，完成大气气溶胶光学厚度反演。本发明专利技术城市大气气溶胶光学厚度反演方法易于实现，反演精度高，可以有效的用于城区气溶胶监测。
【技术实现步骤摘要】
【专利摘要】本专利技术公开一种基于MODIS数据反演城市大气气溶胶光学厚度的方法。该方法包括：1)根据6S大气辐射传输方程，建立数学模型计算四种基本气溶胶沙尘性、水溶性、海洋性以及煤烟性的体积比浓度，得到自定义的气溶胶类型；2)根据MODIS数据中2.1μm通道的地表反射率与红、蓝通道的地表反射率的关系，得出红、蓝通道的地表反射率的估算值；3)根据6S大气辐射传输方程，结合步骤1中的气溶胶类型，分别建立红通道下气溶胶光学厚度查找表和蓝通道下气溶胶光学厚度查找表；4)对查找表中的数据进行三次样条插值，得出指定地表反射率下的气溶胶光学厚度，完成大气气溶胶光学厚度反演。本专利技术城市大气气溶胶光学厚度反演方法易于实现，反演精度高，可以有效的用于城区气溶胶监测。【专利说明】一种基于MODIS数据反演城市大气气溶胶光学厚度的方法
本专利技术涉及大气遥感领域的气溶胶特性研究，尤其是涉及一种基于MODIS数据反演城市大气气溶胶光学厚度的方法。
技术介绍
近年来，由于城市交通污染、市政建设、周边地区工业气体排放和自然界生态环境的严重破坏等人为因素影响，城市空气质量受到很大影响，空气污染严重。为减少污染源的影响，对城市大气环境质量的监测变得尤为重要，气溶胶粒子因为其在大气监测中的重要作用受到广泛的关注和研究。大气气溶胶由不同相态的物质组成，虽然含量很少，但其可以通过吸收和散射影响太阳辐射，进而影响地-气系统的辐射收支平衡，对大气中发生的许多物理化学过程都有重要影响，因此其是气候变化模拟和环境遥感中的重要因子。气溶胶光学厚度是气溶胶最重要的参数之一，是表征大气混浊度的重要物理量。探测气溶胶光学厚度可以采用地基探测方法和卫星遥感方法。地基探测可以比较准确地测量气溶胶信息，但获取的只是空间点上的数据，不能反映大区域气溶胶时空分布。遥感技术具有覆盖面积广、信息获取方便、快捷等特点，相对于地面观测可以更高效地获取大气气溶胶信息。当前利用Terra/Aqua卫星上搭载的MODIS传感器观测的高分辨率卫星遥感资料反演气溶胶光学厚度可以克服地基探测的不足，为人们实时了解大区域范围内的气溶胶变化提供了可能。不少学者针对卫星遥感反演气溶胶光学厚度这一问题做过相关研究。例如中国专利公告CN 102636143A (公告日期:2012年8月15日)公开了一种气溶胶光学厚度遥感反演方法。该方法依据纯像元指数提取遥感图像中的纯像元，生成纯像元掩膜，进而确定纯像元地表反射率，最后根据查找表将遥感观测的辐射亮度反演为气溶胶光学厚度。该方法虽然不受混合像元影响，能有效的用于城区的气溶胶监测，但是它选用6S大气辐射传输模型中的既定的气溶胶类型，把整个地球大气气溶胶只划分为沙尘性以及城市性显然非常粗糙。
技术实现思路
本专利技术的目的在于提供一种基于MODIS数据的反演城市大气气溶胶光学厚度的方法，通过建立数学模型来确定最适当的大气气溶胶类型，由6S大气福射传输模型建立气溶胶光学厚度的查找表，并通过对查找表中数据的三次样条插值得出反演结果。为实现上述目的，本专利技术采用以下步骤:I)根据6S大气辐射传输方程，建立数学模型计算四种基本气溶胶沙尘性、水溶性、海洋性以及煤烟性的体积比浓度，得到自定义的气溶胶类型。具体为，令沙尘性气溶胶体积比浓度为V1，水溶性气溶胶体积比浓度为V2，海洋性气溶胶体积比浓度为V3，煤烟性气溶胶体积比浓度为％，其中v4 = 1-(V^vfv3);根据待反演城市的特征,分别确定νρν2、ν3、V4的临界值，设V1的最小值为vlmin，V1的最大值为Vlmax, V2的最小值为V2min, V2的最大值为V2max, V3的最小值为v3min，V3的最大值为V3max ；M0DIS数据中卫星观测的表观反射率在I波段(660nm)下为p ^，将V1从Vlniin到Vlniax进行递增迭代,将V2从V2niin到V2niax进行递增迭代,将V3从V3min到V3max进行递增迭代,并根据V4的表达式得到V4的值,将每一组Vp V2> V3> V4的值带入6S大气辐射传输方程分别计算在660nm下的卫星反射率P1，计算P ^与P i的差值的平方，此数据用ε表示,将ε与对应的VpVyVyV4的值都记录下来；当VpVyV3都达到I时迭代完毕。从记录的所有数据中取出ε最小的200万组数据，取这200万组数据对应的V1的平均数作为自定义的气溶胶类型中沙尘性气溶胶体积浓度；取这200万组数据对应的V2的平均数作为自定义的气溶胶类型中水溶性气溶胶体积浓度；取这200万组数据对应的V3的平均数作为自定义的气溶胶类型中海洋性气溶胶体积浓度；取这200万组数据对应的V4的平均数作为自定义的气溶胶类型中煤烟性气溶胶体积浓度。2)根据MODIS数据中2.1 μ m通道的地表反射率与红、蓝通道的地表反射率的关系，得出红、蓝通道的地表反射率的估算值。3)根据6S大气辐射传输方程，结合步骤I中的气溶胶类型，分别建立红通道下气溶胶光学厚度查找表和蓝通道下气溶胶光学厚度查找表；具体的，查找表中，气溶胶光学厚度从0.00到2.50之间迭代，地表反射率从0.00到0.10之间迭代。4)对查找表中的数据进行三次样条插值，得出指定地表反射率下的气溶胶光学厚度，完成大气气溶胶光学厚度反演。所述的步骤I)中Vp v2、V3递增迭代的步长均为0.001。所述的步骤I)中通过6S大气辐射传输方程计算在660nm下的卫星反射率P i时，方程中气溶胶光学厚度取为太阳光度计观测的气溶胶光学厚度。所述的步骤3)所建的两个查找表中，气溶胶光学厚度的迭代步长为0.05，地表反射率的迭代步长为0.01。本专利技术城市大气气溶胶光学厚度反演方法易于实现，反演精度高，可以有效的用于城区气溶胶监测。【专利附图】【附图说明】图1是本专利技术的工作流程示意图。图2是采用杭州地区2013年10月的MODIS LlB数据对本专利技术的方法进行操作的操作结果与杭州站太阳光度计的测量结果的对比图。【具体实施方式】如图1所示，本专利技术包含以下步骤:I)根据6S大气辐射传输方程，建立数学模型计算四种基本气溶胶沙尘性、水溶性、海洋性以及煤烟性的体积比浓度，得到自定义的气溶胶类型。具体为，令沙尘性气溶胶体积比浓度为V1，水溶性气溶胶体积比浓度为V2，海洋性气溶胶体积比浓度为V3，煤烟性气溶胶体积比浓度为％，其中v4 = 1-(V^vfv3);根据待反演城市的特征,分别确定νρν2、ν3、V4的临界值，设V1的最小值为vlmin，V1的最大值为Vlmax, V2的最小值为V2min, V2的最大值为V2max, V3的最小值为v3min，V3的最大值为V3max ；M0DIS数据中卫星观测的表观反射率在I波段(660nm)下为p ^，将V1从Vlniin到Vlniax进行递增迭代,将V2从V2niin到V2niax进行递增迭代,将V3从V3min到V3max进行递增迭代,并根据V4的表达式得到V4的值,将每一组Vp V2> V3> V4的值带入6S大气辐射传输方程分别计算在660nm下的卫星反射率P1，计算P ^与P i的差值的平方，此数据用ε表示,将ε与对应的VpVyVyV4的值都记录下来；当VpVyV3都达到I时迭代本文档来自技高网...
【技术保护点】
一种基于MODIS数据反演城市大气气溶胶光学厚度的方法，其特征在于包含以下步骤：1)根据6S大气辐射传输方程，建立数学模型计算四种基本气溶胶沙尘性、水溶性、海洋性以及煤烟性的体积比浓度，得到自定义的气溶胶类型。具体为，令沙尘性气溶胶体积比浓度为v1，水溶性气溶胶体积比浓度为v2，海洋性气溶胶体积比浓度为v3，煤烟性气溶胶体积比浓度为v4，其中v4＝1‑(v1+v2+v3)；根据待反演城市的特征，分别确定v1、v2、v3、v4的临界值，设v1的最小值为v1min，v1的最大值为v1max，v2的最小值为v2min，v2的最大值为v2max，v3的最小值为v3min，v3的最大值为v3max；MODIS数据中卫星观测的表观反射率在1波段(660nm)下为ρ0，将v1从v1min到v1max进行递增迭代，将v2从v2min到v2max进行递增迭代，将v3从v3min到v3max进行递增迭代，并根据v4的表达式得到v4的值，将每一组v1、v2、v3、v4的值带入6S大气辐射传输方程分别计算在660nm下的卫星反射率ρ1，计算ρ0与ρ1的差值的平方，此数据用ε表示，将ε与对应的v1、v2、v3、v4的值都记录下来；当v1、v2、v3都达到1时迭代完毕。从记录的所有数据中取出ε最小的200万组数据，取这200万组数据对应的v1的平均数作为自定义的气溶胶类型中沙尘性气溶胶体积浓度；取这200万组数据对应的v2的平均数作为自定义的气溶胶类型中水溶性气溶胶体积浓度；取这200万组数据对应的v3的平均数作为自定义的气溶胶类型中海洋性气溶胶体积浓度；取这200万组数据对应的v4的平均数作为自定义的气溶胶类型中煤烟性气溶胶体积浓度。2)根据MODIS数据中2.1μm通道的地表反射率与红、蓝通道的地表反射率的关系，得出红、蓝通道的地表反射率的估算值。3)根据6S大气辐射传输方程，结合步骤1中的气溶胶类型，分别建立红通道下气溶胶光学厚度查找表和蓝通道下气溶胶光学厚度查找表；具体的，查找表中，气溶胶光学厚度从0.00到2.50之间迭代，地表反射率从0.00到0.10之间迭代。4)对查找表中的数据进行三次样条插值，得出指定地表反射率下的气溶胶光学厚度，完成大气气溶胶光学厚度反演。...
【技术特征摘要】
【专利技术属性】
技术研发人员：范娇，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。