基于空域感知的块效应检测方法、装置、设备和存储介质与流程

1.本技术涉及视屏及图像的压缩技术领域，特别是涉及一种基于空域感知的块效应检测方法、装置、设备和存储介质。


背景技术：

2.在视频及图像的压缩过程中，随着压缩比率及压缩参数的变化，解码后的视频或图像会产生不同程度的质量损失，即有损压缩现象。压缩算法利用人眼对于图像中某些频率信息的不敏感即在有损压缩过程中的一些信息丢失不会引起感知的特性，在损失信息在人们可接受范围之内的条件下，算法就可以通过损失一些原始信息以换取更大的压缩比。其中影响视频或图像直观观看质量的重要噪声就是压缩噪声，即视频压缩过程中产生的块效应现象。
3.块效应的检测方法根据原理可以分为基于空域的块效应检测和基于变换域的检测，基于空域的块效应评价算法先定位块效应发生的位置，然后通过分析空间域块边界处的灰度变化来估计块效应严重程度。基于变换域的块效应评价算法的核心思想是在图像进行块编码时，图像块的尺寸是不变的，因而块效应会在复原后的图像上周期性的出现，因此借助dct、离散傅里叶变换(discrete fourier transform，dft)、dwt等周期性的变换，可以实现对图像块效应的评价。


技术实现要素：

4.基于此，有必要针对上述技术问题，提供一种能够提高检测效率及检测精准度的基于空域感知的块效应检测方法、装置、设备和存储介质。
5.一种基于空域感知的块效应检测方法，所述方法包括：
6.获取待检测的图像，所述图像为被压缩后再进行解压后的受损图像；
7.对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
8.根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
9.根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
10.根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
11.通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
12.在其中一实施例中，对所述灰度图像进行划分时，各所述单元块的尺寸按照对图像进行压缩时分块的尺寸划分。
13.在其中一实施例中，所述根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差包括：
14.其中，计算所述水平方向的像素差采用以下公式：
15.dh＝t(x,y 1)-t(x,y),y＝1,2,...,n-1
16.在上式中，用t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，dh为所述水平方向上的像素差。
17.在其中一实施例中，根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度包括：
18.根据所述水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，再根据所述水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内部的块效应强度；
19.根据所述垂直方向上的像素差进行计算，得到垂直方向上各单元块之间的块效应强度，再根据所述垂直方向上的像素差以及各单元块之间的块效应强度进行计算得到垂直方向上单元块内部的块效应强度；
20.根据水平方向上单元块内部的块效应强度以及垂直方向上单元块内部的块效应强度进行计算得到所述单元块内部的平均块效应强度；
21.根据水平方向上各单元块之间的块效应强度以及垂直方向上各单元块之间的块效应强度进行计算得到所述单元块之间的平均块效应强度。
22.在其中一实施例中，根据所述水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，采用以下公式：
[0023][0024]
在上式中，表式向上取整，p为常数，bs为单元块的大小。
[0025]
在其中一实施例中，根据所述水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内的块效应强度，采用以下公式：
[0026][0027]
在上式中，p为常数，bs为单元块的大小。
[0028]
在其中一实施例中，根据所述预处理后的图像中各像素点的亮度值进行计算，得到各所述单元块内的方差值采用以下公式：
[0029][0030]
在上式中，t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，bs为单元块的大小。
[0031]
一种基于空域感知的块效应检测装置，所述装置包括：
[0032]
图像获取模块，用于获取待检测的图像，所述图像为被压缩后再进行解压后的受
损图像；
[0033]
图像预处理模块，用于对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
[0034]
像素差计算模块，用于根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0035]
平均块效应强度计算模块，用于根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，用于得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0036]
方差值计算模块，用于根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
[0037]
块效应检测模块，用于通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
[0038]
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
[0039]
获取待检测的图像，所述图像为被压缩后再进行解压后的受损图像；
[0040]
对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
[0041]
根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0042]
根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0043]
根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
[0044]
通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
[0045]
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
[0046]
获取待检测的图像，所述图像为被压缩后再进行解压后的受损图像；
[0047]
对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
[0048]
根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0049]
根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0050]
根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
[0051]
通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效
应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
[0052]
上述基于空域感知的块效应检测方法、装置、设备和存储介质，通过对图像进行预处理后，先计算图像水平方向上以及垂直方向上的像素差，再根据这两个方向上的像素差进行计算分别得到图像中单元块内部的平均块效应强度以及单元块之间的平均块效应强度，并且还计算单元块内部的方差作为块效应的补充，最后根据单元块内部及之间的平均块效应强度以及块内方差通过非线性函数进行融合，并根据最后的结果对图像的块效应进行检测，降低块效应的检测难度的同时提高检测效率及精准度。
附图说明
[0053]
图1为一个实施例中基于空域感知的块效应检测方法的流程示意图；
[0054]
图2为一个实施例中基于空域感知的块效应检测方法的实施步骤框架示意图；
[0055]
图3为一个实施例中基于空域感知的块效应检测方法装置的结构框图；
[0056]
图4为一个实施例中计算机设备的内部结构图。
具体实施方式
[0057]
为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
[0058]
如图1所示，提供了一种基于空域感知的块效应检测方法，包括以下步骤：
[0059]
步骤s100，获取待检测的图像，图像为被压缩后再进行解压后的受损图像；
[0060]
步骤s110，对图像进行预处理，预处理包括将图像进行灰度转换后得到对应的灰度图像，再将灰度图像划分为多个单元块；
[0061]
步骤s120，根据预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0062]
步骤s130，根据水平方向上的像素差以及垂直方向上的像素差进行计算，得到图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0063]
步骤s140，根据预处理后的图像中各像素点的亮度值进行计算，得到单元块内部的方差值；
[0064]
步骤s150，通过非线性函数对单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对图像的块效应进行检测。
[0065]
在步骤s100中，待进行块效应检测的图像是基于dct(离散余弦变换)进行压缩后再进行恢复后的受损图像。而该图像可以是一张单独的图像，或者是视屏中的任意一帧图像。
[0066]
在步骤s110中，在对灰度图像进行划分时，各单元块的尺寸按照对图像进行压缩时分块的尺寸划分。在采用dct压缩技术时，一般以8的倍数对图像进行分块计算，也就是说各图像块中均包括有8*8或者16*16的像素点，而发生块效应也大多出现在这些图像块的边界上，故将灰度图像按照相同的方式进行划分。最终灰度图像中的单元块的数量根据图像
的尺寸大小决定，而各单元块的尺寸一般为8*8或者16*16。
[0067]
在步骤s120中，根据预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差采用以下公式进行计算，由于水平方向上与垂直方向上的像素差计算方式类似，在这里给出水平方向的像素差计算公式作为例子进行说明，公式为：
[0068]dh
＝t(x,y 1)-t(x,y),y＝1,2,...,n-1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0069]
在公式(1)中，用t(x,y)表示尺寸为m*n的图像中各像素点的亮度值，而dh表示水平方向上的像素差。
[0070]
由于在基于8
×
8的dct的编码技术在对dct变换后的系数进行量化时会造成量化误差，主要表现为块内的块效应和块间的块效应。块内的块效应主要是由于丢失了高频的dct系数造成，而块间的块效应主要是由于块边界的不连续性造成的。所以，在步骤s130中，根据水平方向上的像素差以及垂直方向上的像素差进行计算，得到图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度，具体计算包括：
[0071]
根据水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，再根据水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内部的块效应强度；
[0072]
根据垂直方向上的像素差进行计算，得到垂直方向上各单元块之间的块效应强度，再根据垂直方向上的像素差以及各单元块之间的块效应强度进行计算得到垂直方向上单元块内部的块效应强度；
[0073]
根据水平方向上单元块内部的块效应强度以及垂直方向上单元块内部的块效应强度进行计算得到单元块内部的平均块效应强度；
[0074]
根据水平方向上各单元块之间的块效应强度以及垂直方向上各单元块之间的块效应强度进行计算得到单元块之间的平均块效应强度。
[0075]
同样的，由于计算水平方向上的单元块内部的块效应强度及单元块之间的块效应强度的方法与计算垂直方向上的是类似的，所以在本文中以计算水平方向上的块效应强度为例进行说明。
[0076]
具体的，根据水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，采用以下公式：
[0077][0078]
在公式(2)中，表式向上取整，p为常数，bs为单元块的大小。
[0079]
进一步的，根据水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内的块效应强度，采用以下公式：
[0080][0081]
在公式(3)中，p为常数，bs为单元块的大小。
[0082]
接着，再通过计算水平和垂直方向的块效应强度平均值作为整幅图像对应的评价指标。具体的，也就是根据水平方向上单元块内部的块效应强度以及垂直方向上单元块内部的块效应强度进行计算得到单元块内部的平均块效应强度。根据水平方向上各单元块之
间的块效应强度以及垂直方向上各单元块之间的块效应强度进行计算得到单元块之间的平均块效应强度。
[0083]
其中，计算单元块之间的平均块效应强度采用以下公式：
[0084]
o＝(oh ov)/2
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0085]
在公式(4)中，oh为水平方向上各单元块之间的块效应强度，ov为垂直方向上各单元块之间的块效应强度。
[0086]
其中，计算单元块内部的平均块效应强度采用以下公式：
[0087]
i＝(ih iv)/2
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0088]
在公式(5)中，ih为水平方向上单元块内部的块效应强度，iv为垂直方向上单元块内部的块效应强度。
[0089]
在步骤s140中，根据预处理后的图像中各像素点的亮度值进行计算，得到各单元块内的方差值采用以下公式：
[0090][0091][0092]
在公式(6)中，t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，bs为单元块的大小。
[0093]
在步骤s150中，通过非线性函数对单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算采用以下公式：
[0094][0095]
在公式(7)中，a、α1、α2均为常数。
[0096]
最后根据计算结果与预设的阈值对图像的块效应进行检测。
[0097]
在利用上述方法进行实际检测时，还可按照图2的算法框架流程图进行实施。
[0098]
上述基于空域感知的块效应检测方法中，通过对图像进行预处理后，先计算图像水平方向上以及垂直方向上的像素差，再根据这两个方向上的像素差进行计算分别得到图像中单元块内部的平均块效应强度以及单元块之间的平均块效应强度，并且还计算单元块内部的方差作为块效应的补充，最后根据单元块内部及之间的平均块效应强度以及块内方差通过非线性函数进行融合，并根据最后的结果对图像的块效应进行检测。本方法提出的块效应评价散发相对于主流的全参考块效应评价算法、非特定失真块效应评价或者特定失真的块效应评价算法，算法结果在与人类mos值的单调一致性有了很大的提高；相对主流的基于滤波、分块、频域变换思想的噪声评价算法，算法的稳定性和准确性上都有极大的提高。并且在本方法中还独创性的提供了像素差、单元块间及单元块内部的块效应强度的计算公式以及单元块内部的方差值计算公式，进一步简化的块效应检测的步骤。
[0099]
应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行
完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
[0100]
在一个实施例中，如图3所示，提供了一种基于空域感知的块效应检测装置，包括：图像获取模块200、图像预处理模块210、像素差计算模块220、平均块效应强度计算模块230、方差值计算模块240以及块效应检测模块250，其中：
[0101]
图像获取模块200，用于获取待检测的图像，所述图像为被压缩后再进行解压后的受损图像；
[0102]
图像预处理模块210，用于对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
[0103]
像素差计算模块220，用于根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0104]
平均块效应强度计算模块230，用于根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，用于得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0105]
方差值计算模块240，用于根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
[0106]
块效应检测模块250，用于通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
[0107]
关于基于空域感知的块效应检测装置的具体限定可以参见上文中对于基于空域感知的块效应检测方法的限定，在此不再赘述。上述基于空域感知的块效应检测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0108]
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图4所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于空域感知的块效应检测方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0109]
本领域技术人员可以理解，图4中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0110]
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
[0111]
获取待检测的图像，所述图像为被压缩后再进行解压后的受损图像；
[0112]
对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
[0113]
根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0114]
根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0115]
根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
[0116]
通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
[0117]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0118]
对所述灰度图像进行划分时，各所述单元块的尺寸按照对图像进行压缩时分块的尺寸划分。
[0119]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0120]
所述根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差包括：
[0121]
其中，计算所述水平方向的像素差采用以下公式：
[0122]dh
＝t(x,y 1)-t(x,y),y＝1,2,...,n-1
[0123]
在上式中，用t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，dh为所述水平方向上的像素差。
[0124]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0125]
根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度包括：
[0126]
根据所述水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，再根据所述水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内部的块效应强度；
[0127]
根据所述垂直方向上的像素差进行计算，得到垂直方向上各单元块之间的块效应强度，再根据所述垂直方向上的像素差以及各单元块之间的块效应强度进行计算得到垂直方向上单元块内部的块效应强度；
[0128]
根据水平方向上单元块内部的块效应强度以及垂直方向上单元块内部的块效应强度进行计算得到所述单元块内部的平均块效应强度；
[0129]
根据水平方向上各单元块之间的块效应强度以及垂直方向上各单元块之间的块效应强度进行计算得到所述单元块之间的平均块效应强度。
[0130]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0131]
根据所述水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，采用以下公式：
[0132]
[0133]
在上式中，表式向上取整，p为常数，bs为单元块的大小。
[0134]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0135]
根据所述水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内的块效应强度，采用以下公式：
[0136][0137]
在上式中，p为常数，bs为单元块的大小。
[0138]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0139]
根据所述预处理后的图像中各像素点的亮度值进行计算，得到各所述单元块内的方差值采用以下公式：
[0140][0141]
在上式中，t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，bs为单元块的大小。
[0142]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
[0143]
获取待检测的图像，所述图像为被压缩后再进行解压后的受损图像；
[0144]
对所述图像进行预处理，所述预处理包括将所述图像进行灰度转换后得到对应的灰度图像，再将所述灰度图像划分为多个单元块；
[0145]
根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差；
[0146]
根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应强度以及各单元块之间的平均块效应强度；
[0147]
根据所述预处理后的图像中各像素点的亮度值进行计算，得到所述单元块内部的方差值；
[0148]
通过非线性函数对所述单元块内部的平均块效应强度、各单元块之间的平均块效应强度以及单元块内部的方差值进行融合计算，并根据计算结果与预设阈值进行比较，以对所述图像的块效应进行检测。
[0149]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤对所述灰度图像进行划分时，各所述单元块的尺寸按照对图像进行压缩时分块的尺寸划分。
[0150]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤所述根据所述预处理后的图像进行计算，分别得到水平方向上的像素差以及垂直方向上的像素差包括：
[0151]
其中，计算所述水平方向的像素差采用以下公式：
[0152]dh
＝t(x,y 1)-t(x,y),y＝1,2,...,n-1
[0153]
在上式中，用t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，dh为所述水平方向上的像素差。
[0154]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤根据所述水平方向上的像素差以及垂直方向上的像素差进行计算，得到所述图像中单元块内部的平均块效应
强度以及各单元块之间的平均块效应强度包括：
[0155]
根据所述水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，再根据所述水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内部的块效应强度；
[0156]
根据所述垂直方向上的像素差进行计算，得到垂直方向上各单元块之间的块效应强度，再根据所述垂直方向上的像素差以及各单元块之间的块效应强度进行计算得到垂直方向上单元块内部的块效应强度；
[0157]
根据水平方向上单元块内部的块效应强度以及垂直方向上单元块内部的块效应强度进行计算得到所述单元块内部的平均块效应强度；
[0158]
根据水平方向上各单元块之间的块效应强度以及垂直方向上各单元块之间的块效应强度进行计算得到所述单元块之间的平均块效应强度。
[0159]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤根据所述水平方向上的像素差进行计算，得到水平方向上各单元块之间的块效应强度，采用以下公式：
[0160][0161]
在上式中，表式向上取整，p为常数，bs为单元块的大小。
[0162]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤根据所述水平方向上的像素差以及各单元块之间的块效应强度进行计算得到水平方向上单元块内的块效应强度，采用以下公式：
[0163][0164]
在上式中，p为常数，bs为单元块的大小。
[0165]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0166]
根据所述预处理后的图像中各像素点的亮度值进行计算，得到各所述单元块内的方差值采用以下公式：
[0167][0168]
在上式中，t(x,y)表示尺寸为m*n的所述图像中各像素点的亮度值，bs为单元块的大小。
[0169]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强
型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。
[0170]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0171]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。