图像与多媒体技术・Image&Multimedia Technology 基于Harris角点检测的数字图像反取证技术研究 文/高飞 王正光 fspecial(’gaussian’,[5 1],1.5) Harri S角点检9/4就属于基于 即: 加无关角点,破坏图像中原有角点的位置。 模版的检测算法。通过识别出图 像的角点就可以大致描绘出图像 的轮廓。而且通过比对和匹配图 像的角点还可以对复制一粘贴伪 造区域进行检测和定位。本文通 过对Harri S角点检测算法进行系 W(x,y)=[0.1201;O.2339;0.2921;0.2339; 0.120l J。 Harris角点所做的高斯滤波也可看成是某 个脉冲和图像卷积的结果。所以使用逆滤波器 可以降低高斯滤波对于图像的作用。不过这种 令 逆滤波器是提前运用的,即在脉冲对图像进行 AI=0.1201 蜒+0.2339"目+o 2921’碍+0.2339 40.1201‘鹾 统的分析和实验研究旨在寻找一 种方法在不改变图像大致原貌的 A3 l如1・研+0.2339*耐,eO.2921・研m2339 mlⅪ1 职 基础上隐藏,篡改,增加图像角点, 使Harri s角点检测算法无法正常 A2 0.1201 IxlIyl+0.2339 Ix2Iy3+O.2921 I 识别出图像应有的角点。 3Iy3+0.2339 Ix4Iy5+0.1201 Ix5Iy5 则 【关键词】Harti S角点检测数字图像反取证 M=[ A I ] :^ 一 + 技术 其中^和 分别是M的两个特征值。通过计算 R的值即可确定窗口是否包含角点。如果一个 像素点的R值大于某个闽值并且大于周围相 角点是图像的一种极其重要的特征。准确 邻8个像素点的R值,即可确认该点为角点。 地找到一幅图像的角点有助于理解和分析该图 所以说要想处理一张图片使得Harris角 像。同时角点检测又是图像取证的一种基本方 点检测算法无法正确地检测出角点即改变像素 法。角点检测可分为三类,分别是基于图像灰 的R值即可。 度的角点检测、基于二值图像的角点检测、基 由于Harris算法在做检测之前先进行了高 于轮廓的角点检测。而基于图像灰度的检测又 斯滤波,所以各个点之间并不是单独存在的, 分为梯度、模版和模板梯度三种类型。Harris 而是互相之间有关联的,改变一个像素的值会 角点检测就属于基于模版的检测算法。通过 引发周围其他像素R值的改变。这样,经过 识别出图像的角点就可以大致描绘出图像的轮 处理后再进行角点计算所得的结果是未知的, 廓。而且通过比对和匹配图像的角点还可以对 无法预测的。因为每幅图片像素之间的关系并 复制一粘贴伪造区域进行检测和定位。 不相同,如果只是改变其中一个像素,处理后 本文通过对Harris角点检测算法进行系 的结果依赖于图片本身。Harris检测算法做的 统的分析和实验研究旨在寻找一种方法在不改 高斯滤波既有利又有弊,好的一方面是处理了 变图像大致原貌的基础上隐藏,篡改,增加图 噪声点的问题,使得图像中的噪声点无法被当 像角点,使Harris角点检测算法无法正常识别 成角点检测出来。但是不好的一方面是引入了 出图像应有的角点。 相关性,使得各个像素之间关联性增强。怎样 1 Har ri s算法 解除这种相关性成了Harris角点反检测的重要 Harris角点检测算法的本质是用一个窗口 步骤。 在图像上沿着周围八个方向移动,然后通过计 2基于Har r j s角点检测的反取证算法及 算公式: 逆滤波器分析 M= )[ zx 2 矿Ixy] 所谓的基于Harris角点检测的反取证算法 是指对一幅图像进行处理,使得在用Harris算 求得窗口中是否存在角点,其中w(x,y) 法对这幅图像进行角点检测时所检测出的角点 是一个高斯平滑滤波器,普遍性地,在使用 多为伪角点,而真实的角点则被隐藏。基于这 Matlab计算图像角点时,该高斯滤波函数取: 个目的,可以对图像中的角点进行隐藏,并增 104・电子技术与软件工程Electronic Technology&Software Engineering 干扰之前所做的预先的逆滤波。换成另外一种 角度来看,提前做的处理可以当成噪声脉冲, 那么Harris角点检测所做的高斯滤波就成为了 噪声脉冲的逆滤波。提前经过处理的图像在经 过进行高斯滤波之后,就又变回了原图像,这 样Harris在寻找角点时各个像素之间的关联就 可以降低,在这个基础上更改某个像素的值并 不会对于其他的点产生影响。 3代码实现 3.1逆滤波器的代码 forl=1:k ofrO=一3:3 for02=一3:3 Pl=sum(dOub1 e(current(C(1,2)+O一 2:C(1,2)+o+2,c(1,1)+o2)). h2); P=sum(pl(:)); current(c(1,2)十o,c(1,1)+o2)= uint8(p); . end end end 3.2 Harri S角点检测的Matlab代码 ofri=2:height-1 ofrj=2:width一1 if R(ij)>0.01 Rmax&&R(ij)>R(i一 1J一1)&&R(ij)>R(i—Ij)&&R(ij)>R(i一1J+1) &&R(iJ)>R(i,j一1)&&R(ij)>R(ij+1)&& R(ij)>R(i+lj一1)&&R(ij)>R(i+ij)&&ROj) >R(i+1j+1) result(i,j) 1; cnt=cnt+l; end; mage&Multimedia Technology・图像与多媒体技术 表l:1000张图片处理后统计数据 次数 l 原角点数 l06-38 处理后角点数 102.41479 消除角点数 89.45337 增加角点数 85.48338 未处理角点数 16.931404 未处理角点数百分比 17.8431340909091% 2 1O6-38 90.410503 90.41264 74.43837 l5.9721329 16.7026954476807% 3 106-38 79.046087 90.92818 63.58949 15.45659164 16.1315910l5118800% 4 106.38 73.213290 91.22722 58.05573 15.15755627 15.7519466882416OO% 5 1O6_38 65.826366 91.45659 50.898l7 14.92818864 l5.49139943844490O% 6 106_38 60.637727 91.65273 45.90568 14.73204716 l5.1562O59395248O0% 7 1O6.38 56.82958 1 91.77920 42.22400 14.60557342 15.131228833693300% 8 lO6-38 53.803858 91.90139 39.32047 14.48338692 l4.89099510248l100% 9 l06.38 51.357984 91.96141 36.9346l 14.42336549 14.893O844228695OO% l0 106-38 49.658092 91.95069 35.3 1082 14.34726688 14.7934744708423O0% l5 106.38 45.136120 92.08574 30.92390 14.2122l865 14.566559676O259O0% 20 106-38 43.790996 92.13290 29.62593 14.16505895 l4.526758812O950OO% 50 106_38 43.033226 92.241l5 28.88960 14.14362272 14.477791317494600% 1O0 l06-38 43.025723 92.24115 28.88210 l4.14362272 14.4777913 l7494600% 500 1O6.38 43.032154 92.16184 28.89603 14.13612O04 14.477897904967600% end; end; 入分析意图寻找一种反取证算法对图像中的角 点进行隐藏和修改,使得基于Harris角点检测 的图像取证方法失效。 通过对Harris算法进行深入研究可以了解 Harti S角点检测方法[J】.计算机应用与 软件,2008(12):242-243. 4反取证算法综分析 因为Harris算法中寻找角点是用每个像 素的R值与其周围像素的R值做比较,所以 只更改其中一个像素值很可能会影响其他相邻 像素的R值的比较结果。所以对一幅图像进 行多次反取证算法的处理可以隐藏尽可能多的 角点。 [2]赵万金,龚声蓉,刘纯平,沈项军.一种 自适应的Harri S角点检测算法【J].计算 机工程,2008(1 0):212-214. 角点检测算法的基本原理和作用,从而找出一 种合适并且高效的反取证算法。反取证的实现 和应用可以对篡改,抄袭图片的常用手段进行 了解和认识并且在了解这些手段的前提下对于 图像篡改等情况进行防范和分析。本文中所提 出的反取证算法有着简单,对于图像的依赖性 【3】张黎,马燕,李顺宝.基于细节增强的自 适应Harri S角点检测算法[J】.计算机应 用与软件,2O1 2(04):253-255. 【4】刘睿,王锋,陈卫东,孔李军.基于小波 变换多尺度Harri S角点检测算法[J].微 计算机信息,2009(1 8):244—245. 表1是对1000张不同图片经过实验后得 较强,不同图像的处理效果不一样,无法滤除 到的平均统计数据。从中可以发现对图像进行 8次反取证算法处理是较为妥当的 掉图像中的噪声点,会使一些噪声被当作角点 检测出来,对于复杂的图片处理效果不明显。 作者单位 1.南开大学滨海学院天津市 300270 5结论 参考文献 本文通过对Harris角点检测算法进行深 [1]王德超,2.天津师范大学天津市300387 涂亚庆,张铁楠.一种多尺度 Electronic Technology&Software Engineering电子技术与软件工程・105