发明名称 基于参数自优化多群粒子群优化的视频超分辨重建方法
摘要 本发明公开了一种基于参数自优化的多子群粒子群优化的视频超分辨重建方法,主要解决已有技术中运动估计误差造成重建质量较低,无法较好地完成目标高速运动场景重建、计算复杂度高的问题。其实现过程是:通过获取同一场景的若干低分辨率视频序列,将动态场景内的点表示为空间和时间的三维坐标,建立超分辨率前后的视频序列之间的约束条件,统一考虑时域和空域重建问题,选择恰当的目标函数,采用参数自优化的多群优化方法,把超分辨率问题转化成为满足约束条件的优化问题进行求解,求解结果即为最佳高分辨率重建结果。本发明具有重建视频质量高,不受运动估计影响,智能算法自适应性好,计算量小的优点,可应用于高分辨率视频的获取。
申请公布号 CN102355589B 申请公布日期 2013.09.25
申请号 CN201110327875.9 申请日期 2011.10.26
申请人 西安电子科技大学 发明人 闫允一;郭宝龙;胡颖颖;吴宪祥;朱娟娟;孙伟
分类号 H04N9/888(2006.01)I;G06N3/00(2006.01)I 主分类号 H04N9/888(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于参数自优化多群粒子群优化的视频超分辨重建方法,包括如下步骤:(1)采集同一场景中的若干低分辨率视频序列,将所有元素按照x-y-t的辞书顺序排列成向量<img file="FDA00003001032200011.GIF" wi="57" he="60" />(2)建立视频超分辨重建过程中的约束矩阵:<maths num="0001"><![CDATA[<math><mrow><mi>M</mi><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>=</mo><mover><mi>L</mi><mo>&RightArrow;</mo></mover><mo>+</mo><mover><mi>N</mi><mo>&RightArrow;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow></math>]]></maths>其中M是一个与空时模糊核函数对应的稀疏Toeplitz矩阵,向量<img file="FDA00003001032200013.GIF" wi="54" he="61" />由为待求解高分辨率序列中的所有元素按照x-y-t的辞书式顺序排列而成的向量,<img file="FDA00003001032200014.GIF" wi="59" he="73" />为噪声向量;(3)设超分辨率重建过程中的目标函数为f,其具体形式为:<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mo>=</mo><msup><mrow><mrow><mo>|</mo><mo>|</mo><mi>M</mi><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>-</mo><mover><mi>L</mi><mo>&RightArrow;</mo></mover></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>Sm</mi><msub><mrow><mo>(</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><mi>Dt</mi><msub><mrow><mo>(</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>2</mn><mo>)</mo></mrow><mo></mo></mrow></math>]]></maths>其中:<maths num="0003"><![CDATA[<math><mrow><mi>Sm</mi><msub><mrow><mo>(</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>&lambda;</mi><mi>x</mi></msub><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>W</mi><mi>x</mi></msub><msub><mi>G</mi><mi>x</mi></msub><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&lambda;</mi><mi>y</mi></msub><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>W</mi><mi>y</mi></msub><msub><mi>G</mi><mi>y</mi></msub><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&lambda;</mi><mi>t</mi></msub><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>W</mi><mi>t</mi></msub><msub><mi>G</mi><mi>t</mi></msub><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>表示图像光滑程度,<maths num="0004"><![CDATA[<math><mrow><mi>Dt</mi><msub><mrow><mo>(</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>D</mi><mi>x</mi></msub><mrow><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo></mrow></mrow><mn>2</mn></msup><mo>+</mo><mo>|</mo><mo>|</mo><msup><mrow><msub><mi>D</mi><mi>y</mi></msub><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>&epsiv;</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow></math>]]></maths>表示图像的细节保持程度,式中,λ<sub>x</sub>表示x方向上的图像平滑系数,W<sub>x</sub>为x方向上权重矩阵,G<sub>x</sub>表示沿着x方向的二阶差分算子,λ<sub>y</sub>表示y方向上的图像平滑系数,W<sub>y</sub>为y方向上权重矩阵,G<sub>y</sub>表示沿着y方向的二阶差分算子,λ<sub>t</sub>表示时域t方向上的图像平滑系数,W<sub>t</sub>为时域t方向上权重矩阵,G<sub>t</sub>表示沿着t方向的二阶差分算子,D<sub>x</sub>表示沿x方向的梯度算子,D<sub>y</sub>表示沿y方向的梯度算子,ε是一个正的常数,当时域比空域更平滑,则令λ<sub>x</sub>=λ<sub>y</sub><λ<sub>t</sub>,当空域比时域更平滑,则令λ<sub>x</sub>=λ<sub>y</sub>>λ<sub>t</sub>;(4)以式1)为约束条件,以式2)为优化的目标函数,高分辨视频序列的向量<img file="FDA00003001032200018.GIF" wi="50" he="57" />为令式2)取得最小值的解,其求解过程如下:(4a)设定求解种群的粒子个数Z<sub>s</sub>,每个求解种群的粒子为一个候选的高分辨率视频序列向量<img file="FDA00003001032200021.GIF" wi="92" he="63" />初始化求解种群,即将低分辨视频序列在空域或时域进行线性插值;(4b)设定参数配置种群的粒子个数Z<sub>c</sub>,每个配置种群内的C-粒子为求解种群优化过程所需的参数,并按如下步骤初始化参数配置种群:(4b1)将K-子群进化参数ωk,ck<sub>1</sub>,ck<sub>2</sub>,r-子群进化参数ωr,cr<sub>1</sub>,cr<sub>2</sub>,K-子群的繁殖率ρ<sub>k</sub>,这7个参数的取值范围设为[0,4];(4b2)将r-子群的繁殖率ρ<sub>r</sub>的取值范围设为[0,20];(4b3)将r-子群所占比重P<sub>r</sub>的取值范围设为[0,0.5];(4b4)将K-子群的速度收缩因子比T<sub>k</sub>,r-子群的速度收缩因子比T<sub>r</sub>的取值范围设为[0,1];(4c)计算求解种群中的每个粒子的对应的目标函数值,按函数值从大到小排列求解粒子,将前Z<sub>s</sub>*P<sub>r</sub>,0<P<sub>r</sub><1个粒子划为r-子群,其余的Z<sub>s</sub>*(1-P<sub>r</sub>)个粒子为K-子群;(4d)在每个配置种群粒子的配置下,将r-子群和K-子群分别按照r-策略和K-策略进行优化:(4d1)设K-子群的进化参数与r-子群的进化参数条件为:<maths num="0005"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>&omega;k</mi><mo>></mo><mi>&omega;r</mi></mtd></mtr><mtr><mtd><msub><mi>ck</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>cr</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>ck</mi><mn>2</mn></msub><mo>&lt;</mo><msub><mi>cr</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>k</mi></msub><mo>&lt;</mo><msub><mi>T</mi><mi>r</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>其中,ωk是K-子群的参数惯性因子,ck<sub>1</sub>是K-子群的自身认知因子,ck<sub>2</sub>是K-子群的社会认知因子,T<sub>k</sub>是K-子群的速度收缩因子比;ωr是r-子群的参数惯性因子,cr<sub>1</sub>是r-子群的自身认知因子,cr<sub>2</sub>是r-子群的社会认知因子,T<sub>r</sub>是r-子群的速度收缩因子比;(4d2)设K-子群每个粒子产生ρ<sub>k</sub>个后代,r-子群每个粒子产生ρ<sub>r</sub>个后代,且满足条件ρ<sub>r</sub>>ρ<sub>k</sub>;(4d3)将K-子群和r-子群的所有后代一起按目标函数值从小到大的排序排列;(4d4)设定违反约束程度阈值Vc,计算各粒子违反约束程度值,若粒子违反约束程度值小于Vc,则该粒子可行,否则不可行;(4d5)在满足违反约束条件阈值条件下,保留前Z<sub>s</sub>个粒子,作为进化后的新粒子,其余粒子舍弃;若满足违反约束条件阈值粒子数少于Z<sub>s</sub>时,将其余粒子按违反约束条件从小到大排序,选取较小部分使后代粒子总数达到Z<sub>s</sub>;(4e)求解当前配置下的种群最佳粒子,并以进化得到的最佳粒子适应度来评价对应的配置粒子,执行配置种群的优化过程;(4f)配置种群优化过程结束后得到的最佳求解粒子,即为最佳重建的高分辨率视频序列的向量<img file="FDA00003001032200031.GIF" wi="72" he="63" />(5)将<img file="FDA00003001032200032.GIF" wi="50" he="56" />重新排列成为三维矩阵,获得重建后的高分辨率视频序列。
地址 710071 陕西省西安市太白南路2号