发明名称 基于3维Log-Gabor变换的视频图像融合性能评价方法
摘要 本发明公开了一种基于3维Log-Gabor变换的视频图像融合性能评价方法,主要解决输入视频有噪声或背景运动情况下,现有技术不能准确评价融合算法性能的问题。其实现步骤是:利用3维Log-Gabor变换对输入视频和融合后视频进行多方向,多尺度分解;用视频图像的3维相位一致性构建时-空相位一致性评价因子;用视频图像的3维Log-Gabor变换系数幅值构建时-空信息提取评价因子;将时-空一致性评价因子和时-空信息提取评价因子结合,构建全局时-空性能评价因子,根据该因子的计算结果评价视频融合算法性能;结合人眼视觉ST-CSF公式、3维梯度结构张量设计局部和全局参数。本发明能准确评价有噪声或背景运动情况下融合算法性能,可用于评价视频图像融合算法性能。
申请公布号 CN102946548B 申请公布日期 2015.02.18
申请号 CN201210493342.2 申请日期 2012.11.27
申请人 西安电子科技大学 发明人 张强;陈闵利;王龙
分类号 H04N17/00(2006.01)I 主分类号 H04N17/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于3维Log‑Gabor变换的视频图像融合性能评价方法,包括如下步骤:(1)在3维频域球坐标系<img file="FDA0000588200840000011.GIF" wi="176" he="73" />中,构建3维Log‑Gabor滤波器:<img file="FDA0000588200840000012.GIF" wi="1233" he="198" />其中,k=1,2,..,N表示第k通道的Log‑Gabor滤波器,N为总的通道数,ρ<sub>k</sub>表示该通道滤波器的中心径向频率,θ<sub>k</sub>和<img file="FDA0000588200840000013.GIF" wi="54" he="61" />分别为相应通道滤波器的方向角和仰角,θ<sub>k</sub>和<img file="FDA0000588200840000014.GIF" wi="57" he="59" />共同决定了滤波器的中心方向,σ<sub>k,ρ</sub>和σ<sub>k,ζ</sub>分别决定该滤波器的径向带宽和角度带宽,<img file="FDA0000588200840000015.GIF" wi="290" he="76" />为球坐标系中某一点<img file="FDA0000588200840000016.GIF" wi="178" he="74" />与该滤波器中心频点<img file="FDA0000588200840000017.GIF" wi="226" he="75" />之间的夹角;(2)采用3维Log‑Gabor滤波器对两个输入视频Va、Vb和两个输入视频融合后的视频Vf分别进行多尺度,多方向分解,得到第一个输入视频Va的复数分解系数<img file="FDA0000588200840000018.GIF" wi="334" he="85" />第二个输入视频Vb的复数分解系数<img file="FDA0000588200840000019.GIF" wi="292" he="86" />和融合后的视频Vf的复数分解系数<img file="FDA00005882008400000110.GIF" wi="322" he="90" />(3)把两个输入视频Va、Vb和两个输入视频融合后视频Vf分别分成M个大小相同且互补重叠的3维时‑空子块R(j),其中j=1,2,...,M;(4)针对第j个子块R(j),构建时‑空一致性融合性能评价因子Q<sub>STC</sub>(Va,Vb,Vf|R<sub>j</sub>):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>STC</mi></msub><mrow><mo>(</mo><mi>Va</mi><mo>,</mo><mi>Vb</mi><mo>,</mo><mi>Vf</mi><mo>|</mo><msub><mi>R</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><msub><mi>&omega;</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>Z</mi><mi>fa</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>Z</mi><mi>fb</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&omega;</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&omega;</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mtd><mtd><msub><mi>Z</mi><mi>ab</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>></mo><mo>=</mo><mi>Th</mi><mn>1</mn></mtd></mtr><mtr><mtd><mi>max</mi><mrow><mo>(</mo><msub><mi>Z</mi><mi>fa</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>Z</mi><mi>fb</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mtd><mtd><msub><mi>Z</mi><mi>ab</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>&lt;</mo><mrow><mi>Th</mi><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00005882008400000111.GIF" wi="1433" he="242" /></maths>其中,ω<sub>a</sub>(j)和ω<sub>b</sub>(j)分别为输入视频Va和输入视频Vb相应子块R(j)的权值,Z<sub>fa</sub>(j)是PC<sub>f</sub>(x,y,t)和PC<sub>a</sub>(x,y,t)在子块R(j)内的去均值归一化互相关系数,Z<sub>fb</sub>(j)是PC<sub>f</sub>(x,y,t)和PC<sub>b</sub>(x,y,t)在子块R(j)内的去均值归一化互相关系数,Z<sub>ab</sub>(j)是PC<sub>a</sub>(j)和PC<sub>b</sub>(x,y,t)在子块R(j)内的去均值归一化互相关系数,PC<sub>a</sub>(x,y,t)、PC<sub>b</sub>(x,y,t)和PC<sub>f</sub>(x,y,t)分别为两个输入视频Va、Vb和两个输入视频融合后视频Vf的时‑空相位一致性特征,Th1为用于判断输入视频图像在该区域内相似度关系的阈值,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>Th</mi><mn>1</mn><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>Z</mi><mi>ab</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000588200840000021.GIF" wi="404" he="148" /></maths>(5)针对第j个子块R(j),构建时‑空信息提取融合性能评价因子Q<sub>STE</sub>(Va,Vb,Vf|R<sub>j</sub>):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>STE</mi></msub><mrow><mo>(</mo><mi>Va</mi><mo>,</mo><mi>Vb</mi><mo>,</mo><mi>Vf</mi><mo>|</mo><msub><mi>R</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&chi;</mi><mi>k</mi></msub><msubsup><mi>Q</mi><mi>STE</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>Va</mi><mo>,</mo><mi>Vb</mi><mo>,</mo><mi>Vf</mi><mo>|</mo><msub><mi>R</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&chi;</mi><mi>k</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000588200840000022.GIF" wi="980" he="275" /></maths>其中,N表示滤波器通道总数,χ<sub>k</sub>是第k通道下局部区域R(j)对应的权系数,<img file="FDA0000588200840000023.GIF" wi="412" he="86" />是第k通道下子块R(j)的时‑空信息提取融合性能评价因子,“|”为参与运算各个视频与指定运算区域的分隔线;(6)按照如下方式对时‑空一致性融合性能评价因子Q<sub>STC</sub>(Va,Vb,Vf|R<sub>j</sub>)和时‑空信息提取融合性能评价因子Q<sub>STE</sub>(Va,Vb,Vf|R<sub>j</sub>)进行组合,构建第j个子块R(j)的局部时‑空融合性能评价因子:<img file="FDA0000588200840000024.GIF" wi="1508" he="106" />其中,<img file="FDA0000588200840000026.GIF" wi="70" he="80" />表示一种非负操作,当操作数为正数时,输出结果为原数,否则为0;参数δ取值范围为[0,1];(7)构建全局时‑空融合性能评价因子Q<sub>ST_CE</sub>(Va,Vb,Vf)为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>ST</mi><mo>_</mo><mi>CE</mi></mrow></msub><mrow><mo>(</mo><mi>Va</mi><mo>,</mo><mi>Vb</mi><mo>,</mo><mi>Vf</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>&gamma;</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><msub><mi>Q</mi><mrow><mi>ST</mi><mo>_</mo><mi>CE</mi></mrow></msub><mrow><mo>(</mo><mi>Va</mi><mo>,</mo><mi>Vb</mi><mo>,</mo><mi>Vf</mi><mo>|</mo><msub><mi>R</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mi>&gamma;</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000588200840000025.GIF" wi="1054" he="286" /></maths>其中,γ(j)为第j个子块R(j)的全局权系数,γ(j)=max(ω<sub>a</sub>(j),ω<sub>b</sub>(j)),式中,ω<sub>a</sub>(j)和ω<sub>b</sub>(j)为构建局部时‑空一致性客观评价因子Q<sub>STC</sub>(Va,Vb,Vf|R<sub>j</sub>)过程中所采用的局部权系数;(8)根据全局时‑空融合性能评价因子Q<sub>ST_CE</sub>(Va,Vb,Vf)的计算结果,判断视频融合算法性能,计算结果越大则表明视频融合性能越好。
地址 710071 陕西省西安市太白南路2号