基于压缩量化与邻域空间LBP 算子的图像哈希算法.docx

1、基于压缩量化与邻域空间LBP 算子的图像哈希算法基于压缩量化与邻域空间LBP 算子的图像哈希算法 基于压缩量化与邻域空间LBP 算子的图像哈希算法 王彦超 （平顶山教育学院，平顶山 467000） 摘要：目的 为了解决哈希算法的感知鲁棒性与伪造检测能力不高的问题，提出基于特征压缩机制与邻域空间局部二值模式的紧凑图像哈希算法。方法 首先利用2D线性插值技术，对输入图像进行预处理；嵌入Ring分割技术，将其变为二次图像；再利用Gabor滤波技术对其完成过滤；考虑到图像的颜色特征与其内在的空间关系，基于局部二值模式LBP设计邻域空间LBP算子，提取滤波图像的特征；构建特征压缩量化准则，输出紧凑的哈希

2、二值数组；迭代Logistic映射，输出随机序列，通过量化每个序列值输出密钥流，通过构建动态引擎设计分段异加密模型，实现紧凑哈希序列的加密，获取图像哈希；最后计算原始哈希序列与待检测哈希序列的Hamming距离，实现图像信息的安全认证。结果 与已有的哈希生成机制相比，文中算法所输出的哈希序列更紧凑，对旋转、伽马校正等篡改操作具有更好的感知鲁棒。结论 所提哈希技术具备较高的安全性，在包装图标检索、信息检测等领域具有较好的价值。 关键词：图像哈希；特征压缩量化规则；邻域空间局部二值模式；分段异加密；Gabor滤波；决策阈值 强大的图像编辑工具能够随意篡改图像而且无明显的痕迹，对图像信息造成了巨大威

3、胁12。为了能够对图像的真伪实行准确检测，研究人员设计了图像内容认证技术，其中之一就是图像哈希技术34。如张秋余等5提出了基于非下采样轮廓波变换的彩色图像感知哈希算法。此技术忽略了图像的颜色特征与其内在的空间关系，且其输出的哈希序列维数较高。Yong Soo等6提出了基于层次直方图的图像哈希算法。此哈希技术主要是利用了图像的直方图分布特性，导致其对旋转篡改的识别能力较弱。余俊伟等7提出了基于局部不变矩和DWT特征矩阵的图像哈希算法。由于旋转操作会导致子块与其旋转版本会存在较大的内容差异，且单纯利用图像的 V分量来生产哈希序列，忽略了图像的颜色特征。为了提高哈希算法的感知鲁棒性与紧凑度，文中提出

4、基于特征压缩机制与邻域空间局部二值模式的紧凑图像哈希算法。 1 文中图像哈希算法 所提的紧凑图像哈希算法过程见图1。 图1 文中图像哈希认证算法过程Fig.1 The propsoed image Hash authentication algorithm process 1.1 基于2D线性插值的图像预处理 首先，引入2D线性插值技术8，固定哈希尺寸： 式中：x0, x1为初始起点；y0, y1为x0, x1的插值结果；x0zx1, y0zy1。 再利用式（1）与式（2），联合卷积掩模3来实现初始图像的预处理。令TG（i,j）为卷积掩模中（i,j）处的元素，其模型为： 式中：为卷积掩模的标准

5、差； 另外，文中将图像RGB变为YCbCr空间9： 式中：R, G, B为输入图像的红、绿、蓝三分量；Y, Cb, Cr为亮度、蓝色与红色浓度偏移量。以图2a为例，利用式（13）对其插值后，输出结果见图2b。 图2 图像预处理结果Fig.2 Image preprocessing results 1.2 基于Ring分割的二次图像输出 文中基于Ring分割技术10输出二次图像。令f0（x,y）的尺寸为 mm；n为环形数量；Rk是像素值集合。通过计算每个像素与图像中心之间的距离以及计算环形半径来对输入图像实施合理的Ring分割，其示意见图3。令为第k（k=1,2,3n）个环形的半径，那么，对于m

6、m的图像f0（x,y），其令rn的计算函数为： 式中：为向下取整运算。 图3 Ring分割过程Fig.3 Ring segmentation process 为了确定其他环形半径，首先计算内接圆A的面积以及每个环形的平均面积uA见式（6）。 则最内环形区域的面积r1为： 对于其他环形区域 rk（k=1,2,3n?1）的面积，其计算模型为： 接下来，需估算其他位置像素与图像中心间的距离。令（xc,yc）代表图像中心位置，当m为偶数时，xc=yc=m/2+0.5；若 m 为奇数，则 xc=yc=（m+1）/2+0.5。则任意像素与中心像素的距离d为： 根据式（5）、式（7）、式（8）以及式（9），

7、获取距离d后，文中将这些像素值变为n个集合： 随后，将 Rk（k=1,2,3n）中的元素按照从小到大的顺序重排，生成新的矢量uk，确保uk与旋转操作无关。再利用线性插值技术，将uk映射为矢量vk，其尺寸为uA1。那么，对vk进行重排，生成二次图像V。 以图2c为例，对其进行逆时针旋转32，得到其旋转版本，见图4a；并利用Ring分割对其分环处理，结果见图4b，最终生成的二次图像见4c。 图4 二次图像的生成Fig.4 Secondary image generation 1.3 基于Gabor滤波器的二次图像平滑处理 利用Gabor滤波技术来平滑二次图像，其模型为11： 式中：x,y为Gabo

8、r滤波的窗口尺寸；x, y分别为x,y的标准差；为正弦波的波长；j为正弦波的方向。 式中：k为Gabor滤波的方向数目。 利用式（13）处理二次图像，结果见图 4d。由图 4d可知，Gabor滤波图像可显著改善哈希序列对噪声、光照亮度干扰的鲁棒性。 1.4 基于领域空间LBP算子的特征提取 文中考虑图像的颜色特征与其内在的空间关系，并基于局部二值模式12，设计了邻域空间LBP算子。传统的LBP仅依靠单个像素区域内的特征，其模型为12： 式中：S（x）为中心和领域点的灰度差值；gc, gi分别为中心像素及其邻域像素点；LR,P为半径为R且存在P个邻域点的LBP算子。这种LBP算子没有考虑图像色彩

9、与空间结构关，为此，文中设计邻域空间LBP算子。 首先，依据式（4），将RGB变为YCbCr空间；并利用YCbCr中任意2个分量p1, p2来估算2个领域与中心像素的偏差 式中： Ip 1（g i ） 为分量 p1中的第 i个领域值；Ip 1 （g c ）,I p 2（g c ） 分别为分量p1, p2的中心领域值。 通过对比D p1 i ,D p2i ，得到新的灰度差值： 为了与式（17）对应，文中将式（19）进行分割，得到子灰度差： 随后，根据式（16）对 s1（p1ip2i）, s2（p1ip2i）完成编码： 再结合 F1（p1p2）, F2（p1p2），获取色彩 LBP算子 LBP（p

10、1p2）。将空间结构关系引入到 LBP（p1p2）中，得到邻域空间 LBP机制。为了充分利用相邻区域的空间结构特性，文中考虑图像中2个相邻的像素区域： 式中：为2个LBP中心像素点的距离；r为LBP像素区的半径；c为左侧LBP中心位置。 由于 LBP提取的特征维数较高。为了降低特征维数，文中在任意一个LBP区域择取4个相邻像素点，见图5。 图 5 不同（r, ）下的 P(r, )Fig.5 P(r, ) with different (r, ) 通过考虑色彩R, G, B, Cr这4个分量，联合式（22），得到邻域空间LBP算子见式（23）。 利用文中设计的邻域空间 LBP算子提取 Gabor

11、滤波图像的特征，输出特征矢量 FV =v1 , v 2 , v3 .v N 。以图2a为目标，利用式（24）对其处理，获取LBP特征图像，见图6。由图6可知，目标的鲁棒特征被充分提取出来，无背景干扰特征。 1.5 基于特征压缩量化准则的低维哈希序列形成 文中构建了压缩量化准则，以降低 Fv=v1, v2,v3vN的维数，输出紧凑哈希 Fv=v1, v2, v3vN。其步骤如下所述。 图6 邻域空间LBP算子的输出结果Fig.6 Output results of LBP operator in neighborhood space 1）引入最邻近域思想13，对Fv=v1, v2, v3vN中的

12、每个元素完成相似性连接；若b1,b2为FV的相邻两元素，且|b1，b2|2时，将b1认为是b2的邻近点，再将二者进行连接。通过对Fv=v1, v2, v3vN中的每个元素完成该过程 则可输出邻接图。 2）计算邻接图中任意2个邻近点的相似度权重： 3）根据式（24）计算的权重，将Fv=v1, v2, v3vN完成映射，对其维数进行压缩： 式中：为实数。 得到式（25）的特征矢量后，将代入 F =T（F ） VV 中，可输出压缩特征矢量FV。随后，对低维特征矢量FV进行量化，输出紧凑哈希序列B=（B1, B2, B3BL）。 式中：t为阈值，其计算函数见式（27）。 式中：N为哈希长度；i为矢量v

13、i中所有元素的方差。 1.6 基于分段异加密模型与 Hamming距离的哈希认证 为了改善哈希序列的安全度，文中设计哈希加密机制。删除哈希数组B=（B1, B2, B3BN）中的B1元素，再统计其他元素值之和： 再根据Sum值，估算初值E0： 同时，引入Logistic映射14，设置初值x0, 进行迭代N次，生成混沌数组x1, x2, x3xN。 并对x1, x2, x3xN进行量化，输出混沌密钥流ki。 根据E0与ki，对B=（B1, B2, B3BN）中的首个元素完成加密。 再构建2个引擎q1, q2来对哈希序列中其他元素进行加密： 式中：q11, i?1, q2i+1, N均为动态引擎；

14、N为哈希序列长度。 利用q1, q2，设计了哈希加密机制： 式中：Bi为加密哈希值。 反复执行上述哈希加密过程，当i=N?1时结束。对于哈希BN，根据引擎q2对其进行加密： 执行整个加密机制后，可输出最终的哈希序列H=H1, H2, H3HN。 完成哈希加密后，再对其实施认证。若初始图像对应的哈希序列，待检测图像的 哈 希 序 列 是 。文中利用Hamming距离15D来计算H0, H1的相似度： 式中：为异或运算。 根据式（37），若D值低于用户阈值W时，认为接收到的图像的内容为真实；反之，其内容被篡改。 2 实验结果与分析 借助 UCID数据库16对所提哈希认证算法完成安全性测试，同时为了

15、体现所提哈希技术的优异性，将文献6，18视为对照组，并利用ROC曲线17对3种哈希技术进行评估。关键实验参数为：x0=1.5, =3.2,n=8, k=6。 2.1 阈值W的优化 在UCID图像库16中选择1500幅图像视为样本，将如下内容操作赋予每一幅图像：椒盐噪声分别为0.01, 0.03, 0.05, 0.1；亮度调整因子为0.6, 0.8, 1.2,2.6；伽玛校正因子为0.1, 0.4, 0.6, 1；JPEG压缩因子为 0.6, 0.8, 1.4, 1.6；旋转角度为 10, 30, 60, 110；缩放尺度为0.5, 0.9, 1.2 1.5。 不同图像的Hamming距离以及频

16、数分布情况见图 7。根据图 7中数据可知，就相同图像而言，当Hamming距离小于0.53时，其频数变化较大。对于不同的图像对，其Hamming距离大于0.53，频数变化剧烈，因此，在文中哈希技术中，取阈值 W=0.53进行鲁棒性与安全性实验测试。 图7 阈值的优化测试结果Fig.7 Threshold optimization test results 2.2 哈希算法的鲁棒性测试 鲁棒性是哈希序列的重要评估指标6，为此，从UCID图像库中任意选择4幅图像，见图8ad，并上述的篡改内容对其进行处理；同时，利用式（37）输出其Hamming距离，测试数据见图9a9f。根据图9中显示的Hammi

17、ng距离，所提哈希技术的Hamming距离D均要小于0.53。 2.3 哈希算法的敏感性测试 若图像在网络传输中遭遇攻击使其数据被篡改，则攻击前后的图像哈希序列时存在较大差异3。文中对图像施加6种攻击，见图10。同时估算二者的Hamming距离，图10bf的Hamming距离分别为0.543, 0.591,0.637, 0.564, 0.602。利用文中哈希技术对其进行检测后，它们所对应的Hamming距离均超过阈值0.53，这显示文中哈希算法对上述 6种攻击具有强烈的敏感性，可对图像内容信息的真伪进行准确认证。 2.4 哈希算法的安全性测试 图8 测试图像Fig.8 Test image 图

18、9 文中哈希技术的感知鲁棒性测试Fig.9 Perceptual robustness test of the proposed Hash technique 图10 6种数字操作对应的图像Fig.10 Images corresponding to six digital operations 理想的哈希算法是具有的较高的安全性7，因此文中对2400组错误密钥进行测试，统计其对应的Hamming距离，结果见图11。由统计结果可知，对于这2400组错误密钥，所提哈希算法产生的Hamming距离值均超过了阈值0.53。这显示了所提哈希技术具备理想的安全性。 2.5 不同哈希算法的鲁棒性分析 在U

19、CID图像库16中任意选择800幅图像作为测试样本，利用文中算法、文献6，18对其进行处理，得到ROC曲线见图12。根据ROC曲线可知，所提哈希算法对应的ROC特性最佳，尤其是面对旋转与伽马校正攻击，在 PFPR=0.2时，所提哈希技术对应的 PTPR=0.983。文献6，18的感知鲁棒性均要低于所提哈希技术，其中，文献6对旋转攻击的识别能力最弱，文献18的鲁棒性虽然要优于文献6，但是仍然不理想，在 PFPR=0.2时，二者的 PTPR分别为0.861, 0.925。 图11 所提哈希算法的安全性测试Fig.11 Security test of the proposed Hash algor

20、ithm 图12 3种算法的ROC曲线测试Fig.12 ROC curve test of three algorithms 3 结语 提出了一种新的紧凑图像哈希算法来压缩哈希序列长度与提高其感知鲁棒性。联合 Ring分割与Gabor滤波技术，将输入图像转变为抗旋转的二次图像，并基于经典的 LBP算子，考虑图像的颜色特征与其内在的空间关系，设计了邻域空间 LBP算子，增强其特征描述能力，有效提取鲁棒特征。另外，通过设计相应的加密函数来进一步改善哈希序列的安全性。实验结果显示所提哈希算法的认证能力更好，具有理想的ROC特性曲线。 参考文献： 1?AVU?O?LU , KA?AR S, PEHLI

21、VAN I. Secure Image Encryption Algorithm Design Using A Novel Chaos Based S-BoxJ. Chaos, Solitons and Fractals,2017, 95(16): 92101. 2郭静博, 孙琼琼. 改进的引力模型耦合明文像素相关交叉机制的图像加密算法J. 包装工程, 2016,37(13): 165172.GUO Jing-bo, SUN Qiong-qiong. Image Encryption Algorithm Based on Improved Gravity Model Coupled with

22、Plain Pixel Correlation Crossing MechanismJ.Packaging Engineering, 2016, 37(13): 165172. 3王彦超, 郭静博, 周丽宴. 基于数据投影降维机制与对称局部二值模式的紧凑图像哈希算法J. 激光与光电子学进展, 2017, 54(2): 416.WANG Yan-chao, GUO Jing-bo, ZHOU Li-yan. Compact Image Hashing Algorithm Based on Data Projection Dimensionality Reduction Mechanism and

23、 Symmetric Local Two Value ModelJ. Laser & Optoelectronics Progress, 2017, 54(2): 416. 4SUN Rui, ZENG Wen-jun. Secure and Robust Image Hashing Via Compressive SensingJ. Multimedia Tools and Applications, 2014, 70(3): 16511665. 5张秋余, 杨子, 李凯. 基于非下采样轮廓波变换的彩色图像感知哈希算法J. 兰州理工大学学报(自然科学版), 2016, 42(5): 9510

24、1.ZHANG Qiu-yu, YANG Zi, LI Kai. Color Image Perception Hash Algorithm Based on Non-sampling Contourlet TransformJ. Journal of Lanzhou University of Technology(Natural Science Edition), 2016, 42(5):95101. 6YONG S C, PARK J H. Image Hash Generation Method Using Hierarchical HistogramJ. Multimedia Too

25、ls and Applications, 2015, 61(1): 181194. 7余俊伟. 基于局部不变矩和DWT特征矩阵的图像哈希算法D. 广西: 广西师范大学, 2015.YU Jun-wei. Image Hashing Algorithm Based on Local Moment Invariants and DWT Feature MatrixD. Guangxi: Guangxi Normal University, 2015. 8ROS?ONIEC S. An Example of Two-Dimensional Interpolation Using a Linear C

26、ombination of Bicubic B-SplinesJ. International Journal of Electronics and Telecommunications, 2012, 45(6): 109118. 9方旺盛, 李玉南, 张蓉. 一种视觉模型引导的小波域彩色图像盲水印算法J. 计算机应用研究, 2011,28(7): 27192722.FANG Wang-sheng, LI Yu-nan, ZHANG Rong. A Blind Watermarking Algorithm for Color Image Guided by Visual ModelJ. Com

27、puter Application Research, 2011, 28(7): 27192722. 10TANG Zhen-jun, ZHANG Xian-quan. Robust Perceptual Image Hashing Based on Ring Partition and NMFJ. IEEE Transactions on Knowledge & Data Engineering, 2014, 26(3): 711724. 11王延年, 牛飞婷, 刘婷. 基于Gabor滤波的指纹图像增强方法J. 激光杂志, 2016, 11(1): 117119.WANG Yan-nian,

28、 NIU Fei-ting, LIU Ting. Fingerprint Image Enhancement Method Based on Gabor FilterJ.Journal of Laser, 2016, 11(1): 117119. 12李春利, 沈鲁娟. 基于改进LBP算子的纹理图像分类方法J. 计算机工程与设计, 2016, 37(1): 232236.LI Chun-li, SHEN Lu-juan. Texture Image Classification Method Based on Improved LBP OperatorJ.Computer Engineerin

29、g and Design, 2016, 37(1):232236. 13LIU Yang, LI Mei-dong, BI Xiao-ru. An Improved Location Difference of Multiple Distances Based Nearest Neighbors Searching AlgorithmJ. Optik-International Journal for Light and Electron Optics,2016, 127(22): 1083810843. 14胡春杰, 陈晓, 陈霞. 基于改进广义 Arnold映射的多混沌图像加密算法J. 包

30、装工程, 2017, 38(3):144149.HU Chun-jie, CHEN Xiao, CHEN Xia. Multi Chaotic Image Encryption Algorithm Based on Improved Generalized Arnold MapJ. Packaging Engineering,2017, 38(3): 144149. 15付海燕. 基于图像哈希的大规模图像检索方法研究D. 辽宁: 大连理工大学, 2015.FU Hai-yan. Research on Image Retrieval Based on Image HashD. Liaoning

31、: Dalian University of Technology, 2015. 16SCHAEFER G, STICH M. UCID-an Uncompressed Color Image DatabaseM. Proceedings of SPIE, Storage and Retrieval Methods and Applications for Multi-media, 2004. 17夏蕾, 周冰. 像素预测误差耦合似然映射的移动-复制图像伪造检测算法J. 量子电子学报, 2016, 33(2):153161.XIA Lei, ZHOU Bing. A motion Detect

32、ion Algorithm Based on Pixel Prediction Error Coupled Likelihood MappingJ. Quantum Electronics Journal, 2016,33(2): 153161. 18刘凯. 基于压缩感知的图像哈希算法D. 广西: 广西师范大学, 2016.LIU Kai. Image Hashing Algorithm Based on Compressed SensingD. Guangxi: Guangxi Normal University, 2016. Image Hashing Algorithm Based on Compression Quantization and Neighborhood Space LBP Operator WANG Yan-chao(Pingdingshan Institute o