发明名称 一种基于图像信息熵的交通参数自动提取方法
摘要 本发明公开了一种基于图像信息熵的交通参数自动提取方法。首先,从彩色视频序列中提取图像帧,并利用混合高斯背景建模方法获得背景图像。其次,利用图像纹理信息对检测区域图像进行预处理。然后,利用图像信息熵与自适应阈值相结合的方法,进行车辆存在检测。最后,计算车辆数、交通流、车道占有率和空间平均速率三个交通参数。实施结果表明,本发明提出的方法解决了运动车辆阴影以及环境光线变化对交通参数提取造成干扰的问题,有效地提高了交通参数提取的准确性。
申请公布号 CN101499214B 申请公布日期 2011.05.11
申请号 CN200910105577.8 申请日期 2009.02.27
申请人 北京大学深圳研究生院 发明人 施行;王一言;邹月娴
分类号 G08G1/01(2006.01)I;G06T7/00(2006.01)I 主分类号 G08G1/01(2006.01)I
代理机构 代理人
主权项 1.一种基于图像信息熵的交通参数自动提取方法,其特征在于,所述方法包括:1)从彩色视频序列中读取一帧原始图像F<sub>n</sub>,并采用混合高斯背景建模算法提取背景图像B<sub>n</sub>;2)分别从当前图像帧F<sub>n</sub>与背景图像B<sub>n</sub>中提取相应检测区域图像FR<sub>n</sub>和BR<sub>n</sub>;3)将图像FR<sub>n</sub>和BR<sub>n</sub>由RGB色彩空间转换为YC<sub>r</sub>C<sub>b</sub>色彩空间,并提取亮度空间图像FRY<sub>n</sub>和BRY<sub>n</sub>;4)将图像FRY<sub>n</sub>和BRY<sub>n</sub>做纹理均方差运算,并将该结果进行二值化处理,从而得到二值图DB<sub>n</sub>;纹理均方差运算的方法为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>d</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><mn>2</mn><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mn>2</mn><mi>N</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>M</mi></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>N</mi></mrow></munderover><mo>[</mo><msub><mi>R</mi><mrow><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>,</mo><mi>FRY</mi></mrow></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mo>-</mo><msub><mi>R</mi><mrow><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>,</mo><mi>BRY</mi></mrow></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msup><mo>]</mo><mn>2</mn></msup><mo>,</mo><mfenced open='' close=''><mtable><mtr><mtd><mi>M</mi><mo>&le;</mo><mi>x</mi><mo>&le;</mo><mi>X</mi><mo>-</mo><mi>M</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>N</mi><mo>&le;</mo><mi>y</mi><mo>&le;</mo><mi>Y</mi><mo>-</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,X,Y为图像尺寸,R<sub>(x,y)</sub>(u,v)的计算方法为:<maths num="0003"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>2</mn><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mn>2</mn><mi>N</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mn>2</mn><mi>M</mi><mo>+</mo><mn>1</mn><mo>-</mo><mi>u</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>2</mn><mi>N</mi><mo>+</mo><mn>1</mn><mo>-</mo><mi>v</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>&times;</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>M</mi><mo>-</mo><mi>u</mi></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mi>v</mi></mrow></msubsup><mi>I</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mi>I</mi><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>u</mi><mo>,</mo><mi>n</mi><mo>+</mo><mi>v</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>M</mi></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>N</mi></mrow></msubsup><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mfenced open='' close=''><mtable><mtr><mtd><mn>0</mn><mo>&le;</mo><mi>u</mi><mo>&le;</mo><mn>2</mn><mi>M</mi></mtd></mtr><mtr><mtd><mn>0</mn><mo>&le;</mo><mi>v</mi><mo>&le;</mo><mn>2</mn><mi>N</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,u,v分别为m,n方向上的坐标,2M+1,2N+1为图像块的尺寸,I(m,n)为坐标点(m,n)处像素点的像素值;二值图DB<sub>n</sub>的计算方法为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>DB</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>d</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>></mo><msub><mi>&tau;</mi><mi>T</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>d</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&tau;</mi><mi>T</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mfenced open='' close=''><mtable><mtr><mtd><mi>M</mi><mo>&le;</mo><mi>x</mi><mo>&le;</mo><mi>X</mi><mo>-</mo><mi>M</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>N</mi><mo>&le;</mo><mi>y</mi><mo>&le;</mo><mi>Y</mi><mo>-</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,τ<sub>T</sub>是设定阀值;5)计算检测区域内的车辆呈现度W<sub>n</sub>以及自适应阈值τ<sub>n</sub>;车辆呈现度W<sub>n</sub>的计算方法为:W<sub>n</sub>=P<sub>V</sub>×H其中,P<sub>v</sub>的计算方法为:<maths num="0006"><![CDATA[<math><mrow><msub><mi>P</mi><mi>v</mi></msub><mo>=</mo><mfrac><msub><mi>K</mi><mi>ab</mi></msub><mi>K</mi></mfrac></mrow></math>]]></maths>其中,K<sub>ab</sub>为动态图块的个数,动态图块为含有动态点的图块,动态点定义为二值图DB<sub>n</sub>中像素值为1的像素点,K为检测区域所分成的图块个数;H的计算方法为:<maths num="0007"><![CDATA[<math><mrow><mi>H</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>p</mi><mi>i</mi></msub><mi>exp</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FSB00000392557800018.GIF" wi="253" he="184" />N<sub>i</sub>为每个图块中动态点的数目;自适应阈值τ<sub>n</sub>的计算方法为:<maths num="0008"><![CDATA[<math><mrow><msub><mi>&tau;</mi><mi>n</mi></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>K</mi><mi>ab</mi></msub></mfrac><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>K</mi><mi>ab</mi></msub></msubsup><msub><mi>p</mi><mi>i</mi></msub><mo>+</mo><msqrt><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>K</mi><mi>ab</mi></msub></msubsup><msup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><mfrac><mn>1</mn><msub><mi>K</mi><mi>ab</mi></msub></mfrac><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>K</mi><mi>ab</mi></msub></msubsup><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>K</mi><mi>ab</mi></msub></mrow></math>]]></maths>6)比较数值W<sub>n</sub>与τ<sub>n</sub>,进行判断,检测区域内若存在车辆,则将当前帧标识符I<sub>n</sub>置为1;若检测区域内不存在车辆,则将当前帧标识符I<sub>n</sub>置为0;W<sub>n</sub>与τ<sub>n</sub>的比较方法为:<maths num="0009"><![CDATA[<math><mrow><msub><mi>I</mi><mi>n</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>W</mi><mi>n</mi></msub><mo>&GreaterEqual;</mo><msub><mi>&tau;</mi><mi>n</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>W</mi><mi>n</mi></msub><mo>&lt;</mo><msub><mi>&tau;</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>7)从视频序列中新读入一帧原始图像F<sub>n+1</sub>,重复执行步骤2至步骤6,最终得到标识符结果为I<sub>n+1</sub>;8)通过比较I<sub>n+1</sub>与I<sub>n</sub>,判定车辆数Q<sub>car</sub>是否进行累加,并相继计算交通流q,车道占有率O,空间平均速率v<sub>s</sub>三项交通参数;9)车辆数Q<sub>car</sub>的计算方法为:<img file="FSB00000392557800023.GIF" wi="690" he="140" />其中I<sub>n-1</sub>为前一帧标识符,I<sub>n</sub>为当前帧标识符;10)交通流q的计算方法为:<maths num="0010"><![CDATA[<math><mrow><mi>q</mi><mo>=</mo><mfrac><msub><mi>Q</mi><mi>car</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mfrac></mrow></math>]]></maths>其中,T<sub>d</sub>为观测时间;11)车道占有率O的计算方法为:<maths num="0011"><![CDATA[<math><mrow><mi>O</mi><mo>=</mo><mfrac><msub><mi>T</mi><mi>p</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><msub><mi>T</mi><mi>p</mi></msub><mo>=</mo><mfrac><msub><mi>N</mi><mi>f</mi></msub><mi>Fr</mi></mfrac></mrow></math>]]></maths>其中,T<sub>p</sub>为车辆通过检测区域的时间累计值,N<sub>f</sub>为在观测时间T<sub>d</sub>内标识符I<sub>n</sub>=1的帧数目,Fr为帧率;12)空间平均速率v<sub>s</sub>的计算方法为:<maths num="0013"><![CDATA[<math><mrow><msub><mi>v</mi><mi>s</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><mi>car</mi></msub><mrow><mo>(</mo><mi>L</mi><mo>+</mo><mi>d</mi><mo>)</mo></mrow></mrow><msub><mi>T</mi><mi>p</mi></msub></mfrac></mrow></math>]]></maths>其中,L为检测区域的长度,d为通过车辆的平均长度。
地址 518055 广东省深圳市南山区西丽镇丽水路大学城北京大学校区