发明名称 基于全方位视觉传感器的火灾监控装置
摘要 一种基于全方位计算机视觉的火灾监控方法与装置,主要是利用全方位计算机视觉传感器对现场进行监视,将摄得的连续全方位图像输入计算机,不断进行图像处理和分析,通过火灾的简单的模式分类、火焰面积变化特征判断、火焰形体变化特征判断、火焰颜色特征判断、火焰闪动规律特征判断、火焰整体运动特征判断以及火灾强度特征判断得到一个综合判断结果,根据该计算综合判断量化值进行不同的处理,降低火灾的误报率,本发明提供了一种全新的更快速、更准确、更可靠的多参数智能火灾监控方法与装置。
申请公布号 CN100538757C 申请公布日期 2009.09.09
申请号 CN200510061876.8 申请日期 2005.12.07
申请人 浙江工业大学 发明人 汤一平;金顺敬;叶永杰;邓飞;顾小凯
分类号 G08B17/00(2006.01)I;G08B17/12(2006.01)I 主分类号 G08B17/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 代理人 王 兵;袁木棋
主权项 1、一种基于全方位视觉传感器的火灾监控装置,该火灾监控装置包括微处理器、用于监视现场的视频传感器、用于与外界通信的通信模块,所述的微处理器包括:图像数据读取模块,用于读取从视频传感器传过来的视频图像信息;文件存储模块,用于将视频传感器采集的数据存储到存储器中;现场实时播放模块,用于连接外部的显示设备,将监视现场画面实时播放;所述的视频传感器的输出与微处理器通信连接,其特征在于:所述的视频传感器是全方位视觉传感器,所述的视觉传感器包括用于反射监控领域中物体的外凸反射镜面、透明圆柱体、摄像头,所述的外凸反射镜面朝下,所述的透明圆柱体支撑外凸反射镜面,用于拍摄外凸反射镜面上成像体的摄像头位于透明圆柱体的内部,摄像头位于外凸反射镜面的虚焦点上;所述的微处理器还包括:传感器标定模块,用于对全方位视觉传感器的参数进行标定,建立空间的实物图像与所获得的视频图像的对应关系;颜色模型转换模块,用于将彩色图像中的每一个像素点的颜色从RGB颜色空间转换到(Cr,Cb)空间颜色模型;图像展开处理模块,用于将采集的圆形视频图像展开为全景柱状图,根据圆形全方位图像上的一个点(x<sup>*</sup>,y<sup>*</sup>)和矩形柱状全景图上的一个点(x<sup>**</sup>,y<sup>**</sup>)的对应关系,建立(x<sup>*</sup>,y<sup>*</sup>)与(x<sup>**</sup>,y<sup>**</sup>)的映射矩阵,式(1)所示:<img file="C200510061876C00032.GIF" wi="757" he="107" />上式中,<img file="C200510061876C00033.GIF" wi="92" he="92" />是映射矩阵,<img file="C200510061876C00034.GIF" wi="207" he="84" />是圆形全方位图像上的像素矩阵,<img file="C200510061876C00035.GIF" wi="281" he="104" />是矩形柱状全景图上的像素矩阵;运动对象检测模块,用于将所获得的当前帧现场视频图像与一个相对比较稳定的基准参考图像进行差值运算,图像相减的计算公式如式(2)表示:f<sub>d</sub>(X,t<sub>0</sub>,t<sub>i</sub>)=f(X,t<sub>i</sub>)-f(X,t<sub>0</sub>)  (2)上式中,f<sub>d</sub>(X,t<sub>0</sub>,t<sub>i</sub>)是实时拍摄到图像与基准参考图像间进行图像相减的结果;f(X,t<sub>i</sub>)是实时拍摄到图像;f(X,t<sub>0</sub>)是基准参考图像;并将当前图像中与相邻K帧的图像相减计算公式如(3)所示:f<sub>d</sub>(X,t<sub>i-k</sub>,t<sub>i</sub>)=f(X,t<sub>i</sub>)-f(X,t<sub>i-k</sub>)    (3)上式中,f<sub>d</sub>(X,t<sub>i-k</sub>,t<sub>i</sub>)是实时拍摄到图像与相邻K帧图像间进行图像相减的结果;f(X,t<sub>i-k</sub>)是相邻K帧时的图像;如f<sub>d</sub>(X,t<sub>0</sub>,t<sub>i</sub>)≥阈值、f<sub>d</sub>(X,t<sub>i-k</sub>,t<sub>i</sub>)≥阈值成立时,判定为可疑火焰对象;如f<sub>d</sub>(X,t<sub>0</sub>,t<sub>i</sub>)≥阈值、f<sub>d</sub>(X,t<sub>i-k</sub>,t<sub>i</sub>)&lt;阈值,判定静止对象,并用式(4)来更新替换基准参考图像:<maths num="0001"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>X</mi><mo>,</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&DoubleLeftArrow;</mo><mi>f</mi><mrow><mo>(</mo><mi>X</mi><mo>,</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>-</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>如f<sub>d</sub>(X,t<sub>0</sub>,t<sub>i</sub>)&lt;阈值,判定为静止对象;连通区域计算模块,用于对当前图像进行标记,像素灰度为0的小区表示此小区无可疑火焰,像素灰度为1则表示此小区有可疑火焰,计算当前图像中的像素是否与当前像素周围相邻的某一个点的像素相等,如灰度相等判断为具有连通性,将所有具有连通性的像素作为一个连通区域;模式分类模块,用于判定为可疑火焰后,对每个连通区域求出其面积Si,并依照如下:1)若S<sub>i</sub>&lt;阈值1,则该变化区域为噪声点;2)若S<sub>i</sub>&gt;阈值2,则该变化区域为大面积的红外光线变化;3)若阈值1&lt;S<sub>i</sub>&lt;阈值2,则该变化区域为可疑的火焰区域;判断为可疑的火焰区域后,计算相邻帧变化图像的相似度ε<sub>i</sub>,如公式(5)所示,<maths num="0002"><![CDATA[<math><mrow><msub><mi>&epsiv;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mrow><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>&Omega;</mi></mrow></munder><msub><mi>b</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&cap;</mo><msub><mi>b</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><munder><mi>&Sigma;</mi><mrow><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>&Omega;</mi></mrow></munder><msub><mi>b</mi><mi>i</mi></msub><mrow><mo></mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Union;</mo></mrow><msub><mi>b</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,b<sub>i</sub>(x,y)为上一帧中可疑的火焰区域,b<sub>i+1</sub>(x,y)为当前帧中可疑的火焰区域;根据上述计算结果,模式分类为:1)若ε<sub>i</sub>≤阈值1,则图像模式为快速运动的亮点;2)若ε<sub>i</sub>≥阈值2,则图像模式为固定红外发光区域;3)若阈值1&lt;ε<sub>i</sub>&lt;阈值2,则图像模式为火焰;火焰颜色特征判断模块,用于在(Cr,Cb)空间的分布模型,计算发光源是否落在火焰图像在(Cr,Cb)空间的分布模型内,计算公式由式(6)所示:<maths num="0003"><![CDATA[<math><mrow><msub><mi>W</mi><mi>firecolor</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo></mo><mo>{</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo></mo><mo>[</mo><mi>A</mi><mo>*</mo><msup><mrow><mo>(</mo><mi>Cr</mi><mo>-</mo><mover><mi>Cr</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>2</mn><mi>B</mi><mo>*</mo><mrow><mo>(</mo><mi>Cr</mi><mo>-</mo><mover><mi>Cr</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>*</mo><mrow><mo>(</mo><mi>Cb</mi><mo>-</mo><mover><mi>Cb</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>+</mo><mi>C</mi><mo>*</mo><msup><mrow><mo>(</mo><mi>Cb</mi><mo>-</mo><mover><mi>Cb</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo></mo></mrow><mo>}</mo><mo></mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(6)中,W<sub>firecolor</sub>是颜色特征量,<img file="C200510061876C00043.GIF" wi="68" he="82" />,<img file="C200510061876C00044.GIF" wi="64" he="71" />是火焰点Cr、Cb的样本标准均值,A、B、C分别是由样本标准差和均值计算出来的系数;<img file="C200510061876C00043.GIF" wi="68" he="82" />=144.6;<img file="C200510061876C00044.GIF" wi="64" he="71" />=117.5;A=3.7*10<sup>-3</sup>;B=4.1*10<sup>-3</sup>;C=4.5*10<sup>-3</sup>;火焰面积变化特征判断模块,用于根据每帧图像的发光面积Si进行递推计算,求在下一帧图像的发光面积的递推值<img file="C200510061876C00051.GIF" wi="120" he="58" />计算公式由式(7)给出;<maths num="0004"><![CDATA[<math><mrow><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>*</mo><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><mi>k</mi><mo>*</mo><msub><mi>S</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="C200510061876C00053.GIF" wi="97" he="58" />是下一帧图像的发光面积的递推平均值,<img file="C200510061876C00054.GIF" wi="77" he="57" />是目前帧图像的发光面积的递推平均值,Si是目前帧发光面积的计算值,K为系数,小于1,得到下式(8)<maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mtext>t</mtext><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msubsup></mrow><mn>3</mn></mfrac><mo>></mo><mfrac><mrow><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msubsup><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>i</mi><mrow><mtext>t</mtext><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>4</mn><mo>)</mo></mrow></mrow></msubsup></mrow><mn>3</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>如果上述(8)成立,火焰面积扩展W<sub>fire area</sub>量化值取为1;如上式(8)不成立,火焰面积扩展W<sub>fire area</sub>量化值取为0;火焰闪动规律特征判断模块,用于在上述(8)成立时,将发光面积的计算值Si与目前帧图像的发光面积的递推平均值<img file="C200510061876C0005090610QIETU.GIF" wi="88" he="63" />求差计算,如两者的差的正负发生变化就进行记数,计算在一定时间段中其变化的次数作为火焰闪动的变化频率f<sub>fenquncy</sub>,并将该f<sub>fenquncy</sub>与所设定的阈值f<sub>c</sub>相比较,当f<sub>fenquncy</sub>≥f<sub>c</sub>时,火焰闪动W<sub>fire bicker</sub>量化值取为1;当f<sub>fenquncy</sub>&lt;f<sub>c</sub>时,火焰闪动W<sub>fire bicker</sub>量化值0;火灾强度判断模块,用于在火焰面积扩展W<sub>fire area</sub>量化值为1和火焰颜色特征W fire color的计算值大于0.5时,计算公式如式(9)所示;<maths num="0006"><![CDATA[<math><mrow><msub><mi>W</mi><mi>fire indensity</mi></msub><mo>=</mo><mfrac><msub><mi>S</mi><mi>i</mi></msub><mi>&Sigma;S</mi></mfrac><mo>*</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式(9)中,∑S是整个监控面积,W<sub>fire indensity</sub>是火灾强度值;火灾发生综合判断模块,用于根据火焰模式、火焰颜色特征值、火焰面积、火焰闪动规律综合判断是否有火灾发生,其计算式如(10)所示:W<sub>fire alarm</sub>=K<sub>fire pattern</sub>×ε<sub>i</sub>+K<sub>fire color</sub>×W<sub>fire color</sub>+K<sub>fire area</sub>×W<sub>fire area</sub>   (10)+K<sub>fire bicker</sub>×W<sub>fire bicker</sub>+K<sub>fire indensity</sub>×W<sub>fire indensity</sub>式中:K<sub>fire pattern</sub>为火焰模式的加权系数;      K<sub>fire color</sub>为火焰颜色特征的加权系数;      K<sub>fire area</sub>为火焰面积变化的加权系数;      K<sub>fire bicker</sub>为火焰闪动的加权系数;      K<sub>fire indensity</sub>为火灾强度的加权系数;如K<sub>alarm</sub>≤W<sub>fire alarm</sub>,判定为火灾告警,通过通信模块通知管理人员。
地址 310014浙江省杭州市下城区朝晖六区