发明名称 一种基于I1I2I3颜色空间的火焰检测方法
摘要 一种基于I1I2I3颜色空间的火焰检测方法,主要通过分辨图像像素的颜色分布确定是否属于火焰区域。本发明提出的方法中把图像分块处理,逐块计算每个子块像素在I1I2I3颜色空间的三个均值;进一步地利用多个规则,即等式组决定该三个均值是否属于火焰颜色。该方法简单易行,具有高实时性。本方法不涉及视频序列,因此与镜头中可能存在的运动无关。因此,本方法不仅可以检出中图片中是否存在火焰区域,而且保证了检测效果的稳定性,满足实时性的要求。适用于基于计算机视觉的森林着火检测,大型场馆着火检测等领域。
申请公布号 CN104933739A 申请公布日期 2015.09.23
申请号 CN201510424599.6 申请日期 2015.07.17
申请人 山东大学 发明人 刘兆广;张行愚;杨阳
分类号 G06T7/40(2006.01)I 主分类号 G06T7/40(2006.01)I
代理机构 济南金迪知识产权代理有限公司 37219 代理人 吕利敏
主权项 一种基于I1I2I3颜色空间的火焰检测方法,其特征在于,所述检测方法包括步骤如下:1)通过公式(I),将目标图像中每个像素值从RGB颜色空间转换到I1I2I3颜色空间:<img file="FDA0000762241170000011.GIF" wi="1162" he="510" />2)将所述目标图像分给成若干个4x4的子块,对每个子块分配一个是否为火焰的标记,形成标记矩阵为flag(i,j),其中i,j为子块在flag矩阵中的横纵坐标;3)针对每个子块,分别计算I1,I2,I3三个分量的均值:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>m</mi><msup><mn>1</mn><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>4</mn><mo>*</mo><mi>i</mi></mrow><mrow><mn>4</mn><mo>*</mo><mi>i</mi><mo>+</mo><mn>3</mn></mrow></munderover><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>y</mi><mo>=</mo><mn>4</mn><mo>*</mo><mi>j</mi></mrow><mrow><mn>4</mn><mo>*</mo><mi>j</mi><mo>+</mo><mn>3</mn></mrow></munderover><mi>I</mi><mn>1</mn><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000762241170000012.GIF" wi="827" he="203" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>m</mi><msup><mn>2</mn><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>4</mn><mo>*</mo><mi>i</mi></mrow><mrow><mn>4</mn><mo>*</mo><mi>i</mi><mo>+</mo><mn>3</mn></mrow></munderover><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>y</mi><mo>=</mo><mn>4</mn><mo>*</mo><mi>j</mi></mrow><mrow><mn>4</mn><mo>*</mo><mi>j</mi><mo>+</mo><mn>3</mn></mrow></munderover><mi>I</mi><mn>2</mn><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mo>(</mo><mrow><mi>I</mi><mi>I</mi></mrow><mo>)</mo></mrow>]]></math><img file="FDA0000762241170000013.GIF" wi="1577" he="204" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>m</mi><msup><mn>3</mn><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>4</mn><mo>*</mo><mi>i</mi></mrow><mrow><mn>4</mn><mo>*</mo><mi>i</mi><mo>+</mo><mn>3</mn></mrow></munderover><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>y</mi><mo>=</mo><mn>4</mn><mo>*</mo><mi>j</mi></mrow><mrow><mn>4</mn><mo>*</mo><mi>j</mi><mo>+</mo><mn>3</mn></mrow></munderover><mi>I</mi><mn>3</mn><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000762241170000014.GIF" wi="845" he="201" /></maths>4)均值归一化,对步骤3)中的三个分量的均值归一化处理,处理后的数据范围为0‑1,方法如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>m</mi><mn>1</mn><mo>=</mo><mfrac><mrow><mi>m</mi><msup><mn>1</mn><mo>&prime;</mo></msup></mrow><mn>255</mn></mfrac></mrow>]]></math><img file="FDA0000762241170000015.GIF" wi="308" he="202" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>m</mi><mn>2</mn><mo>=</mo><mfrac><mrow><mi>m</mi><msup><mn>2</mn><mo>&prime;</mo></msup></mrow><mn>255</mn></mfrac></mrow>]]></math><img file="FDA0000762241170000016.GIF" wi="323" he="193" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>m</mi><mn>3</mn><mo>=</mo><mfrac><mrow><mi>m</mi><msup><mn>3</mn><mo>&prime;</mo></msup></mrow><mn>255</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mrow><mi>I</mi><mi>I</mi><mi>I</mi></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000762241170000017.GIF" wi="1062" he="227" /></maths>5)预先在I1I2I3颜色空间统计出火焰对应区域,当所述m1、m2、m3的值分别位于所述火焰对应区域时,则认为该4x4子块为火焰区域,具体判别算法满足下述规则:a)m1和m2的值满足等式(IV)或者等式(V)<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mrow><mo>(</mo><mn>1.3085</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>-</mo><mn>0.618</mn><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mrow><mo>(</mo><mo>-</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>1.6</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mrow><mo>(</mo><mn>5.4</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>-</mo><mn>3.846</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mrow><mi>I</mi><mi>V</mi></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000762241170000021.GIF" wi="1219" he="233" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mrow><mo>(</mo><mn>1.3085</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>-</mo><mn>0.618</mn><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mrow><mo>(</mo><mo>-</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>1.6</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mrow><mo>(</mo><mo>-</mo><mn>5.875</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>5.0612</mn><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000762241170000022.GIF" wi="1239" he="235" /></maths>b)m1和m3的值满足等式(VI)<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mn>2.8</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>-</mo><mn>1.044</mn><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mn>1.7857</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>-</mo><mn>0.5064</mn><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mo>-</mo><mn>4.6667</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>3.8167</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mn>0.5</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>0.51</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mn>0.8667</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>0.3047</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mn>7.5</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>-</mo><mn>2.415</mn><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mo>-</mo><mn>3.3333</mn><mo>*</mo><mi>m</mi><mn>1</mn><mo>+</mo><mn>1.5933</mn><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mrow><mi>V</mi><mi>I</mi></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000762241170000023.GIF" wi="1259" he="532" /></maths>c)m2和m3的满足等式(VII)<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mrow><mn>1.4375</mn><mo>*</mo><mi>m</mi><mn>2</mn><mo>-</mo><mn>0.0656</mn></mrow><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mrow><mo>-</mo><mn>1.2727</mn><mo>*</mo><mi>m</mi><mn>2</mn><mo>+</mo><mn>1.6418</mn></mrow><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mrow><mo>-</mo><mn>0.2857</mn><mo>*</mo><mi>m</mi><mn>2</mn><mo>+</mo><mn>1.1286</mn></mrow><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mrow><mn>1.1429</mn><mo>*</mo><mi>m</mi><mn>2</mn><mo>+</mo><mn>0.4857</mn></mrow><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>2</mn><mo>-</mo><mn>0.38</mn><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mn>3</mn><mo>-</mo><mrow><mo>(</mo><mrow><mo>-</mo><mn>1.4444</mn><mo>*</mo><mi>m</mi><mn>2</mn><mo>+</mo><mn>1.2889</mn></mrow><mo>)</mo></mrow><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mrow><mi>V</mi><mi>I</mi><mi>I</mi></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000762241170000024.GIF" wi="1354" he="456" /></maths>6)针对flag矩阵,使用3x3的窗口,从上到下,从左到右,步长为1,滑过所述flag矩阵,当4x4子块的均值m1、m2、m3满足步骤5)所述的等式,则将该4x4子块标记为1;如果所述3x3的窗口中超过半数的窗口标记为1,则判断位于3x3的窗口中心的子块所在的4x4子块所对应的图像区域为火焰区域,否则为非火焰区域。
地址 250199 山东省济南市历城区山大南路27号
您可能感兴趣的专利