首页 > 专利 > 杭州电子科技大学 > 一种基于城市区域网格自适应的PM2.5浓度推测方法专利详情

一种基于城市区域网格自适应的PM2.5浓度推测方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2016-03-15
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2016-08-03
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2018-08-07
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2036-03-15
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201610146147.0 申请日 2016-03-15
公开/公告号 CN105740643B 公开/公告日 2018-08-07
授权日 2018-08-07 预估到期日 2036-03-15
申请年 2016年 公开/公告年 2018年
缴费截止日
分类号 G06F17/18G06N3/02 主分类号 G06F17/18
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 0
引用专利数量 0 被引证专利数量 0
非专利引证
引用专利 被引证专利
专利权维持 2 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 戴国骏、郭鸿杰、张桦、吴以凡、仇建 第一发明人 戴国骏
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
浙江杭州金通专利事务所有限公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
王佳健
摘要
本发明涉及一种基于城市区域网格自适应的PM2.5浓度推测方法。本发明首先将城市区域网格化,利用出租车采集原始PM2.5浓度数据,再制定标准,将网格分辨率标准化,利用提取到的网格特征离线训练网格分辨率细化等级和网格特征间的对应关系,实时推测城市浓度时,首先利用提取到的实时网格特征和训练模型,计算出网格所要调整的分辨率,重新划分网格,最后利用线性回归模型推测出每个网格的PM2.5浓度数据,得到城市区域PM2.5浓度数据分布图。本发明方法系统可扩展性强,精确度高,计算量小,为居民出行活动提供参考,并帮助管理部门寻找污染源,改善城市环境质量。
  • 摘要附图
    一种基于城市区域网格自适应的PM2.5浓度推测方法
  • 说明书附图:图1
    一种基于城市区域网格自适应的PM2.5浓度推测方法
  • 说明书附图:图2
    一种基于城市区域网格自适应的PM2.5浓度推测方法
  • 说明书附图:图3
    一种基于城市区域网格自适应的PM2.5浓度推测方法
  • 说明书附图:图4
    一种基于城市区域网格自适应的PM2.5浓度推测方法
  • 说明书附图:图5
    一种基于城市区域网格自适应的PM2.5浓度推测方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2018-08-07 授权
2 2016-08-03 实质审查的生效 IPC(主分类): G06F 19/00 专利申请号: 201610146147.0 申请日: 2016.03.15
3 2016-07-06 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种基于城市区域网格自适应的PM2.5浓度推测方法,其特征在于该方法的具体步骤是:
步骤(1)、城市区域网格化,并采集城市PM2.5原始数据;
具体是将城市区域划分成500m*500m大小的初始分辨率网格,车顶安装PM2.5采集设备的出租车随机行驶于城市道路,采集城市区域PM2.5浓度数据;
步骤(2)、网格分辨率标准化;
定义四种不同分辨率的网格,分别为1000m*1000m,500m*500m,250m*250m,125m*125m,并定义四种相应的网格分辨率细化等级,分别为-1,0,1,2;制定网格分辨率细化等级计算,如下式所示:
MDV=(|xij-xi-1j|+|xij-xi+1j|+|xij-xij-1|+|xij-xij-1|)/4
其中xij表示城市区域网格中第i行,第j列的网格PM2.5浓度,MDV表示中心网格PM2.5浓度和其四周网格PM2.5浓度的平均差值,L表示网格分辨率细化等级,其中:
3
等级-1:当一个网格和其周围网格的PM2.5浓度都小于75ug/m 时,合并这四个相邻的网格为1000m*1000m分辨率网格;
等级1:当一个网格的PM2.5浓度大于115ug/m3,并且该网格的PM2.5浓度与四周网格的PM2.5浓度平均差值在5ug/m3和10ug/m3之间时,将该网格划分成四个相同的,分辨率为250m*
250m的网格;
等级2:当一个网格的PM2.5浓度大于115ug/m3,并且该网格的PM2.5浓度与四周网格的PM2.5浓度平均差值大于10ug/m3时,将该网格划分成十六个相同的,分辨率为125m*125m的网格;
等级0:其他网格保持其初始500m*500m分辨率;
步骤(3)、网格特征提取;
定义六种网格特征:
温度特征Ft:该特征表示初始分辨率下各个网格区域的气温大小;
湿度特征Fh:该特征表示初始分辨率下各个网格区域的湿度大小;
天气特征Fw:该特征表示初始分辨率下各个网格区域的天气情况,它分为晴,阴,小雨,大雨,雪;
风力特征Fw_p:该特征表示初始分辨率下各个网格区域的风力大小;
交通状况特征Ftr:该特征表示初始分辨率下各个网格区域的交通状况,它分为畅通,缓行,拥挤,严重拥挤;
位置特征Fl:该特征表示初始分辨率下各个网格区域的地理位置;
步骤(4)、离线训练;
利用人工神经网络模型训练网格分辨率细化等级与网格特征之间的对应关系,得到网格分辨率细化模型;
步骤(5)、重新划分网格;
依据步骤(4)训练获得的网格分辨率细化模型和步骤(3)实时提取的网格特征,推测各个初始分辨率的网格需要细化的等级,自适应调整网格大小,重新划分整个城市区域的网格;
步骤(6)、在线推测;
根据步骤(1)实时采集的PM2.5浓度原始数据,利用线性回归模型,在自适应分辨率网格下推测无原始数据网格的PM2.5浓度;具体为:首先计算得到回归系数a,b,如下式所示;
其中X为区域内所有网格的温度,湿度,风力,天气,交通状况,位置特征集,X为n×6矩阵,XT为矩阵X的转置矩阵,XT为6×n, 为所有网格的特征平均值矩阵, 为1×6矩阵,y为所有网格的PM2.5浓度值,y为n×1矩阵,为所有网格的PM2.5浓度平均值,n为城市区域网格个数,得到系数b为6×1矩阵,系数a为常数;
然后建立城市任意网格PM2.5浓度值与网格特征之间的线性回归模型;
Y=xb+a
其中Y为任意网格的PM2.5浓度值,x为对应网格的网格特征集;
最后,利用此回归模型实时推测城市区域PM2.5浓度分布。
说明书

技术领域

[0001] 本发明涉及传感器网络技术、计算机应用技术和数据挖掘技术,是一种基于城市区域网格自适应的PM2.5浓度推测方法。

背景技术

[0002] PM2.5被认为是对公众健康和环境造成最大影响的污染物,特别是像北京这样的发展中国家大城市更是深受PM2.5的危害。由于复杂的城市结构和多样的城市功能区,传统的监测站监测PM2.5的方法根本无法告诉城市居民实时准确的PM2.5浓度。
[0003] 近年来,许多学者试图提供更加细粒度的城市污染物浓度分布,他们将PM2.5采集装置安装在公交车或者出租车顶,让交通工具实时自主的采集城市PM2.5数据,结合机器学习,数据挖掘等技术推测出交通工具未覆盖区域的PM2.5浓度,最后给出城市区域细粒度,固定分辨率的PM2.5浓度分布图。这种方法能提供精确的城市PM2.5浓度分布,但是其将城市区域划分成 100m*100m或者200m*200m的过小的固定分辨率网格增加了算法计算量,降低了系统扩展性,不利于推广到大范围PM2.5浓度监测。
[0004] 总的来说,城市区域网格自适应调整下推测城市PM2.5浓度十分重要,过大的网格导致无法接受的推测误差,而过小的网格导致过多的计算量。

发明内容

[0005] 本发明针对现有城市网格划分不足和传统监测PM2.5浓度的缺点,结合传感器网络,数据挖掘等技术,提出了一种基于城市区域网格自适应的 PM2.5浓度推测方法。
[0006] 本发明主要由以下几个步骤构成:1、城市区域网格化,并采集城市 PM2.5原始数据2、网格分辨率标准化3、网格特征提取4、离线训练5、重新划分网格6、在线推测。
[0007] 本发明方法的具体步骤是:
[0008] 步骤(1)、城市区域网格化,并采集城市PM2.5原始数据。具体是将城市区域划分成500m*500m大小的初始分辨率网格,车顶安装PM2.5采集设备的出租车随机行驶于城市道路,采集城市区域PM2.5浓度数据。
[0009] 步骤(2)、网格分辨率标准化。本发明方法定义了4种不同分辨率的网格,分别为1000m*1000m,500m*500m,250m*250m,125m*125m,并定义了4种相应的网格分辨率细化等级,分别为-1,0,1,2。根据PM2.5国家标准和PM2.5浓度变化对人体健康的影响,本发明方法制定了网格分辨率细化等级计算方法如下式所示:
[0010] MDV=(|xij-xi-1j|+|xij-xi+1j|+|xij-xij-1|+|xij-xij-1|)/4
[0011]
[0012] 其中xij表示城市区域网格中第i行,第j列的网格PM2.5浓度,MDV 表示中心网格PM2.5浓度和其四周网格PM2.5浓度的平均差值,L表示网格分辨率细化等级。以下对该公式做进一步说明:
[0013] 等级-1:当一个网格和其周围网格的PM2.5浓度都小于75ug/m3(空气质量等级为“良”)时,合并这4个相邻的网格为1000m*1000m分辨率网格。
[0014] 等级1:当一个网格的PM2.5浓度大于115ug/m3(空气质量等级为“中度污染”),并且该网格的PM2.5浓度与四周网格的PM2.5浓度平均差值在 5ug/m3和10ug/m3之间时,将该网格划分成4个相同的,分辨率为 250m*250m的网格。
[0015] 等级2:当一个网格的PM2.5浓度大于115ug/m3(空气质量等级为“中度污染”),并且该网格的PM2.5浓度与四周网格的PM2.5浓度平均差值大于10ug/m3时,将该网格划分成16个相同的,分辨率为125m*125m的网格。
[0016] 等级0:其他网格保持其初始500m*500m分辨率。
[0017] 步骤(3)、网格特征提取。PM2.5浓度受温度,湿度,风力,交通状况等特征因素影响,据此本发明方法定义了6种网格特征:
[0018] 温度特征(Ft):该特征表示初始分辨率下各个网格区域的气温大小。
[0019] 湿度特征(Fh):该特征表示初始分辨率下各个网格区域的湿度大小。
[0020] 天气特征(Fw):该特征表示初始分辨率下各个网格区域的天气情况,它分为晴,阴,小雨,大雨,雪。
[0021] 风力特征(Fw_p):该特征表示初始分辨率下各个网格区域的风力大小。
[0022] 交通状况特征(Ftr):该特征表示初始分辨率下各个网格区域的交通状况,它分为畅通,缓行,拥挤,严重拥挤。
[0023] 位置特征(Fl):该特征表示初始分辨率下各个网格区域的地理位置。
[0024] 步骤(4)、离线训练。本发明方法利用人工神经网络(Artificial Neutral Network,ANN)模型训练网格分辨率细化等级与网格特征之间的对应关系,得到网格分辨率细化模型。
[0025] 步骤(5)、重新划分网格。依据步骤(4)训练获得的网格分辨率细化模型和步骤(3)实时提取的网格特征,推测各个初始分辨率的网格需要细化的等级,自适应调整网格大小,重新划分整个城市区域的网格。
[0026] 步骤(6)、在线推测。根据步骤(1)实时采集的PM2.5浓度原始数据,利用线性回归模型(Linear Regression,LR)在自适应分辨率网格下推测无原始数据网格的PM2.5浓度。具体为:首先计算得到回归系数a,b,如下式所示;
[0027]
[0028] 其中X为区域内所有网格的温度,湿度,风力,天气,交通状况,位置特征集,X为n×6矩阵,XT为矩阵X的转置矩阵,XT为6×n, 为所有网格的特征平均值矩阵,为1×6矩阵,y为所有网格的PM2.5浓度值,y为n×1 矩阵,为所有网格的PM2.5浓度平均值,n为城市区域网格个数,得到系数b为6×1矩阵,系数a为常数。
[0029] 然后建立城市任意网格PM2.5浓度值与网格特征之间的线性回归模型;
[0030] Y=xb+a
[0031] 其中Y为任意网格的PM2.5浓度值,x为对应网格的网格特征集;
[0032] 最后,利用此回归模型实时推测城市区域PM2.5浓度分布。
[0033] 本发明的有益效果是:
[0034] (1)、本发明方法系统可扩展性强,可以通过调整初始分辨率被运用到大范围的城市空气质量监测;
[0035] (2)、本发明方法自适应调整网格大小,确保了推测精度前提下,大大降低了计算复杂度;
[0036] (3)、本发明方法细粒度城市区域PM2.5浓度推测方法可以告知城市居民实时确切的空气质量情况,为居民出行活动提供参考,并帮助管理部门寻找污染源,改善城市环境质量。

实施方案

[0042] 以下结合附图对本发明作进一步说明。
[0043] 如图1所示,一种基于城市区域网格自适应的PM2.5浓度推测方法具体步骤是:
[0044] 步骤(1)、城市区域网格化,并采集城市PM2.5原始数据。如图2所示,将城市区域划分成500m*500m大小的初始分辨率网格,车顶安装PM2.5采集设备的出租车随机行驶于城市道路,采集城市区域PM2.5浓度数据,并将实时采集的数据通过通用分组无线业务(General Packet Radio Service, GPRS)传输至服务器,图中标圆圈网格代表车辆覆盖的网格,即有PM2.5原始数据的网格。
[0045] 步骤(2)、网格分辨率标准化。本发明方法定义了4种不同分辨率的网格,分别为1000m*1000m,500m*500m,250m*250m,125m*125m,并定义了4种相应的网格分辨率细化等级,分别为-1,0,1,2。根据PM2.5国家标准和PM2.5浓度变化对人体健康的影响,本发明方法制定了网格分辨率细化等级计算方法如下式所示:
[0046] MDV=(|xij-xi-1j|+|xij-xi+1j|+|xij-xij-1|+|xij-xij-1|)/4
[0047]
[0048] 其中xij表示城市区域网格中第i行,第j列的网格PM2.5浓度,MDV 表示中心网格PM2.5浓度和其四周网格PM2.5浓度的平均差值,L表示网格分辨率细化等级。以下对该公式做进一步说明:
[0049] 等级-1:当一个网格和其周围网格的PM2.5浓度都小于75ug/m3(空气质量等级为“良”)时,合并这4个相邻的网格为1000m*1000m分辨率网格。
[0050] 等级1:当一个网格的PM2.5浓度大于115ug/m3(空气质量等级为“中度污染”),并且该网格的PM2.5浓度与四周网格的PM2.5浓度平均差值在 5ug/m3和10ug/m3之间时,将该网格划分成4个相同的,分辨率为 250m*250m的网格。
[0051] 等级2:当一个网格的PM2.5浓度大于115ug/m3(空气质量等级为“中度污染”),并且该网格的PM2.5浓度与四周网格的PM2.5浓度平均差值大于10ug/m3时,将该网格划分成16个相同的,分辨率为125m*125m的网格。
[0052] 等级0:其他网格保持其初始500m*500m分辨率。
[0053] 步骤(3)、网格特征提取。PM2.5浓度受温度,湿度,风力,交通状况等特征因素影响,据此本发明方法定义了6种网格特征:
[0054] 温度特征(Ft):该特征表示初始分辨率下各个网格区域的气温大小。
[0055] 湿度特征(Fh):该特征表示初始分辨率下各个网格区域的湿度大小。
[0056] 天气特征(Fw):该特征表示初始分辨率下各个网格区域的天气情况,它分为晴,阴,小雨,大雨,雪。
[0057] 风力特征(Fw_p):该特征表示初始分辨率下各个网格区域的风力大小。
[0058] 交通状况特征(Ftr):该特征表示初始分辨率下各个网格区域的交通状况,它分为畅通,缓行,拥挤,严重拥挤。
[0059] 位置特征(Fl):该特征表示初始分辨率下各个网格区域的地理位置。
[0060] 步骤(4)、离线训练。本发明方法利用人工神经网络(Artificial Neutral Network,ANN)模型训练网格分辨率细化等级与网格特征之间的对应关系,得到网格分辨率细化模型。如图3所示,本发明方法人工神经网络包括三层:输入层,隐藏层,输出层,其中输入层有8个节点,输出层有1 个节点, 和Lk分别表示网格k的温度特征,湿度特征,天气特征,风力特征,交通状况特征,位置特征和网格分辨率细化等级。
[0061] 步骤(5)、重新划分网格。依据步骤(4)训练获得的网格分辨率细化模型和步骤(3)实时提取的网格特征,推测各个初始分辨率的网格需要细化的等级,自适应调整网格大小,重新划分整个城市区域的网格。如图4 所示,将初始分辨率城市区域的每个网格特征提取出来输入步骤(4)获得的离线训练模型,输出每个网格相应的细化等级,依据其细化等级,自适应调整网格大小,得到自适应分辨率城市区域网格。
[0062] 步骤(6)、在线推测。根据步骤(1)实时采集的PM2.5浓度原始数据,利用线性回归模型(Linear Regression,LR)在自适应分辨率网格下推测无原始数据网格的PM2.5浓度。具体为:首先计算得到回归系数a,b,如下式所示;
[0063]
[0064] 其中X为区域内所有网格的温度,湿度,风力,天气,交通状况,位置特征集,X为n×6矩阵,XT为矩阵X的转置矩阵,XT为6×n, 为所有网格的特征平均值矩阵, 为1×6矩阵,y为所有网格的PM2.5浓度值,y为n×1 矩阵,为所有网格的PM2.5浓度平均值,n为城市区域网格个数,得到系数b为6×1矩阵,系数a为常数。
[0065] 然后建立城市任意网格PM2.5浓度值与网格特征之间的线性回归模型。
[0066] Y=xb+a
[0067] 其中Y为任意网格的PM2.5浓度值,x为对应网格的网格特征集。
[0068] 最后,利用此回归模型实时推测城市区域PM2.5浓度分布,城市区域 PM2.5浓度分布如图5所示。

附图说明

[0037] 图1示出本发明方法的流程图;
[0038] 图2示出城市区域PM2.5原始数据采集场景图;
[0039] 图3示出本发明方法离线训练图;
[0040] 图4示出本发明方法重新划分网格过程图;
[0041] 图5示出网格自适应分辨率下,城市区域PM2.5网格浓度推测结果图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号