发明名称 一种结合视觉分析的图像清晰度衡量方法
摘要 本发明公开了一种结合视觉分析的图像清晰度衡量方法,本发明根据视觉传输模型实现分频率分析,获得了分频率分析后的调制图像f;然后通过获得了分频率分析后的调制图像f;本发明提出的方法结合了视觉传输模型,更为符合主观要求,实现主客观的一致性评价;此外,利用局部强度差异来衡量图像清晰度,也符合清晰度变化的客观要求,整个方法运行高效而稳定。
申请公布号 CN105118060A 申请公布日期 2015.12.02
申请号 CN201510511802.3 申请日期 2015.08.19
申请人 杭州电子科技大学 发明人 赵巨峰;高秀敏;逯鑫淼;张辉朝
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 杭州君度专利代理事务所(特殊普通合伙) 33240 代理人 杜军
主权项 一种结合视觉分析的图像清晰度衡量方法,其特征在于,包括如下步骤:(1)、根据视觉传输模型实现分频率分析对于待评价的原始观测图像g,空间域像素坐标用(i,j)表示,并且其在傅里叶频率域的表达为G,频率域像素坐标用(u,v)表示;使用的视觉传输模型为:A(ρ)=(0.2+22.995ρ)e<sup>‑0.92ρ</sup>,其中ρ为归一化空间频率,量纲为周/像素;利用视觉传输模型实现对图像的分频率分析,分析过程如下:利用A(ρ)=(0.2+22.995ρ)e<sup>‑0.92ρ</sup>对G进行分频率分析,首先将A(ρ)表达为与G等尺寸的图像;假设G的像素尺寸为(2M+1)×(2N+1),则其截止空间频率为<img file="FDA0000784431070000011.GIF" wi="293" he="124" />周/像素宽度,该数值即对角线的一半;由于模型A(ρ)中的空间频率ρ是归一化的,对于G,其归一化因子就是<img file="FDA0000784431070000012.GIF" wi="314" he="124" />于是对ρ进行反归一化,则G的空间频率<img file="FDA0000784431070000013.GIF" wi="504" he="150" />于是可以获得对应的视觉传输模型A(ρ<sub>0</sub>);于是,对于G中的任意像素位置(u,v),空间频率<img file="FDA0000784431070000014.GIF" wi="318" he="90" />于是根据A(ρ<sub>0</sub>)可构造(2M+1)×(2N+1)尺寸的频率分析图A<sub>map</sub>;根据获取的A<sub>map</sub>,对G进行分频率分析调制,求取调制图F=GA<sub>map</sub>,即G与A<sub>map</sub>之间像素对像素相乘;通过傅里叶反变换,获取在空域的调制图像f;于是,获得了分频率分析后的调制图像f;(2)、局部强度差异衡量图像清晰度对步骤1中的分频率分析后的调制图像f进行清晰度的考察;在局部区域内,本发明定义图像f在(i,j)处的局部强度差异Dif记录如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>D</mi><mi>i</mi><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>=</mo><msub><mi>&Omega;</mi><mrow><mi>c</mi><mi>e</mi><mi>n</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>,</mo><mo>&ForAll;</mo><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>&Element;</mo><mi>&Omega;</mi></mrow></munder><mo>&lsqb;</mo><mi>w</mi><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><mi>D</mi><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000784431070000015.GIF" wi="1600" he="198" /></maths>在上式中,Ω是图像f中一个局部窗口区域,其尺寸为S×S,S通常为奇数,对于这么一个窗口内的像素,中心像素为(i,j)即Ω<sub>center</sub>,而(x,y)为该区域内的其余任意像素位置,而上式中的w表示空间权重,定义为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>D</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>i</mi><mi>j</mi><mo>,</mo><mi>x</mi><mi>y</mi><mo>)</mo></mrow><mo>/</mo><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></msup></mrow>]]></math><img file="FDA0000784431070000021.GIF" wi="560" he="149" /></maths>其中D<sub>s</sub>(ij,xy)表示像素(x,y)与像素(i,j)的空间距离(用像素间隔表示),而σ<sup>2</sup>是伸缩系数;而D(f<sub>ij</sub>,f<sub>xy</sub>)为像素间的灰度差异,其定义为D(f<sub>xy</sub>,f<sub>ij</sub>)=|f<sub>xy</sub>‑f<sub>ij</sub>|;随着窗口的移动,(i,j)也逐个移动,于是获得了局部强度差异图Dif,其与原始图像g等尺寸;对Dif进行归一化;最后定义图像清晰度为局部强度差异图Dif的均值:Clearness=mean(Dif)。
地址 310018 浙江省杭州市下沙高教园区2号大街