主权项 |
一种基于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>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>x</mi><mo>·</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</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>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>y</mi><mo>·</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</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>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</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>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</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>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</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>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</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>α</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时,则判断该区域背景残余区域,不进行标记。 |