发明名称 一种基于鲁棒主成分分析的监控图像背景建模方法
摘要 本发明公开了一种基于鲁棒主成分分析的监控图像背景建模方法,该方法基于稀疏与低秩矩阵分解理论,以鲁棒主成分分析(RPCA)作为基础,采用截断核范数取代传统的核范数逼近矩阵低秩约束,在增广拉格朗日乘子法框架下,使用一种收敛更快的非精确增广拉格朗日乘子法(Inexact augmented Lagrange multiplier,IALM),直接从监控图像中分离出前景物体和背景模型。该方法恢复出的低秩矩阵即为背景图像矩阵,而稀疏大噪声矩阵则是前景物体位置矩阵。本文方法在动态纹理背景、光照渐变和雾霾天气等较复杂场景中均能够较精确地检测出前景物体,且恢复出的背景矩阵更为低秩,能够更加简洁高效地解决背景建模这一实际问题。
申请公布号 CN106056607A 申请公布日期 2016.10.26
申请号 CN201610373816.8 申请日期 2016.05.30
申请人 天津城建大学 发明人 张艳;刘树东;阮磊;袁兰兰;张亚飞
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 重庆百润洪知识产权代理有限公司 50219 代理人 刘立春
主权项 一种基于鲁棒主成分分析的监控图像背景建模方法,其特征在于:把每一帧图像的像素值按列首尾相接组成一个m维列向量,那么连续的n帧图像表示为m×n维的矩阵D,具有极大相似性的背景部分用待恢复的低秩矩阵A表示,而分布范围很小的前景部分表示为稀疏矩阵E,且D=A+E,在鲁棒主成分分析基础上采用截断核范数取代传统的核范数逼近矩阵低秩约束,截断核范数定义式为<maths num="0001"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><mi>A</mi><mo>|</mo><msub><mo>|</mo><mrow><mi>T</mi><mo>,</mo><mo>*</mo></mrow></msub><mo>=</mo><mo>|</mo><mo>|</mo><mi>A</mi><mo>|</mo><msub><mo>|</mo><mo>*</mo></msub><mo>-</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><msup><mi>XX</mi><mi>T</mi></msup><mo>=</mo><mi>I</mi><mo>,</mo><msup><mi>YY</mi><mi>T</mi></msup><mo>=</mo><mi>I</mi></mrow></munder><mi>T</mi><mi>r</mi><mrow><mo>(</mo><msup><mi>XAZ</mi><mi>T</mi></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001004078300000011.GIF" wi="1396" he="94" /></maths>其中<img file="FDA0001004078300000012.GIF" wi="89" he="69" />表示A的截断核范数,||A||<sub>*</sub>表示A的核范数,Tr(XAZ<sup>T</sup>)为矩阵XAZ<sup>T</sup>的迹,X、Z为任意的矩阵,X∈R<sup>r×m</sup>,Z∈R<sup>r×n</sup>,于是形成式(2)所示的优化问题:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><munder><mi>min</mi><mrow><mi>A</mi><mo>,</mo><mi>E</mi></mrow></munder></mtd><mtd><mrow><mo>|</mo><mo>|</mo><mi>A</mi><mo>|</mo><msub><mo>|</mo><mo>*</mo></msub><mo>-</mo><munder><mi>max</mi><mrow><msup><mi>XX</mi><mi>T</mi></msup><mo>=</mo><mi>I</mi><mo>,</mo><msup><mi>YY</mi><mi>T</mi></msup><mo>=</mo><mi>I</mi></mrow></munder><mi>T</mi><mi>r</mi><mrow><mo>(</mo><msup><mi>XAZ</mi><mi>T</mi></msup><mo>)</mo></mrow><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>E</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub></mrow></mtd><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mi>D</mi><mo>=</mo><mi>A</mi><mo>+</mo><mi>E</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001004078300000013.GIF" wi="1621" he="102" /></maths>其中||E||<sub>1</sub>表示矩阵E的l<sub>1</sub>范数,λ用于平衡A的低秩程度和E的稀疏程度;为求解问题(2),设计了一种迭代方法,第l次迭代为Step1:计算A<sub>l</sub>的奇异值分解U∑V<sup>T</sup>,其中U={u<sub>1</sub>,u<sub>2</sub>…u<sub>m</sub>}∈R<sup>m×m</sup>,V={v<sub>1</sub>,v<sub>2</sub>…v<sub>n</sub>}∈R<sup>n×n</sup>,<img file="FDA0001004078300000014.GIF" wi="487" he="143" />和σ={σ<sub>1</sub>,σ<sub>2</sub>…σ<sub>n</sub>},得到U<sub>r</sub>=(u<sub>1</sub>,u<sub>2</sub>…u<sub>r</sub>)<sup>T</sup>和V<sub>r</sub>=(v<sub>1</sub>,v<sub>2</sub>…v<sub>r</sub>)<sup>T</sup>,r是矩阵A的秩,式(2)可转换为如式(3)所示的凸优化问题:<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><munder><mi>min</mi><mrow><mi>A</mi><mo>,</mo><mi>E</mi></mrow></munder></mtd><mtd><mrow><mo>|</mo><mo>|</mo><mi>A</mi><mo>|</mo><msub><mo>|</mo><mo>*</mo></msub><mo>-</mo><mi>T</mi><mi>r</mi><mrow><mo>(</mo><msubsup><mi>U</mi><mi>r</mi><mi>T</mi></msubsup><msub><mi>AV</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>E</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub></mrow></mtd><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mi>D</mi><mo>=</mo><mi>A</mi><mo>+</mo><mi>E</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001004078300000015.GIF" wi="1517" he="103" /></maths>Step2:采用增广拉格朗日乘子法求解式(3)所示的凸优化问题,迭代直至收敛到这个子问题的最优解,恢复出A和E,实现背景和前景的分离。
地址 300384 天津市西青区津静路26号