发明名称 基于离散余弦变换和稀疏表示的无参考图像质量评价方法
摘要 本发明公开了一种基于离散余弦变换和稀疏表示的无参考图像质量评价方法,主要解决现有技术对无参考图像质量评价不准确的问题。其实现步骤是:输入一副灰度图像,对其进行离散余弦变换并提取自然场景统计特征;提取一系列不同失真类型和不同内容的图像的自然场景统计特征,结合平均主观差异分数构建原始特征字典;对原始特征字典进行聚类,根据测试图像特征与原始特征字典中各类的近似程度自适应选择原子组成稀疏表示字典;利用稀疏表示在特征空间张成测试图像特征并计算稀疏表示系数,结合稀疏表示字典中的主观评价值进行线性加权求和,得到图像质量测度。本发明与主观评价结果具有较好的一致性,适用于对各种失真类型图像的质量评价。
申请公布号 CN104376565B 申请公布日期 2017.03.29
申请号 CN201410695579.8 申请日期 2014.11.26
申请人 西安电子科技大学 发明人 张小华;焦李成;温阳;王爽;钟桦;田小林;朱虎明
分类号 G06T7/00(2017.01)I 主分类号 G06T7/00(2017.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于离散余弦变换和稀疏表示的无参考图像质量评价方法,包括如下步骤:(1)读入一幅灰度图像I,对其进行离散余弦变换,提取一系列与主观感知相关联的自然场景统计特征f;(2)组建训练图像的原始特征字典;(2a)重复步骤(1),提取n幅训练图像的自然场景统计特征,并构成特征矩阵F,F=[f<sub>1</sub>,f<sub>2</sub>,…,f<sub>i</sub>,…,f<sub>n</sub>],其中f<sub>i</sub>为第i幅训练图像的自然场景统计特征,i=1,2,…,n;(2b)整合n幅训练图像的平均主观差异分数,并构成质量向量M,M=[m<sub>1</sub>,m<sub>2</sub>,...,m<sub>i</sub>,…,m<sub>n</sub>],其中m<sub>i</sub>为第i幅训练图像的平均主观差异分数,i=1,2,…,n;(2c)将质量向量M和特征矩阵F进行对应结合,构建原始特征字典D:<img file="FDA0001193411560000011.GIF" wi="670" he="143" />其中<img file="FDA0001193411560000012.GIF" wi="475" he="143" />(2d)用K‑means算法将原始特征字典D聚为H类,H类中第k类聚类中心为C<sub>k</sub>,其中<img file="FDA0001193411560000013.GIF" wi="267" he="143" />Fc<sub>k</sub>为特征聚类中心,Mc<sub>k</sub>为质量聚类中心,k=1,2…H;(3)读入一幅测试图像I',重复步骤(1),提取测试图像I'的自然场景统计特征f';(4)计算测试图像I'的自然场景统计特征f'与原始特征字典D中第k类的特征聚类中心Fc<sub>k</sub>的欧式距离dis<sub>k</sub>,得到测试图像I'与原始特征字典D中第k类的近似程度P<sub>k</sub>,其中<img file="FDA0001193411560000014.GIF" wi="699" he="151" />(5)根据近似程度P<sub>k</sub>,从原始特征字典D的第k类中选择与特征聚类中心Fc<sub>k</sub>的欧式距离最小的前N<sub>k</sub>个样本作为测试图像I'的稀疏表示原子<img file="FDA0001193411560000015.GIF" wi="99" he="78" />其中,p=1,2,…,N<sub>k</sub>,k=1,2,…,H H类原子共同构成针对测试图像I'的稀疏表示字典D':<img file="FDA0001193411560000016.GIF" wi="734" he="143" />其中<img file="FDA0001193411560000017.GIF" wi="875" he="158" /><img file="FDA0001193411560000018.GIF" wi="60" he="71" />为从原始特征字典D中选择的第k类中的第p个样本的自然场景统计特征,<img file="FDA0001193411560000019.GIF" wi="62" he="70" />为第k类中的第p个样本的平均主观差异分数,p=1,2,…,N<sub>k</sub>,N<sub>k</sub>=δ*P<sub>k</sub>,k=1,2,…,H,δ为正常数;(6)根据稀疏表示的方法求解测试图像I'的自然场景统计特征f'在特征矩阵F'下的稀疏表示系数:<img file="FDA0001193411560000021.GIF" wi="692" he="109" />其中argmin表示将使目标函数λ·||α||<sub>1</sub>+||f'‑D′α||<sub>2</sub>最小的α=[α<sup>1</sup>,α<sup>2</sup>,…,α<sup>k</sup>,…,α<sup>H</sup>]<sup>T</sup>赋值给α<sup>*</sup>,<img file="DA00011934115633804.GIF" wi="683" he="71" /><img file="FDA0001193411560000022.GIF" wi="591" he="71" />R<sup>L</sup>表示L维实数空间,<img file="FDA0001193411560000023.GIF" wi="513" he="135" />表示向量β=[β<sub>1</sub>,β<sub>2</sub>,…,β<sub>l</sub>,…,β<sub>L</sub>]<sup>T</sup>的1范数,<img file="FDA0001193411560000024.GIF" wi="289" he="140" />表示向量β的2范数,λ是用于平衡保真项和正则项的正常数,T代表转置操作;(7)根据所构建的稀疏表示计算测试图像I'的质量:<img file="FDA0001193411560000025.GIF" wi="637" he="270" />其中<img file="FDA0001193411560000026.GIF" wi="62" he="78" />为稀疏表示字典D'中第k类第p个原子<img file="FDA0001193411560000027.GIF" wi="56" he="78" />的平均主观差异分数,α<sub>k,p</sub>表示测试图像I'的自然场景统计特征f'在特征矩阵F'中的第k类第p个原子<img file="FDA0001193411560000028.GIF" wi="59" he="75" />的稀疏表示系数,Q为测试图像I'的最终质量测度。
地址 710071 陕西省西安市太白南路2号