CN111079848A - 一种基于图像的空气质量等级评估方法 - Google Patents
一种基于图像的空气质量等级评估方法 Download PDFInfo
- Publication number
- CN111079848A CN111079848A CN201911333157.5A CN201911333157A CN111079848A CN 111079848 A CN111079848 A CN 111079848A CN 201911333157 A CN201911333157 A CN 201911333157A CN 111079848 A CN111079848 A CN 111079848A
- Authority
- CN
- China
- Prior art keywords
- image
- air quality
- channel
- pixel
- ratio
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/46—Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
- G06V10/462—Salient features, e.g. scale invariant feature transforms [SIFT]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
一种基于图像的空气质量等级评估方法,涉及机器学习相关技术,尤其涉及图像处理技术。目的是为解决现有的空气质量检测成本昂贵且不能覆盖城市中每一个角落的问题。本发明首先对图像进行HSI空间的转换,然后对图像的饱和度通道S和亮度通道I做DCT变换将它们转换到频域,最后提取出图像的边缘像素比、高亮度像素比、低饱和度像素比及低频率像素比作为特征,将其输入到SVM中进行分类,从而构建出空气质量等级的评估模型。该方法所构建的空气质量等级评估模型可以较准确的估计图像采集地的空气质量状况,将空气质量检测变得低成本、细粒度,便于人们随时随地的监督周边的空气污染情况。且该方法更适用于对空气湿度适中、图像内容丰富且天空区域占比小于二分之一的图像进行评估。
Description
技术领域
一种基于图像的空气质量等级评估方法,涉及机器学习相关技术,尤其涉及图像处理技术。
背景技术
近些年来,随着我国工业化进程的不断推进,国家经济飞速发展的同时,环境污染也在不断加剧,绝大多数城市都出现了不同层次的污染。空气污染不仅会对日常生活中的交通起居等造成影响,同时还会危害我们的身体健康,因此,空气污染的日常监测和治理是十分重要的。目前我国获取空气质量的方法主要是在城市所设立监测站中通过各种精密的仪器来检测并计算污染物浓度,再将它换算成空气质量等级或指数,最后发布至气象局供人们查询。但是这种方法的所需仪器的成本较高,且它的检测方式粗粒度的,不能够覆盖城市的每个角落。最近几年,有很多学者提出了将图像处理技术与机器学习方法相结合来检测空气质量的方法。但是受到所选取的特征及机器学习方法的限制,目前所提出的方法都存在着如适用性不佳或图像采集地点受限等缺点。
发明内容
本发明的目的是为解决传统与现有的空气质量检测方法的不足,让空气质量检测变得低成本、细粒度,便于人们随时随地的监督周边的空气污染情况而提供的一种基于图像的空气质量等级评估方法。
该基于图像的空气质量等级评估方法包括以下步骤:
步骤一、从网络上收集具有空气质量等级标签及相关污染物信息的图像,将这些图像做为模型构建所用的空气质量图像库;
步骤二、将图像从RGB空间转移到HSI空间下,再对空气质量图像库中不同空气质量等级的图像进行探究与分析;
步骤三、由步骤二可知,在HSI空间下,空气污染会对图像的饱和度通道和亮度度通道产生较大影响,而对色度通道所产生的影响并不明显,最终选择将边缘梯度、亮度、饱和度以及频域的能量信息作为依据来评估图像的空气质量等级;
步骤四、运用图像处理技术对选定的图像特征进行提取,具体包括:
步骤四一、在HSI空间下的I通道中利用Canny算子来对图像进行边缘检测,计算边缘点数量与图像总像素数量的比值;
步骤四二、图像在HSI空间中的I通道代表了图像的亮度通道,对图像中亮度像素值大于A的像素点进行统计,计算它与图像总像素数量的比值;
步骤四三、图像在HSI空间中的S通道代表了图像的饱和度通道,对图像中饱和度像素值小于B的像素点进行统计,计算它与图像总像素数量的比值;
步骤四四、使用DCT变换将上述图像中的H、S、I分量分别转换到频域空间,并计算S和I通道的频域图像中的低频信息数量与图像总像素点数量的比值;
步骤五、将步骤四中所提取的高亮度像素比、低饱和度像素比、I通道的边缘像素比、S和I通道的低频率像素比共五维特征输入到SVM(支撑向量机)中进行分类,构建出基于图像的空气质量等级评估模型。
进一步地,所述方法还包括:
步骤六、对模型进行参数优化并进行性能与适用性测试。
进一步地,步骤四二中的A=128。
进一步地,步骤四三中的B=20。
上述方法首先将图像转换到HSI空间下,再提取图像的边缘像素比、高亮度像素比、低饱和度像素比及低频率像素比作为特征,然后结合SVM分类器构建了空气质量等级评估模型。通过相关的性能及适用性实验进行测试,实验结果表明,本方法所构建的空气质量等级评估模型可以较准确的估计图像采集地的空气质量状况。该方法有效的解决了现有的空气质量检测成本昂贵且不能覆盖城市中每一个角落的问题,将空气质量检测变得低成本、细粒度。但该方法也存在一定的局限性,更适用于对空气湿度适中、图像内容丰富且天空区域占比小于二分之一的图像进行评估。
附图说明
图1为本发明所述的基于图像的空气质量等级评估方法的流程图;
图2为步骤二中用于探究与分析的不同等级的空气质量图像,(a)、(b)、(c)、(d)分别为空气质量为优、轻度污染、重度污染和严重污染的图像。
图3为步骤四一中对不同空气质量等级图像的I通道进行边缘检测的图像示例,(a)、(b)、(c)、(d)分别为空气质量为优、轻度污染、重度污染和严重污染的图像。
图4为步骤四二中的不同空气质量等级图像的亮度分布直方图示例,(a)、(b)、(c)、(d)分别为空气质量为优、轻度污染、重度污染和严重污染的图像。
图5为步骤四三中的不同空气质量等级图像的饱和度分布直方图示例,(a)、(b)、(c)、(d)分别为空气质量为优、轻度污染、重度污染和严重污染的图像。
图6为步骤六适用性测试中误差过大的部分图像示例,(a)的真实空气质量等级为良,但评估空气质量等级为严重污染,原因是图像中没有天空区域;(b)的真实空气质量等级为优,但评估空气质量等级为严重污染,原因是图像中只有天空区域;(c)的真实空气质量等级为良,但评估空气质量等级为重度污染,原因是空气湿度为80%;(d)的真实空气质量等级为良,但评估空气质量等级为严重污染,原因是空气湿度为82%。
具体实施方式:如图1所示,本实施方式所述的基于图像的空气质量等级评估方法包含以下步骤:
步骤一、由于基于图像评估空气质量等级的研究没有公开的图像库,所以我们从网络上选取了100幅具有空气质量等级标签及相关污染物信息的图像,并将这些图像做为模型构建所用的空气质量图像库。
步骤二、由于HSI空间更符合人眼的视觉特性,且它的分量相关性远小于RGB空间,所以我们先将空气质量图像库中的图像从RGB空间转移到HSI空间,转换公式如下:
然后在HSI空间下对不同空气质量等级的图像进行探究与分析。
步骤三、由步骤二可知,在HSI空间下,空气污染会对图像的饱和度通道和亮度度通道产生较大影响,而对色度通道所产生的影响并不明显,所以空气污染的确会对图像质量造成影响。同时由于空气污染物中的悬浮颗粒物的存在,使得从物体表面反射的光线遭受到了严重的散射,从而使一部分光线无法到达成像设备中,导致图像质量降低,会出现清晰度下降、边缘模糊、亮度变化不均和饱和度下降等现象,严重破坏了图像的视觉效果。从图像的频域角度来看,由于空气污染会造成图像的边界模糊、灰度分布集中和细节丢失等现象,这也使得图像的频域会出现高频和低频信息丢失较为严重而中频部分保留较多的现象。因此,我们选择将边缘梯度、亮度、饱和度以及频域的能量信息作为依据来评估图像的空气质量等级。
步骤四、运用图像处理技术对选定的图像特征进行提取,具体包括:
步骤四一、传统的边缘检测方法侧重于灰度图像的边缘检测,而通过灰度图像的检测方法只能检测到彩色图像90%的边缘信息,但仍然有剩下的10%的边缘信息需要利用颜色信息才能得到。在彩色图像的HSI颜色空间中,彩色图像的边缘是由色度、饱和度和亮度的不连续性造成的,而亮度分量I几乎包含了彩色图像所有边缘信息。所以这里选择在HSI空间下的I(亮度)通道中利用Canny算子来对图像进行边缘检测。
Canny边缘检测的流程如下:首先使用高斯滤波器来平滑图像,滤除噪声,其中高斯滤波器核的生成方程式如下:
其中,x和y分别为模板中心位置点的横、纵坐标,而标准差σ的大小决定了高斯函数的宽度;然后计算每个像素点的梯度强度与方向,而像素点的梯度强度G与梯度方向θ的表达式如下:
其中Gx和Gy分别代表横向与纵向边缘检测的图像灰度值;再利用非极大值抑制来消除边缘检测所带来的杂散响应,它的主要流程为:将当前像素的梯度强度与沿正负梯度方向上的两个像素进行比较,如果当前像素的梯度强度大于另外两个像素的梯度强度,则该像素点保留为边缘点,否则该像素点将被抑制;最后利用双阈值检测来确定真实与潜在的边缘点并将它们连接起来。
通过边缘检测可知,随着空气质量等级的下降,所计算出的边缘点的数量越来越少。最后,计算边缘点数量与图像总像素数量的比值,将其作为特征;
步骤四二、图像在HSI空间中的I通道代表了图像的亮度通道,当空气质量较好时图像的亮度分布更均匀,而随着空气质量等级的下降,整体亮度越来越向高亮度区域集中。对图像中亮度像素值大于A的像素点进行统计。再计算它与图像总像素数量的比值作为特征。随着空气质量等级的下降它的数值越来越大。为了取得更好的分类效果,A的优选值为128;
步骤四三、图像在HSI空间中的S通道代表了图像的饱和度通道,当空气质量较好时图像整体的饱和度较高,而随着空气质量等级的下降,整体饱和度越来越低。对图像中饱和度像素值小于B的像素点进行统计。再计算它与图像总像素数量的比值作为特征。随着空气质量等级的下降它的数值越来越大。为了取得更好的分类效果,B的优选值为20;
步骤四四、使用DCT变换将上述图像中的H、S、I分量分别转换到频域空间,DCT变换的公式如下:
其中,f(i,j)为原始图像中某一点的数值,而原始图像的大小为N*N(一般取N=8),这里若原始图像的长与宽不相等时,需要将其补齐后再进行变换;
经过DCT变换后,会使图像中能量和系数较大的像素点集中在左上角,而能量和系数较小的像素点集中在右下角,即代表目标图像中的轮廓及灰度分布特征的信息集中体现在左上角,而目标图像的细节特征集中体现在右下角。与图像在空间域的表现相似,经DCT变换后的频域图像同样在饱和度(S)通道和亮度(I)通道受到空气污染的影响较大,而在色度(S)通道的影响并不明显,且由于图像的高频信息的系数较小,提取存在困难。所以只计算S和I通道的频域图像中的低频信息数量与图像总像素点数量的比值作为特征。随着空气质量等级的降低,它的数值变得越来越小。
步骤五、将步骤四中所提取的高亮度像素比、低饱和度像素比、I通道的边缘像素比、S和I通道的低频率像素比共五维特征输入到SVM(支撑向量机)中,支撑向量机进行分类的核心思想是将低维空间的线性不可分问题转换至高维空间使其线性可分,在高维特征空间内通过一个超平面来区分不同的空气质量等级。在实际应用时,需要恰当地选择核函数,由于图像分类问题是线性不可分的,所以本方法选用的是RBF径向基内核。整体流程为:首先将图像集按照步骤四中的特征提取流程进行特征提取,然后将特征与其所对应的空气质量等级标签一同输入到支撑向量机中进行分类,最后得到空气质量等级的评估模型。
步骤六、在上述构建模型的过程中,还需要对模型进行参数优化,是利用图像集通过网格搜索与交叉验证相结合的方法来找到最优参数的。最后要进行性能与适用性的测试。测试结果表明,该方法更适用于对空气湿度适中、首要污染物为颗粒物、图像内容丰富且天空区域占比小于二分之一的图像进行评估。
Claims (4)
1.一种基于图像的空气质量等级评估方法,其特征在于,包括以下步骤:
步骤一、从网络上收集具有空气质量等级标签及相关污染物信息的图像,将这些图像做为模型构建所用的空气质量图像库;
步骤二、将图像从RGB空间转移到HSI空间下,再对空气质量图像库中不同空气质量等级的图像进行探究与分析;
步骤三、由步骤二可知,在HSI空间下,空气污染会对图像的饱和度通道和亮度度通道产生较大影响,而对色度通道所产生的影响并不明显,最终选择将边缘梯度、亮度、饱和度以及频域的能量信息作为依据来评估图像的空气质量等级;
步骤四、运用图像处理技术对选定的图像特征进行提取,具体包括:
步骤四一、在HSI空间下的I通道中利用Canny算子来对图像进行边缘检测,计算边缘点数量与图像总像素数量的比值;
步骤四二、图像在HSI空间中的I通道代表了图像的亮度通道,对图像中亮度像素值大于A的像素点进行统计,计算它与图像总像素数量的比值;
步骤四三、图像在HSI空间中的S通道代表了图像的饱和度通道,对图像中饱和度像素值小于B的像素点进行统计,计算它与图像总像素数量的比值;
步骤四四、使用DCT变换将上述图像中的H、S、I分量分别转换到频域空间,并计算S和I通道的频域图像中的低频信息数量与图像总像素点数量的比值;
步骤五、将步骤四中所提取的高亮度像素比、低饱和度像素比、I通道的边缘像素比、S和I通道的低频率像素比共五维特征输入到SVM(支撑向量机)中进行分类,构建出基于图像的空气质量等级评估模型。
2.根据权利要求1所述的方法,其特征在于,还包括:
步骤六、对模型进行参数优化并进行性能与适用性测试。
3.根据权利要求1所述的方法,其特征在于,步骤四二中的A=128。
4.根据权利要求1所述的方法,其特征在于,步骤四三中的B=20。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911333157.5A CN111079848B (zh) | 2019-12-23 | 2019-12-23 | 一种基于图像的空气质量等级评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911333157.5A CN111079848B (zh) | 2019-12-23 | 2019-12-23 | 一种基于图像的空气质量等级评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111079848A true CN111079848A (zh) | 2020-04-28 |
CN111079848B CN111079848B (zh) | 2021-02-12 |
Family
ID=70316651
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911333157.5A Expired - Fee Related CN111079848B (zh) | 2019-12-23 | 2019-12-23 | 一种基于图像的空气质量等级评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111079848B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115359431A (zh) * | 2022-10-21 | 2022-11-18 | 山东省生态环境监测中心 | 一种大气环境污染源污染程度评估方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104517126A (zh) * | 2014-12-26 | 2015-04-15 | 北京邮电大学 | 一种基于图像分析的空气质量评估方法 |
CN105445204A (zh) * | 2015-12-07 | 2016-03-30 | 太原理工大学 | 一种基于镜头蓝色光波分析的空气质量等级判别方法 |
US20160370333A1 (en) * | 2015-06-22 | 2016-12-22 | International Business Machines Corporation | Generating fine resolution air pollution estimates |
CN106651036A (zh) * | 2016-12-26 | 2017-05-10 | 东莞理工学院 | 空气质量预报系统 |
CN108052980A (zh) * | 2017-12-21 | 2018-05-18 | 南京理工大学 | 基于图像的空气质量等级检测方法 |
-
2019
- 2019-12-23 CN CN201911333157.5A patent/CN111079848B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104517126A (zh) * | 2014-12-26 | 2015-04-15 | 北京邮电大学 | 一种基于图像分析的空气质量评估方法 |
US20160370333A1 (en) * | 2015-06-22 | 2016-12-22 | International Business Machines Corporation | Generating fine resolution air pollution estimates |
CN105445204A (zh) * | 2015-12-07 | 2016-03-30 | 太原理工大学 | 一种基于镜头蓝色光波分析的空气质量等级判别方法 |
CN106651036A (zh) * | 2016-12-26 | 2017-05-10 | 东莞理工学院 | 空气质量预报系统 |
CN108052980A (zh) * | 2017-12-21 | 2018-05-18 | 南京理工大学 | 基于图像的空气质量等级检测方法 |
Non-Patent Citations (1)
Title |
---|
邓益侬 等: "改进YOLOv3网络在图像中评价空气质量", 《计算机工程与应用》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115359431A (zh) * | 2022-10-21 | 2022-11-18 | 山东省生态环境监测中心 | 一种大气环境污染源污染程度评估方法及系统 |
CN115359431B (zh) * | 2022-10-21 | 2023-02-28 | 山东省生态环境监测中心 | 一种大气环境污染源污染程度评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111079848B (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108052980B (zh) | 基于图像的空气质量等级检测方法 | |
CN103186904B (zh) | 图片轮廓提取方法及装置 | |
CN106557764B (zh) | 一种基于二进制编码字符水尺和图像处理的水位识别方法 | |
CN105787937B (zh) | 一种基于osm的高分辨遥感影像道路变化检测方法 | |
CN103048329B (zh) | 一种基于主动轮廓模型的路面裂缝检测方法 | |
CN103235952B (zh) | 一种基于高分辨率遥感影像的城市绿度空间的度量方法 | |
CN108009591A (zh) | 一种基于深度学习的接触网关键部件识别方法 | |
CN105512666A (zh) | 一种基于视频的河道垃圾识别方法 | |
CN109376740A (zh) | 一种基于视频的水尺读数检测方法 | |
CN107169953A (zh) | 基于hog特征的桥梁混凝土表面裂缝检测方法 | |
CN113313677B (zh) | 一种卷绕锂电池x光图像的质量检测方法 | |
CN102982534B (zh) | 基于弦线切线法的Canny边缘检测双阈值获取方法 | |
CN107480607B (zh) | 一种智能录播系统中站立人脸检测定位的方法 | |
Li et al. | Research on vehicle license plate location based on neural networks | |
CN112149543B (zh) | 一种基于计算机视觉的建筑扬尘识别系统与方法 | |
CN109118548A (zh) | 一种综合智能水质识别方法 | |
CN105354847A (zh) | 基于滑动比较窗口自适应分割的水果表面缺陷检测方法 | |
CN108510496B (zh) | 基于图像dct域的svd分解的模糊检测方法 | |
CN109961070A (zh) | 一种输电线路智能化图像监测中区分雾体浓度的方法 | |
CN104268590A (zh) | 基于互补性组合特征与多相回归的盲图像质量评价方法 | |
CN104700405A (zh) | 一种前景检测方法和系统 | |
CN108665468B (zh) | 一种提取直线塔绝缘子串的装置及方法 | |
CN111079848B (zh) | 一种基于图像的空气质量等级评估方法 | |
CN103295238B (zh) | 安卓平台上基于roi运动检测的视频实时定位方法 | |
CN113469097B (zh) | 一种基于ssd网络的水面漂浮物多相机实时检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210212 Termination date: 20211223 |
|
CF01 | Termination of patent right due to non-payment of annual fee |