发明名称 一种低照度环境下对地监控的运动目标检测方法
摘要 本发明公开了一种低照度环境下对地监控的运动目标检测方法,属于数字图像处理领域;首先,采用稀疏表示理论中字典学的方法求得字典;然后,采用正交匹配追踪方法求得所监控场景的背景图像和当前帧图像的稀疏系数;最后,根据稀疏系数的分布不同判断运动目标区域。本发明的方法对于低照度引起的大噪声有很好的鲁棒性,保证了对地监控视频在低照度环境下对运动目标检测的准确性和适应性。
申请公布号 CN103870847B 申请公布日期 2015.04.08
申请号 CN201410074755.6 申请日期 2014.03.03
申请人 中国人民解放军国防科学技术大学 发明人 刘煜;张茂军;王炜;徐玮;张政;肖华欣
分类号 G06K9/66(2006.01)I;G06T7/20(2006.01)I 主分类号 G06K9/66(2006.01)I
代理机构 长沙正奇专利事务所有限责任公司 43113 代理人 卢宏
主权项 一种低照度环境下对地监控的运动目标检测方法,其特征在于,包括如下步骤:1)从任意的对地监控视频中选取M帧图像,将所述M帧图像分为N个包含了<img file="FDA0000606096510000011.GIF" wi="208" he="86" />个像素的图像块,所述n为像素个数;将所有的图像块矢量化后,得到字典的训练样本<img file="FDA0000606096510000012.GIF" wi="275" he="106" />2)从训练样本Y中随机选取k个向量作为初始字典D<sub>0</sub>∈R<sup>n×k</sup>;3)令算法循环次数J=1;4)采用正交匹配追踪方法求得y<sub>i</sub>的稀疏系数α<sub>i</sub>,使得对于任意的i=1,2,…N,都满足如下约束:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mi>min</mi><msub><mi>&alpha;</mi><mi>i</mi></msub></munder><msubsup><mrow><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>D</mi><mrow><mi>J</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup></mrow>]]></math><img file="FDA0000606096510000013.GIF" wi="423" he="125" /></maths>  s.t||α<sub>i</sub>||<sub>0</sub><λ其中,λ为正则化参数,那么对于字典D<sub>J‑1</sub>=[d<sub>1</sub>,d<sub>2</sub>,…d<sub>k</sub>],其稀疏系数矩阵A为[α<sub>1</sub>,α<sub>2</sub>,…α<sub>N</sub>],5)利用稀疏系数矩阵A对字典D<sub>J‑1</sub>=[d<sub>1</sub>,d<sub>2</sub>,…d<sub>k</sub>]进行更新:5.1)令j=1;5.2)令<img file="FDA0000606096510000017.GIF" wi="675" he="94" />其中<img file="FDA0000606096510000018.GIF" wi="69" he="76" />表示稀疏系数矩阵A中原子d<sub>j</sub>对应的第j行;5.3)令Ω<sub>j</sub>为N×|ω<sub>j</sub>|大小的矩阵,该矩阵在(ω<sub>j</sub>(i),i)处的值为1,其他点为0;5.3)对于原子d<sub>j</sub>,其稀疏表示误差为<img file="FDA0000606096510000014.GIF" wi="426" he="140" />其中d<sub>e</sub>表示字典中第e个原子,<img file="FDA0000606096510000015.GIF" wi="70" he="87" />表示稀疏系数矩阵A中原子d<sub>e</sub>对应的第e行,其去掉零输入后的系数表示误差为<img file="FDA0000606096510000016.GIF" wi="281" he="92" />5.4)对去掉零输入后的系数表示误差<img file="FDA0000606096510000021.GIF" wi="80" he="91" />进行奇异值分解,得到<img file="FDA0000606096510000022.GIF" wi="317" he="92" />则更新后的<img file="FDA0000606096510000023.GIF" wi="64" he="107" />为矩阵U的第一列,更新后的<img file="FDA0000606096510000024.GIF" wi="79" he="91" />为矩阵V的第一列与Δ(1,1)的乘积;5.5)令j=j+1,直到j=N,就完成对字典D<sub>J‑1</sub>的更新;6)得到更新后的字典D<sub>J</sub>;7)令算法循环次数J=J+1,重复步骤4)~6),直到算法循环次数J达到预设的次数,即得到一个字典D;8)对于任意监控测试视频,假设所监控场景的背景图像为I<sub>B</sub>,对其进行分块,每个图像块的大小为<img file="FDA0000606096510000025.GIF" wi="202" he="92" />个像素,共计M个图像块,采用步骤4)中的正交匹配追踪方法求得每个图像块上的稀疏系数X<sub>B</sub>=[x<sub>1</sub>,x<sub>2</sub>,…x<sub>M</sub>];同理,对于任意的测试图像I<sub>T</sub>,其稀疏系数为X<sub>T</sub>=[x′<sub>1</sub>,x′<sub>2</sub>,…x′<sub>M</sub>];9)利用下式计算所述测试图像I<sub>T</sub>第l个测试图像块的前景检测值I<sub>F</sub>(l):<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>I</mi><mrow><mi>F</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>0</mn></msub><mo>-</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>0</mn></msub><mo>|</mo></mtd></mtr><mtr><mtd><msub><mi>I</mi><mrow><mi>F</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>-</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>|</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000606096510000026.GIF" wi="717" he="252" /></maths>10)判断前景检测值I<sub>F1</sub>(l)和I<sub>F2</sub>(l)是否大于预设阈值T<sub>1</sub>和T<sub>2</sub>,若二者均大于阈值,则认为第l个测试图像块属于前景;否则,认为第l个测试图像块属于背景,从而得到了低照度运动检测结果。
地址 410073 湖南省长沙市开福区德雅路109号