CN110222633A - 基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法 - Google Patents
基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法 Download PDFInfo
- Publication number
- CN110222633A CN110222633A CN201910481037.3A CN201910481037A CN110222633A CN 110222633 A CN110222633 A CN 110222633A CN 201910481037 A CN201910481037 A CN 201910481037A CN 110222633 A CN110222633 A CN 110222633A
- Authority
- CN
- China
- Prior art keywords
- image
- color
- moment
- representing
- sliding window
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000008569 process Effects 0.000 title claims abstract description 24
- 239000002910 solid waste Substances 0.000 title claims abstract description 20
- 238000000605 extraction Methods 0.000 claims abstract description 31
- 238000002485 combustion reaction Methods 0.000 claims abstract description 30
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 238000000513 principal component analysis Methods 0.000 claims abstract description 13
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000004056 waste incineration Methods 0.000 claims description 14
- 239000010813 municipal solid waste Substances 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 6
- 238000012706 support-vector machine Methods 0.000 claims description 6
- 238000002834 transmittance Methods 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 5
- 230000001629 suppression Effects 0.000 abstract 1
- 238000009826 distribution Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012843 least square support vector machine Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000004939 coking Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000003344 environmental pollutant Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 231100000719 pollutant Toxicity 0.000 description 2
- 239000000779 smoke Substances 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 239000002956 ash Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000010881 fly ash Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- 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/56—Extraction of image or video features relating to colour
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
- Incineration Of Waste (AREA)
Abstract
本发明涉及基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法。国内城市固体废物(MSW)的组分复杂,领域专家通常依据经验观测火焰视频图像识别焚烧炉燃烧状态进而调整MSW焚烧(MSWI)过程操作参数,难以维持稳定的运行工况。首先,对焚烧火焰图像进行去雾和去噪预处理,提高图像清晰度;接着,将图像转换到适合视觉系统识别的HSV空间后采用滑窗分块提取颜色矩特征,采用主成分分析(PCA)提取潜在特征以消除高维颜色矩特征间的共线性;最后,以提取的相互独立的潜在特征为输入,采用最小二乘支持向量机(LSSVM)算法构建燃烧工况识别模型。基于国内某厂的实际焚烧图像仿真验证了所提方法的有效性。
Description
技术领域
本发明属于城市固体废物焚烧领域。
背景技术
城市固体废物(MSW)的产生数量随工业化、城市化的发展以及人民生活水平的日益提高已达到8%的全球年增长率[1][2]。MSW焚烧(MSWI)发电技术因其减量效果好、处理速度快和能源利用高等优点已经逐渐成为我国解决“垃圾围城”问题的主要方式[3]。
MSWI过程的固废焚烧状态与焚烧装置和蒸汽发电装置的安全运行、全流程的运行优化化控制、污染物生成量及排放量密切相关。目前,发展中国家的MSWI 仍然存在许多亟待解决的问题[4],其中最为突出的问题是污染物排放不达标[5][6],其与燃烧过程的运行状态不稳定密切相关。我国MSW的成分特殊、热值波动大,国外引进成套设备难以自动运行,多由领域专家摸索经验后手动控制运行,易造成炉膛内结焦、积灰、腐蚀等问题,甚至引起炉膛爆炸[7]。因此, MSWI过程的优化运行需要准确识别焚烧状态和及时调整操作参数,以确保焚烧的充分性,提高运行的经济性和环保性。国内某MSWI发电厂的工艺流程如图1 所示。
由图1可知,我国MSWI过程对焚烧状态的识别主要依靠领域专家观测焚烧炉内的火焰视频图像,凭经验识别焚烧状态进而调整操作参数,存在随意性和差异性,难以维持MSWI过程运行在最佳燃烧工况。因此,非常有必要研究基于机器视觉的MSWI燃烧工况识别技术。
目前,识别MSWI过程的燃烧工况主要依据固废焚烧图像中火焰燃烬线分布的不同位置进行,机器视觉在MSWI火焰状态识别中的研究还不够成熟。针对复杂多变的MSWI环境,如何构造可靠、鲁棒的焚烧工况识别模型还是一个开放性难题。MSWI过程工况识别存在诸多难点,如焚烧过程中炉膛内固废灰烬的影响、图像传输中的噪声干扰等因素导致火焰边缘模糊,MSW成分及其热值的大范围波动性、燃烧区域划分的特定性等因素需要依据焚烧对象特性提取特征。此外,如何构建分类器模型对非典型工况进行准确识别也是需要深入的问题。
针对焚烧图像预处理,文献[8]提出的基于暗通道先验的去雾算法有效提高了图像能见度;文献[9]通过对比实验,验证了中值滤波算法能够有效抑制图像传输噪声;文献[10]指出,中值滤波算法能够有效去除脉冲噪声和保留图像边缘;但上述预处理算法在MSWI过程火焰图像中的应用还未见报道。
针对焚烧图像特征提取,文献[11]依据据火焰的颜色特点和空间分布特点,提取火焰平均亮度、火焰图像锋面中心点位置、火焰锋面轴向位移和火焰颜色分量比等4个特征;文献[12]提取焚烧火焰图像有效区平均灰度及其方差、火焰高温区面积、着火面积、火焰中心水平偏移距离等7个特征;文献[13]提取了感兴趣区域的色彩特征、全局形态特征和局部形态特征。上述研究所提取的特征重点在于火焰的大小、位置和亮度,其提取过程较为繁琐。通常,MSWI过程的焚烧图像的直观特征是颜色。目前常用的颜色特征有颜色直方图、颜色不变量、颜色矩、颜色熵、颜色聚合向量和颜色相关图,其中颜色矩具有简单有效、处理效率高等特点。针对MSWI火焰图像在不同燃烧区域的特点,提取不同区域颜色矩特征的研究还未见报道。针对焚烧图像的高维特征,文献[12]运用粗糙集理论进行维数约简,但该方法主要是进行特征的删减,未能有效去除特征间的多重共线性。文献[14]采用主成分分析法(PCA)处理高维图像特征获得相互独立的潜在变量。
针对焚烧工况识别,文献[12]采用BP神经网络模型识别MSWI燃烧状态,但方法要求数据之间存在明显差异性;文献[15]的研究表明,与其他方法相比,支持向量机(SVM)和最小二乘支持向量机(LSSVM)具有较高的分类准确率,能够很好的解决小样本、非线性问题。相比SVM,LSSVM具有较快的求解速度快,求解所需计算资源少。
综上,本文提出基于火焰图像颜色特征的MSWI过程燃烧工况识别方法。首先,对焚烧火焰图像进行去雾和去噪预处理,提高图像清晰度;接着,将图像转换到适合视觉系统识别的HSV空间后采用滑窗分块提取颜色矩特征,采用 PCA提取潜在特征以消除高维颜色矩特征间的共线性;最后,以提取的相互独立的潜在特征为输入,采用LSSVM算法构建燃烧工况识别模型。基于国内某厂的实际焚烧图像验证了所提方法的有效性。
发明内容
MSW由专用运输车收集后运至卸料大厅,倾倒至密封的存放池内后由人工操控的吊斗将MSW放入焚烧炉进料斗内,再由给料机将其推至炉排炉。在焚烧炉的炉排内依次经历干燥、点燃、燃烧和烧尽四个阶段。相应的燃烧工况可以分为:(1)燃烬线前移;(2)燃烬线正常;(3)燃烬线后移。通过图像中火焰燃烬线的分布位置能够识别工况。不同工况的描述和相应操作为:(1)燃烬线前移,即靠近推料器下方,易出现结焦阻碍下料口等情况,此时需要增加上炉排或者推料器速度,降低一次风量;(2)燃烬线正常,火焰集中在燃烧线附近,火焰明亮,垃圾充分燃烧,燃烧状态良好;(3)燃烬线后移,固废在燃烬区焚烧,造成垃圾焚烧不充分,需要减慢燃烬炉排推进速度,增加一次风量。当火焰燃烬线正常时是最佳燃烧状态,符合工业生产要求。
根据上述描述可知,当前主要依靠人工方式观察焚烧炉内固废焚烧视频,通过观察火焰燃烬线的分布位置判断工况,但是该方式易受操作者经验,状态和工作态度的影响。
随着图像处理技术的不断发展和成熟,本文将机器视觉技术应用到固废焚烧中,提出了由图像预处理、基于滑窗分块的颜色矩特征提取和燃烧工况识别3 个模块组成的MSWI过程燃烧工况识别方法,如图2所示。
图2中,表示原始图像集合,In(s,t)表示第n幅图像,N表示全部图像数量,即建模样本数量,(s,t)表示像素点的空间坐标;表示经过中值滤波后得到的图像。表示经预处理之后得到的全部火焰图像; Xcolor∈RN×P表示提取颜色矩后特征集合,P表示颜色矩特征的维数;Zcolor∈RN×M表示对颜色矩特征进行潜在变量提取后的特征集合,M表示潜在变量的维数;表示燃烧工况的识别结果。
上述模块的功能是:
(1)图像预处理模块:对图像进行去雾和去噪处理。
(2)基于滑窗分块的颜色矩特征提取模块:对图像进行颜色空间转换,利用滑窗分块并提取图像的颜色矩特征,将不同滑块的颜色矩特征进行串行融合后进行潜在特征提取和选择。
(3)燃烧工况识别模块:基于选择的潜在特征构建燃烧工况识别模型。
固废焚烧图像的噪声主要来源于采集和传输过程,其造成图像质量差进而影响后续燃烧工况的识别。在实际焚烧图像采集过程中,固废焚烧内部环境复杂,如较大的光线强弱变化以及焚烧过程中伴有的飞灰和烟雾,会导致火焰图像的清晰度变差,进而在图像生成和传输过程中增加了图像噪声的随机性。
采用暗通道先验去雾算法还原了固废焚烧火焰图像的颜色和能见度。以固废焚烧原始图像数据集中的第n张图片In(s,t)为例,其为Len*Wid*3大小的 RGB图像,各颜色空间下像素点总数NP=Len*Wid,其中Len和Wid分别表示图像的长度和宽度。基于暗通道先验的图像去雾算法可描述如下。
首先,定义局部块Ω=k×l,以Ω大小的模板对有雾图像In(s,t)在R、G、B 三个通道里做最小值滤波得到
接着,计算大气光成分A,步骤如下首先在暗通道图中选择图像像素点个数 (Np)千分之一个像素值最高的点Np_num,如下所示:
并记录这些点对应的坐标(s,t);再根据这些点的坐标(s,t)分别在原图像In的三个通道里找到对应像素点求和分别得到r_sum、g_sum和b_sum,统一表示为(r_sum,g_sum,b_sum);最后得到大气光成分A=[Ar Ag Ab],Ar、Ag和Ab分别表示R、G、B三个通道下对应A值的大小,即Ar=r_sum/NP;Ag=g_sum/NP; Ab=b_sum/NP。
然后,基于以下公式计算透射率
其中λ是深度系数,该参数的设定可以为远处的景物保留一定量的雾,使图像更为自然,一般视实际应用而定,在本文中λ取值是0.95。
利用如下公式得到去雾后的图像Jn(s,t),
当透射率接近0时,直接衰减项也接近0,这样会导致复原图像产生噪声,考虑到原始焚烧图像的暗度,需对设定一个阈值 r0。r0的典型值是0.1。则式(3)改写为:
其中,max(·)表示取最大值函数。
采用非线性中值滤波算法对图像进行去噪,能够在图像去噪的同时保留火焰边缘。假定S表示大小为mwin*nwin的矩形窗口,中值滤波器是指mwin*nwin的窗口中将像素值从小到大排序,将中间像素值赋给对应模板中心的像素,即
其中,表示去噪声后的第n幅图像。
基于滑窗分块的颜色矩特征提取包括颜色空间转换、图像滑窗分块、颜色矩特征提取、特征串行组合和潜在特征提取等共5个模块,分别描述如下。
(1)颜色空间转换
由于采集到的固废焚烧图像具有明显的亮度变化,不易与背景进行区分,故将图像从RGB颜色空间转换到符合人类视觉系统的HSV颜色空间,过程如下所示,
其中,表示转到HSV颜色空间的第n幅图像。
(2)图像滑窗分块
以为例,采用滑动窗口进行图像分块特征提取。采用无重叠像素的滑窗,令G=msli*nsli为滑窗大小,其中G为单个滑窗包含的像素点数,msli和nsli分别表示滑窗的长和宽,令d2表示每个颜色空间下进行滑窗的次数,其计算如下:
其中滑窗大小与图像大小成比例。
(3)颜色矩特征提取
选择3个颜色空间的一阶矩、二阶矩和三阶矩作为图像特征。以第n幅图像的第q次滑窗块为例,各颜色空间下的一阶矩二阶矩和三阶矩如下所示:
其中,和分别表示H、S和V三个颜色空间下的一阶矩;和分别表示H、S和V三个颜色空间下的二阶矩;和分别表示H、S和V三个颜色空间下的三阶矩;表示大小为G的滑窗内第g个像素点值。
(4)特征串行组合
全部样本的颜色矩特征Xcolor表示为,
其中,和分别表示第n张图片的第q个滑窗对应的一阶矩、二阶矩和三阶矩。由上式可知,每幅图像提取的特征维数为P=d2×3×3,故 Xcolor∈RN×P。
(5)潜在特征提取
采用PCA对潜在特征进行提取,将颜色矩特征从高维空间映射到低维空间。此处将累计方差贡献率的阈值记为θPCA,该阈值用于确定主成分的个数,即获得潜在特征变量的数量。一般实际应用中,要求θPCA的取值大于等于85%,因此本文选取θPCA为90%。所提取与选择的潜在特征数量M,则获得的输入特征数据集表示为,
其中,zn表示第n个样本,zm表示第m个变量。
构建基于LS-SVM的燃烧工况识别模型。
LS-SVM的优化问题描述为:
其中,w表示权重,c表示正则化项,ξn表示误差变量,保证一定的容错率, yn表示第n个样本对应的类别,Ψ(zn)表示非线性映射,将输入数据从低维映射到高维,b表示偏差,n=1,2,…,N,N表示全部图像数量。
用拉格朗日法求解上述优化问题,构造如下方程为:
式中αn为拉格朗日乘子,将式(17)分别对w,b,ξn和αn求导,并令等式为零可得到:
LS_SVM优化问题转化为求解线性方程:
其中K(zi,zj)为支持向量机的核函数,i=1,2,...,N,j=1,2,...,N。本文中,核函数选择径向基函数:
其中zi和zj表示图像样本,σ表示核函数的宽度。
则LS-SVM的分类决策函数表达式为:
其中sgn(·)表示符号函数,表示图像的分类结果。
附图说明
图1国内某MSWI发电厂工艺流程
图2基于火焰图像颜色特征的MSWI过程燃烧工况识别
图3固废焚烧状态1
图4固废焚烧状态2
图5固废焚烧状态3
图6工况1的去雾图像
图7工况2的去雾图像
图8工况3的去雾图像
图9工况1的去噪图像
图10工况2的去噪图像
图11工况3的去噪图像
图12工况1的颜色空间转换结果
图13工况2的颜色空间转换结果
图14工况3的颜色空间转换结果
图15滑窗块操作后的工况1
图16滑窗块操作后的工况2
图17滑窗块操作后的工况3
图18不同工况下特征串行后的分布
图19 PCA提取的潜在特征选择曲线
具体实施方式
本文中的建模数据参考于北京某MSWI焚烧企业。以间隔一分钟的速度对焚烧火焰视频进行截图。样本数量N=270,其中训练样本180个,测试样本90 个,图片大小为1436*507。炉膛内左右两侧均安装有摄像机,采集到的图像分为左右两部分。结合现场专家经验,将焚烧工况分为三种,火焰分别集中位于燃一段、燃烧线和燃烬段。三种工况分别如图3-图5所示。
首先采用暗通道先验去雾算法对图像进行去雾,参数设定值:局部块Ω=9*9, 阈值r0=0.1,λ=0.95。不同工况下的处理结果如图6-图8所示。
由图6-8可知,经去雾处理后,焚烧图片中的烟雾明显减少,火焰颜色明亮且易于与背景分离,但是仍然存在随机噪声,需进一步对图像进行去噪。
采用S=5*5模板对图像进行去噪处理。不同工况下的处理结果如图9-图11 所示。
由图9-11可知,图像经中值滤波去噪处理后,噪点明显减少,提高图像质量的同时保留了边缘细节。
(1)颜色空间转换结果
三种工况下对图像进行颜色空间转换结果。如图12-图14所示。
根据上述图像可知,与背景相比,火焰的亮度和颜色较为突出,因此火焰图像由RGB转换到HSV颜色空间之后,其明亮程度分布就能很好的体现出来,在不同工况下,分布也不相同。
(2)滑窗分块结果
本次实验中采用的滑窗大小为287*101,将图像分为25块,如图15-图17 所示。
由图15-图17可知,经过滑窗块操作得到不同工况下差别化的焚烧区域,其滑窗块分别位于焚烧火焰图像的燃一段,燃烧线和燃烬段。
利用上述滑窗块提取的颜色矩是表征焚烧工况的关键特征。将所有滑窗块提取到的特征串行组合后得到225维特征,不同工况下这些特征的分布如图18 所示。
从图18中可以看出,采用225维的数据集在不同工况下的分布差异不明显,不利于分类器进行识别,因此需要对颜色矩特征进行潜在变量提取。
(3)潜在特征提取结果
利用PCA进行潜在变量提取结果如表1所示。
表1不同维数下的累计方差贡献率
PCA潜在特征选择曲线如图19所示。
结合表1和图19,本文选择阈值θPCA=0.9,最终确定的输入样本的维数大小为27。
基于不同颜色特征的分类器模型比较
本实验中采LSSVMlabv1_8_R2009b_R2011a工具箱进行仿真,对三种特征提取方法分别运行5次得到的结果如表2所示。
表2基于不同颜色特征的分类器模型比较结果
由表可知:第1种方法不进行滑窗操作,对9维的数据集仿真得到平均准确率是50.44%,最高准确率是52.22%,对应的c和σ分别取1.55和4.28;第2种方法引入滑窗操作,输入数据为225维,仿真结果平均准确率是64%,最高准确率是68.89%;对应的c和σ分别取6.58和60.17;第3种方法对225 数据进行潜在特征提取,得到27维的数据,仿真结果平均准确率是71.67%,最高准确率是75%,对应的c和σ分别取17.34和13.24。
实验结果表明,第3种方法进行滑窗操作和潜在变量提取之后得到的数据无论是平均准确率还是最高准确率都要优于前两种方法,所提方法具有有效性。
本文提出了基于火焰图像颜色特征提取的MSWI过程燃烧工况识别方法。该方法的贡献表现在:(1)基于暗通道的图像去雾算法提高了焚烧火焰图像的清晰度;(2)基于滑窗分块在HSV空间提取高维颜色矩特征和基于PCA 再提取潜在特征消除高维颜色矩特征间的共线性的策略能够提取有效特征和降低特征维数。基于北京某固废焚烧发电企业的工业过程数据,仿真验证了所提方法的有效性和可行性。
。
Claims (1)
1.基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法,其特征在于:
首先定义如下:表示原始图像集合,In(s,t)表示第nth幅图像,N表示全部图像数量,即建模样本数量,(s,t)表示像素点的空间坐标;表示经过中值滤波后得到的图像;表示经预处理之后得到的全部火焰图像;Xcolor∈RN×P表示提取颜色矩后特征集合,P表示颜色矩特征的维数;Zcolor∈RN×M表示对颜色矩特征进行潜在变量提取后的特征集合,M表示潜在变量的维数;表示燃烧工况的识别结果;
包括以下模块:
(1)图像预处理模块:对图像进行去雾和去噪处理;
(2)基于滑窗分块的颜色矩特征提取模块:对图像进行颜色空间转换,利用滑窗分块并提取图像的颜色矩特征,将不同滑块的颜色矩特征进行串行融合后进行潜在特征提取和选择;
(3)燃烧工况识别模块:基于选择的潜在特征构建燃烧工况识别模型;
采用暗通道先验去雾算法还原了固废焚烧火焰图像的颜色和能见度;以固废焚烧原始图像数据集中的第n张图片In(s,t)为例,其为Len*Wid*3大小的RGB图像,各颜色空间下像素点总数NP=Len*Wid,其中Len和Wid分别表示图像的长度和宽度;
基于暗通道先验的图像去雾算法描述如下;
首先,定义局部块Ω=k×l,以Ω大小的模板对有雾图像In(s,t)在R、G、B三个通道里做最小值滤波得到
接着,计算大气光成分A,步骤如下:首先在暗通道图中选择图像像素点个数(Np)千分之一个像素值最高的点Np_num,如下所示:
并记录这些点对应的坐标(s,t);再根据这些点的坐标(s,t)分别在原图像In的三个通道里找到对应像素点求和分别得到r_sum、g_sum和b_sum,统一表示为(r_sum,g_sum,b_sum);最后,得到大气光成分A=[Ar Ag Ab],Ar、Ag和Ab分别表示R、G、B三个通道下对应A值的大小,即Ar=r_sum/NP;Ag=g_sum/NP;Ab=b_sum/NP;
然后,基于以下公式计算透射率
其中λ是深度系数,该参数的设定可以为远处的景物保留一定量的雾,使图像更为自然,一般视实际应用而定,在本文中λ取值0.95;,;
利用如下公式得到去雾后的图像Jn(s,t),
当透射率接近0时,直接衰减项也接近0,这样会导致复原图像产生噪声,考虑到原始焚烧图像的暗度,需对设定一个阈值r0,r0的典型值是0.1;则式(3)改写为:
其中,max(·)表示取最大值函数;
采用非线性中值滤波算法对图像进行去噪,能够在图像去噪的同时保留火焰边缘;假定S表示大小为mwin*nwin的矩形窗口,中值滤波器是指mwin*nwin的窗口中将像素值从小到大排序,将中间像素值赋给对应模板中心的像素,即
其中,表示去噪声后的第nth幅图像;
基于滑窗分块的颜色矩特征提取包括颜色空间转换、图像滑窗分块、颜色矩特征提取、特征串行组合和潜在特征提取等共5个模块,分别描述如下;
(1)颜色空间转换
将图像从RGB颜色空间转换到符合人类视觉系统的HSV颜色空间,过程如下所示,
其中,表示转到HSV颜色空间的第nth幅图像;
(2)图像滑窗分块
以为例,采用滑动窗口进行图像分块特征提取;采用无重叠像素的滑窗,令G=msli*nsli为滑窗大小,其中G为单个滑窗包含的像素点数,msli和nsli分别表示滑窗的长和宽,令d2表示每个颜色空间下进行滑窗的次数,其计算如下:
其中滑窗大小与图像大小成比例;
(3)颜色矩特征提取
选择3个颜色空间的一阶矩、二阶矩和三阶矩作为图像特征;以第n幅图像的第q次滑窗块为例,各颜色空间下的一阶矩二阶矩和三阶矩如下所示:
其中,和分别表示H、S和V三个颜色空间下的一阶矩;和分别表示H、S和V三个颜色空间下的二阶矩; 和分别表示H、S和V三个颜色空间下的三阶矩;表示大小为G的滑窗内第g个像素点值;
(4)特征串行组合
全部样本的颜色矩特征Xcolor表示为,
其中,和分别表示第n张图片的第q个滑窗对应的一阶矩、二阶矩和三阶矩;由上式可知,每幅图像提取的特征维数为P=d2×3×3,故Xcolor∈RN×P;
(5)潜在特征提取;
采用PCA对进行潜在特征提取,将颜色矩特征从高维空间映射到低维空间;此处将累计方差贡献率的阈值记为θPCA,该阈值用于确定主成分的个数,即获得潜在特征变量的数量;一般实际应用中,要求θPCA的取值大于等于85%,因此本文中选取θPCA为90%;所提取与选择的潜在特征数量M,则获得的输入特征数据集表示为,
其中,zn表示第n个样本,zm表示第m个变量;
构建基于LS-SVM的燃烧工况识别模型;
LS-SVM的优化问题描述为:
其中,w表示权重,c表示正则化项,ξn表示误差变量,-,yn表示第n个样本对应的类别,Ψ(zn)表示非线性映射,将输入数据从低维映射到高维空间,b表示偏差,n=1,2,…,N,N表示全部图像数量;用拉格朗日法求解上述优化问题,构造如下方程为:
式中αn为拉格朗日乘子,将式(17)分别对w,b,ξn和αn求导,并令等式为零得到:
LS_SVM优化问题转化为求解线性方程:
其中K(zi,zj)为支持向量机的核函数,i=1,2,...,N,j=1,2,...,N;本文中,核函数选择径向基函数:
其中zi和zj表示图像样本,σ表示核函数的宽度;
则LS-SVM的分类决策函数表达式为:
其中sgn(·)表示符号函数,表示图像的分类结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910481037.3A CN110222633B (zh) | 2019-06-04 | 2019-06-04 | 基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910481037.3A CN110222633B (zh) | 2019-06-04 | 2019-06-04 | 基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110222633A true CN110222633A (zh) | 2019-09-10 |
CN110222633B CN110222633B (zh) | 2021-04-27 |
Family
ID=67819233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910481037.3A Active CN110222633B (zh) | 2019-06-04 | 2019-06-04 | 基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110222633B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110566929A (zh) * | 2019-10-17 | 2019-12-13 | 山东鸿源工业智能技术有限公司 | 一种双炉排工业锅炉实时控制方法 |
CN110717495A (zh) * | 2019-09-30 | 2020-01-21 | 北京工业大学 | 基于多尺度颜色矩特征和随机森林的固废焚烧工况识别方法 |
CN111461003A (zh) * | 2020-03-31 | 2020-07-28 | 湖南大学 | 基于视频图像序列特征提取的燃煤工况识别方法 |
CN112016552B (zh) * | 2020-11-02 | 2021-02-12 | 矿冶科技集团有限公司 | 基于泡沫颜色的混合浮选工况识别方法和系统 |
CN113344909A (zh) * | 2021-07-01 | 2021-09-03 | 中国石油大学(北京) | 一种火电锅炉透火焰高温滤镜结焦识别显示方法和装置 |
CN114155383A (zh) * | 2021-09-28 | 2022-03-08 | 北京工业大学 | 基于燃烧火焰颜色矩特征的城市固废焚烧过程风量设定方法 |
CN116089825A (zh) * | 2023-04-07 | 2023-05-09 | 中国环境科学研究院 | 一种基于统计熵的固废指纹特征提取方法 |
CN117173141A (zh) * | 2023-09-11 | 2023-12-05 | 山东博昂信息科技有限公司 | 一种基于火焰图像特征的冶炼观测系统 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2832279A1 (en) * | 2011-04-05 | 2012-10-11 | Mascoma Corporation | Methods for the improvement of product yield and production in a microorganism through the addition of alternate electron acceptors |
CN103092164A (zh) * | 2013-01-04 | 2013-05-08 | 清华大学 | 基于物联网的城市固废收运智能监控装置及方法 |
CN103617414A (zh) * | 2013-11-09 | 2014-03-05 | 中国科学技术大学 | 一种基于最大间距准则的火灾颜色模型及火灾火焰和烟雾识别方法 |
CN105183935A (zh) * | 2015-07-20 | 2015-12-23 | 昆明理工大学 | 一种火焰燃烧状况及稳定性的评价方法 |
CN105740866A (zh) * | 2016-01-22 | 2016-07-06 | 合肥工业大学 | 一种具有仿反馈调节机制的回转窑烧成状态识别方法 |
WO2017018796A1 (ko) * | 2015-07-27 | 2017-02-02 | 부산대학교 산학협력단 | 바이러스 기반 변색 센서에 의한 농작물의 원산지 판별 방법 및 바이러스 기반 변색 센서를 포함하는 농작물의 원산지 판별 장치 |
CN109086589A (zh) * | 2018-08-02 | 2018-12-25 | 东北大学 | 一种结合手势识别的智能终端人脸解锁方法 |
CN109300110A (zh) * | 2018-08-24 | 2019-02-01 | 哈尔滨理工大学 | 一种基于改进颜色模型的森林火灾图像检测方法 |
CN109442474A (zh) * | 2018-11-12 | 2019-03-08 | 西安艾贝尔科技发展有限公司 | 一种气化炉火焰检测装置及检测方法 |
CN109711345A (zh) * | 2018-12-27 | 2019-05-03 | 南京林业大学 | 一种火焰图像识别方法、装置及其存储介质 |
CN109754019A (zh) * | 2019-01-10 | 2019-05-14 | 燕山大学 | 一种连续监测锅炉燃烧状况的方法 |
-
2019
- 2019-06-04 CN CN201910481037.3A patent/CN110222633B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2832279A1 (en) * | 2011-04-05 | 2012-10-11 | Mascoma Corporation | Methods for the improvement of product yield and production in a microorganism through the addition of alternate electron acceptors |
CN103092164A (zh) * | 2013-01-04 | 2013-05-08 | 清华大学 | 基于物联网的城市固废收运智能监控装置及方法 |
CN103617414A (zh) * | 2013-11-09 | 2014-03-05 | 中国科学技术大学 | 一种基于最大间距准则的火灾颜色模型及火灾火焰和烟雾识别方法 |
CN105183935A (zh) * | 2015-07-20 | 2015-12-23 | 昆明理工大学 | 一种火焰燃烧状况及稳定性的评价方法 |
WO2017018796A1 (ko) * | 2015-07-27 | 2017-02-02 | 부산대학교 산학협력단 | 바이러스 기반 변색 센서에 의한 농작물의 원산지 판별 방법 및 바이러스 기반 변색 센서를 포함하는 농작물의 원산지 판별 장치 |
CN105740866A (zh) * | 2016-01-22 | 2016-07-06 | 合肥工业大学 | 一种具有仿反馈调节机制的回转窑烧成状态识别方法 |
CN109086589A (zh) * | 2018-08-02 | 2018-12-25 | 东北大学 | 一种结合手势识别的智能终端人脸解锁方法 |
CN109300110A (zh) * | 2018-08-24 | 2019-02-01 | 哈尔滨理工大学 | 一种基于改进颜色模型的森林火灾图像检测方法 |
CN109442474A (zh) * | 2018-11-12 | 2019-03-08 | 西安艾贝尔科技发展有限公司 | 一种气化炉火焰检测装置及检测方法 |
CN109711345A (zh) * | 2018-12-27 | 2019-05-03 | 南京林业大学 | 一种火焰图像识别方法、装置及其存储介质 |
CN109754019A (zh) * | 2019-01-10 | 2019-05-14 | 燕山大学 | 一种连续监测锅炉燃烧状况的方法 |
Non-Patent Citations (6)
Title |
---|
MONGKOL PATCHARAVONGSIRI 等: ""Effect of rotation rate on rotary bio-drying of dewatered municipal solid waste"", 《ICMCE》 * |
XIUCHENG YANG 等: ""Mapping of urban water bodies from sentinel-2 MSI imageryat 10m resolution via NDWI-based image sharpening"", 《REMOTE SENSING》 * |
刘小虎 等: ""基于改进梯度边缘特征的早期火焰识别方法"", 《消防科学与技术》 * |
崔宁宁: ""炉内火焰图像特征提取关键技术研究"", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
徐欢: ""基于模糊聚类的炉膛火焰分割技术研究"", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
白晓静: ""基于数字成像的单颗粒燃料燃烧特性研究与炉膛火焰监测"", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110717495A (zh) * | 2019-09-30 | 2020-01-21 | 北京工业大学 | 基于多尺度颜色矩特征和随机森林的固废焚烧工况识别方法 |
CN110717495B (zh) * | 2019-09-30 | 2024-01-26 | 北京工业大学 | 基于多尺度颜色矩特征和随机森林的固废焚烧工况识别方法 |
CN110566929A (zh) * | 2019-10-17 | 2019-12-13 | 山东鸿源工业智能技术有限公司 | 一种双炉排工业锅炉实时控制方法 |
CN110566929B (zh) * | 2019-10-17 | 2021-03-26 | 山东鸿源工业智能技术有限公司 | 一种双炉排工业锅炉实时控制方法 |
CN111461003A (zh) * | 2020-03-31 | 2020-07-28 | 湖南大学 | 基于视频图像序列特征提取的燃煤工况识别方法 |
CN112016552B (zh) * | 2020-11-02 | 2021-02-12 | 矿冶科技集团有限公司 | 基于泡沫颜色的混合浮选工况识别方法和系统 |
CN113344909A (zh) * | 2021-07-01 | 2021-09-03 | 中国石油大学(北京) | 一种火电锅炉透火焰高温滤镜结焦识别显示方法和装置 |
CN113344909B (zh) * | 2021-07-01 | 2023-12-08 | 中国石油大学(北京) | 一种火电锅炉透火焰高温滤镜结焦识别显示方法和装置 |
CN114155383A (zh) * | 2021-09-28 | 2022-03-08 | 北京工业大学 | 基于燃烧火焰颜色矩特征的城市固废焚烧过程风量设定方法 |
CN114155383B (zh) * | 2021-09-28 | 2024-05-31 | 北京工业大学 | 基于燃烧火焰颜色矩特征的城市固废焚烧过程风量设定方法 |
CN116089825A (zh) * | 2023-04-07 | 2023-05-09 | 中国环境科学研究院 | 一种基于统计熵的固废指纹特征提取方法 |
CN117173141A (zh) * | 2023-09-11 | 2023-12-05 | 山东博昂信息科技有限公司 | 一种基于火焰图像特征的冶炼观测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110222633B (zh) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110222633B (zh) | 基于火焰图像颜色特征提取的城市固废焚烧过程燃烧工况识别方法 | |
CN110717495B (zh) | 基于多尺度颜色矩特征和随机森林的固废焚烧工况识别方法 | |
CN109101878B (zh) | 一种用于秸秆燃值估计的图像分析系统及图像分析方法 | |
CN103778418A (zh) | 一种输电线路杆塔图像监测系统的山火图像识别方法 | |
CN111814750A (zh) | 基于深度学习目标检测和图像识别的智能垃圾分类方法及系统 | |
Tang et al. | Prediction of the ash content of flotation concentrate based on froth image processing and BP neural network modeling | |
CN116543241B (zh) | 泄露气体云团的检测方法、装置、存储介质和电子设备 | |
CN106289070A (zh) | 测量不规则物体长宽的方法 | |
CN107644431A (zh) | 一种低照度情况下视频监控的有效浓缩方法 | |
Subramanian et al. | Learning forest wildfire dynamics from satellite images using reinforcement learning | |
Zhao et al. | Learning class-specific spectral patterns to improve deep learning-based scene-level fire smoke detection from multi-spectral satellite imagery | |
CN116524397A (zh) | 一种管道缺陷智能检测分析方法 | |
CN114973000B (zh) | 一种垃圾焚烧炉的燃烧优化控制方法 | |
CN116385368A (zh) | 一种基于生成对抗网络的光伏电池缺陷检测数据集增广方法 | |
CN114155383B (zh) | 基于燃烧火焰颜色矩特征的城市固废焚烧过程风量设定方法 | |
CN115761459A (zh) | 一种面向桥隧表观病害识别的多场景自适应方法 | |
CN113902680B (zh) | 引入注意力机制的多尺度特征融合ssd的螺栓松动检测方法 | |
Duan et al. | Recognition of combustion condition in MSWI process based on multi-scale color moment features and random forest | |
Jiang et al. | Deep learning of qinling forest fire anomaly detection based on genetic algorithm optimization | |
Zhang | Computer image processing and neural network technology for thermal energy diagnosis of boiler plants | |
Gao et al. | A framework to identify guano on photovoltaic modules in offshore floating photovoltaic power plants | |
CN102622760A (zh) | 一种新图像分割算法 | |
Zhang et al. | Recognition of Combustion Conditions in MSWI Process Using Convolutional Neural Network | |
Hui et al. | Discussion and analyze on image fusion technology | |
CN114518737B (zh) | 面向mswi过程多模态数据驱动预测模型的仿实时验证平台 |
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 |