[发明专利]JPEG-缩放-JPEG图像初次量化步长估计方法

权利要求书

1.JPEG-缩放-JPEG图像初次量化步长估计方法，其特征在于，该方法包括以下步骤：

选取一组候选缩放因子并对待测图像进行逆缩放操作，对逆缩放操作后的图像DC通道DCT系数的PSD直方图进行伽马校正和归一化处理，其中，具有最显著峰的PSD直方图对应的缩放因子作为估计的缩放因子；

利用估计的缩放因子对待测图像进行逆缩放操作，得到重构的初次压缩图像；

根据重构的初次压缩图像DC通道DCT系数PSD直方图中最显著峰的位置信息，利用几何拟合计算DC通道的初次量化步长；

根据AC通道的DCT系数直方图的分布特性，利用滤波算法和极大似然估计计算AC通道的初次量化步长。

2.根据权利要求1所述的JPEG-缩放-JPEG图像初次量化步长估计方法，其特征在于，所述对逆缩放操作后的图像DC通道DCT系数的PSD直方图进行伽马校正和归一化处理，包括：

其中，为伽马校正后的PSD直方图的第e个数值；为待测图像的PSD直方图的第e个数值；max(f_psd)为待测图像的PSD直方图中的最大数值；γ为伽马值；为伽马校正和归一化后的PSD直方图的第e个数值；N为伽马校正后的PSD直方图中的数量。

3.根据权利要求1所述的JPEG-缩放-JPEG图像初次量化步长估计方法，其特征在于，所述DC通道的初次量化步长的计算公式为：

其中，Q为初次量化步长；round()为四舍五入函数；fp^*为最显著峰对应的频率；K为频谱宽度；bias为超参数。

4.根据权利要求1所述的JPEG-缩放-JPEG图像初次量化步长估计方法，其特征在于，所述根据AC通道的DCT系数直方图的分布特性，利用滤波算法和极大似然估计计算AC通道的初次量化步长，包括：

对DCT系数直方图进行滤波；滤波的方法为：

输入：H_Origin：重构的初次压缩图像AC通道DCT系数直方图；

pks_sorted：直方图的峰值大小(未排序)；

locs_sorted：直方图的峰值频率位置(未排序)；

pks_unsorted：直方图的峰值大小(按大小从大到小排序)；

locs_unsorted：直方图的峰值频率位置(按大小从大到小排序)；

index_valley：直方图的波谷的频率坐标；

index_peak：直方图的波峰的频率坐标；

输出：H_filted：过滤后的AC通道DCT系数直方图

begin

for i＝1to length(index_peak)-1do

index_{prex_peak}＝locs_unsorted[index_peak[i]]；//获得前置峰值坐标

index_{suff_peak}＝locs_unsorted[index_peak[i+1]]；//获得后置峰值坐标

//从前置峰值坐标到后置峰值坐标遍历

for j＝index_{prex_peak}to index_{suff_peak}do

//获得前置峰值坐标到后置峰值之间的最小波谷值的位置

index_{true_cursor}＝0；

//获得前置峰值坐标到后置峰值之间的最小波谷值的数值

value_{temp_min}＝pks_sorted[1]；

if index_{cursor_valley}in[index_{prex_p}eak，index_{suff_peak}]then

//获取真实的最小波谷值

if H_Origin[locs_unsorted[m]]＜＝value_{temp_min}then

//更新最小波谷值大小

value_{temp_min}＝H_Origin[locs_unsprted[m]]；

//更新最小波谷值位置

index_{true_cursor}＝locs_unsorted[m]；

for n in locs_unsorteddo

//从峰值位置到最小波谷值中的数值全部加至峰值位置

sum_temp＝sum_temp+H_Origin[n]；

H_Origin[n]＝0；

for n in locs_unsorteddo；

//从最小波谷值到峰值位置中的数值全部加至峰值位置

if nindex_{true_cursor}n＜index_{suff_peak}then

sum_temp＝sum_temp+H_Origin[n]；

H_Origin[n]＝0；

H_filted＝H_Origin；//获滤波后的AC通道DCT系数直方图

end。

5.根据权利要求1所述的JPEG-缩放-JPEG图像初次量化步长估计方法，其特征在于，利用极大似然估计计算AC通道的初次量化步长，包括：

其中，likeh为AC通道的初次量化步长；max为取最大值函数；D为过滤后的AC通道的DCT系数直方图的值；λ是步长；1(·)为判断函数。