发明名称 一种基于NMI特征的静止目标消除方法
摘要 本发明属于光电产品应用技术领域,具体涉及一种基于NMI(Normalized Moment of Inertia)特征的静止目标消除方法。本算法首先利用多高斯背景建模方法和三帧累积帧差分的方法获得运动目标前景区域,记为Fgb和Fgd。然后分别对Fgb和Fgd进行前景区域标记。本算法采用图像的NMI(Normalized Moment of Inertia)不变特性,即图像的归一化转动惯量作为特征匹配参数,NMI特征匹配正确率高、计算量小、速度快,而且具有抗灰度畸变性及对几何畸变有较好的保持不变性。
申请公布号 CN104463902A 申请公布日期 2015.03.25
申请号 CN201310452893.9 申请日期 2013.09.25
申请人 北京环境特性研究所 发明人 杨文佳;柴智;李亚鹏;李香祯;王楠
分类号 G06T7/20(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 一种基于NMI特征的静止目标消除方法,其特征在于:包括以下步骤:1)在拍摄设备静止状态下采集当前视频序列图像;2)采用多高斯混合模型算法检测运动目标,获取运动目标前景图Fgb;3)采用三帧累积差分算法检测运动目标,获取运动目标前景图Fgd;4)计算Fgb图中前景区域的质心(x<sub>o</sub>,y<sub>o</sub>),设区域大小为(M,N),f(x,y)为在图像位置(x,y)处的值,前景为1,背景为0:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>x</mi><mi>o</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>x</mi><mo>&CenterDot;</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac></mtd><mtd><msub><mi>y</mi><mi>o</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>y</mi><mo>&CenterDot;</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000387289780000011.GIF" wi="1613" he="284" /></maths>5)计算围绕Fgb图中质心(x<sub>o</sub>,y<sub>o</sub>)的转动惯量记为<img file="FDA0000387289780000015.GIF" wi="159" he="70" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>J</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>o</mi></msub><mo>,</mo><msub><mi>y</mi><mi>o</mi></msub><mo>)</mo></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>|</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>x</mi><mi>o</mi></msub><mo>,</mo><msub><mi>y</mi><mi>o</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000387289780000012.GIF" wi="1576" he="157" /></maths>6)计算Fgb图中归一化后转动惯量;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>NMI</mi><mo>=</mo><mfrac><msqrt><msub><mi>J</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>o</mi></msub><mo>,</mo><msub><mi>y</mi><mi>o</mi></msub><mo>)</mo></mrow></msub></msqrt><mi>m</mi></mfrac><mo>=</mo><mfrac><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></msqrt><mrow><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000387289780000013.GIF" wi="1121" he="303" /></maths>7)采用相同的计算方法对三帧累积差分得到的前景区域Fgd进行归一化转动惯量的计算,得到NMI2;8)利用NMI1和NMI2计算α值:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mrow><mo>|</mo><mi>NMI</mi><mn>1</mn><mo>-</mo><mi>NMI</mi><mn>2</mn><mo>|</mo></mrow><mrow><mi>NMI</mi><mn>1</mn></mrow></mfrac></mrow>]]></math><img file="FDA0000387289780000014.GIF" wi="413" he="141" /></maths>9)当α值小于阈值Th时,则判断该区域背景残余区域,不进行标记。
地址 100076 北京市海淀区永定路50号