发明名称 智能移动终端抗噪声干扰运动目标提取方法
摘要 本发明涉及智能移动终端抗噪声干扰运动目标提取方法,经设定图像空间转换的优化转换参数,原始图像转换处理得到灰度图像后,选择高斯滤波对预处理图像滤波,并选取三幅相同帧间隔的连续预处理图像做差分处理,以克服单独使用三帧差分法提取运动目标质量的不足,二值化处理所得差分图像后,对二值化图像做逻辑“或”以及逻辑“与”处理,以得到潜在运动目标区域的二值化图像,利用逻辑“与”自动摒弃新暴露出来的像素点,消除“鬼影”现象,然后建立背景模型并实时更新,同时基于高斯模型检测潜在运动目标区域内的运动目标,避免了基于帧间差分方法容易引入的“空洞”问题,最后根据95%置信区间理论,对运动前景提取得到运动目标。
申请公布号 CN105405153A 申请公布日期 2016.03.16
申请号 CN201510717145.8 申请日期 2015.10.29
申请人 宁波大学 发明人 郑紫微;熊欧
分类号 G06T7/20(2006.01)I;G06T7/40(2006.01)I;G06T5/10(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 宁波诚源专利事务所有限公司 33102 代理人 刘凤钦
主权项 智能移动终端抗噪声干扰运动目标提取方法,其特征在于,包括以下步骤:(1)对智能移动终端获取的原始图像S(i,j)进行灰度化处理,实现颜色空间转换,得到灰度图像Y(i,j);其中,原始图像S(i,j)转换为灰度图像Y(i,j)采用如下转换公式:Y(i,j)=0.257×r(i,j)+0.504×g(i,j)+0.239×b(i,j);其中,r(i,j)、g(i,j)、b(i,j)分别对应原始图像S(i,j)的颜色空间三维坐标中的红、绿、蓝颜色分量,Y(i,j)是该像素点对应的灰度值;(2)对灰度图像Y进行预处理,得到预处理图像I,并采用高斯滤波对预处理图像I进行处理,高斯函数为均值是零的二维高斯函数g(i,j),其中,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msup><mi>i</mi><mn>2</mn></msup><mo>+</mo><msup><mi>j</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000833140230000011.GIF" wi="518" he="199" /></maths>其使用的图像模板算子为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mn>1</mn><mn>273</mn></mfrac><mo>&times;</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>4</mn></mtd><mtd><mn>7</mn></mtd><mtd><mn>4</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>4</mn></mtd><mtd><mn>16</mn></mtd><mtd><mn>26</mn></mtd><mtd><mn>16</mn></mtd><mtd><mn>4</mn></mtd></mtr><mtr><mtd><mn>7</mn></mtd><mtd><mn>26</mn></mtd><mtd><mn>41</mn></mtd><mtd><mn>26</mn></mtd><mtd><mn>7</mn></mtd></mtr><mtr><mtd><mn>4</mn></mtd><mtd><mn>16</mn></mtd><mtd><mn>26</mn></mtd><mtd><mn>16</mn></mtd><mtd><mn>4</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>4</mn></mtd><mtd><mn>7</mn></mtd><mtd><mn>4</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000833140230000012.GIF" wi="591" he="375" /></maths>(3)对连续间隔m帧的三帧预处理图像I<sub>k‑m</sub>(i,j)、I<sub>k</sub>(i,j)和I<sub>k+m</sub>(i,j)做差分,得到两个差分图像D<sub>(k,k‑m)</sub>(i,j)和D<sub>(k+m,k)</sub>(i,j);其中,D<sub>(k,k‑m)</sub>(i,j)=|I<sub>k</sub>(i,j)‑I<sub>k‑m</sub>(i,j)|,D<sub>(k+m,k)</sub>(i,j)=|I<sub>k+m</sub>(i,j)‑I<sub>k</sub>(i,j)|;其中,I<sub>k‑m</sub>(i,j)、I<sub>k</sub>(i,j)和I<sub>k+m</sub>(i,j)分别表示预处理序列图像I(i,j)中的某三帧预处理图像,m表示预处理序列图像I(i,j)之间相邻间隔的帧数,m∈Z,且m∈[1,5];(4)对所得两个差分图像D<sub>(k,k‑m)</sub>(i,j)和D<sub>(k+m,k)</sub>(i,j)分别做二值化处理,得到对应的二值化图像R<sub>(k,k‑m)</sub>(i,j)和R<sub>(k+m,k)</sub>(i,j);其中,二值化处理准则如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>255</mn></mtd><mtd><mrow><msub><mi>D</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&gt;</mo><mi>T</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>e</mi><mi>l</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>255</mn></mtd><mtd><mrow><msub><mi>D</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&gt;</mo><mi>T</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>e</mi><mi>l</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000833140230000013.GIF" wi="1678" he="167" /></maths>其中,T表示阈值,m表示各序列图像I(i,j)之间相邻间隔的帧数;(5)对所得到的两幅相邻的二值化图像R<sub>(k,k‑m)</sub>(i,j)和R<sub>(k+m,k)</sub>(i,j)做逻辑“或”操作和逻辑“与”操作,分别得到对应运算后的联合二值化“或”图像Bor<sub>k</sub>(i,j)和联合二值化“与”图像Band<sub>k</sub>(i,j),并将联合二值化“或”图像Bor<sub>k</sub>(i,j)和联合二值化“与”图像Band<sub>k</sub>(i,j)进行逻辑“与”运算,得到潜在运动目标区域的二值化图像B<sub>k</sub>(i,j);其中,联合二值化“或”图像为<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Bor</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>255</mn></mtd><mtd><mrow><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&cup;</mo><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&NotEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&cup;</mo><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000833140230000021.GIF" wi="1046" he="167" /></maths>联合二值化“与”图像为<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>Band</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>255</mn></mtd><mtd><mrow><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&cap;</mo><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&NotEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&cap;</mo><msub><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000833140230000022.GIF" wi="1086" he="174" /></maths>潜在运动目标区域的二值化图像<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>B</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>255</mn></mtd><mtd><mrow><msub><mi>Bor</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&cap;</mo><msub><mi>Band</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&NotEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msub><mi>Bor</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&cap;</mo><msub><mi>Band</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000833140230000023.GIF" wi="951" he="158" /></maths>(6)对潜在运动目标区域内进行基于高斯模型的运动目标检测,对高斯分布的参数均值<img file="FDA0000833140230000024.GIF" wi="59" he="85" />以及方差<img file="FDA0000833140230000025.GIF" wi="114" he="79" />赋初值,建立背景模型;其中,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mn>0</mn></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><mi>N</mi><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>,</mo><msup><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mn>0</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><mn>1</mn><mo>/</mo><mi>N</mi><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>&lsqb;</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mn>0</mn></msubsup><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000833140230000026.GIF" wi="982" he="135" /></maths>N表示被用来进行模型初始化的图像的数量,<img file="FDA0000833140230000027.GIF" wi="58" he="78" />表示第k帧图像中(i,j)位置上的像素观察值;(7)对建立的背景模型以设定的更新方式进行实时更新;其中,设定的更新方式为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>&alpha;x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>,</mo><msup><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>&alpha;</mi><msup><mrow><mo>(</mo><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000833140230000028.GIF" wi="1334" he="86" /></maths>其中,α表示当前帧图像中的像素对所对应位置上的背景模型的更新所做出的贡献程度,α=0.0025;(8)计算像素位置上高斯分布的95%置信区间的上限<img file="FDA0000833140230000029.GIF" wi="85" he="75" />和下限<img file="FDA00008331402300000210.GIF" wi="110" he="76" />其中,上限<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>h</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>+</mo><mn>1.96</mn><mo>&times;</mo><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA00008331402300000211.GIF" wi="444" he="79" /></maths>下限<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>l</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>-</mo><mn>1.96</mn><mo>&times;</mo><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mi>k</mi></msubsup><mo>;</mo></mrow>]]></math><img file="FDA00008331402300000212.GIF" wi="437" he="78" /></maths>(9)根据95%置信区间的理论,对运动前景提取得到运动目标<img file="FDA00008331402300000213.GIF" wi="124" he="78" /><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><msubsup><mi>x</mi><mi>l</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>x</mi><mi>h</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mn>255</mn></mtd><mtd><mrow><mi>e</mi><mi>l</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mfenced><mo>.</mo></mrow>]]></math><img file="FDA00008331402300000214.GIF" wi="660" he="166" /></maths>
地址 315211 浙江省宁波市江北区风华路818号