一种基于卫星偏振技术的PM2.5浓度遥感估算方法与流程
本发明涉及遥感技术领域，特别涉及一种基于卫星偏振技术的pm2.5浓度遥感估算方法。
背景技术：
大气细粒子气溶胶通常是指悬浮于大气中的微小粒子，粒径范围可以从0.001μm到几十微米，作为地-气系统的重要组成部分，它主要通过直接辐射强迫和间接辐射强迫影响着气候。大气细粒子气溶胶同时具有显著的环境效应，直径在10μm以下的气溶胶颗粒物(pm10)可到达人类呼吸系统的支气管区；直径小于2.5μm的气溶胶颗粒物(pm2.5)微粒可到达肺泡区，最终导致心血管和哮喘疾病的增加，直接对人类健康造成显著的影响，威胁着人类的生存与社会可持续发展。与粒径较粗的大气颗粒物相比，pm2.5粒径小，富含大量的有毒、有害物质，且在大气中的停留时间长、输送距离远，因而对人体健康和大气环境质量的影响更大。
近年来我国的大气污染问题比较严重。2013年1月中国东部自北向南的大部分区域发生了严重的霾污染。在北京地区，共发生5次空气质量指数(aqi)超过200(重度污染)的灰霾天气。在这5次重度污染过程中，北京的首要污染物均为细颗粒物pm2.5，是造成严重灰霾天气的主要原因，因此，研究pm2.5浓度具有重要的意义。国家环境保护部于2013年1月开始发布了全国主要地面监测站点观测的pm2.5浓度数据。尽管地面pm2.5观测技术能够获得实时、高精度的pm2.5浓度，但是其覆盖范围有限，主要分布在重点城市区域，广大的农村区域站点较少，导致其空间性代表性不足，对于研究全国或重点区域的大气环境质量存在一定的局限性。卫星遥感可以连续、大面积获取动态数据，可以有效弥补地面pm2.5监测站点的不足。因为气溶胶光学厚度与pm2.5浓度有着较强的相关性，所以基于气溶胶光学厚度的遥感估算方法被广泛应用于pm2.5浓度的估算工作中。现有的pm2.5浓度遥感估算方法普遍精度不高，使得pm2.5浓度反演的效果总体不理想。由于大气细粒子气溶胶散射辐射具有强偏振特性，而大多数陆地表面散射辐射具有弱偏振特性且其时空变化较小，因此表观偏振反射率主要来自大气细粒子气溶胶，利用偏振信息可以有效的将大气细粒子气溶胶和地表的贡献区分开来，pm2.5是大气细粒子气溶胶颗粒物，所以，利用卫星偏振技术监测pm2.5浓度具有独特的优势。
技术实现要素：
针对目前卫星遥感估算pm2.5浓度算法普遍精度不高的现状，为提高pm2.5浓度遥感估算方法的空间分辨率、覆盖率和精度，本发明提出了一种基于卫星偏振技术的pm2.5浓度遥感估算方法。
本发明由下述技术方案实现：
一种基于卫星偏振遥感技术的pm2.5浓度估算方法，所述方法包括如下步骤：
基于卫星偏振遥感数据反演得到单个像元的大气细粒子气溶胶光学厚度；
对大气细粒子气溶胶光学厚度进行垂直校正和湿度校正，获取近地面大气细粒子气溶胶的干消光系数；
基于近地面大气细粒子气溶胶的干消光系数和地面观测的pm2.5浓度，构建卫星偏振遥感估算pm2.5浓度模型，得到估算的pm2.5浓度；
选取相应时间、地点的地面pm2.5浓度数据，验证上述卫星偏振遥感估算pm2.5浓度模型的精度。
进一步的，所述基于卫星偏振遥感数据反演得到单个像元的大气细粒子气溶胶光学厚度的步骤包括：
基于大气细粒子气溶胶类型和矢量辐射传输模型建立大气细粒子气溶胶参数查找表；
基于卫星偏振遥感数据、地表偏振模型和大气细粒子气溶胶参数查找表反演得到单个像元的大气细粒子气溶胶光学厚度。
进一步的，所述卫星偏振遥感数据为基于parasol和gf-5/dpc载荷观测的数据，所述地表偏振模型表达式为：
γ＝αcos(-cosθs×cosθv-sinθs×sinθv×cosδφ)
其中，
进一步的，所述基于卫星偏振遥感数据和地表偏振模型反演得到单个像元的大气细粒子气溶胶光学厚度的步骤包括：
将parasol和gf-5/dpc载荷同一个像元多个角度(具体的，可以选择3-12个)真实观测的表观偏振反射率减去所述由地表偏振模型计算得到的地表偏振反射率值，得到卫星观测的大气细粒子气溶胶偏振反射率值；
将所述卫星观测的大气细粒子气溶胶偏振反射率值与通过大气细粒子气溶胶参数查找表模拟的大气细粒子气溶胶偏振反射率进行比较，得到总体误差ε，所述总体误差ε最小所对应的大气细粒子气溶胶光学厚度值即为所述反演得到的单个像元的大气细粒子气溶胶光学厚度。
其中，所述ε的表达式为：
其中，ε为总体误差，λi代表波长，i代表第i个偏振波段，一共有两个，分别是670nm和865nm的偏振波段，j代表第j个观测角度，n代表观测角度总数，
进一步的，所述对大气细粒子气溶胶光学厚度进行垂直校正和湿度校正，获取近地面大气细粒子气溶胶的干消光系数的步骤具体为：
所述垂直校正、湿度校正和干消光系数的计算公式分别如下：
β0＝aod/h
f(rh)＝1/(1-rh)
β1＝β0/f(rh)
其中,aod为大气细粒子气溶胶光学厚度，h为大气边界层高度，f(rh)为吸湿增长因子，rh为大气相对湿度，β1为大气细粒子气溶胶的干消光系数。
进一步的，所述卫星偏振遥感估算pm2.5浓度模型的表达式为：
其中，(μi,vi)为第i个样本点的经纬度坐标，a0(μi,vi)是第i个样本点上的初始回归参数，回归参数是观测点地理位置的函数，ak(μi,vi)是第i个样本点上的第k个回归参数，p是样本点的回归参数总数，εi为第i个样本点上的误差项，pm2.5i为第i个样本点上的pm2.5浓度，βik为第i个样本点上的近地面大气细粒子气溶胶的干消光系数。
进一步的，所述选取相应时间、地点的地面pm2.5浓度数据，验证上述卫星偏振遥感估算pm2.5浓度模型估算的pm2.5浓度精度的步骤包括：
选取卫星偏振遥感估算pm2.5浓度模型估算的pm2.5浓度的结果；
根据经纬度最邻近的原则，选取当天卫星过境前后1小时内的地面pm2.5浓度数据，并对所述地面pm2.5浓度数据进行平均；
验证所述卫星偏振遥感估算pm2.5浓度模型估算的pm2.5浓度精度。
综上所述，本发明提供了一种基于卫星偏振技术的pm2.5浓度遥感估算方法，可以精确的反演某一区域pm2.5浓度，通过最终得到的基于观测点地理位置信息的回归模型，提高了pm2.5浓度遥感估算的空间分辨率、空间覆盖率和精度，同时为多角度偏振载荷在大气环境监测方面的成功应用提供技术关键支撑。
附图说明
图1为本发明的卫星偏振遥感估算pm2.5浓度的方法流程图；
图2为本发明的基于卫星偏振遥感数据反演得到单个像元的大气细粒子气溶胶光学厚度的方法流程图；
图3为具体实施例中的估算pm2.5浓度卫星遥感技术路线图；
图4为具体实施例中的pm2.5浓度卫星遥感专题图；
图5为具体实施例中的基于地面pm2.5监测站点的卫星偏振遥感技术估算pm2.5浓度模型的精度验证结果图；
图6为具体实施例中的基于地面每个城市pm2.5监测站点平均值的卫星偏振遥感技术估算pm2.5浓度模型的精度验证结果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式并参照附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
本发明提供了一种基于卫星偏振技术的pm2.5浓度遥感估算方法。如图1所示，本发明的估算方法包括如下步骤：
步骤s100,基于卫星偏振遥感数据反演得到单个像元的大气细粒子气溶胶光学厚度(aod)。
进一步的，该基于卫星偏振遥感数据反演得到单个像元的大气细粒子气溶胶光学厚度的步骤包括，如图2所示：
步骤s110，基于大气细粒子气溶胶类型和矢量辐射传输模型建立大气细粒子气溶胶参数查找表。
具体的，选择大气细粒子气溶胶类型，利用mie散射理论计算细粒子气溶胶的偏振特性和散射特性，6种细粒子气溶胶类型如表1所示。
表16种细粒子气溶胶类型
其中，rm为气溶胶粒子中值半径，s为标准差，c为粒子浓度；mr为复折射指数的实部，mi为复折射指数的虚部，ssa为单次散射反照率，670nm和865nm为波长。
具体的，通过设定不同的太阳和卫星观测几何参数：21个太阳天顶角、20个观测天顶角、37个相对方位角，以及不同的大气细粒子气溶胶参数:6种细粒子气溶胶类型、12个大气细粒子aod，考虑到观测数据所在的波段，使用rt3矢量辐射传输模型进行计算，得到多组具有偏振信息和强度信息的大气参数组合而成的大气细粒子气溶胶参数查找表，如表2所示。
表2查找表构建
步骤s120，基于卫星偏振遥感数据、地表偏振模型和大气细粒子气溶胶参数查找表反演得到单个像元的大气细粒子气溶胶光学厚度。
具体的，利用地表偏振模型及其α、β经验系数，在地表偏振模型中输入parasol和gf-5/dpc观测的太阳天顶角、观测天顶角、相对方位角以及地表类型ndvi(归一化植被指数)，计算得到parasol和gf-5/dpc载荷观测的地表偏振反射率。地表偏振模型公式如下：
γ＝αcos(-cosθs×cosθv-sinθs×sinθv×cosδφ)公式(2)
其中，
将parasol和gf-5/dpc载荷同一个像元多个(具体的，可以选择3-12个)角度真实观测的表观偏振反射率值减去上述计算得到的地表偏振反射率，得到卫星观测的大气细粒子气溶胶偏振反射率，与通过大气细粒子气溶胶参数查找表模拟的大气细粒子气溶胶偏振反射率进行比较，得到总体误差ε，所述总体误差ε最小所对应的大气细粒子气溶胶光学厚度值即为最终结果，从而得到高空间分辨率的大气细粒子气溶胶aod(6kmⅹ7km或3.3kmⅹ3.3km)。
其中，总体误差ε可由下式计算；
其中，ε为总体误差，λi代表波长，i代表第i个偏振波段，一共有两个，分别是670nm和865nm的偏振波段，j代表第j个观测角度，n代表观测角度总数，
步骤s200，对大气细粒子气溶胶光学厚度进行垂直校正和湿度校正，获取近地面大气细粒子气溶胶的干消光系数。
具体的，利用大气边界层高度和相对湿度数据对大气细粒子aod开展垂直校正和湿度校正，获取近地面大气细粒子气溶胶的干消光系数。垂直校正、湿度校正和干消光系数的计算公式分别如下：
β0＝aod/h公式(4)
f(rh)＝1/(1-rh)公式(5)
β1＝β0/f(rh)公式(6)
其中,aod为大气细粒子气溶胶光学厚度，h为大气边界层高度，f(rh)为吸湿增长因子，rh为大气相对湿度，β1为大气细粒子气溶胶的干消光系数。
步骤s300，基于近地面大气细粒子气溶胶的干消光系数和地面观测的pm2.5浓度，构建卫星偏振遥感估算pm2.5浓度模型，得到估算的pm2.5浓度。
具体的，将pm2.5浓度的空间特性纳入公式(7)模型中，假定线性回归模型中的回归参数是观测点地理位置的函数，构建了基于βk的卫星偏振遥感估算pm2.5浓度模型，模型表达式为：
其中，(μi,vi)为第i个样本点的经纬度坐标，a0(μi,vi)是第i个样本点上的初始回归参数，回归参数是观测点地理位置的函数，ak(μi,vi)是第i个样本点上的第k个回归参数，p是样本点的回归参数总数，εi为第i个样本点上的误差项，pm2.5i为第i个样本点上的pm2.5浓度，βik为第i个样本点上的近地面大气细粒子气溶胶的干消光系数。
步骤s400，选取相应时间、地点的地面pm2.5浓度数据，验证上述卫星偏振遥感估算pm2.5浓度模型估算的pm2.5浓度精度。
具体的，选取卫星偏振遥感估算pm2.5浓度的结果，根据经纬度最邻近的原则，再选取当天卫星过境前后1小时内的地面pm2.5浓度数据，并对所述地面pm2.5浓度数据进行平均，用于验证卫星估算的pm2.5浓度精度，评价卫星偏振遥感估算pm2.5浓度模型的精度。
下面以一个具体的实施例，按照上述步骤对本发明进行进一步的说明。
以京津冀地区为例，具体如图3所示，是基于parasol载荷估算京津冀区域的pm2.5浓度方法流程图。如图4所示，基于卫星偏振遥感估算pm2.5浓度的结果，利用arcmap软件将估算的pm2.5浓度制成专题图。2013年5月19号，京津冀区域北部空气质量总体较好，其中承德市、张家口市、北京市大部分区域pm2.5浓度小于35μg/m3,保定市、廊坊市、石家庄市、天津市、沧州市、衡水市、邢台市的部分区域pm2.5浓度较高。图4有效的显示了京津冀地区pm2.5浓度的空间变化情况，为京津冀区域空气质量监测提供了关键技术支撑。
本发明进一步使用2013年5月19号的地面pm2.5浓度数据进行对比验证。利用本发明方法估算的京津冀区域pm2.5浓度和地面实际测量值的对比验证结果见图5和图6。图5是基于地面每个pm2.5监测站点的卫星偏振遥感估算pm2.5浓度模型的精度验证结果，图5给出了线性拟合的斜率、截距、相关系数(r)和均方根误差(rmse)，结果显示r2为0.75，rmse为15.42，拟合斜率分别为0.65，这些结果显示了本发明方法的精度高。图6是基于地面每个城市pm2.5监测站点平均值的卫星偏振遥感估算pm2.5浓度模型的精度验证结果，由图6可以用得出每个城市的pm2.5浓度估算值与地面观测值基本一致，这些结果证明本发明可以有效的用于所处区域的pm2.5浓度监测。
综上所述，本发明提供了一种基于卫星偏振技术的pm2.5浓度遥感估算方法，可以有效地用于精确估算某一区域pm2.5浓度，其空间分辨率、覆盖率和精度更高，有效弥补了地面pm2.5监测站点的不足，为多角度偏振载荷在大气环境监测方面的成功应用提供技术关键支撑。
应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。
- 还没有人留言评论。精彩留言会获得点赞！