CN106528997B - 一种绘制区域颗粒物小时浓度分布图的方法 - Google Patents

一种绘制区域颗粒物小时浓度分布图的方法 Download PDF

Info

Publication number
CN106528997B
CN106528997B CN201610959276.1A CN201610959276A CN106528997B CN 106528997 B CN106528997 B CN 106528997B CN 201610959276 A CN201610959276 A CN 201610959276A CN 106528997 B CN106528997 B CN 106528997B
Authority
CN
China
Prior art keywords
concentration
matrix
particulate matter
data
particle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610959276.1A
Other languages
English (en)
Other versions
CN106528997A (zh
Inventor
张祺杰
刘召贵
彭建波
栾旭东
林远
方华炳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Skyray Instrument Co Ltd
Original Assignee
Jiangsu Skyray Instrument Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu Skyray Instrument Co Ltd filed Critical Jiangsu Skyray Instrument Co Ltd
Priority to CN201610959276.1A priority Critical patent/CN106528997B/zh
Priority to FR1670765A priority patent/FR3058221A1/fr
Publication of CN106528997A publication Critical patent/CN106528997A/zh
Application granted granted Critical
Publication of CN106528997B publication Critical patent/CN106528997B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06—Investigating concentration of particle suspensions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N2015/0042—Investigating dispersion of solids
    • G01N2015/0046—Investigating dispersion of solids in gas, e.g. smoke

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Image Processing (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及一种绘制区域颗粒物小时浓度分布图的方法；本发明要解决的问题是采集监测区域颗粒物浓度数据并绘制出该区域的颗粒物小时浓度分布图。颗粒物分析仪将监测点颗粒物浓度发送到数据处理平台，监测区域需要安装多个颗粒物分析仪才能很好地将监测区域的浓度分布图绘制出来；本发明通过将监测区域颗粒物分析仪的数量定制为3个或以上，根据分析仪器摆放位置形成原始浓度数据矩阵，将每个颗粒物分析仪分析得到的浓度值放入矩阵，用拉普拉斯插入法根据颗粒物分析仪监测浓度数据产生矩阵中非颗粒物分析仪监测浓度值，再用双三次插值法对浓度矩阵中数据进行平滑，然后对矩阵中各个浓度值着色成对应颜色值，从而绘制出颗粒物浓度区域分布图；其优点在于减少了颗粒物分析仪的数量，并且适当扩大了测量区域，绘制的分布图更直观、更精确。

Description

一种绘制区域颗粒物小时浓度分布图的方法
技术领域
本发明涉及一种绘制区域颗粒物小时浓度分布图的方法，具体涉及区域颗粒物浓度数据采集形成浓度矩阵、拉普拉斯插入法完成浓度矩阵数据、双三次插值法平滑浓度矩阵数据和绘制区域浓度分布图的方法。
背景技术
区域颗粒物小时浓度分布图，是通过在监测区域安装多个颗粒物分析仪，使用经纬度坐标将监测区域形成矩形监测区。要想颗粒物浓度分布图更精确，必须安装足够多的颗粒物分析仪。
每台颗粒物分析仪连续监测被监测点颗粒物浓度，并将每小时监测得到的颗粒物浓度发送到数据处理平台。由于颗粒物分析仪数量不可能平均分布在被监测区域，所以需要使用插值算法通过周边颗粒物分析仪得到的浓度数据计算没有安装颗粒物分析仪地方的颗粒物浓度。
本发明从区域颗粒物浓度数据采集形成浓度矩阵、拉普拉斯插入法完成浓度矩阵数据、双三次插值法平滑浓度矩阵数据和绘制区域浓度分布图各个阶段入手，涉及一种绘制区域颗粒物浓度分布图的方法。
发明内容
为了直观、精确地绘制区域颗粒物小时浓度分布图。本发明内容如下：
1.区域颗粒物浓度数据采集形成浓度矩阵
(1)在监测区域安装不少于3台颗粒物分析仪，各颗粒物分析仪尽量安装在监测区域边界附近且间距尽量大。
(2)各颗粒物分析仪连续监测监测点颗粒物；每个小时将监测到的颗粒物浓度平均值通过无线数据传输单元传输到数据处理平台。
(3)数据处理平台根据各颗粒物分析仪某小时的浓度平均值初始化区域浓度矩阵。数据处理平台根据安装的颗粒物分析仪经纬度坐标，图形点数精度值形成一个覆盖监测区域的矩阵，对矩阵中浓度初始值设置为0，根据各颗粒物分析仪经纬度坐标，将数据处理平台选取的某小时颗粒物浓度值设置到矩阵中。
2.拉普拉斯插入法完成浓度矩阵数据
(1)先将浓度矩阵中初始的0值设置为1.e99。
(2)使用拉普拉斯插入法计算值为1.e99的点的值。
3.双三次插值法平滑浓度矩阵数据
通过双三次插值法重新计算矩阵中从第三行到倒数第二行，第二列到倒数第三列的值。
4.绘制区域浓度分布图
(1)设置浓度区间对应的显示颜色。
(2)新建图片，将图片高度宽度定义为矩阵行列数。
(3)根据浓度矩阵中浓度值，将图片的对应位置像素着色为浓度对应的显示颜色。
附图说明
图一为本发明实施流程图
图二为本发明绘制的效果图
具体实施方式
下面对照图一对本发明实施过程进行详细描述。
步骤1：安装颗粒物分析仪，在监测区域安装4台颗粒物(PM2.5)分析仪，尽量安装在监测区域的四个角上。
步骤2：颗粒物分析仪发送数据到数据处理平台，4台颗粒物分析仪连续监测监测点PM2.5浓度，每小时将PM2.5浓度平均值通过无线数据传输单元(dtu数据传输模块)发送到数据处理平台。
步骤3：形成原始浓度矩阵，先得到监测区域坐标最小经度和最大经度值，再得到监测区域纬度最小值和最大值。从最小经度到最大经度均分为800份，从最小纬度到最大纬度值均分为600份，初始化原始浓度矩阵大小为800*600，根据4台颗粒物分析仪经纬度坐标将其小时平均浓度值设置到矩阵对应位置。
步骤4：拉普拉斯插入法完成浓度矩阵数据，先将原始浓度矩阵中非颗粒物分析仪监测到的浓度值设置为1.e99，再通过拉普拉斯插入法重新计算浓度值为1.e99的位置的浓度。
步骤5：初始化结果矩阵，定义一个与原始浓度矩阵大小相同的实数矩阵。初始值都设置为浓度矩阵数据对应坐标值。
步骤6：使用双三次插值法平滑浓度矩阵数据，从矩阵第3行开始到倒数第2行，对每行第2列开始到倒数第3列。通过以下6部得到第i行第j列位置数据。
①y0＝yij
②y1＝y(i+1)j
③y2＝y(i+1)(j+1)
④y3＝yi(j+1)
⑤t＝j/c,u＝i/r，其中r和c分别为数据总行数和列数
⑥y(i,j)＝(1-t)(1-u)y0+t(1-u)y1+tuy2+(1-t)uy3
步骤7：设置浓度对应的显示颜色，浓度小于等于35颜色为绿色；浓度小于等于75颜色为黄色；浓度小于等于115颜色为橙色；浓度小于等于150颜色为红色；浓度小于等于250颜色为紫色；浓度小于等于500颜色为褐紫色；浓度大于500的颜色为黑色。
步骤8：初始化绘制图片，新建图片，将图片宽度定义为结果浓度矩阵列数，将图片高度设置为结果浓度矩阵行数。
步骤9：为图片各像素着色，根据浓度矩阵中浓度值，将图片的对应位置像素着色为浓度对应的显示颜色。

Claims (4)

1.一种绘制区域颗粒物小时浓度分布图的方法，特征为包括：
区域颗粒物浓度数据采集形成浓度矩阵、拉普拉斯插入法完成浓度矩阵数据、双三次插值法平滑浓度矩阵数据和绘制区域浓度分布图；
其中，所述区域颗粒物浓度数据采集形成浓度矩阵包括以下步骤：
(1)在监测区域安装不少于3台颗粒物分析仪，各颗粒物分析仪安装在监测区域的四个角上；
(2)各颗粒物分析仪连续监测监测点颗粒物；每个小时将监测到的颗粒物浓度平均值通过无线数据传输单元传输到数据处理平台；
(3)形成原始浓度矩阵，先得到监测区域坐标最小经度和最大经度值，再得到监测区域纬度最小值和最大值；从最小经度到最大经度均分为800份，从最小纬度到最大纬度值均分为600份，初始化原始浓度矩阵大小为800*600，根据每台颗粒物分析仪经纬度坐标将其小时平均浓度值设置到矩阵对应位置；
所述双三次插值法平滑浓度矩阵数据包括如下步骤：
使用双三次插值法从矩阵第3行开始到倒数第2行，对每行第2列开始到倒数第3列，通过下6部得到第i行第j列位置数进而得到平滑浓度矩阵数据：
①y0＝yij
②y1＝y(i+1)j
③y2＝y(i+1)(j+1)
④y3＝yi(j+1)
⑤t＝j/c，u＝i/r，其中r和c分别为数据总行数和列数
⑥y(i，j)＝(1-t)(1-u)y0+t(1-u)y1+tuy2+(1-t)uy3。
2.如权利要求1所述的绘制区域颗粒物小时浓度分布图的方法，所述的颗粒物浓度指监测区域安装的颗粒物分析仪每小时监测得到的颗粒物浓度的平均值；颗粒物种类根据颗粒物分析仪的采样切割器种类不同分为PM2.5、PM10、TSP。
3.如权利要求1所述的绘制区域颗粒物小时浓度分布图的方法，所述的拉普拉斯插入法完成浓度矩阵数据特征为：
(1)先将浓度矩阵中初始的0值设置为1.e99；
(2)使用拉普拉斯插入法计算值为1.e99的点的值。
4.如权利要求1所述的绘制区域颗粒物小时浓度分布图的方法，所述的绘制区域浓度分布图特征为：
(1)设置浓度区间对应的显示颜色；
(2)新建图片，将图片像素定义为矩阵行列数；
(3)根据浓度矩阵中浓度值，将图片的对应位置像素着色为浓度对应的显示颜色。
CN201610959276.1A 2016-10-28 2016-10-28 一种绘制区域颗粒物小时浓度分布图的方法 Active CN106528997B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610959276.1A CN106528997B (zh) 2016-10-28 2016-10-28 一种绘制区域颗粒物小时浓度分布图的方法
FR1670765A FR3058221A1 (fr) 2016-10-28 2016-12-18 Methode pour tracer le profil de concentration horaire de matieres particulaires dans une zone.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610959276.1A CN106528997B (zh) 2016-10-28 2016-10-28 一种绘制区域颗粒物小时浓度分布图的方法

Publications (2)

Publication Number Publication Date
CN106528997A CN106528997A (zh) 2017-03-22
CN106528997B true CN106528997B (zh) 2020-08-14

Family

ID=58325838

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610959276.1A Active CN106528997B (zh) 2016-10-28 2016-10-28 一种绘制区域颗粒物小时浓度分布图的方法

Country Status (2)

Country Link
CN (1) CN106528997B (zh)
FR (1) FR3058221A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110058548A (zh) * 2019-04-18 2019-07-26 深圳市赛为智能股份有限公司 降低城市开放空间颗粒物浓度的方法、装置及计算机设备
CN110910480A (zh) * 2019-09-29 2020-03-24 谢国宇 基于颜色模式映射关系的环境监测图像渲染方法
CN111125206B (zh) * 2019-12-26 2020-11-17 中科三清科技有限公司 空气污染物数据的处理方法和装置
CN111912755B (zh) * 2020-08-07 2021-08-10 山东中煤工矿物资集团有限公司 一种矿用粉尘浓度传感器、传感器系统及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101900662A (zh) * 2010-08-02 2010-12-01 山东大学 基于粒子系统的区域污染物浓度可视化方法
CN102254330A (zh) * 2010-07-29 2011-11-23 山东大学 基于图像处理的空气污染数据可视化方法
CN102353622A (zh) * 2011-07-01 2012-02-15 黑龙江科技学院 一种煤矿井下工作面粉尘浓度的监测及测量方法
CN103674794A (zh) * 2013-12-16 2014-03-26 中国科学院遥感与数字地球研究所 遥感监测近地表细颗粒物质量浓度pm2.5的多元回归方法
CN104483248A (zh) * 2014-12-30 2015-04-01 中科宇图天下科技有限公司 一种区域颗粒物立体监测系统和方法
CN105678085A (zh) * 2016-01-12 2016-06-15 环境保护部卫星环境应用中心 一种pm2.5浓度的估算方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012165991A1 (en) * 2011-05-31 2012-12-06 Schlumberger Holdings Limited Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102254330A (zh) * 2010-07-29 2011-11-23 山东大学 基于图像处理的空气污染数据可视化方法
CN101900662A (zh) * 2010-08-02 2010-12-01 山东大学 基于粒子系统的区域污染物浓度可视化方法
CN102353622A (zh) * 2011-07-01 2012-02-15 黑龙江科技学院 一种煤矿井下工作面粉尘浓度的监测及测量方法
CN103674794A (zh) * 2013-12-16 2014-03-26 中国科学院遥感与数字地球研究所 遥感监测近地表细颗粒物质量浓度pm2.5的多元回归方法
CN104483248A (zh) * 2014-12-30 2015-04-01 中科宇图天下科技有限公司 一种区域颗粒物立体监测系统和方法
CN105678085A (zh) * 2016-01-12 2016-06-15 环境保护部卫星环境应用中心 一种pm2.5浓度的估算方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"城市热岛与空气污染物特性分析及其影响参数研究";贺广兴;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20150215(第02期);第B027-8页 *

Also Published As

Publication number Publication date
FR3058221A1 (fr) 2018-05-04
CN106528997A (zh) 2017-03-22

Similar Documents

Publication Publication Date Title
CN106528997B (zh) 一种绘制区域颗粒物小时浓度分布图的方法
CN105157590B (zh) 一种基于三维激光扫描技术的建构物健康监测系统
CN104679822B (zh) 一种地理国情外业调绘方法
US9123174B2 (en) Method and apparatus for displaying a simulated application of at least one coating to a digital image
CN102901444B (zh) 一种基于mp小波滤波的零件尺寸检测方法及其检测系统
CN104200193A (zh) 一种果树产量估测方法及装置
CN103914860A (zh) 晶体条位置查找表生成方法及装置
CN103366045B (zh) 基于格子Boltzmann的流体可视化仿真方法
US11068714B2 (en) Method and apparatus for determining hot zones of human flow based on cartography
JP2018066943A (ja) 地目変化判読支援装置、地目変化判読支援方法、及びプログラム
CN104980730A (zh) 一种基于同心圆理论定位光心的方法
CN117131441B (zh) 夜间光污染监测方法、装置、计算机设备和存储介质
CN103340600A (zh) 监护仪上波形失真处理方法及处理系统
CN111428285A (zh) 一种噪音评估方法、装置及存储介质
CN116681959A (zh) 基于机器学习的锋线识别方法、装置、存储介质及终端
CN107230245B (zh) 一种基于Unity3D着色器显示机房地面温场分布的方法和装置
CN103714552B (zh) 运动阴影消除方法及装置和智能视频分析系统
CN113160340A (zh) 城市色彩定量分析与评价方法、系统及存储介质
CN113686513B (zh) 液体泄露定位方法、装置和电子设备
CN113724179B (zh) 一种空间亮度评价指标的计算方法及装置
CN116844075B (zh) 一种耕地环境判定方法及系统
CN112667715B (zh) 基于航班收益数据的可视化分析方法、系统及电子设备
CN219417227U (zh) 一种植被覆盖度测量系统
CN114581786B (zh) 一种根据地面图像估测建筑面积的方法和装置
CN107168933B (zh) 一种空间影响区域分析方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhang Qijie

Inventor after: Liu Zhaogui

Inventor after: Peng Jianbo

Inventor after: Luan Xudong

Inventor after: Lin Yuan

Inventor after: Fang Huabing

Inventor before: Liu Zhaogui

Inventor before: Peng Jianbo

Inventor before: Luan Xudong

Inventor before: Lin Yuan

Inventor before: Fang Huabing

CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhang Qijie

Inventor after: Liu Zhaogui

Inventor after: Peng Jianbo

Inventor after: Luan Xudong

Inventor after: Lin Yuan

Inventor after: Fang Huabing

Inventor before: Zhang Qijie

Inventor before: Liu Zhaogui

Inventor before: Peng Jianbo

Inventor before: Luan Xudong

Inventor before: Lin Yuan

Inventor before: Fang Huabing

GR01 Patent grant
GR01 Patent grant