地统计插值与卫星遥感联合反演地面PM2.5的方法及系统

摘要

本发明提供一种地统计插值与卫星遥感联合反演地面PM2.5的方法及系统，其中的方法包括：基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；基于构建的卫星AOD反演模型，获取待反演地区的地面PM2.5的第二浓度值；将第一浓度值与第二浓度值进行融合，估算待反演地区的地面PM2.5的最终浓度值；其中，根据第一浓度值和第二浓度值的均方根误差建立权重关系，根据权重关系融合第一浓度值与第二浓度值，以估算待反演地区的地面PM2.5的最终浓度值。通过本发明既可以提高地面监测站点较少情况下的PM2.5的估算精度，又可以弥补AOD数据缺失处PM2.5无法估算的缺点。

权利要求

1.一种地统计插值与卫星遥感联合反演地面PM2.5的方法，包括：

基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；其中，获取所述第一浓度值的过程包括：当各PM2.5站点所监测到的PM2.5观测数据之间存在时空相关性时，根据待反演地区的地面PM2.5数据构建时空变异函数模型并拟合，根据对所述时空变异函数模型的拟合结果，采用时空克里金插值模型估算待反演地区的地面PM2.5的第一浓度值；

基于构建的卫星AOD反演模型，获取待反演地区的地面PM2.5的第二浓度值；其中，所述卫星AOD反演模型包括混合效应模型和地理加权回归模型，获取所述第二浓度值的过程包括：根据待反演地区的地面PM2.5数据、AOD数据和气象数据构建随时间变化的混合效应模型，并根据所述混合效应模型对待反演地区的地面PM2.5进行初步估算，以及根据待反演地区的地理坐标构建地理加权回归模型，根据所述地理加权回归模型获取所述混合效应模型的残差，根据所述混合效应模型的初步估算结果和所述残差估算待反演地区的地面PM2.5的第二浓度值；

根据所述第一浓度值和第二浓度值的均方根误差建立权重关系，当所述第二浓度值与所述第一浓度值的均方根误差差值在预设范围内时，根据所述权重关系融合所述第一浓度值与第二浓度值，以估算待反演地区的地面PM2.5的最终浓度值。

2.如权利要求1所述的地统计插值与卫星遥感联合反演地面PM2.5的方法，在获取待反演地区的地面PM2.5的第一浓度值之前，还包括：根据空间变异函数与时间相关函数获取待反演地区内各PM2.5站点所监测到的PM2.5观测数据的时空相关性。

3.如权利要求1所述的地统计插值与卫星遥感联合反演地面PM2.5的方法，其中，在构建混合效应模型的过程中，将待反演地区的AOD数据、地面PM2.5数据和气象数据重采样到创建的与AOD数据分辨率相当的网格后进行数据匹配；其中，所述网格由若干网格单元构成；所述数据匹配过程包括：构建所述AOD数据像元点的泰森多边形，与所述网格进行空间叠置分析，将AOD数据分配到所述网格中的每个网格单元中；

将所述网格中的某个网格单元内所有PM2.5站点同一天采集到的地面PM2.5数据进行平均后赋值给对应的网格单元；

将高于待反演地区所对应的网格单元的分辨率的气象数据进行平均后赋值给对应的网格单元，将不高于待反演地区所对应的网格单元的分辨率的气象数据，采用距离反比加权的插值方法将所述气象数据插值并重采样到与所述网格单元的分辨率一致。

4.如权利要求3所述的地统计插值与卫星遥感联合反演地面PM2.5的方法，其中，在根据地理加权回归模型获取所述混合效应模型的残差的过程中，根据待反演地区的每个PM2.5站点的位置坐标，获取所述地理加权回归模型的回归参数；

根据所述回归参数、与所述PM2.5站点对应的AOD数据、所述PM2.5站点的位置坐标的截距与斜率获取所述混合效应模型的残差。

5.如权利要求4所述的地统计插值与卫星遥感联合反演地面PM2.5的方法，其中，对所述混合效应模型和所述地理加权回归模型进行交叉验证，根据交叉验证结果调整所述混合效应模型和地理加权回归模型，或确定最终的混合效应模型和地理加权回归模型。

6.如权利要求5所述的地统计插值与卫星遥感联合反演地面PM2.5的方法，其中，将确定的最终的混合效应模型和地理加权回归模型的估算结果相加，估算所述待反演地区的地面PM2.5的第二浓度值。

7.如权利要求1所述的地统计插值与卫星遥感联合反演地面PM2.5的方法，其中，将所述第一浓度值与第二浓度值进行融合的过程中，获取所述第一浓度值、第二浓度值的均方根误差权重，以及所述第二浓度值与所述第一浓度值的均方根误差差值；

当所述差值在预设范围内时，根据所述第一浓度值与所述第一浓度值的均方根误差权重、所述第二浓度值与所述第二浓度值的均方根误差权重对所述第一浓度值与第二浓度值进行融合；

当所述差值不在预设范围内时，不对所述第一浓度值与第二浓度值进行融合。

8.一种地统计插值与卫星遥感联合反演地面PM2.5的系统，包括：

第一浓度值获取单元，用于基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；其中，获取的所述第一浓度值的过程包括：当各PM2.5站点所监测到的PM2.5观测数据之间存在时空相关性时，根据待反演地区的地面PM2.5数据构建时空变异函数模型并拟合，根据对所述时空变异函数模型的拟合结果，采用时空克里金插值模型估算待反演地区的地面PM2.5的第一浓度值；

第二浓度值获取单元，用于基于构建的卫星AOD反演模型获取待反演地区的地面PM2.5的第二浓度值；其中，所述卫星AOD反演模型包括混合效应模型和地理加权回归模型，获取所述第二浓度值的过程包括：根据待反演地区的地面PM2.5数据、AOD数据和气象数据构建随时间变化的混合效应模型，并根据混合效应模型对待反演地区的地面PM2.5进行初步估算，以及根据待反演地区的地理坐标构建地理加权回归模型，根据所述地理加权回归模型获取所述混合效应模型的残差，根据所述混合效应模型的初步估算结果和所述残差估算待反演地区的地面PM2.5的第二浓度值；

PM2.5浓度值估算单元，用于根据所述第一浓度值和第二浓度值的均方根误差建立权重关系，当所述第二浓度值与所述第一浓度值的均方根误差差值在预设范围内时，根据所述权重关系融合所述第一浓度值与第二浓度值，估算待反演地区的地面PM2.5的最终浓度值。

9.如权利要求8所述的地统计插值与卫星遥感联合反演地面PM2.5的系统，所述PM2.5浓度值估算单元进一步包括：误差权重获取单元，用于分别获取所述第一浓度值、第二浓度值的均方根误差权重；

差值获取单元，用于获取所述第二浓度值与所述第一浓度值的均方根误差差值；

融合单元，用于根据所述误差权重获取单元所获取的第一浓度值与第一浓度值的均方根误差权重、所述差值获取单元所获取的第二浓度值与第二浓度值的均方根误差权重对所述第一浓度值与第二浓度值进行融合。

说明书

地统计插值与卫星遥感联合反演地面PM2.5的方法及系统

技术领域

[0001] 本发明涉及气溶胶监测技术领域，更为具体地，涉及一种地统计插值与卫星遥感联合反演地面PM2.5的方法及系统。

背景技术

[0002] 随着经济的快速发展，工业活动与机动车尾气等人为排放的有害气体的急剧增加，导致空气质量持续恶化。PM2.5是指空气中空气动力学粒径小于2.5微米的颗粒物。与大粒径颗粒物相比，PM2.5粒径小，富含大量的有毒有害物质且在大气中的停留时间长、输送距离远，因而对人体和大气环境质量的影响很大。大量的流行病学研究证明，PM2.5与哮喘、呼吸道感染、肺癌、心血管疾病等存在一定的关联性。因此，对区域PM2.5的监测成为研究大气环境的关键点之一。

[0003] 目前，国内外对区域PM2.5的估算方法按照估算尺度可以分为两大类，一类是区域尺度的PM2.5浓度估算方法，此类方法主要是利用卫星遥感估算PM2.5浓度分布，另一类是对城市尺度污染物浓度进行估算的方法，此类方法主要基于土地覆盖数据、PM2.5观测站点、路网分布、气象数据等辅助数据，通过地统计插值、土地利用回归模型等方法获取PM2.5浓度的时空分布。

[0004] 卫星遥感很早就被应用于大气污染监测，但真正应用于定量估算颗粒物浓度主要始于MODIS(moderate-resolution Imaging Spectroradiometer，中分辨率成像光谱仪)、MISR(Multi-angle Imaging Spectrometer，多角度成像光谱仪)等可提供较高精度气溶胶数据的传感器升空以后。其中，卫星气溶胶光学厚度(Aerosol Optical Depth，AOD)在国际上被广泛用来估算区域PM2.5浓度。然而，卫星AOD代表气溶胶在整层大气上的垂直消光贡献，PM2.5则为近地面颗粒物干燥条件下的单位体积内的质量浓度，二者之间的关系受气溶胶垂直分布、气象条件、下垫面等影响较大。利用AOD估算PM2.5浓度目前还没有广泛适用的方法；同时，卫星AOD的反演精度、以及区域的云覆盖、高地表反射率、重霾污染等造成的值缺失也会对PM2.5的估算结果带来很多不确定性。

[0005] 地统计插值的方法可以有效的利用已知空间数据估算邻近区域的未知值，广泛应用在城市尺度的空气污染研究中。其中，变异函数分析工具可以有效的对PM2.5的时空变异特征进行分析。克里金插值是最常用的地统计插值方法，它可以有效的将少量的监测站点扩展到整个研究区域。但是，地统计插值方法受监测站点的分布特征以及密度影响较大，当监测站点间隔较远、采样密度降低时，可能无法较好的反映PM2.5的时空分布特征，造成插值结果精度较低。

[0006] 如上所述，目前在对区域PM2.5进行估算的方法中，卫星遥感与地统计插值这两大类方法基本上完全独立，且两类方法各自存在缺点。例如，在利用卫星遥感对区域PM2.5进行估算时，则无法估算出卫星AOD缺失处的PM2.5浓度；在利用地统计插值对区域PM2.5进行估算时，若地面监测站点较少则会降低PM2.5的估算精度。因此有必要将两类方法融合，以研究出一种精度更高的空间插值模型。

发明内容

[0007] 鉴于上述问题，本发明的目的是提供一种地统计插值与卫星遥感联合反演地面PM2.5的方法及系统，以解决现有的估算方法对区域PM2.5的估算精度不高的问题。

[0008] 根据本发明的一个方面，提供一种地统计插值与卫星遥感联合反演地面PM2.5的方法，包括：

[0009] 基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；其中，获取第一浓度值的过程包括：根据待反演地区的地面PM2.5数据构建时空变异函数模型并拟合，根据对时空变异函数模型的拟合结果，采用时空克里金插值模型估算待反演地区的地面PM2.5的第一浓度值；

[0010] 基于构建的卫星AOD反演模型，获取待反演地区的地面PM2.5的第二浓度值；其中，卫星AOD反演模型包括混合效应模型和地理加权回归模型，获取第二浓度值的过程包括：根据待反演地区的地面PM2.5数据、AOD数据和气象数据构建随时间变化的混合效应模型，并根据混合效应模型对待反演地区的地面PM2.5进行初步估算，以及根据待反演地区的地理坐标构建地理加权回归模型，根据地理加权回归模型获取混合效应模型的残差，根据混合效应模型的初步估算结果和残差估算待反演地区的地面PM2.5的第二浓度值；

[0011] 将第一浓度值与第二浓度值进行融合，估算待反演地区的地面PM2.5的最终浓度值；其中，根据第一浓度值和第二浓度值的均方根误差建立权重关系，根据权重关系融合第一浓度值与第二浓度值，以估算待反演地区的地面PM2.5的最终浓度值。

[0012] 另一方面，本发明提供一种地统计插值与卫星遥感联合反演地面PM2.5的系统，包括：

[0013] 第一浓度值获取单元，用于基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；其中，获取第一浓度值的过程包括：根据待反演地区的地面PM2.5数据构建时空变异函数模型并拟合，根据对时空变异函数模型的拟合结果，采用时空克里金插值模型估算待反演地区的地面PM2.5的第一浓度值；

[0014] 第二浓度值获取单元，用于基于构建的卫星AOD反演模型获取待反演地区的地面PM2.5的第二浓度值；其中，卫星AOD反演模型包括混合效应模型和地理加权回归模型，获取第二浓度值的过程包括：根据待反演地区的地面PM2.5数据、AOD数据和气象数据构建随时间变化的混合效应模型，并根据混合效应模型对待反演地区的地面PM2.5进行初步估算，以及根据待反演地区的地理坐标构建地理加权回归模型，根据地理加权回归模型获取混合效应模型的残差，根据混合效应模型的初步估算结果和残差估算待反演地区的地面PM2.5的第二浓度值；

[0015] PM2.5浓度值估算单元，用于将第一浓度值与第二浓度值进行融合，估算待反演地区的地面PM2.5的最终浓度值；其中，根据第一浓度值和第二浓度值的均方根误差建立权重关系，根据权重关系融合第一浓度值与第二浓度值，估算待反演地区的地面PM2.5的最终浓度值。

[0016] 利用上述根据本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的方法及系统，将地统计插值方法中的时空克里金插值模型与卫星遥感方法相融合，从而既可以提高地面监测站点较少情况下的PM2.5的估算精度，又可以弥补AOD数据缺失处PM2.5无法估算的缺点。

[0017] 为了实现上述以及相关目的，本发明的一个或多个方面包括后面将详细说明并在权利要求中特别指出的特征。下面的说明以及附图详细说明了本发明的某些示例性方面。然而，这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外，本发明旨在包括所有这些方面以及它们的等同物。

附图说明

[0018] 通过参考以下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：

[0019] 图1为根据本发明实施例的地统计插值与卫星遥感联合反演地面PM2.5的方法流程示意图；

[0020] 图2为根据本发明实施例的地统计插值与卫星遥感联合反演地面PM2.5的系统的逻辑结构框图。

[0021] 在所有附图中相同的标号指示相似或相应的特征或功能。

具体实施方式

[0022] 以下将结合附图对本发明的具体实施例进行详细描述。

[0023] 针对前述现有的估算方法对区域PM2.5的估算精度不高的问题，本发明通过将地统计插值模型中的时空克里金插值模型与卫星遥感相结合，根据误差权重将两种模型的估算结果进行融合，从而既能保证PM2.5的估算精度，又可以弥补AOD缺失处PM2.5无法估算的缺点。

[0024] 在对本发明进行说明前，先对本发明中涉及的概念和术语进行说明。

[0025] 在本发明中采用美国宇航局Aqua卫星上的中分辨率成像光谱仪(MODerate-resolution Imaging Spectroradiometer，简称MODIS)第六版(C6)卫星遥感AOD数据与PM2.5地面监测站点数据相结合，以估算并反演我国全国尺度的PM2.5浓度分布。

[0026] 为了说明本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的方法，图1示出了根据本发明实施例的地统计插值与卫星遥感联合反演地面PM2.5的方法流程。

[0027] 如图1所示，本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的方法包括：

[0028] S110：基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；其中，获取第一浓度值的过程包括：根据待反演地区的地面PM2.5数据构建时空变异函数模型并拟合，根据对时空变异函数模型的拟合结果，采用时空克里金插值模型估算待反演地区的地面PM2.5的第一浓度值。

[0029] 进一步地，在获取第一浓度值之前，还可以根据搜集到的地面PM2.5站点监测到的PM2.5数据，利用空间变异函数与时间相关函数计算各PM2.5站点所监测到的PM2.5数据的时间相关特性与空间相关特性。如果PM2.5站点所监测到的PM2.5数据之间具有时空相关性，则可采用地统计插值模型中的时空克里金插值模型估算待反演地区的PM2.5的浓度值(即第一浓度值)。

[0030] S120：基于构建的卫星AOD反演模型，获取待反演地区的地面PM2.5的第二浓度值；其中，卫星AOD反演模型包括混合效应模型和地理加权回归模型，获取第二浓度值的过程包括：根据待反演地区的地面PM2.5数据、AOD数据和气象数据构建随时间变化的混合效应模型，并根据混合效应模型对待反演地区的地面PM2.5进行初步估算，以及根据待反演地区的地理坐标构建地理加权回归模型，根据地理加权回归模型获取混合效应模型的残差，根据混合效应模型的初步估算结果和地理加权回归模型所获取的残差估算待反演地区的地面PM2.5的第二浓度值。

[0031] 具体地，在构建混合效应模型的过程中，搜集待反演地区的地面气象数据、AOD数据、地面PM2.5数据，将搜集到的上述数据重采样到创建的与AOD数据分辨率相当的网格后进行数据匹配。其中，上述网格由若干网格单元构成。

[0032] 上述数据匹配过程包括：构建AOD数据像元点的泰森多边形，与网格进行空间叠置分析，将AOD数据分配到网格中的每个网格单元中；将网格中的某个网格单元内所有PM2.5站点同一天采集到的地面PM2.5数据进行平均后赋值给对应的网格单元；将高于待反演地区所对应的网格单元的分辨率的气象数据进行平均后赋值给对应的网格单元，将不高于待反演地区所对应的网格单元的分辨率的气象数据，采用距离反比加权的插值方法将气象数据插值并重采样到与网格单元的分辨率一致。根据上述AOD数据、PM2.5数据和气象数据构建随时间变化的混合效应模型。

[0033] 为了进一步提高PM2.5的估算精度，还需根据待反演地区的地理坐标构建地理加权回归模型，以获取上述混合效应模型的残差。

[0034] 具体地，构建地理加权回归模型的过程包括：根据待反演地区的每个PM2.5站点的位置坐标，获取地理加权回归模型的回归参数；根据回归参数、与PM2.5站点对应的AOD数据、PM2.5站点的位置坐标的截距与斜率获取混合效应模型的残差。

[0035] 为了进一步提高估算精度，对上述混合效应模型和地理加权回归模型进行交叉验证，根据交叉验证结果调整混合效应模型和地理加权回归模型，或确定最终的混合效应模型和地理加权回归模型。根据确定的最终的混合效应模型和地理加权回归模型的估算结果相加，获取待反演地区的地面PM2.5的第二浓度值。

[0036] 具体地，当交叉验证结果显示混合效应模型和地理加权回归模型存在过度拟合时，通过调整AOD数据、气象数据对混合效应模型和地理加权模型进行调整，然后对调整后的混合效应模型和地理加权模型再次进行交叉验证，直到交叉验证结果显示调整后的混合效应模型和地理加权模型的过度拟合在预设的范围内时，根据调整后的混合效应模型和地理加权模型估算待反演地区的地面PM2.5的浓度值(即将确定的最终的混合效应模型和地理加权回归模型的估算结果相加，估算待反演地区的地面PM2.5的第二浓度值)。

[0037] S130：将第一浓度值与第二浓度值进行融合，估算待反演地区的地面PM2.5的最终浓度值；其中，根据第一浓度值和第二浓度值的均方根误差建立权重关系，根据权重关系融合第一浓度值与第二浓度值，以估算待反演地区的地面PM2.5的最终浓度值。

[0038] 具体地，在将第一浓度值与第二浓度值进行融合的过程中，获取第一浓度值、第二浓度值的均方根误差权重，以及第二浓度值与第一浓度值的均方根误差差值；当差值在预设范围内时，根据第一浓度值与第一浓度值的均方根误差权重、第二浓度值与第二浓度值的均方根误差权重对第一浓度值与第二浓度值进行融合；当差值不在预设范围内时，不对第一浓度值与第二浓度值进行融合。

[0039] 下面对本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的方法进行示例性说明。

[0040] 具体地，本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的方法包括如下步骤：

[0041] 1、搜集地面PM2.5监测站点数据，利用空间变异函数与时间相关函数计算其时间相关特性与空间相关特性，如果PM2.5站点数据之间具有较强的时空相关性，则可选择用时空克里金插值方法进行后续计算。空间变异特征计算方法如公式1所示，时间相关特征计算方法如公式2、公式3所示：

[0042]     (公式1)

[0043]     (公式2)

[0044]     (公式3)

[0045] 其中，公式1中γ(h)为区域化变量在si与si+h处的值Z(si)与Z(si+h)的差的方差之半；h为两点间距离，又称为滞后距离，N(h)为距离h之间用来计算变异函数值的样本对数。

[0046] 公式2，3中，n是时间序列变量yn的数目，yn-t是与yn时间间隔t距离的变量， 是时间距离t的自相关函数， 为均值。

[0047] 2、进行实验时空变异函数计算，针对计算结果选择合适的时空变异函数模型，拟合时空变异函数模型，根据时空变异函数模型的拟合结果，并利用时空克里金模型对待反演地区的地面PM2.5进行插值估算(即基于时空克里金插值模型对PM2.5进行估算)。

[0048] 3、搜集地面气象数据，将数据重采样到创建的网格中，提取PM2.5站点对应网格的AOD数据与气象数据，利用混合效应模型拟合能反映AOD数据与PM2.5数据时间变化特征的模型。首先确定PM2.5估算网格的分辨率。搜集的气象数据若分辨率过高，则对网格内的气象数据计算平均值；若气象数据分辨率过低，则利用距离反比加权的插值方法，将气象数据插值并重采样到与格网分辨率一致。建立AOD数据、PM2.5数据、气象数据随时间变化的混合效应模型，对于时间序列较长的数据，选择分季节建模。

[0049] 由于AOD数据与PM2.5之间关系受气象等各种因子影响，每天都存在变化，仅仅使用简单的多元线性回归模型无法很好的反应这种时间差异性。本过程中使用利用带有随机截距和坡度的混合效应模型建立AOD数据与PM2.5数据之间的关系。针对PM2.5数据与AOD数据之间的相关特性，建立的混合效应模型如下所示：

[0050]     (公式4)

[0051] PMs，t表示在第s个站点第t天的PM2.5浓度观测值；AODs，t表示格网中对应第s个观测站点在第t天的AOD值；b0和b0,t分别表示固定截距和随机截距；b1和b1,t分别表示AOD的固定斜率和随机斜率；b2-b5表示相应气象数据的随机斜率。st表示站点的随机效应的影响，εs，t为系统误差。PM_Ms,t表示CMAQ模式模拟的第s个站点第t天的PM2.5浓度值，PBLs,t表示第s个站点第t天的边界层高度，RHs,t表示第s个站点第t天的相对湿度，Winds,t表示第s个站点第t天的二维风速。统计模型中AOD固定效应代表所有研究时间段内AOD对PM2.5的平均影响，而AOD的随机效应解释了AOD数据与PM2.5数据之间的逐日变化性。

[0052] 4、利用地理加权回归模型建立混合效应模型的残差与AOD数据之间随空间变化的模型。在该步骤中，根据每个PM2.5站点位置的坐标，求取地理加权回归模型的回归参数，在该模型中选用高斯核函数法以及交叉确认法选择核函数和带宽。地理加权回归在时间尺度上是单一的，也就是说每次地理加权回归只能针对一个时间截面，因此，针对全国的数据每天分别使用地理加权回归做一次回归。建立的地理加权回归模型如下所示：

[0053] PM_resist＝β0,s+β1,sAODst+εst   (公式5)

[0054] 其中PM_resist表示上一步骤中利用混合效应模型得到的第t天，站点s的残差部分，AODst表示第t天站点s的AOD值，β0,s与β1,s表示特定地理位置的截距与斜率，β基于高斯权重函数计算的每个地理坐标点的回归参数；εst表示网格单元s在第t天的误差项。

[0055] 5、将第4步与第5步的拟合结果相加，即为利用卫星AOD反演模型(即利用混合效应模型和地理加权回归模型)估算得到的PM2.5浓度值。并利用交叉验证的方法对混合效应模型和地理加权回归模型进行交叉验证。至于如何进行验证不在本发明所述步骤之内，在此不再赘述。

[0056] 6、将基于时空克里金插值模型与基于卫星AOD反演模型此两种模型的估算结果进行融合。计算基于时空克里金插值模型与卫星AOD反演模型的均方根误差RMSE，根据两种模型误差建立权重关系，并将两种模型结果进行融合。在该步骤中根据每天交叉验证的RMSE作为权重来将两种模型结果进行融合。RMSESTK表示基于时空克里金插值模型验证的均方根误差，RMSERS表示基于卫星AOD反演模型验证得到的均方根误差。基于时空克里金插值模型与卫星AOD反演模型的权重分别如以下所示：

[0057] WeightSTK＝RMSERS/(RMSERS+RMSESTK)  (公式6)

[0058] WeightRS＝RMSESTK/(RMSERS+RMSESTK)  (公式7)

[0059] Change＝RMSERS-RMSESTK   (公式8)

[0060] 其中，weightSTK表示基于时空克里金插值模型的均方根误差权重；weightRS表示基于卫星AOD反演模型的均方根误差权重；Change表示两种模型的RMSE差值，当Change在-5到+5μg/m3以内时，则认为两种模型估算精度接近，则利用公式9将两个模型的估算结果进行整合。公式9如下所示：

[0061] PMMerge＝PMSTK*WeightSTK+PMRS*WeightRS    (公式9)

[0062] 其中，PMMerge表示估算的待反演地区的地面PM2.5的最终浓度值；PMSTK表示基于时空克里金插值模型估算得到的PM2.5浓度值；PMRS表示基于卫星AOD反演模型估算得到的PM2.5浓度值。

[0063] 具体地，当Change>5μg/m3时，说明在该天，基于时空克里金插值模型的估算结果明显优于卫星AOD反演模型的估算结果，因此两者不进行融合，直接利用基于时空克里金插值模型的估算结果作为最终PM2.5的估算结果。当Change<-5μg/m3时，说明基于卫星AOD反演模型的估算结果明显优于基于时空克里金插值模型的估算结果，但AOD由于云覆盖等原因存在较多缺失值，此时AOD缺失处利用基于时空克里金插值模型的估算结果作为补充，将最终的结果输出成.tiff格式栅格数据。

[0064] 与上述方法相对应，本发明还提供一种地统计插值与卫星遥感联合反演地面PM2.5的系统，图2示出了根据本发明实施例的地统计插值与卫星遥感联合反演地面PM2.5的系统的逻辑结构。

[0065] 如图2所示，本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的系统200包括第一浓度值获取单元210、第二浓度值获取单元220和PM2.5浓度值估算单元230。

[0066] 其中，第一浓度值获取单元210用于基于时空克里金插值模型和构建的时空变异函数模型，获取待反演地区的地面PM2.5的第一浓度值；其中，获取第一浓度值的过程包括：根据待反演地区的地面PM2.5数据构建时空变异函数模型并拟合，根据对时空变异函数模型的拟合结果，采用时空克里金插值模型估算待反演地区的地面PM2.5的第一浓度值。

[0067] 此外，在获取第一浓度值之前，还可以根据搜集到的地面PM2.5站点监测到的PM2.5数据，利用空间变异函数与时间相关函数计算各PM2.5站点所监测到的PM2.5数据的时间相关特性与空间相关特性。如果PM2.5站点所监测到的PM2.5数据之间具有时空相关性，则可采用地统计插值模型中的时空克里金插值模型估算待反演地区的PM2.5的浓度值(即第一浓度值)。

[0068] 第二浓度值获取单元220用于基于构建的卫星AOD反演模型，获取待反演地区的地面PM2.5的第二浓度值；其中，卫星AOD反演模型包括混合效应模型和地理加权回归模型，获取第二浓度值的过程包括：根据待反演地区的地面PM2.5数据、AOD数据和气象数据构建随时间变化的混合效应模型，并根据混合效应模型对待反演地区的地面PM2.5进行初步估算，以及根据待反演地区的地理坐标构建地理加权回归模型，根据地理加权回归模型获取混合效应模型的残差，根据混合效应模型的初步估算结果和地理加权回归模型所获取的残差估算待反演地区的地面PM2.5的第二浓度值。

[0069] 具体地，在构建混合效应模型的过程中，搜集待反演地区的地面气象数据、AOD数据、地面PM2.5数据，将搜集到的上述数据重采样到创建的与AOD数据分辨率相当的网格后进行数据匹配。其中，上述网格由若干网格单元构成。

[0070] 上述数据匹配过程包括：构建AOD数据像元点的泰森多边形，与网格进行空间叠置分析，将AOD数据分配到网格中的每个网格单元中；将网格中的某个网格单元内所有PM2.5站点同一天采集到的地面PM2.5数据进行平均后赋值给对应的网格单元；将高于待反演地区所对应的网格单元的分辨率的气象数据进行平均后赋值给对应的网格单元，将不高于待反演地区所对应的网格单元的分辨率的气象数据，采用距离反比加权的插值方法将气象数据插值并重采样到与网格单元的分辨率一致。根据上述AOD数据、PM2.5数据和气象数据构建随时间变化的混合效应模型。

[0071] 为了进一步提高PM2.5的估算精度，还需根据待反演地区的地理坐标构建地理加权回归模型，以获取上述混合效应模型的残差。

[0072] 具体地，构建地理加权回归模型的过程包括：根据待反演地区的每个PM2.5站点的位置坐标，获取地理加权回归模型的回归参数；根据回归参数、与PM2.5站点对应的AOD数据、PM2.5站点的位置坐标的截距与斜率获取混合效应模型的残差。

[0073] 为了进一步提高估算精度，对上述混合效应模型和地理加权回归模型进行交叉验证，根据交叉验证结果调整混合效应模型和地理加权回归模型，或确定最终的混合效应模型和地理加权回归模型。将确定的最终的混合效应模型的估算结果和地理加权回归模型所获取的残差相加，获取待反演地区的地面PM2.5的第二浓度值。

[0074] 具体地，当交叉验证结果显示混合效应模型和地理加权回归模型存在过度拟合时，通过调整AOD数据、气象数据对混合效应模型和地理加权模型进行调整，然后对调整后的混合效应模型和地理加权模型再次进行交叉验证，直到交叉验证结果显示调整后的混合效应模型和地理加权模型的过度拟合在预设的范围内时，根据调整后的混合效应模型和地理加权模型估算待反演地区的地面PM2.5的浓度值(即根据确定的最终的混合效应模型和地理加权回归模型估算待反演地区的地面PM2.5的第二浓度值)。

[0075] PM2.5浓度值估算单元230用于将第一浓度值与第二浓度值进行融合，估算待反演地区的地面PM2.5的最终浓度值；其中，根据第一浓度值和第二浓度值的均方根误差建立权重关系，根据权重关系融合所述第一浓度值与第二浓度值以估算待反演地区的地面PM2.5的最终浓度值。

[0076] 此外，PM2.5浓度值估算单元230进一步包括误差权重获取单元231、差值获取单元232和融合单元233。其中，误差权重获取单元231用于分别获取第一浓度值、第二浓度值的均方根误差权重；差值获取单元232用于获取第二浓度值与第一浓度值的均方根误差差值；
融合单元233用于根据误差权重获取单元231所获取的第一浓度值与第一浓度值的均方根误差权重、差值获取单元232所获取的第二浓度值与第二浓度值的均方根误差权重对第一浓度值与第二浓度值进行融合。

[0077] 通过上述可知，本发明提供的地统计插值与卫星遥感联合反演地面PM2.5的方法及系统，将时空克里金插值模型与卫星AOD反演模型相结合，根据误差权重将两种模型的估算结果进行融合，既可以保证PM2.5的估算精度，又可以弥补AOD缺失处PM2.5无法估算的缺点。

[0078] 如上参照附图以示例的方式描述了根据本发明的地统计插值与卫星遥感联合反演地面PM2.5的方法及系统。但是，本领域技术人员应当理解，对于上述本发明所提出的地统计插值与卫星遥感联合反演地面PM2.5的方法及系统，还可以在不脱离本发明内容的基础上做出各种改进。因此，本发明的保护范围应当由所附的权利要求书的内容确定。