JPEG压缩域增强的参数估计方法、装置、介质及计算设备

jpeg压缩域增强的参数估计方法、装置、介质及计算设备
技术领域
1.本发明涉及图像检测和处理技术领域，尤其涉及一种jpeg压缩域增强的参数估计方法、装置、存储介质及计算设备。


背景技术：

2.jpeg(joint photographic experts group)即联合图像专家组，是用于连续色调静态图像压缩的一种标准，文件后缀名为.jpg或.jpeg，是最常用的图像文件格式。
3.jpeg压缩域增强是一种增强jpeg图像对比度的处理操作。在将两幅jpeg图像进行合成时，如果两者的对比度存在差异，则容易留下明显的合成痕迹。一种减轻合成痕迹的常用做法是分别调整两幅图像的对比度使之匹配，这会导致合成图像中的某个局部的增强参数与其它局部的增强参数不一致。此外，如果两幅jpeg图像的压缩参数本身就不同，也会导致合成图像中某个局部的压缩参数与其它局部的压缩参数不一致。因此，通过估计可疑图像各个局部的增强参数和压缩参数，并进行一致性比对，有助于判明该可疑图像是否为合成图像，为防范合成图像用于非法目的提供技术保障。
4.jpeg压缩域增强图像可以看作是一幅经过jpeg压缩后再经过增强的图像，因此其参数的估计包含两个方面：增强参数的估计和压缩参数的估计。目前发现的一些图像参数估计方案主要为估计jpeg压缩位图的压缩参数而设计，由于建模分析时未考虑压缩域增强的影响，导致此类方法若直接用于估计压缩域增强图像的压缩参数，往往效果很差甚至无法实施，并且缺乏估计增强参数的能力。
5.目前能同时估计增强参数和压缩参数的方法尚未见报道，为估计jpeg压缩位图的压缩参数而设计的方法并不适用于同时估计jpeg压缩域增强图像的增强参数和压缩参数。


技术实现要素：

6.鉴于现有技术存在的不足，本发明提供了一种jpeg压缩域增强的参数估计方法、装置、存储介质及计算设备，可以用于对给定的一幅经过jpeg压缩域增强的位图，同时精确地估计其jpeg压缩域增强图像的增强参数和压缩参数，并可进一步用于图像篡改检测。
7.为了实现上述的目的，本发明采用了如下的技术方案：
8.一种jpeg压缩域增强的参数估计方法，包括：
9.对输入的jpeg压缩域增强图像的像素矩阵分块，挑选出适于参数估计的可用块；
10.对每一个可用块进行离散余弦变换，得到对应的系数矩阵dk；
11.对系数矩阵dk中的每一个空间频率(i,j)，计算对应的分频似然函数l
ij
；
12.根据分频似然函数l
ij
，计算总似然函数l；
13.根据总似然函数l计算增强参数估计值
14.根据增强参数估计值和分频似然函数l
ij
计算压缩参数估计值
15.作为其中一种实施方式，所述对对系数矩阵dk中的每一个空间频率(i,j)，计算对
应的分频似然函数l
ij
的步骤，包括：
16.遍历增强参数的候选值集合s中的每一个增强参数候选值λ
t
，对系数矩阵dk的每个系数进行反向增强操作，得到反向增强系数矩阵；
17.根据反向增强系数矩阵，计算每个空间频率(i,j)对应的分频似然函数l
ij
。
18.作为其中一种实施方式，所述对输入的jpeg压缩域增强图像的像素矩阵分块，挑选出适于参数估计的可用块的步骤，包括：
19.对输入的jpeg压缩域增强图像的像素矩阵进行8
×
8分块；
20.排除分块中的奇异块，所述奇异块的最大灰度值为255，或最小灰度值为0，或最大灰度值等于最小灰度值。
21.作为其中一种实施方式，所述根据反向增强系数矩阵，计算每个空间频率(i,j)对应的分频似然函数l
ij
的步骤，包括：
22.根据拉普拉斯概率质量函数p
l
和高斯截断概率密度函数g，计算出概率分布函数p；
23.根据反向增强系数矩阵和概率分布函数p，计算出分频似然函数l
ij
。
24.作为其中一种实施方式，所述分频似然函数l
ij
满足：
[0025][0026]
其中，k为分块索引，1≤k≤k，为系数矩阵dk中空间频率(i,j)处的系数，0≤i,j≤7；q
ij
为空间频率(i,j)处的整数量化步长，取值范围为1≤q
ij
≤qb；概率分布函数p满足：表示整数集；
[0027]
其中，拉普拉斯概率质量函数p
l
满足：
[0028][0029]
其中，高斯截断概率密度函数g满足：
[0030][0031]
σ2＝1/12，b＝4，
[0032]
作为其中一种实施方式，所述总似然函数l满足：
[0033][0034]
所述增强参数估计值满足：
[0035][0036]
s为增强参数的候选值集合，s＝{λ1,λ2,
…
,λ
t
,
…
,λ
t
}；
[0037]
所述根据增强参数估计值和分频似然函数l
ij
计算压缩参数估计值包括：
[0038]
将增强参数估计值代入每一个空间频率(i,j)对应的分频似然函数l
ij
中，通过下式得到压缩参数估计值
[0039][0040]
本发明的另一目的在于提供一种jpeg压缩域增强的参数估计装置，包括：
[0041]
筛分模块，用于对输入的jpeg压缩域增强图像的像素矩阵分块，挑选出适于参数估计的可用块；
[0042]
变换模块，用于对每一个可用块进行离散余弦变换，得到对应的系数矩阵dk；
[0043]
分频似然函数计算模块，用于对系数矩阵dk中的每一个空间频率(i,j)，计算对应的分频似然函数l
ij
；
[0044]
总似然函数计算模块，用于根据分频似然函数l
ij
，计算总似然函数l；
[0045]
第一估算模块，用于根据总似然函数l计算增强参数估计值
[0046]
第二估算模块，用于根据增强参数估计值和分频似然函数l
ij
计算压缩参数估计值
[0047]
作为其中一种实施方式，所述分频似然函数计算模块包括：
[0048]
变换模块，用于对每一个可用块进行离散余弦变换，得到对应的系数矩阵dk；
[0049]
反向增强模块，用于遍历增强参数的候选值集合s中的每一个增强参数候选值λ
t
，对系数矩阵dk的每个系数进行反向增强操作，得到反向增强系数矩阵；
[0050]
运算模块，用于根据反向增强系数矩阵，计算每个空间频率(i,j)对应的分频似然函数l
ij
。
[0051]
本发明的又一目的在于提供一种存储介质，所述存储介质内存储有多条指令，所述指令适于由处理器加载并执行所述jpeg压缩域增强的参数估计方法的步骤。
[0052]
本发明的又一目的在于提供一种计算设备，包括所述存储介质和适于实现各指令的处理器。
[0053]
本发明针对压缩域增强的特点，通过计算相应的分频似然函数和总似然函数，根据总似然函数的得到增强参数估计值，根据分频似然函数得到压缩参数估计值，能够同时精准估计两类参数。
附图说明
[0054]
图1为本发明实施例的一种jpeg压缩域增强的参数估计方法的流程示意图；
[0055]
图2为本发明实施例的一种jpeg压缩域增强的参数估计装置的结构框图。
具体实施方式
[0056]
在本发明中，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本技术的描述中，“多个”的含义是至少两个，例如两
个，三个等，除非另有明确具体的限定。
[0057]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
[0058]
参阅图1，本发明实施例提供了一种jpeg压缩域增强的参数估计方法，包括：
[0059]
s01、对输入的jpeg压缩域增强图像的像素矩阵分块，挑选出适于参数估计的可用块，排除奇异块。
[0060]
在步骤s01中，对输入的jpeg压缩域增强图像像素矩阵分块，挑选出适于增强检测的可用块的步骤，具体可以包括：
[0061]
s011、对输入的jpeg压缩域增强图像的像素矩阵进行8
×
8分块。
[0062]
例如，可以记输入的jpeg压缩域增强图像像素矩阵为i，对输入的jpeg压缩域增强图像像素矩阵i进行8
×
8分块。
[0063]
s012、排除分块中的奇异块，奇异块的最大灰度值为255，或最小灰度值为0，或最大灰度值等于最小灰度值。
[0064]
考虑到如果一个分块中的最大灰度值为255，或最小灰度值为0，或最大灰度值等于最小灰度值，则该分块会对增强检测会产生不良影响，因此此类奇异块需要弃之不用。排除所有奇异块后，剩下的分块称为可用块，可用块数量记为k。
[0065]
s02、对每一个可用块进行离散余弦变换，得到对应的系数矩阵dk。
[0066]
其中的上标k表示分块索引，1≤k≤k。记系数矩阵dk中空间频率(i,j)处(0≤i,j≤7)的系数为
[0067]
s03、对系数矩阵dk中的每一个空间频率(i,j)，计算对应的分频似然函数l
ij
。
[0068]
计算对应的分频似然函数l
ij
前，还需要实现对系数矩阵dk进行处理，其计算过程具体包括：
[0069]
s031、遍历增强参数的候选值集合s中的每一个增强参数候选值λ
t
，对系数矩阵dk的每个系数进行反向增强操作，得到反向增强系数矩阵。
[0070]
这里，设定增强参数的候选值集合s＝{λ1,λ2,
…
,λ
t
,
…
,λ
t
}。遍历每一个候选值λ
t
，对系数进行反向增强操作：
[0071][0072]
其中的f为增强算子，表示用增强参数λ
t
进行反向增强后的dct(discrete cosine transform，离散余弦变换)系数。
[0073]
s032、根据反向增强系数矩阵，计算每个空间频率(i,j)对应的分频似然函数l
ij
：
[0074][0075]
其中，式(2)中的q
ij
为空间频率(i,j)处的整数量化步长，且1≤q
ij
≤qb，最大取值qb可供设定；log表示对数运算；p为概率分布函数，定义为：
[0076]
[0077]
其中，式(3)中的表示整数集，p
l
为拉普拉斯概率质量函数，g为高斯截断概率密度函数。
[0078]
因此，计算每个空间频率(i,j)对应的分频似然函数l
ij
，首先需要根据拉普拉斯概率质量函数p
l
和高斯截断概率密度函数g，计算出概率分布函数p；然后才根据反向增强系数矩阵和概率分布函数p，计算出分频似然函数l
ij
。
[0079]
式(3)中，拉普拉斯概率质量函数p
l
定义为：
[0080][0081]
其中，式(4)中，ρ
ij
为尺度参数，可通过进行计算，exp表示以自然常数e＝2.71828
…
为幂的指数函数。
[0082]
式(3)中，高斯截断概率密度函数g定义为：
[0083][0084]
其中，式(5)中的σ2＝1/12，b＝4，z的取值为π表示圆周率。
[0085]
s04、根据分频似然函数l
ij
，计算总似然函数l。
[0086]
在分频似然函数l
ij
式(2)的基础上，进一步计算总似然函数l：
[0087][0088]
总似然函数l(λ
t
)是空间频率集合f中所有频率对应的分频似然函数l
ij
的最大值之和，空间频率集f为可设参数。
[0089]
s05、根据总似然函数l计算增强参数估计值
[0090]
增强参数可根据下式估计：
[0091][0092]
s为增强参数的候选值集合，s＝{λ1,λ2,
…
,λ
t
,
…
,λ
t
}。
[0093]
s06、根据增强参数估计值和分频似然函数l
ij
计算压缩参数估计值
[0094]
压缩参数估计值的计算具体包括：将增强参数估计值代入每一个空间频率(i,j)对应的分频似然函数l
ij
中，通过下式得到压缩参数估计值
[0095][0096]
在其他实施方式中，上述步骤s05中的通过最大化l(λ
t
)来估计增强参数，也可以变更为：
[0097]
[0098]
上述步骤s06中的通过最大化来估计压缩参数，也可以变更为：
[0099][0100]
针对压缩域增强的特点，本实施例通过上述一系列分频似然函数和总似然函数的构造和计算过程，通过总似然函数的最大值估计增强参数，通过分频似然函数的最大值估计压缩参数，能够同时实现jpeg压缩域增强图像增强参数和压缩参数的准确估计。
[0101]
本实施例还对上述的jpeg压缩域增强图像参数的估算方法进行了实验测试，具体情况如下：
[0102]
(1)测试图库的选择：采用ucid图像库，共包含1338幅图像，对ucid库中的每一幅图像，进行质量因子qf为50、60、70、80、90的jpeg压缩，然后对所得的jpeg图像进行增强参数λ为0.7、0.8、0.9、1.0、1.1、1.2、1.3、1.4、1.5、1.6、1.8、2.0、2.2、2.4、2.6、2.8、3.0的jpeg压缩域增强，得到经jpeg压缩域增强的图像；
[0103]
(2)参数设定：设定增强参数的候选值集合s＝{0.6,0.61,0.62,
…
,3.29,3.30}，共包含271个候选值；qb设定为25；空间频率集合f设定为包含前9个交流频率，即f＝{(0,1),(0,2),(0,3),(1,0),(1,1),(1,2),(2,0),(2,1),(3,0)}。
[0104]
(3)测试结果：不同qf情况下的增强参数和压缩参数(只统计前9个交流频率)的估计准确率如下：
[0105][0106]
可以看出，对于不同的jpeg压缩质量因子qf，本实施例的jpeg压缩域增强的参数估计方法都能实现精准的增强参数估计和压缩参数估计，验证了本实施例估计方法的有效性。
[0107]
如图2所示，本发明还提供了一种jpeg压缩域增强的参数估计装置，包括筛分模块1、变换模块2、分频似然函数计算模块3、总似然函数计算模块4、第一估算模块5以及第二估算模块6，筛分模块1用于对输入的jpeg压缩域增强图像的像素矩阵分块，挑选出适于参数估计的可用块；变换模块2用于对每一个可用块进行离散余弦变换，得到对应的系数矩阵dk；分频似然函数计算模块3用于对系数矩阵dk中的每一个空间频率(i,j)，计算对应的分频似然函数l
ij
；总似然函数计算模块4用于根据分频似然函数l
ij
，计算总似然函数l；第一估算模块5用于根据总似然函数l计算增强参数估计值第二估算模块6用于根据增强参数估计值和分频似然函数l
ij
计算压缩参数估计值
[0108]
具体地，筛分模块1可包括划分模块11和筛选模块12，划分模块11用于对输入的
jpeg压缩域增强图像的像素矩阵进行8
×
8分块，筛选模块12用于排除分块中的奇异块，奇异块的最大灰度值为255，或最小灰度值为0，或最大灰度值等于最小灰度值。
[0109]
具体地，分频似然函数计算模块3可包括变换模块31、反向增强模块32和运算模块33，变换模块31用于对每一个可用块进行离散余弦变换，得到对应的系数矩阵dk；反向增强模块32用于遍历增强参数的候选值集合s中的每一个增强参数候选值λ
t
，对系数矩阵dk的每个系数进行反向增强操作，得到反向增强系数矩阵；运算模块33用于根据反向增强系数矩阵，计算每个空间频率(i,j)对应的分频似然函数l
ij
。
[0110]
而运算模块33具体还可以包括第一子运算模块331和第二子运算模块332，第一子运算模块331用于根据拉普拉斯概率质量函数p
l
和高斯截断概率密度函数g计算出概率分布函数p，第二子运算模块332用于根据反向增强系数矩阵和概率分布函数p，计算出分频似然函数l
ij
。
[0111]
本发明的又一目的在于提供一种存储介质，存储介质内存储有多条指令，指令适于由处理器加载并执行jpeg压缩域增强的参数估计方法的步骤。
[0112]
本发明的又一目的在于提供一种计算设备，包括存储介质和适于实现各指令的处理器。
[0113]
另外，本发明还提供了一种存储介质及计算设备，该存储介质内存储有多条指令，该指令适于由处理器加载并执行上述的jpeg压缩域增强的参数估计方法的步骤，该存储介质为计算设备的一部分。处理器在一些实施例中可以是中央处理器(central processing unit，cpu)、控制器、微控制器、微处理器、或其他数据处理芯片。该处理器通常用于控制计算设备的总体操作。本实施例中，该处理器用于运行存储介质中存储的程序代码或者处理数据。
[0114]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台终端(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
[0115]
本发明针对压缩域增强的特点，通过计算相应的分频似然函数和总似然函数，根据总似然函数的得到增强参数估计值，根据分频似然函数得到压缩参数估计值，能够同时精准估计两类参数。
[0116]
以上所述仅是本技术的具体实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本技术原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本技术的保护范围。