一种扬尘源监测方法

摘要

本申请公开了一种扬尘源监测方法，包括：对监测区域的中分辨率卫星影像数据进行处理，获得空间分辨率为10米的多个波段数据产品；在所述多个波段建立空间分辨率为10米包含扬尘源卫星影像的样本集；以所述样本集作为特征向量输入支持向量机，确定支持向量机算法的参数；将所述多个波段数据产品输入所述支持向量机算法，输出各像元为扬尘源的置信度；选择置信度在阈值以上的全部像元。本申请可实现及时准确跟踪、高可靠度显示和监测扬尘源的变化情况。

权利要求

1.一种扬尘源监测方法，其特征在于，包括以下步骤：对监测区域的中分辨率卫星影像数据进行处理，获得空间分辨率为10米的多个波段数据产品；所述中分辨率卫星的时间分辨率小于或等于5天；

在所述多个波段建立空间分辨率为10米包含扬尘源卫星影像的样本集；

所述样本集包括扬尘源、植被、水面、裸地、建筑和道路 这六种地表类型对应的卫星影像数据；样本集按6:4分配为训练集和测试集，即随机挑选样本集中的每一类的60％作为训练集，用于训练分类模型；每一类的剩余40％留作测试集，用于测试模型精度；每类地表类型对应的样本数量不少于1000个；

以所述样本集作为特征向量输入支持向量机，采用GridSearch算法进行10‑折交叉验证方法确定支持向量机算法的参数；所述支持向量机算法是RBF核函数，所述参数是RBF核函数中的C参数和G参数；

将所述多个波段数据产品输入所述支持向量机算法，输出各像元为扬尘源的置信度；

选择置信度在阈值以上的全部像元。

2.根据权利要求1所述的一种扬尘源监测方法，其特征在于，还包括以下步骤：将所述样本集分为训练集和测试集；

将所述训练集的数据输入支持向量机，确定支持向量机算法的参数；

用所述测试集的数据，通过混淆矩阵评价方法确定支持向量机算法的最佳参数。

3.根据权利要求1所述的一种扬尘源监测方法，其特征在于，还包括以下步骤：对所述影像数据进行预处理，处理方式包括选取无云或少云区域的数据，以及几何校正、大气校正、裁剪、镶嵌中的至少一种。

4.根据权利要求1所述的一种扬尘源监测方法，其特征在于，所述中分辨率卫星的空间分辨率为10米～60米。

5.根据权利要求1所述的一种扬尘源监测方法，其特征在于，所述多个波段是波长为

0.490μm、0.560μm、0.665μm、0.705μm、0.740μm、0.783μm、0.842μm、0.865μm、1.610μm、2.190μm的波段。

6.根据权利要求1所述的一种扬尘源监测方法，其特征在于，还包括以下步骤，对波长为0.705μm、0.740μm、0.783μm、0.865μm、1.610μm、2.190μm的波段对应的所述影像数据进行重采样。

说明书

一种扬尘源监测方法

技术领域

[0001] 本申请涉及遥感领域，尤其涉及一种扬尘源监测的方法。

背景技术

[0002] 扬尘是指地表松散颗粒物质在自然力或人力作用下进入到环境空气中形成的一定粒径范围的空气颗粒物。扬尘源是指在自然力或人力作用下各种不经过排气筒、无组织、
无规则排放地表松散颗粒物质的颗粒物排放源。扬尘源主要包括以下几类：土壤扬尘源，道
路扬尘源，施工扬尘源，堆场扬尘源。城市范围内扬尘源中施工扬尘源和堆场扬尘源的分布
最为广泛，也是扬尘管控的关键。

[0003] 利用卫星数据进行扬尘源识别及监测，一方面应满足数据空间分辨率以及时间分辨率的要求，另一方面还应满足卫星数据的可获得性、经济性。

[0004] 现有技术多采用高分辨率遥感卫星，结合人工目视解译提取城市范围内扬尘源的空间分布信息，其优点在于因为卫星影像数据精度高，目视解译结果准确度高，但由于数据
时间分辨率低，无法以月、季度的频次进行提取，无法满足现有的环境监管需求，且目视解
译耗时长，成本高，受人为主观因素的一定影响。中分辨率遥感卫星数据的时间分辨率可以
满足月、季度频次的扬尘源提取工作，但基于中分辨率卫星数据的提取目前多以单一相关
指数(如裸地指数)为基础，通过阈值判断提取扬尘源，这种方法虽然满足环境监管频次要
求，但是准确率较低，误差较大,不具有普遍适用性。也就是说，扬尘源的变化周期很快，仅
凭高分辨率卫星影像无法满足监管需求。而中分辨率卫星识别出的数据类别单一，只能给
出有无扬尘源污染，缺乏带有置信度的高可靠性数据。

[0005] 因此，对于扬尘源监测，需解决的问题是卫星数据及时准确跟踪扬尘源的变化，同时识别结果需带有置信度的数据，提升监测的准确性。

发明内容

[0006] 本发明提出一种扬尘源监测方法，解决现有技术无法及时准确跟踪扬尘源变化，且对扬尘源污染缺乏置信度等定量评价的问题。

[0007] 本申请实施例提供一种扬尘源监测方法，包括以下步骤：

[0008] 对监测区域的中分辨率卫星影像数据进行处理，获得空间分辨率为10米的多个波段数据产品；

[0009] 在所述多个波段建立空间分辨率为10米抱包含扬尘源卫星影像的样本集；

[0010] 以所述样本集作为特征向量输入支持向量机，确定支持向量机算法的参数；

[0011] 将所述多个波段数据产品输入所述支持向量机算法，输出各像元为扬尘源的置信度；

[0012] 选择置信度在阈值以上的全部像元。

[0013] 作为本申请进一步优化的实施例，还包括以下步骤：

[0014] 将所述样本集分为训练集和测试集；

[0015] 将所述训练集的数据输入支持向量机，确定支持向量机算法的参数；

[0016] 用所述测试集的数据，通过混淆矩阵评价方法确定支持向量机算法的最佳参数。

[0017] 作为本申请进一步优化的实施例，还包括以下步骤：对所述影像数据进行预处理，处理方式包括选取无云或少云区域的数据，以及几何校正、大气校正、裁剪、镶嵌中的至少
一种。

[0018] 优选地，所述样本集还包括植被、水面、裸地、建筑和道路的卫星影像中至少一种。

[0019] 优选地，所述支持向量机算法是RBF核函数，所述参数是RBF核函数中的C参数和G参数。

[0020] 优选地，采用GridSearch算法进行10‑折交叉验证方法确定所述支持向量机算法的参数。

[0021] 优选地，所述中分辨率卫星的时间分辨率小于或等于5天。

[0022] 优选地，所述中分辨率卫星的空间分辨率为10米～60米。

[0023] 优选地，所述多个波段是波长为0.490μm、0.560μm、0.665μm、0.705μm、0.740μm、0.783μm、0.842μm、0.865μm、1.610μm、2.190μm的波段。

[0024] 优选地，还包括以下步骤，对波长为0.705μm、0.740μm、0.783μm、0.865μm、1.610μm、2.190μm的波段对应的所述影像数据进行重采样。

[0025] 本申请的至少一个实施例具有以下有益效果：本申请提出基于中分辨率卫星的扬尘源提取方法，其空间分辨率大小满足扬尘源的识别要求，可以及时跟踪扬尘源的变化，提
升监测的准确性。也就是说，在满足监管需求的频率要求的前提下，采用支持向量机模型自
动分类提取扬尘源的方法，得到扬尘源的空间分布信息，准确度高于现有基于中分辨率卫
星数据的其他扬尘源提取技术。当采用哨兵二号系列卫星数据时，该卫星具有高达5天的时
间分辨率和10米的空间分辨率。作为数据源，以一月一次的监测频率，既可以保证去除大部
分云层覆盖的干扰，也可以保证结果的全面准确性，且提取目标不具有移动性和时相变化，
即使研究区域内一个月内无完整一景无云影像，多景合成的无云/少云影像也可以作为数
据进行使用。

[0026] 同时，如果仅用目视解译的方法对中分辨率数据进行解译，并不能够准确快速的提取扬尘源；本申请首次采用支持向量机的分类置信度输出作为扬尘源目视解译的主要判
断依据，根据高置信度点群的分布密度情况，可以简单快速的定位扬尘源的分布范围，提取
扬尘源的空间分布信息。

附图说明

[0027] 此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：

[0028] 图1为扬尘源监测方法的实施例流程图；

[0029] 图2为训练和校对支持向量机算法最优参数的实施例流程图；

[0030] 图3为样本集中的多个波段数据产品组成情况示意图；

[0031] 图4为对监测区域进行扬尘源识别的结果示意图。

具体实施方式

[0032] 为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本申请一
部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做
出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

[0033] 以下结合附图，详细说明本申请各实施例提供的技术方案。

[0034] 图1为扬尘源监测方法的实施例流程图。本申请提出的扬尘源监测方法，至少包括以下步骤：

[0035] 步骤10：对监测区域的中分辨率卫星影像数据进行处理，获得空间分辨率为10米的多个波段数据产品

[0036] 在步骤10中，对监测区域的中分辨率卫星影像数据进行处理，获得N个像元的多个波段数据产品，所述像元空间分辨率为10米。监测区域是以所述监测区域内全部大气国控
站点周边指定距离缓冲区的最大覆盖范围集合确定的，再将所述区域转换为矢量文件。所
述矢量文件为二维的矢量面，是由一系列相连的x，y坐标对定义的闭合形状。

[0037] 例如，以监测区域内全部大气国控站点的周边5km做为缓冲区，所有站点周边5km缓冲区的覆盖范围叠加则为扬尘源监测区域。

[0038] 优选地，所述中分辨率卫星的时间分辨率小于或等于5天。

[0039] 优选地，所述中分辨率卫星的空间分辨率为10米～60米。

[0040] 例如，所述中分辨率卫星是哨兵二号系列卫星，哨兵二号系列卫星可实现时间分辨率小于或等于5天且各波段空间分辨率分别为10m、20m或60m。

[0041] 哨兵二号系列卫星共13个波段，其中：空间分辨率为10m的波段有4个，分别是0.490μm、0.560μm、0.665μm、0.842μm的波段；空间分辨率为20m的波段有6个，分别是0.705μ
m、0.740μm、0.783μm、0.865μm、1.610μm、2.190μm；空间分辨率为60m的波段有3个，分别是
0.443μm、0.945μm、1.375μm。

[0042] 表1、哨兵二号系列卫星波段表

[0043]

[0044] 需要说明的是，本申请不采纳空间分辨率为60m的波段数据，只采纳空间分辨率为10m和20m的10个波段数据。

[0045] 在步骤10中，所述卫星影像数据包含空间分辨率为10m和20m的10个波段数据。其中，空间分辨率为10m的4个波段，即波长为0.490μm、0.560μm、0.665μm、0.842μm的波段；空
间分辨率为20m的6个波段，分别是波长为0.705μm、0.740μm、0.783μm、0.865μm、1.610μm、
2.190μm的波段。

[0046] 优选地，对波长为0.705μm、0.740μm、0.783μm、0.865μm、1.610μm、2.190μm的波段对应的所述影像数据进行重采样处理，获得空间分辨率为10米的数据。

[0047] 优选地，对所述影像数据进行预处理，处理方式包括选取无云或少云区域的数据，以及几何校正、大气校正、裁剪、镶嵌中的至少一种。

[0048] 在步骤10中，对所述影像数据按照无云或少云区域进行选取，再按照所述矢量文件对所述影像数据剪裁或镶嵌，所述裁剪是按照矢量文件选取所述监测区域内的所述影像
数据，所述镶嵌是将所述监测区域边缘穿过所述同一像元内部的所述影像数据拼接为一个
像元。按照所述影像数据类型进行几何和大气校正处理。

[0049] 例如，所述影像数据为经过几何校正和大气校正的哨兵2号卫星L2A级数据产品。如果L2A级数据暂时没有分发，则下载哨兵2号卫星L1C级数据产品，并使用Sen2Cor进行几
何精校正和大气校正处理，得到哨兵2号卫星L2A级数据产品。

[0050] 步骤20：在所述多个波段建立空间分辨率为10米包含扬尘源卫星影像的样本集；

[0051] 例如，所述样本集是对指定地表类型进行高清影像采样，所述多个波段与步骤10中的多个波段相同，得到每个指定地表类型每个像元对应的多个波段数据。

[0052] 优选地，在步骤20中，所述指定地表类别还包括分别植被、水面、裸地、建筑、道路中的至少一种。

[0053] 步骤30：以所述样本集作为特征向量输入支持向量机，确定支持向量机算法的参数；

[0054] 需要说明的是，支持向量机(Support Vector Machine)是一种结构风险最小的分类器，其基本原理是统计学习理论的VC维和结构风险最小化原理，依据有限的训练样本，寻
求模型复杂度(训练样本的精度)和学习能力(新样本的识别能力)的这种方法通过把原始
的特征向量映射到一个高维空间，在高维空间构建了最大边缘间隔的决策平面，也叫超平
面，超平面两边有2个平行的边界平面，这2个平面式距离正负样本最近的平面，SVM要使得
超平面和这2个边界平面的就最大化，距离越大，分类器的误差就越小，精度就越高。

[0055] 需要说明的是，支持向量机的参数选择可以采用RBF,Linear,Poly,Sigmoid，precomputed或者自定义的核函数。

[0056] 优选地，所述支持向量机算法是RBF核函数，所述参数是RBF核函数中的C参数和G参数。

[0057] 本申请采用径向基函数(RBF)核函数算法作为所述支持向量机的核函数，该算法需要优化的参数有C参数和G参数；C参数为惩罚系数，C太小，容易发生过拟合，C太大，拟合
的结果和实际结果相差太大；G(或gamma)参数是RBF函数的参数之一，G用于调整模型的复
杂度，G值越小模型复杂度越低，反之模型复杂度越大。G的取值直接影响到能否找到一个能
将多累数据分开的最优超平面

[0058] 优选的，采用GridSearch算法进行10‑折交叉验证方法确定所述支持向量机算法的最优参数。

[0059] 例如，本申请通过使用GridSearch算法进行10‑折交叉验证方法确定所述支持向量机算法即径向基函数(RBF)核函数算法的参数，即：输入样本集，调整C参数和G参数直到
支持向量机算法输出的地表类别为扬尘源。

[0060] 步骤40：将所述多个波段数据产品输入所述支持向量机算法，输出各像元为扬尘源的置信度；

[0061] 在步骤40中，将所述多个波段数据产品输入所述支持向量机函数，输出结果为每个像元是否归属扬尘源及对应的置信度。

[0062] 在步骤40中，置信度是指输出结果中，每个像元被分类为某一类地表类别的可能性。该像元对应某一类地表类别的置信度越高，说明该像元被分为这一类地表类别的可能
性越大。支持向量机中的置信度，指的是所述数据产品每个像元到分割超平面的欧式距离。

[0063] 将待分类的多波段影像数据集作为输入，使用训练好的支持向量机模型进行分类。因为支持向量机是一种软分类器，所以选择每个待分类像元在每一个类别的置信度作
为输出，得到分类结果。

[0064] 例如，所述支持向量机函数处理结果为所述检测区域内，每个空间分辨率为10米的像元是否为扬尘源及其对应的置信度分类图。每个像元的置信度从0到1以颜色逐渐变深
显示。

[0065] 步骤50：选择置信度在阈值以上的全部像元。

[0066] 在步骤50中，选取所述支持向量机函数输出地表类别为扬尘源，且置信度在阈值以上的全部像元，勾出所述全部像元区域，即为扬尘源分布图。

[0067] 例如，阈值设定为50％，则勾出所述区域中扬尘源置信度在50％以上的全部像元区域，即为扬尘源分布图。阈值设定为60％，则勾出所述区域中扬尘源置信度在60％以上的
全部像元区域，即为扬尘源分布图。

[0068] 图2为训练和校对支持向量机算法最优参数的实施例流程图。本申请提出的扬尘源监测方法，还包括以下步骤：

[0069] 步骤31：将所述样本集分为训练集和测试集；

[0070] 在步骤31中，按6:4的比例将所述样本集分为训练集和测试集。

[0071] 优选的，在步骤31中，样本集除了扬尘源外，还覆盖植被、水面、裸地、建筑、道路中的至少一种地表类型。

[0072] 例如，样本集覆盖扬尘源、植被、水面、裸地、建筑、道路这六种地表类型对应的卫星影像数据。

[0073] 需要说明的是，所述样本集同样包括步骤10中的10个波段数据。

[0074] 例如：样本集包含有六类不同地表类别的卫星影像数据集，每个样本点均包含对应的10个波段数据。将样本集按6:4分配为训练集和测试集，即随机挑选样本集中的每一类
的60％作为训练集，用于训练分类模型；每一类的剩余40％留作测试集，用于测试模型精
度。，每类地表类型对应的样本数量不少于1000个。

[0075] 步骤32：将所述训练集的数据输入支持向量机，确定支持向量机算法的参数；

[0076] 在步骤32中，将所述测试集作为支持向量机算法即径向基函数(RBF)核函数算法的输入，采用GridSearch算法进行10‑折交叉验证方法确定所述支持向量机算法的参数即C
和G，调整C和G直到支持向量机算法输出的地表类别结果与输入的指定地表类别一致。

[0077] 步骤33：用所述测试集的数据，通过混淆矩阵评价方法确定支持向量机算法的最佳参数。

[0078] 在步骤33中，将所述测试集输入步骤32的得到的所述径向基函数(RBF)核函数算法中，通过混淆矩阵评价方法确定最佳所述C和G参数。

[0079] 例如，使用混淆矩阵分析所述径向基函数(RBF)的精度，设定模型分类精度不低于95％，如果分类模型精度低于95％，重新调整C和gamma的备选值，重复上一步，直到模型测
试精度达到95％以上。

[0080] 图3为样本集中的多个波段数据产品组成情况示意图。

[0081] 所述样本集包括N个样本，每个样本对应一个像元，每个样本包括该像元所包含的特征向量以及所对应的类别标签。所述特征向量为每个像元对应的空间分辨率为10米的多
个波段数据，所述多个波段与图1中所述的波段相同，优选为10个。所述类别标签为每个像
元对应的已知所述地表类别，例如是否扬尘源。

[0082] 例如，将监测区域100划分为空间分辨率为10米的M个像元，从M个像元中选取已知地表类别为扬尘源101的N个像元110，以及已知地表类别为非扬尘源120的M‑N个像元102。
所述非扬尘源可以进一步按照植被、水面、裸地、建筑、道路进行分类。

[0083] 优选地，每类地表类别的像元数量不小于1000个。

[0084] 与步骤10的多波段影像数据相对应，对20米空间分辨率的6个波段(波段5 6 7 11 12和8a)和10米空间分辨率的4个波段(波段2 3 4 8)，用已知地表类别的10米分辨率的10
个波段数据，每个像元均用10个波段的影像数据组成该像元的特征向量，如下表所示。

[0085] 表2、样本数据集示例

[0086]像元序号 标签 特征1 特征2 … 特征10
1 扬尘源 0.1038 0.1104 … 0.1269
2 非扬尘源 0.0812 0.0996 … 0.0336
3 扬尘源 0.1084 0.114 … 0.141
… … … … … …

[0087] 图4为对监测区域进行扬尘源识别的结果示意图。

[0088] 对所述监测区域进行扬尘源识别，选择输出结果地表类别为扬尘源、且标注对应置信度的全部像元区域。

[0089] 还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包
括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要
素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要
素的过程、方法、商品或者设备中还存在另外的相同要素。

[0090] 以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同
替换、改进等，均应包含在本申请的权利要求范围之内。