[发明专利]一种基于梯度关联性的无参考图像质量评价方法有效

专利信息
申请号: 201410284237.7 申请日: 2014-06-23
公开(公告)号: CN104023230A 公开(公告)日: 2014-09-03
发明(设计)人: 刘利雄;化毅;赵清杰;黄华 申请(专利权)人: 北京理工大学
主分类号: H04N17/00 分类号: H04N17/00
代理公司: 暂无信息 代理人: 暂无信息
地址: 100081 北京市*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明提出了一种基于梯度关联性的无参考图像质量评价方法,属于计算机图像分析领域。本方法首先求出图像梯度当中易受失真影响的三种子性质,分别为图像梯度幅值性质,图像梯度方向变化性质和图像梯度幅值变化性质。将这三种性质进行M×M的图像分块,在每一个图像块中求出上述三种性质的统计方差,并使用整幅图像中所有的图像块的统计方差的平均数作为其图像特征,最后使用支持向量机与图像质量评价中两步框架相结合的方法求得图像质量。本发明方法拥有时间复杂度小、主观一致性高、图像特征维数低通用性好等优点,可应用与小型计算设备或与图像质量相关的应用中,具有良好的实用价值。
搜索关键词: 一种 基于 梯度 关联性 参考 图像 质量 评价 方法
【主权项】:
一种基于梯度关联性的无参考图像质量评价方法,其特征在于包括以下步骤:步骤一、对输入的失真图像进行特征提取;首先,对每一幅图像求取其在图像梯度方面的三种不同的子性质,即梯度幅值性质GM、梯度方向变化性质CO和梯度幅值变化性质CM,其中,三种梯度性质分别定义如下:<mrow><mi>GM</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>Gx</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>Gy</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>CO(i,j)=orientation(i,j)‑orientationavg(i,j) (2)<mrow><mi>CM</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>Gx</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Gx</mi><mi>avg</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>Gy</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Gy</mi><mi>avg</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>其中,<mrow><mi>orientation</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>arc </mi><mi>tan</mi><mrow><mo>(</mo><mfrac><mrow><mi>Gy</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mi>Gx</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><msub><mi>orientation</mi><mi>avg</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>arc </mi><mi>tan</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>Gy</mi><mi>avg</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>Gx</mi><mi>avg</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><msub><mi>Gx</mi><mi>avg</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mi>i</mi><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mi>j</mi><mi>N</mi></munderover><mi>Gx</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>N</mi></mrow></mfrac></mrow><mrow><msub><mi>Gy</mi><mi>avg</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mi>i</mi><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mi>j</mi><mi>N</mi></munderover><mi>Gy</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mi>M</mi><mo>&times;</mo><mi>N</mi></mrow></mfrac></mrow>其中,orientation(i,j)表示图像梯度的方向,Gx和Gy为离散的数字图像在X、Y两个正交方向的导数,M、N是为了描述区域内的变化所设定的窗口的大小;并且在求取梯度Gx和Gy中,使用sobel梯度计算算子,并且将其由一组正交方向扩展到两组;其中包含0度和90度的一组正交方向,以及一组‑45度和45度的正交方向;正交方向0度和90度的sobel算子为:<mrow><msup><mi>S</mi><mi>y</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>2</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow><mrow><msup><mi>S</mi><mi>x</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mo>-</mo><mn>2</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>正交方向‑45度和45度的sobel算子为:<mrow><msup><mi>S</mi><mi>y</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow><mrow><msup><mi>S</mi><mi>x</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>2</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>其中上标为y的算子为计算Gy的算子,其中上标为x的算子为计算Gx的算子;然后,将上述所求得的图像梯度性质图像进行分块处理,在每一块梯度性质上求取统计方差,在此基础上将该方差以log函数规范后,以计算平均数的形式将每一块梯度所求取的统计方差融合为整体图像梯度性质的方差,具体如下:将每一种图像梯度性质分割为128×128大小的图像梯度性质图像块,并将每一个图像梯度性质图像块求取其统计方差,即第n块图像块的方差dn;计算公式如式4:dn=∑(h(x)‑E(h(x)))2 (4)其中,h(x)=pdf(θ)其中,θ为每种性质图像中的参数,pdf为该参数的统计分布,h(x)表示对θ进行量化统计后的统计概率表示,E(h(x))表示h(x)的期望;接着将每个求出的方差dn以log函数规范化,并使用平均聚集的方法融合图像梯度性质的图像块所求出的统计方差作为最终的特征f,如式5;<mrow><mi>f</mi><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mo>&times;</mo><mi>log</mi><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>d</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>最后,将输入的失真图像进行降采样,变为第二尺度的图像,并重复上述过程,最终得到一组12维特征向量Feature:Feature=[fGM,fCO,fCM×2orientation×2scale] (6)步骤二、进行特征映射,将经步骤一得到的12维特征向量作为最终的图像质量特征,据此建立图像特征和图像分数之间的映射关系;首先,将图像库分为训练集和测试集;所述测试集用来建立图像特征与图像质量之间的映射关系,测试集用来测试建立的映射关系的机能;使用支持向量机的方法,将测试集中图像的图像特征进行失真分类模型和对应的每个失真分类的质量评价模型的训练;然后,基于无参考图像质量评价中两步框架的流程,在测试集中进行分数的预测测试,即使用失真分类模型将受测图像的失真进行分类,再在失真分类当中使用质量评价模型来预测受测图像质量,从而得到受测图像质量分数,进而可利用现有算法性能标准对其进行评估。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京理工大学,未经北京理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410284237.7/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top