CN108053396A - 一种多失真图像质量的无参考评价方法 - Google Patents

一种多失真图像质量的无参考评价方法 Download PDF

Info

Publication number
CN108053396A
CN108053396A CN201711346145.7A CN201711346145A CN108053396A CN 108053396 A CN108053396 A CN 108053396A CN 201711346145 A CN201711346145 A CN 201711346145A CN 108053396 A CN108053396 A CN 108053396A
Authority
CN
China
Prior art keywords
mrow
msub
image
mfrac
msubsup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201711346145.7A
Other languages
English (en)
Other versions
CN108053396B (zh
Inventor
周玉
李雷达
卢兆林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201711346145.7A priority Critical patent/CN108053396B/zh
Publication of CN108053396A publication Critical patent/CN108053396A/zh
Application granted granted Critical
Publication of CN108053396B publication Critical patent/CN108053396B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/90—Dynamic range modification of images or parts thereof
    • G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20021—Dividing image into blocks, subimages or windows
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30168—Image quality inspection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提出一种多失真图像质量的无参考评价方法，该方法是一种基于双阶结构失真和非局部统计的评价方法，包括步骤：首先对多失真图像进行下采样得到多个分辨率的图像。提取每个分辨率图像的双阶结构特征用来评价图像中的结构失真。所有的结构特征构成了第一组质量评价特征。然后，提取每个分辨率图像的非局部相似性统计特征，将这些特征作为第二组质量评价特征。最后，将图像的两组特征作为输入，利用随机森林训练出多失真图像的质量评价模型。本发明的性能明显优于现有的无参考图像质量评价方法，且具有很好的跨库性能，可扩展性强。

Description

一种多失真图像质量的无参考评价方法
技术领域
本发明涉及图像质量评价领域，尤其是一种多失真图像质量的无参考评价方法。
背景技术
图像质量评价在图像处理领域以及实际应用中都有非常广泛的应用[1]。目前，已经有大量的图像质量评价方法。图像质量评价方法可以分为主观图像质量评价方法和客观图像质量评价方法。主观质量评价方法是指通过人眼来判断图像的质量。这种方法精确度高，但是费时、费力，成本高，尤其是在有海量图像需要进行质量评价时，具有较大的缺陷。客观图像质量评价方法是指通过设计数学模型来评价图像的质量分数。客观质量评价方法克服了主观图像质量评价方法的缺陷。客观图像质量评价方法可以分为全参考、半参考和无参考三种方式。全参考方法必须在获得整个原始无失真图像的情况下才能对失真图像进行质量评价。半参考的方法仅需要部分原始无失真图像的信息即可对失真图像进行质量评价。在实际生活中，一般不可获得原始无失真图像，所以全参考和半参考方法均具有较大的局限性。而无参考方法可以在没有原始无失真图像的情况下，即可对失真图像进行质量评价。因此，无参考方法具有最为广泛的应用和最大的可行性。另外，在所有现有的质量评价方法中，大部分方法是针对包含单一类型失真的图像进行质量评价。例如，针对模糊图像的质量评价方法[2]以及针对JPEG压缩图像的质量评价方法[3]。但是，现实中，我们看到的图像都要经过获取、压缩、传输和存储等过程。这些过程会在图像中引入多种类型的失真[4]，例如图像获取的过程会带来模糊、压缩过程会带来块效应失真，还有一些传输失真等等。因此，现有的为单一失真设计的图像质量评价方法的性能以及应用范围具有一定的限制。
目前，只有三种质量评价方法是专门为多失真图像设计的。方法[5]中提出了一种包含五个步骤的多失真图像质量评价方法。对于一幅多失真图像，首先利用噪声估计方法判断图像中是否包含噪声。如果这幅图像中含有噪声，则先对图像进行去噪。然后，分别利用评价图像清晰度和块效应失真的方法对去噪后的图像进行评价，从而得到图像的清晰度分数和块效应分数；如果图像中没有噪声，则不需要去噪，而是在失真图像上直接求得清晰度分数和块效应分数。最后，将两个分数结合，从而得到最终的多失真图像质量分数。方法[6]在方法[5]的基础上提出了一种包含六个步骤的多失真图像质量评价方法。首先，利用噪声估计方法计算图像中的噪声分数作为分数1。然后，对图像进行去噪。再用评价图像清晰度和块效应失真的方法对去噪后的图像进行评价，计算出图像的清晰度分数和块效应分数。这两个分数被分别作为分数2和分数3。再计算图像中的模糊失真和块效应共同作用下的质量分数，作为分数4。最后，将得到的四个质量分数进行整合得到最终的多失真图像质量分数。方法[4]提出了一种基于结构失真度量的多失真图像质量评价方法。首先，利用低通滤波和下采样的方法得到多个分辨率的图像。然后计算图像的梯度图像，并在梯度图像上计算局部二值化图像，以此来提取图像中的结构信息。将梯度图像的幅值图像作为二值化图像的权值图像对二值化图像进行加权。最后，计算加权后的二值化图像的直方图，并将直方图的求和结果作为失真图像的质量分数。
尽管这三个为多失真图像设计的质量评价方法相比于为单一失真图像设计的方法在性能上有一定的提高，但是，他们仍然存在一些缺陷。尤其对于方法[5][6]，他们只考虑了多失真图像中三种类型的失真，即噪声，模糊和块效应。但是，现实中图像的失真类型不只这三种失真类型，所以，这两种方法的性能和应用范围具有一定的局限性。对于方法[4],它仅仅考虑图像中主要的结构失真。图像中细小的结构，例如图像的纹理，已经被证明同样会影响人眼对于图像的感知[7]。所以，该方法仍然具有提升空间。
基于以上分析，亟需设计一种更加有效的无参考的多失真图像质量评价方法。
[1]W.S.Lin,C.-C.J.Kuo,“Perceptual visual quality metrics:A survey,”J.Vis.Commun.Image Represent.,vol.22,no.4,pp.297-312,May 2011.
[2]Q.B.Sang,H.X.Qi,X.J.Wu,C.F.Li,and A.C.Bovik“No-reference imageblur index based on singular value curve”J.Vis.Commun.Image Represent.,vol.25,no.7,pp.
1625-1630,Oct.2014.
[3]L.D.Li,Y.Zhou,J.J.Wu,W.S.Lin,andH.L.Li“GridSAR:Grid strength andregularity for robust evaluation ofblocking artifacts in JPEG images”J.Vis.Commun.Image Represent.,vol.30,pp.153-163,Jul.2015.
[4]Q.H.Li,W.S.Lin,andY.M.Fang“No-reference quality assessment formultiply distorted imaegs in gradient domain,”IEEE Signal Process.Lett.,vol.23,pp.541-545,Apr.2016.
[5]K.Gu,G.T.Zhai,M.Liu,X.K.Yang,W.J.Zhang,X.H.Sun,W.H.Chen,and Y.Zuo“FISBLIM:A five-step blindmetric for quality assessment ofmultiply distortedimages,”IEEE Workshop on Signal Processing Systems,pp.241-246,Oct.2013.
[6]K.Gu,G.T.Zhai,X.K.Yang,and W.J.Zhang“Hybridno-reference qualitymetric for singly and multiply distorted images,”IEEE Trans.Broadcasting,vol.60,pp.555-567,Sep.2014.
[7]I.Bellamine,and H.Tairi“Motion detection using color structure-texture image decomposition,”Intelligent Systems andcomputerVision(ISCV),2015.
发明内容
发明目的：为克服现有技术的不足，提供一种更加有效的多失真图像质量评价方法，本发明提出一种多失真图像质量的无参考评价方法。
技术方案：本法明提出的技术方案为：
一种多失真图像质量的无参考评价方法，包括步骤：
(1)获取M副原始多失真图像，对每一幅原始多失真图像执行步骤(1-1)至(1-3)
(1-1)记任意一幅原始多失真图像为I0，对图像I0分别进行n次下采样，得到n幅下采样图像，并记Ii为第i个下采样图像，i∈[1,2,…,n]；将图像I0至In归入图像I0的样本集合；
(1-2)对样本集合中每一幅图像进行双阶结构特征提取，获取每一幅图像的一阶结构失真特征和二阶结构失真特征；
(1-3)对样本集合中每一幅图像进行非局部相似性统计特征提取，获取每一幅图像的两个非局部相似性统计特征；
(1-4)步骤(1-2)和(1-3)提取出的所有特征参数共同构成图像I0的特征集合；
(2)将M幅原始多失真图像的特征集合合并形成一个训练参数集合，采用随机森林法从训练参数集合的所有特征参数中学习出一个多失真图像质量评价模型；
(3)通过多失真图像质量评价模型评价待测多失真图像。
进一步的，所述一阶结构特征包括结构信息量特征和结构能量特征；其中，结构信息量特征的计算步骤包括：
(2-1)求出图像Ii的梯度图像gi，i＝0,1,…,n：
▽hi(x,y)＝[-1 0 1]*Ii
▽vi(x,y)＝[-1 0 1]T*Ii
式中，gi(x,y)表示梯度图像gi中像素点(x,y)处的像素值，▽hi(x,y)和▽vi(x,y)分别表示Ii的水平方向梯度和垂直方向梯度；*和T分别表示卷积操作符和转置操作符；
(2-2)对于梯度图像gi中的每一个像素(x,y)，首先确定一个以(x,y)为中心，且大小为W×W的窗口，W表示窗口宽度；计算出窗口中除了像素(x,y)以外的其他所有像素点的平均像素值Ai(x,y)：
(2-3)定义一个函数fi，fi中的元素表示为fi(p,q)，fi(p,q)的值为：满足gi(x,y)＝p且Ai(x,y)＝q的像素点的个数；计算fi(p,q)的概率为：
式中，a1、a2分别为I0的长和宽；
(2-4)计算梯度图像gi的二维熵为：
计算出的二维熵为Ii的结构信息量特征；
结构能量特征的计算步骤包括：
(2-5)将梯度图像gi分为b×b大小的不重叠的图像块，对每个图像块进行奇异值分解，任意一个图像块Bz的奇异值分解结果为：
Bz＝UzSzVz T,z＝1,2,...,R
式中，R表示梯度图像gi中图像块总数，Uz和Vz均为大小为b×b的方阵，并分别满足：Uz TUz＝E和Vz TVz＝E；式中，E为单位矩阵；Sz为大小为b×b的方阵，其对角线上的数据为奇异值，其他位置的数据为0；Bz的平均奇异值为Sz中对角线上数据的平均值，记为Yz；
(2-6)对Yz进行归一化，得到：
式中，Yz′表示Yz归一化值，表示Bz的方差；
(2-7)计算gi的结构能量特征为:
进一步的，所述二阶结构特征提取步骤包括：
(3-1)计算图像Ii的二值化图像LBPi：
式中，表示LBPi中像素点(x,y)处的像素值，s和r分别表示Ii中像素点(x,y)的邻域采样点个数和邻域采样半径；表示Ii(x,y)邻域采样点w的像素值；H()为中间函数；
(3-2)提取二值化图像LBPi的结构信息量特征Qi3和结构能量特征Qi4，Qi3和Qi4即为图像Ii的二阶结构特征。
进一步的，所述非局部相似性统计特征提取的方法为：
(4-1)将Ii分为大小为b×b的相互重叠的图像块，定义第k个图像块为Bk，确定一个以Bk为中心的大小为O×O的搜索窗L；L中和Bk相同大小的搜索块Bl以步长为1的速度按照从左向右从上到下的顺序依次运动，l为搜索块的索引；
(4-2)定义Bl和Bk的相似性为两个块之间的欧氏距离：
(4-3)计算出L中所有搜索块与Bk的相似性，选取相似性最高的v个搜索块；
(4-4)将选取出的v个搜索块和Bk组成一个三维数组，然后对三维数组进行正交三维变换，得到三维变换的小波系数，记小波系数集合为Zit，t＝[1,2,…,E],E表示图像块的总数；
(4-5)对每个小波系数集合Zit执行步骤S1至S4：
S1：计算Zit的一阶绝对矩J1：
式中，z为随机变量，θ表示伽马函数，γit表示Zit的形状参数，y2为一中间参数，
式中，σit表示Zit中所有系数的标准差，求得所有系数的方差为σit 2；
令得到将代入J1的计算公式，得到：
S2：计算Zit的二阶矩为：
J2＝σit 2
S3：令计算出γit：
式中，zj表示集合Zit中第j个小波系数，h表示Zit中小波系数的个数；
(4-6)计算图像Ii中所有图像块的σit 2和γit，根据计算出的σit 2和γit计算：
Qi5和Qi6即为图像Ii的两个非局部相似性统计特征。
有益效果：与现有技术相比，本发明具有以下优势：
1、与现有的无参考图像质量评价方法和多失真图像质量评价方法该评价方法相比，本发明的性能有明显提高；
2、与现有的多失真方法和通用无参质量评价方法相比，本发明的跨库性能最好、可扩展性最强。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
图1所示为本发明的原理流程图，本发明主要包括以下步骤：
1、获取M副原始多失真图像，对每一幅原始多失真图像执行步骤(1-1)至(1-3)
(1-1)下采样：记任意一幅原始多失真图像为I0，对图像I0分别进行n次下采样，得到n幅下采样图像，并记Ii为第i个下采样图像，i∈[1,2,…,n]；将图像I0至In归入图像I0的样本集合；
(1-2)对样本集合中每一幅图像进行双阶结构特征提取，获取每一幅图像的一阶结构失真特征和二阶结构失真特征；
(1-3)对样本集合中每一幅图像进行非局部相似性统计特征提取，获取每一幅图像的两个非局部相似性统计特征；
(1-4)步骤(1-2)和(1-3)提取出的所有特征参数共同构成图像I0的特征集合；
(3)将M副原始多失真图像的特征集合合并形成一个训练参数集合，采用随机森林法从训练参数集合的所有特征参数中学习出一个多失真图像质量评价模型；
(4)通过多失真图像质量评价模型评价待测多失真图像。
下面通过一个具体实施例进一步说明本发明：
步骤一：对于一幅大小为a1×a2的多失真图像I0，通过MATLAB中的内置函数“Imresize”进行4次下采样，得到4幅下采样图像，每次下采样的参数为0.5，即每次下采样后，图像的长和宽分别减小一半。所以，加上原始分辨率的图像，一共可以获得5个不同分辨率的图像，记为I0,I1,I2,I3,I4。
步骤二：一阶结构失真的特征提取：
利用梯度来提取图像的一阶结构：
(1)求每一个多分辨率图像Ii的梯度图像gi：
▽hi(x,y)＝[-1 0 1]*Ii
▽vi(x,y)＝[-1 0 1]T*Ii
式中，gi(x,y)表示梯度图像gi中像素点(x,y)处的像素值，▽hi(x,y)和▽vi(x,y)分别表示Ii的水平方向梯度和垂直方向梯度，*和T分别表示卷积操作符和转置操作符。
然后，从结构信息量和结构能量两个方面来表征图像的结构失真。其中，利用梯度图像的二维熵来计算图像的结构信息量；利用梯度图像的能量来计算图像的结构能量。
梯度图像二维熵：对于梯度图像gi中的每一个像素(x,y)，首先确定一个以(x,y)为中心，且大小为W×W的窗口，W表示窗口宽度；计算出窗口中除了像素(x,y)以外的其他所有像素点的平均像素值Ai(x,y)：
定义一个函数fi，fi中的元素表示为fi(p,q)，fi(p,q)的值为：满足gi(x,y)＝p且Ai(x,y)＝q的像素点的个数；则fi(p,q)的概率为：
式中，a1、a2分别为I0的长和宽；
计算梯度图像gi的二维熵为：
计算出的二维熵为Ii的结构信息量特征。
本实施例中，W＝3。对于每个分辨率图像可以得到一个二维熵特征，因此，5个分辨率图像一共可以得到5个结构信息量特征。
梯度图像gi的能量：首先将梯度图像gi分为b×b大小的不重叠的图像块。对每一个图像块Bz进行奇异值分解，则Bz可以表示为：
Bz＝UzSzVz T,z＝1,2,...,R
R表示梯度图像gi中图像块总数，Uz和Vz均为大小为b×b的方阵，并分别满足：Uz TUz＝E和Vz TVz＝E；式中，E为单位矩阵；Sz为大小为b×b的方阵，其对角线上的数据为奇异值，其他位置的数据为0；Bz的平均奇异值为Sz中对角线上数据的平均值，记为Yz。
为了减小图像内容对结构能量的影响，利用Bz的方差对Yz进行归一化：
式中，Yz′表示Yz归一化值，σz 2表示块Bz的方差。
对图像中每个块Bz归一化后的奇异值求平均，并将结果作为Ii的结构能量特征:
对于每个分辨率图像可以得到1个结构能量特征，一共5个分辨率图像，所以共得到5个结构能量特征。
综上，一共得到10个一阶结构特征，包括5个结构信息量特征和5个结构能量特征。
(2)二阶结构失真的特征提取：
二阶结构代表图像中细小的结构，例如纹理。局部二值化是一种局部描述法，能够有效的提取图像的纹理。局部二值化目前已经广泛应用于纹理分类中。局部二值化图像的计算方法为：
计算每一个多分辨率图像Ii的二值化图像LBPi：
式中，表示LBPi中像素点(x,y)处的像素值，s和r分别表示Ii中像素点(x,y)的邻域采样点个数和邻域采样半径；表示Ii(x,y)邻域采样点w的像素值；H()为中间函数；
和一阶结构上提取特征的方法类似。在每个分辨率图像的二值化图像上提取1个结构信息量特征和1个结构能量特征，分别记为Qi3和Qi4。将它们作为描述多失真图像二阶结构的特征。因为有5个分辨率图像，所以可以计算出10个二阶结构特征。
步骤三：非局部相似性统计特征提取：
非局部相似性统计是指自然图像中的一个块和它邻域内的相似性块的正交变换系数满足广义高斯分布。具体计算方式为：
将每个分辨率图像Ii分为大小为b×b的相互重叠的图像块，以第k个块Bk为例，为了寻找Bk邻域内的相似性块，首先确定一个以Bk为中心的大小为O×O的搜索窗L。L中和Bk相同大小的搜索块Bl以步长为1的速度按照从左向右从上到下的顺序依次运动，其中，l为搜索块的索引。Bl和Bk的相似性定义为两个块之间的欧氏距离：
Dl,k的数值越小，代表块Bl和Bk的相似性越大。按照这种方式，计算出L中所有搜索块和Bk的相似性。然后，选取其中最为相似的v个块。将这v个块和Bk组成一个三维数组。然后对三维数组进行正交三维变换，从而得到三维变换的小波系数。这些小波系数满足广义高斯分布。将这些小波系数构成的集合记为Zit，t＝[1,2,…,E],E表示图像块的总的个数。
对每个小波系数集合Zit执行步骤S1至S3：
S1：计算Zit的一阶绝对矩：
式中，z为随机变量，θ表示伽马函数，γit表示Zit的形状参数，y2为一中间参数，
式中，σit表示Zit中所有系数的标准差，进一步可以求得所有系数的方差σit 2。
令得到将代入J1的计算公式，得到：
S2：计算Zit的二阶矩为：
J2＝σit 2
S3：令 计算出γit：
式中，zj表示集合Zit中第j个小波系数，h表示Zit中小波系数的个数。
2)对于每个系数集合Zit，t＝[1,2,…,E]，可以得到σit 2，γit两个特征。计算图像Ii中所有图像块的σit 2和γit，根据计算出的σit 2和γit计算：
Qi5和Qi6即为图像Ii的两个非局部相似性统计特征。
每个分辨率图像得到两个非局部相似性统计特征，一共是5个分辨率图像，所以可以得到10个非局部相似性统计特征。加上前面计算出的20个结构失真特征，图像I0共有30个特征。
步骤四：质量评价模型的建立，具体为：
选取多张多失真图像，将这些图像的所有特征作为输入，利用随机森林的方法训练出质量评价模型。然后利用该质量评价模型去测试需要测试的多失真图像的质量分数。
实验结果以及性能：
本实施例的实验部分是在3个公开的多失真图像数据集上进行的。即MDID2013库、MLIVE库和MDID库。MDID2013库提供了324幅多失真图像，这些图像包含三种类型的失真，即高斯模糊、JPEG压缩和白噪声。MLIVE库提供了450幅多失真图像，这些图像包含两种组合的失真类型，即高斯模糊和JPEG压缩，以及高斯模糊和白噪声。MDID库提供了1600幅多失真图像。这些图像中包含五种传统类型的失真，即高斯模糊、JPEG压缩、白噪声、JPEG2000压缩和对比度失真。在实验的操作过程中，为了测试本发明方法在待测试的图像数据集上的应用性能，首先将数据集中的图像随机分为80％和20％，其中，80％的图像用作质量评价模型的建立，20％的图像用作模型的性能测试。为了避免偶然性的发生，这个过程循环1000次，4个性能指标的中值被作为最终的性能参数。四个性能指标参数为：1)相关系数(Pearsonlinear correlation coefficient,PLCC)为非线性回归后的线性相关系数；2)均方误差根(Root mean squared error，RMSE)为非线性回归后的标准差；3)Kendall等级相关系数(Kendall’s Rank Correlation Coefficient，KRCC)；4)Spearman相关系数(Spearmanrank order correlation coefficient，SRCC)。PLCC/SRCC/KRCC数值越大，RMSE数值越小，说明算法性能越好。PLCC和RMSE用来评价方法的预测准确性，而SRCC和KRCC用来评价方法的预测单调性。
首先，我们将本发明的性能和现有的多失真图像质量评价方法在三个数据库中进行比较，比较结果如表1所示。为了方便观看，最好的性能已经加粗显示。
表1本发明方法和现有多失真图像质量评价算法的性能对比
由上表，我们可以获得两个方面的信息：(1)对于MDID2013库，本发明的方法和现有的所有多失真图像质量评价方法相比具有明显的优势，即PLCC/SRCC/KRCC的数值明显高于所有方法，而RMSE最低。(2)对于MLIVE库和MDID库，在现有的三种方法中，方法[4]性能最好，但是它的性能仍然低于本发明的性能。综上，本发明方法的性能在三个库上均具有最好的性能。
为了进一步验证本发明的性能，我们将本发明和通用无参考图像质量评价方法进行比较，比较结果如表2所示。通用无参考图像质量评价算法是指无需知道图像的失真类型，即可对图像进行质量评价的算法。
表2本发明方法和通用无参考图像质量评价算法的性能比较表
由表2的数据可知，本发明的性能明显优于通用图像质量评价算法。具体表现为最好的预测准确性和单调性。
对于基于训练的质量评价算法，跨库性能，即可扩展性，是这些方法的重要评价指标。跨库实验，是指用一个图像库中的所有图像提取的特征训练出质量评价模型，然后利用该模型来预测另一个数据库中的所有图像的质量。基于此，对所有基于训练的质量评价方法进行可扩展性的验证实验，实验结果如表3所示。表3中，最好的两种性能被加粗显示。
表3基于训练的算法的跨库性能对比表
由表3的实验结果可知，本发明在所有跨库方式中的性能均能够位列前二。尤其是在MLIVE库上进行训练或者是在MDID2013库进行训练时，本发明能够获得所有方法中最好的跨库性能。在MDID库进行训练，MLIVE库进行测试时，本发明也具有最好的跨库性能。在MDID库训练，MDID2013库进行测试时，本发明的跨库性能略低于方法[4]。但是方法[4]在其他所有跨库方式下，表现的性能均比本发明的方法差。总体来说，本发明的方法具有最强的可扩展性。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.一种多失真图像质量的无参考评价方法，其特征在于，包括步骤：
(1)获取M副原始多失真图像，对每一幅原始多失真图像执行步骤(1-1)至(1-3)
(1-1)记任意一幅原始多失真图像为I0，对图像I0分别进行n次下采样，得到n幅下采样图像，并记Ii为第i个下采样图像，i∈[1,2,…,n]；将图像I0至In归入图像I0的样本集合；
(1-2)对样本集合中每一幅图像进行双阶结构特征提取，获取每一幅图像的一阶结构失真特征和二阶结构失真特征；
(1-3)对样本集合中每一幅图像进行非局部相似性统计特征提取，获取每一幅图像的两个非局部相似性统计特征；
(1-4)步骤(1-2)和(1-3)提取出的所有特征参数共同构成图像I0的特征集合；
(2)将M幅原始多失真图像的特征集合合并形成一个训练参数集合，采用随机森林法从训练参数集合的所有特征参数中学习出一个多失真图像质量评价模型；
(3)通过多失真图像质量评价模型评价待测多失真图像。
2.根据权利要求1所述的一种多失真图像质量的无参考评价方法，其特征在于，所述一阶结构特征包括结构信息量特征和结构能量特征；其中，结构信息量特征的计算步骤包括：
(2-1)求出图像Ii的梯度图像gi，i＝0,1,…,n：
<mrow> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <mo>&amp;dtri;</mo> <msub> <mi>h</mi> <mi>i</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mo>&amp;dtri;</mo> <msub> <mi>v</mi> <mi>i</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>n</mi> </mrow>
<mrow> <mo>&amp;dtri;</mo> <msub> <mi>h</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>&amp;lsqb;</mo> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> <mo>&amp;rsqb;</mo> <mo>*</mo> <msub> <mi>I</mi> <mi>i</mi> </msub> </mrow>
<mrow> <mo>&amp;dtri;</mo> <msub> <mi>v</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>*</mo> <msub> <mi>I</mi> <mi>i</mi> </msub> </mrow>
式中，gi(x,y)表示梯度图像gi中像素点(x,y)处的像素值，和分别表示Ii的水平方向梯度和垂直方向梯度；*和T分别表示卷积操作符和转置操作符；
(2-2)对于梯度图像gi中的每一个像素(x,y)，首先确定一个以(x,y)为中心，且大小为W×W的窗口，W表示窗口宽度；计算出窗口中除了像素(x,y)以外的其他所有像素点的平均像素值Ai(x,y)：
<mrow> <msub> <mi>A</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>l</mi> <mo>=</mo> <mi>y</mi> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <mi>W</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> </mrow> <mrow> <mi>y</mi> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mi>W</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> </mrow> </msubsup> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mi>x</mi> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <mi>W</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> </mrow> <mrow> <mi>x</mi> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mi>W</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> </mrow> </msubsup> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mi>l</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <msup> <mi>W</mi> <mn>2</mn> </msup> </mfrac> </mrow>
(2-3)定义一个函数fi，fi中的元素表示为fi(p,q)，fi(p,q)的值为：满足gi(x,y)＝p且Ai(x,y)＝q的像素点的个数；计算fi(p,q)的概率为：
<mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>,</mo> <mi>q</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>,</mo> <mi>q</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <mo>&amp;times;</mo> <msub> <mi>a</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mrow>
式中，a1、a2分别为I0的长和宽；
(2-4)计算梯度图像gi的二维熵为：
计算出的二维熵为Ii的结构信息量特征；
结构能量特征的计算步骤包括：
(2-5)将梯度图像gi分为b×b大小的不重叠的图像块，对每个图像块进行奇异值分解，任意一个图像块Bz的奇异值分解结果为：
Bz＝UzSzVz T,z＝1,2,...,R
式中，R表示梯度图像gi中图像块总数，Uz和Vz均为大小为b×b的方阵，并分别满足：Uz TUz＝E和Vz TVz＝E；式中，E为单位矩阵；Sz为大小为b×b的方阵，其对角线上的数据为奇异值，其他位置的数据为0；Bz的平均奇异值为Sz中对角线上数据的平均值，记为Yz；
(2-6)对Yz进行归一化，得到：
<mrow> <msubsup> <mi>Y</mi> <mi>z</mi> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <msub> <mi>Y</mi> <mi>z</mi> </msub> <mrow> <msup> <msub> <mi>&amp;sigma;</mi> <mi>z</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
式中，Y′z表示Yz归一化值，表示Bz的方差；
(2-7)计算gi的结构能量特征为:
<mrow> <msub> <mi>Q</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>z</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>R</mi> </msubsup> <msubsup> <mi>Y</mi> <mi>z</mi> <mo>&amp;prime;</mo> </msubsup> </mrow> <mi>R</mi> </mfrac> <mo>.</mo> </mrow>
3.根据权利要求2所述的一种多失真图像质量的无参考评价方法，其特征在于，所述二阶结构特征提取步骤包括：
(3-1)计算图像Ii的二值化图像LBPi：
<mrow> <msubsup> <mi>LBP</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>i</mi> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>w</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>s</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>H</mi> <mrow> <mo>(</mo> <msubsup> <mi>I</mi> <mi>i</mi> <mi>w</mi> </msubsup> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>I</mi> <mi>i</mi> </msub> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mn>2</mn> <mi>w</mi> </msup> </mrow>
<mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>a</mi> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>a</mi> <mo>&lt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中，表示LBPi中像素点(x,y)处的像素值，s和r分别表示Ii中像素点(x,y)的邻域采样点个数和邻域采样半径；表示Ii(x,y)邻域采样点w的像素值；H()为中间函数；
(3-2)提取二值化图像LBPi的结构信息量特征Qi3和结构能量特征Qi4，Qi3和Qi4即为图像Ii的二阶结构特征。
4.根据权利要求3所述的一种多失真图像质量的无参考评价方法，其特征在于，所述非局部相似性统计特征提取的方法为：
(4-1)将Ii分为大小为b×b的相互重叠的图像块，定义第k个图像块为Bk，确定一个以Bk为中心的大小为O×O的搜索窗L；L中和Bk相同大小的搜索块Bl以步长为1的速度按照从左向右从上到下的顺序依次运动，l为搜索块的索引；
(4-2)定义Bl和Bk的相似性为两个块之间的欧氏距离：
<mrow> <msub> <mi>D</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>y</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>b</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>x</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>b</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mi>l</mi> </msub> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>B</mi> <mi>k</mi> </msub> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
(4-3)计算出L中所有搜索块与Bk的相似性，选取相似性最高的v个搜索块；
(4-4)将选取出的v个搜索块和Bk组成一个三维数组，然后对三维数组进行正交三维变换，得到三维变换的小波系数，记小波系数集合为Zit，t＝[1,2,…,E],E表示图像块的总数；
(4-5)对每个小波系数集合Zit执行步骤S1至S4：
S1：计算Zit的一阶绝对矩J1：
<mrow> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </msubsup> <mo>|</mo> <mi>z</mi> <mo>|</mo> <mfrac> <mrow> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> <mrow> <mn>2</mn> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mo>|</mo> <msub> <mi>zy</mi> <mn>2</mn> </msub> <msup> <mo>|</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> </msup> </mrow> </msup> <mi>d</mi> <mi>z</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> <mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </msubsup> <mo>|</mo> <mi>z</mi> <mo>|</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mo>|</mo> <msub> <mi>zy</mi> <mn>2</mn> </msub> <msup> <mo>|</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> </msup> </mrow> </msup> <mi>d</mi> <mi>z</mi> </mrow>
式中，z为随机变量，θ表示伽马函数，γit表示Zit的形状参数，y2为一中间参数，
式中，σit表示Zit中所有系数的标准差，求得所有系数的方差为σit 2；
令得到将代入J1的计算公式，得到：
<mrow> <msub> <mi>J</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>&amp;sigma;</mi> <mfrac> <mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <msqrt> <mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </msqrt> </mfrac> </mrow>
S2：计算Zit的二阶矩为：
J2＝σit 2
S3：令计算出γit：
<mrow> <mfrac> <msubsup> <mi>J</mi> <mn>1</mn> <mn>2</mn> </msubsup> <msub> <mi>J</mi> <mn>2</mn> </msub> </mfrac> <mo>=</mo> <mfrac> <mrow> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mn>2</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>/</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>=</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mrow> <mfrac> <mn>1</mn> <mi>h</mi> </mfrac> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </msubsup> <mo>|</mo> <msub> <mi>z</mi> <mi>j</mi> </msub> <mo>|</mo> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mfrac> <mn>1</mn> <mi>h</mi> </mfrac> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>h</mi> </msubsup> <mo>|</mo> <msub> <mi>z</mi> <mi>j</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
式中，zj表示集合Zit中第j个小波系数，h表示Zit中小波系数的个数；
(4-6)计算图像Ii中所有图像块的σit 2和γit，根据计算出的σit 2和γit计算：
<mrow> <msub> <mi>Q</mi> <mrow> <mi>i</mi> <mn>5</mn> </mrow> </msub> <mo>=</mo> <mover> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mn>2</mn> </msubsup> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>E</mi> </munderover> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> <mn>2</mn> </msubsup> </mrow> <mi>E</mi> </mfrac> </mrow>
<mrow> <msub> <mi>Q</mi> <mrow> <mi>i</mi> <mn>6</mn> </mrow> </msub> <mo>=</mo> <mover> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>E</mi> </munderover> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>t</mi> </mrow> </msub> </mrow> <mi>E</mi> </mfrac> </mrow>
Qi5和Qi6即为图像Ii的两个非局部相似性统计特征。
CN201711346145.7A 2017-12-15 2017-12-15 一种多失真图像质量的无参考评价方法 Active CN108053396B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711346145.7A CN108053396B (zh) 2017-12-15 2017-12-15 一种多失真图像质量的无参考评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711346145.7A CN108053396B (zh) 2017-12-15 2017-12-15 一种多失真图像质量的无参考评价方法

Publications (2)

Publication Number Publication Date
CN108053396A true CN108053396A (zh) 2018-05-18
CN108053396B CN108053396B (zh) 2021-09-14

Family

ID=62133123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711346145.7A Active CN108053396B (zh) 2017-12-15 2017-12-15 一种多失真图像质量的无参考评价方法

Country Status (1)

Country Link
CN (1) CN108053396B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109035217A (zh) * 2018-07-06 2018-12-18 同济大学 一种基于图像二维熵的土体裂隙定量评价方法
CN109325550A (zh) * 2018-11-02 2019-02-12 武汉大学 基于图像熵的无参考图像质量评价方法
CN109410104A (zh) * 2018-11-30 2019-03-01 宁波永新光学股份有限公司 基于机器学习的显微全无线互动教学传输及辅助上课系统
CN109816646A (zh) * 2019-01-21 2019-05-28 武汉大学 一种基于退化决策逻辑的无参考图像质量评价方法
CN110415223A (zh) * 2019-07-17 2019-11-05 西安邮电大学 一种无参考的拼接图像质量评价方法及系统
CN110807757A (zh) * 2019-08-14 2020-02-18 腾讯科技（深圳）有限公司 基于人工智能的图像质量评估方法、装置及计算机设备
CN111612741A (zh) * 2020-04-22 2020-09-01 杭州电子科技大学 一种基于失真识别的精确无参考图像质量评价方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104616310A (zh) * 2015-02-12 2015-05-13 国家电网公司 一种图像质量的评估方法和装置
CN105528776A (zh) * 2015-08-07 2016-04-27 上海仙梦软件技术有限公司 针对jpeg图像格式的显著性细节保持的质量评价方法
CN107147906A (zh) * 2017-06-12 2017-09-08 中国矿业大学 一种虚拟视角合成视频质量的无参考评价方法
CN107343196A (zh) * 2017-07-18 2017-11-10 天津大学 一种混合失真无参考图像质量评价方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104616310A (zh) * 2015-02-12 2015-05-13 国家电网公司 一种图像质量的评估方法和装置
CN105528776A (zh) * 2015-08-07 2016-04-27 上海仙梦软件技术有限公司 针对jpeg图像格式的显著性细节保持的质量评价方法
CN107147906A (zh) * 2017-06-12 2017-09-08 中国矿业大学 一种虚拟视角合成视频质量的无参考评价方法
CN107343196A (zh) * 2017-07-18 2017-11-10 天津大学 一种混合失真无参考图像质量评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YU ZHOU等: ""Blind Quality Index for Multiply Distorted Images"", 《IEEE TRANSACTIONS ON MULTIMEDIA》 *
王志明: "无参考图像质量评价综述", 《自动化学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109035217A (zh) * 2018-07-06 2018-12-18 同济大学 一种基于图像二维熵的土体裂隙定量评价方法
CN109325550A (zh) * 2018-11-02 2019-02-12 武汉大学 基于图像熵的无参考图像质量评价方法
CN109410104A (zh) * 2018-11-30 2019-03-01 宁波永新光学股份有限公司 基于机器学习的显微全无线互动教学传输及辅助上课系统
CN109816646A (zh) * 2019-01-21 2019-05-28 武汉大学 一种基于退化决策逻辑的无参考图像质量评价方法
CN109816646B (zh) * 2019-01-21 2022-08-30 武汉大学 一种基于退化决策逻辑的无参考图像质量评价方法
CN110415223A (zh) * 2019-07-17 2019-11-05 西安邮电大学 一种无参考的拼接图像质量评价方法及系统
CN110807757A (zh) * 2019-08-14 2020-02-18 腾讯科技（深圳）有限公司 基于人工智能的图像质量评估方法、装置及计算机设备
CN110807757B (zh) * 2019-08-14 2023-07-25 腾讯科技（深圳）有限公司 基于人工智能的图像质量评估方法、装置及计算机设备
CN111612741A (zh) * 2020-04-22 2020-09-01 杭州电子科技大学 一种基于失真识别的精确无参考图像质量评价方法
CN111612741B (zh) * 2020-04-22 2023-11-03 杭州电子科技大学 一种基于失真识别的精确无参考图像质量评价方法

Also Published As

Publication number Publication date
CN108053396B (zh) 2021-09-14

Similar Documents

Publication Publication Date Title
CN108053396A (zh) 一种多失真图像质量的无参考评价方法
CN108428227B (zh) 基于全卷积神经网络的无参考图像质量评价方法
CN109325550B (zh) 基于图像熵的无参考图像质量评价方法
CN105208374B (zh) 一种基于深度学习的无参考图像质量客观评价方法
CN108052980B (zh) 基于图像的空气质量等级检测方法
CN101562675B (zh) 基于Contourlet变换的无参考型图像质量评测方法
CN109255358B (zh) 一种基于视觉显著性和深度图的3d图像质量评价方法
CN104376565B (zh) 基于离散余弦变换和稀疏表示的无参考图像质量评价方法
CN102945552A (zh) 基于自然场景统计中稀疏表示的无参考图像质量评价方法
CN106709958A (zh) 一种基于灰度梯度和颜色直方图的图像质量评价方法
CN109872305B (zh) 一种基于质量图生成网络的无参考立体图像质量评价方法
CN106462771A (zh) 一种3d图像的显著性检测方法
CN104036493B (zh) 一种基于多重分形谱的无参考图像质量评价方法
CN104299232B (zh) 一种基于自适应窗方向波域和改进fcm的sar图像分割方法
CN109816646B (zh) 一种基于退化决策逻辑的无参考图像质量评价方法
CN110070539A (zh) 基于信息熵的图像质量评价方法
CN105160667A (zh) 联合梯度和拉普拉斯-高斯信号的盲图像质量评价方法
CN114187261B (zh) 一种基于多维注意力机制的无参考立体图像质量评价方法
CN109344860A (zh) 一种基于lbp的无参考图像质量评价方法
CN109429051A (zh) 基于多视图特征学习的无参考立体视频质量客观评价方法
CN112950596A (zh) 基于多区域多层次的色调映射全向图像质量评价方法
CN104574424B (zh) 基于多分辨率dct边缘梯度统计的无参照图像模糊度评价方法
CN107146216A (zh) 一种基于梯度自相似度的无参考图像客观质量评价方法
CN109685772A (zh) 一种基于配准失真表示的无参照立体图像质量评估方法
Tang et al. Training-free referenceless camera image blur assessment via hypercomplex singular value decomposition

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant