一种基于图像结构相似度的烟雾浓度测量方法和系统与流程

本发明属于烟雾浓度测量领域，更具体地，涉及一种基于图像结构相似度的烟雾浓度测量方法和系统。
背景技术：
：现有技术基于图像结构相似度来检测烟雾，相对于传统的烟雾检测方法具有更强的抗干扰能力，但是它们只能进行视频烟雾定性检测，无法完成视频烟雾浓度的定量测量。现有的烟雾浓度视频测量方法如“基于烟雾浓度的火灾视频检测方法”基于视频图像灰度计算烟雾浓。这种方法的主要依据为浅白色烟雾的存在使背景图像的灰度值增加。实际计算中，利用烟雾化方程，通过灰度不同的相邻背景区域间的烟雾浓度与穿透率近似相等建立方程，求解大气光与烟雾浓度。该方法相比于传统的烟雾浓度测量方法具有结构简单、智能高效等优势，但仍存在不足之处：首先实际生活中，由于烟雾产生的原因不同，生成的烟雾颜色也不同，例如黑烟、白烟、蓝烟等。而白烟以外的其他烟雾是不会使背景图像的灰度值增加的，所以该方法对黑烟、蓝烟等并无很好地检测效果。其次，由于图像的最大灰度值为255，对于灰度值已经很高的背景图像，烟雾的存在使背景图像灰度值提高的空间有限，因此该方法利用灰度值较大的背景图像来计算烟雾浓度是很不准确的。由此可见，现有烟雾浓度测量技术受到烟雾颜色和背景图像灰度的限制，且测量的准确率较低。技术实现要素：针对现有技术的以上缺陷或改进需求，本发明提供了一种基于图像结构相似度的烟雾浓度测量方法和系统，其目的在于提供一种适用于任意背景图像、任意烟雾颜色的方法，由此解决现有烟雾浓度测量技术受到烟雾颜色和背景图像灰度的限制，且测量的准确率较低的技术问题。为实现上述目的，按照本发明的一个方面，提供了一种基于图像结构相似度的烟雾浓度测量方法，包括：(1)采集目标场景的烟雾图像和背景图像，所述背景图像为目标场景未产生烟雾时的图像，计算烟雾图像与背景图像之间的图像结构相似度；(2)提取背景图像的图像特征，利用图像特征描述图像复杂度；利用图像复杂度和图像结构相似度得到烟雾图像的烟雾浓度。进一步的，烟雾浓度m的公式为：m＝-c×(ssim-ssim0)，其中，c表示图像复杂度，ssim表示图像结构相似度，ssim0为预设值。进一步的，预设值为1。进一步的，图像特征包括：角二阶矩asm、对比度con、逆差分矩idm、熵ent和相关性cor。进一步的，图像复杂度c为：c＝con+ent-asm-idm-cor。按照本发明的另一方面，提供了一种基于图像结构相似度的烟雾浓度测量系统，包括：第一模块，用于采集目标场景的烟雾图像和背景图像，所述背景图像为目标场景未产生烟雾时的图像，计算烟雾图像与背景图像之间的图像结构相似度；第二模块，用于提取背景图像的图像特征，利用图像特征描述图像复杂度；利用图像复杂度和图像结构相似度得到烟雾图像的烟雾浓度。进一步的，烟雾浓度m的公式为：m＝-c×(ssim-ssim0)，其中，c表示图像复杂度，ssim表示图像结构相似度，ssim0为预设值。进一步的，预设值为1。进一步的，图像特征包括：角二阶矩asm、对比度con、逆差分矩idm、熵ent和相关性cor。进一步的，图像复杂度c为：c＝con+ent-asm-idm-cor。总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：(1)本发明利用图像复杂度和图像结构相似度得到烟雾图像的烟雾浓度，不受到烟雾颜色和背景图像灰度的限制，且测量的准确率高，应用范围广。(2)优选的，在计算烟雾图像的烟雾浓度时预设值为1，在保证测量结果准确率高的同时，提高了计算效率。(3)优选的，本发明选用角二阶矩asm、对比度con、逆差分矩idm、熵ent和相关性cor，这利用5种图像特征的加权来描述图像复杂度，准确率高。附图说明图1是本发明实施例提供的一种基于图像结构相似度的烟雾浓度测量方法的流程图；图2是本发明实施例的背景图像；图3是本发明实施例中烟密计测量值与结构相似度值的曲线关系图；图4是本发明实施例中烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图5是本发明实施例中三条曲线关系图；图6(a)是本发明实施例中的第一幅梅花图像；图6(b)是本发明实施例中的第二幅梅花图像；图6(c)是本发明实施例中的第三幅梅花图像；图6(d)是本发明实施例中的第四幅梅花图像；图7(a)～7(i)是本发明实施例中的9幅不同图像复杂度的背景图像；图8(a)是本发明实施例烟雾图像中3根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(b)是本发明实施例烟雾图像中4根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(c)是本发明实施例烟雾图像中5根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(d)是本发明实施例烟雾图像中6根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(e)是本发明实施例烟雾图像中7根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(f)是本发明实施例烟雾图像中8根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(a)是本发明实施例中第一背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(b)是本发明实施例中第二背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(c)是本发明实施例中第三背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(d)是本发明实施例中第四背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(e)是本发明实施例中第五背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(f)是本发明实施例中第六背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图。具体实施方式为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。如图1所示，一种基于图像结构相似度的烟雾浓度测量方法，包括：(1)采集目标场景的烟雾图像和背景图像，所述背景图像为目标场景未产生烟雾时的图像，计算烟雾图像与背景图像之间的图像结构相似度；(2)提取背景图像的图像特征，利用图像特征描述图像复杂度；利用图像复杂度和图像结构相似度得到烟雾图像的烟雾浓度。由于图像中烟雾的存在会导致图像失真，图像结构相似度下降，且在一定范围内，烟雾浓度越大，图像失真越严重，图像结构相似度下降越多，因此图像结构相似度与烟雾浓度间存在一定的定量关系。为了寻找这种定量关系，我们录制一段视频，背景图像如图2所示，背景图像表示烟雾图像未产生烟雾前的图像。实验装置有摄像头、实验箱、风扇、烟密计等。其中，摄像头位置固定在实验箱左侧的小孔中，用于录制烟雾浓度随时间变化的视频，分辨率为640×480，帧率为16fps。实验箱大小为0.6m×0.6m，左侧和上侧各有一个圆形小孔，左侧的小孔用于放置摄像头，上侧的小孔用于释放实验箱内的烟雾，背景图像贴于实验箱右侧，也就是录制的烟雾视频的背景。为了使实验箱中的烟雾均匀扩散，分别在右上角与左下角放置一个24v风扇。实验箱底部是一个木质抽屉，抽屉里放置一个铁制托盘，烟雾就是通过点燃抽屉托盘里的棉绳产生的。抽屉上方分别放置一个抽气机与一个烟密计。点燃棉绳产生的烟雾就是通过抽气机抽入到烟箱中的。烟密计用来测量激光发射器与接收器之间的烟雾浓度。总体来说，整个实验过程可以简单地作如下描述：首先点燃抽屉内的棉绳产生烟雾，抽气机将烟雾抽入到实验箱中，风扇使进入到实验箱中的烟雾扩散均匀，烟密计记录下实验箱中烟雾浓度的变化情况，整个过程用摄像头录制成视频。图像结构相似度从图像组成的角度出发，包含亮度、对比度和结构三个因素。给定两幅图像x，y，它们之间的图像结构相似度ssim(x，y)定义为：其中，亮度比较函数：对比度比较函数：结构比较函数：其中，μx为图像x的灰度值均值，μy为图像y的灰度值均值，σx为图像x的无偏估计标准差，σy为图像y的无偏估计标准差，σxy为两幅图像的协方差，c1、c2和c3通常为很小的常数，加入这些常数的目的是为了避免造成分母为0而造成的不稳定现象。一般取：c1＝(k1l)2，c2＝(k2l)2，k1，k2＜＜1，l为图像像素的取值范围。计算视频从第二帧开始的每一帧图像与背景图像(视频第一帧)间的结构相似度，同时利用烟密计记录烟雾浓度在整个视频持续时间内的变化情况。将结构相似度与烟密计测量的烟雾浓度两组数据绘制在同一幅图中，如图3所示。利用matlab中的polyfit函数对图3中的两条曲线进行一阶多项式拟合，得到根据图像结构相似度计算烟雾浓度的转换系数k＝-3.4802，根据转换系数对结构相似度曲线进行f(x)＝-3.4802x变换，得到的曲线将与烟密计测量得到的烟雾浓度曲线接近，结果如图4所示。如图4可以看出，该背景图像下，通过图像结构相似度计算的烟雾浓度很好地吻合了烟密计测量的烟雾浓度，这证明了该方法的可行性。但是，不同背景图像的转换系数不同。当背景图像相对复杂时，即使被一定浓度的烟雾遮挡，仍可以看到背景图像的部分结构信息，此时的图像与背景图像间的结构相似度仍保持较大值，也就是说，复杂的背景图像在被烟雾遮挡后，结构相似度的变化相对缓慢；相反，当背景图像相对简单时，被同样浓度的烟雾遮挡后，背景图像的结构信息损失严重，此时的图像与背景图像间的结构相似度较低，也就是说，简单的背景区域在被烟雾遮挡后，结构相似度的变化相对迅速。如图5所示：假设y1＝x为烟雾浓度的变化曲线，由于背景图像的复杂度越大，在被烟雾遮挡后，与原始背景图像间的结构相似度曲线变化缓慢，因此假设y2＝-x为复杂度较大的背景图像的结构相似度变化曲线，y3＝-2x为复杂度较小的背景图像的结构相似度变化曲线。显然，在通过y2＝-x曲线得到y1＝x曲线时，转换系数为k1＝-1；在通过y3＝-2x曲线得到y1＝x曲线时，转换系数为k2＝-1/2。显然，|k1|＞|k2|。也就是说，转换系数k与背景图像的复杂度相关，且背景图像的复杂度越大，转换系数k的绝对值越大。因此为了统一转换系数，需要对图像的复杂度进行定量描述。这里基于图像纹理特征描述图像复杂度。采用灰度共生矩阵法提取图像纹理特征，在此基础上提取其中五种特征参数：(1)角二阶矩asm：图像的灰度共生矩阵的各个元素值均相等时的asm值小，表示图像的纹理细，能量小；相反，图像的灰度共生矩阵的各个元素值相差较大时的asm值大，表示图像的纹理粗，能量大。(2)对比度con：对比度大的图像纹理深，图像清晰。对比度小的图像纹理浅，图像模糊。(3)逆差分矩idm：逆差分矩大的图像局部纹理的变化少，逆差分矩小的图像局部纹理的变化多。(4)熵ent：熵值大的图像纹理复杂，熵值小的图像纹理简单。(5)相关性cor：相关性的值大代表行或列的相似程度高，表现为图像的复杂度小。相关性的值小代表行或列的相似程度低，表现为图像的复杂度大。为了说明这5种特征参数是如何影响图像复杂度的，这里举例说明：选取4幅大小相同的梅花图像，如图6(a)～6(d)所示。显然，图6(a)～6(d)的图像复杂度依次增加。分别对图6(a)～6(d)所示的四幅图像提取角二阶矩、对比度、逆差分矩、熵、相关性5种特征参数并求平均值，结果如表1所示：表1由表1可以看出，除表中椭圆标记处外，随着图像越来越复杂，角二阶矩的值越来越小，对比度的值越来越大，逆差分矩的值越来越小，熵的值越来越大，相关性的值越来越小。表中三处椭圆标记说明由一种特征参数来衡量图像复杂度会有少量误差，所以本发明用5种特征参数的加权来描述图像复杂度，准确率高，由此，对5种特征参数引入权值，熵、对比度的权值为1，角二阶矩、逆差分矩、相关性的权值为-1，并进行加权求和，作为图像复杂度的值。因此，定义图像复杂度c的计算公式为：c＝con+ent-asm-idm-cor。为了找到转换系数与图像复杂度的关系，这里另外选取九幅不同复杂度的背景图像，如图7(a)～7(i)所示。按照上述方法分别计算转换系数k与图像复杂度c并加以对比，利用计算以图像复杂度的负值代替转换系数的误差值e，所得结果如表2所示。表27(a)7(b)7(c)7(d)7(e)7(f)7(g)7(h)7(i)k＝-1.36k＝-3.18k＝-6.10k＝-1.3k＝-2.4k＝-5.9k＝-3.4k＝-3.6k＝-2.2c＝1.3c＝2.9c＝5.9c＝1.4c＝2.5c＝5.3c＝3.6c＝3.9c＝2.4e＝4％e＝8％e＝3％e＝7％e＝4％e＝10％e＝6％e＝8％e＝9％表2基本上验证了背景图像的复杂度越大，转换系数k的绝对值越大，且对于大部分背景图像来说，以图像复杂度的负值代替转换系数的误差值控制在10％以内。因此，可以确定转换系数k为图像复杂度的负值。也就是说，在任意背景图像条件下，只要计算出该背景图像的复杂度，就可以根据结构相似度计算出烟雾浓度，计算公式为：m＝-c×(ssim-ssim0)，其中m为代表烟雾浓度，单位db/m。c代表图像复杂度，单位db/m。ssim代表烟雾图像与背景图像之间的结构相似度，ssim0代表结构相似度曲线的第一个值，约为1。为了验证本发明方法的正确性和有效性，我们分别在不同烟雾浓度、不同背景条件下进行了若干实验。在不同烟雾浓度下，利用烟密计测量烟雾浓度值，根据图像结构相似度计算出的烟雾浓度值，得到烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图。图8(a)是本发明实施例烟雾图像中3根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(b)是本发明实施例烟雾图像中4根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(c)是本发明实施例烟雾图像中5根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(d)是本发明实施例烟雾图像中6根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(e)是本发明实施例烟雾图像中7根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图8(f)是本发明实施例烟雾图像中8根棉绳同时起烟时的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；由此可见，在不同的烟雾浓度下利用本发明方法计算得到的烟雾浓度与烟雾计实际测量得到的烟雾浓度均吻合，证明了本发明方法的正确性和有效性。图9(a)是本发明实施例中第一背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(b)是本发明实施例中第二背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(c)是本发明实施例中第三背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(d)是本发明实施例中第四背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(e)是本发明实施例中第五背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图；图9(f)是本发明实施例中第六背景图像的烟密计测量值与根据结构相似度计算出的烟雾浓度值的曲线关系图。由此可见，在不同的背景图像下利用本发明方法计算得到的烟雾浓度与烟雾计实际测量得到的烟雾浓度均吻合，证明了本发明方法的正确性和有效性。本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。当前第1页12