CN108109141A - 基于超像素多特征匹配的同源局部复制检测方法 - Google Patents

基于超像素多特征匹配的同源局部复制检测方法 Download PDF

Info

Publication number
CN108109141A
CN108109141A CN201711365742.4A CN201711365742A CN108109141A CN 108109141 A CN108109141 A CN 108109141A CN 201711365742 A CN201711365742 A CN 201711365742A CN 108109141 A CN108109141 A CN 108109141A
Authority
CN
China
Prior art keywords
super
pixel
piecemeal
image
color
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201711365742.4A
Other languages
English (en)
Other versions
CN108109141B (zh
Inventor
杨红颖
牛影
牛盼盼
王向阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Normal University
Original Assignee
Liaoning Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Liaoning Normal University filed Critical Liaoning Normal University
Priority to CN201711365742.4A priority Critical patent/CN108109141B/zh
Publication of CN108109141A publication Critical patent/CN108109141A/zh
Application granted granted Critical
Publication of CN108109141B publication Critical patent/CN108109141B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G06T3/147Transformations for image registration, e.g. adjusting or mapping for alignment of images using affine transformations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/40Analysis of texture
    • G06T7/41Analysis of texture based on statistical description of texture
    • G06T7/45Analysis of texture based on statistical description of texture using co-occurrence matrix computation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image quality inspection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于超像素多特征匹配的同源局部复制篡改检测方法,首先对待检测图像进行高斯平滑滤波预处理,使用BEMD算法求取图像纹理程度比例,自适应初始化超像素分块个数;其次,利用基于SLIC超像素分割算法获取图像超像素分块,使用查色表颜色量化技术和纹理矩分析方法得到超像素分块的特征信息;然后,利用Rg2NN算法和BBF算法进行超像素分块的特征匹配;最后,利用SIFT特征点提取算法、RANSAC方法、ZNCC算法、形态学方法等进行后处理。

Description

基于超像素多特征匹配的同源局部复制检测方法
技术领域
本发明属于数字图像取证技术领域,涉及针对复制粘贴的图像篡改检测方法,特别涉及一种基于超像素多特征匹配的同源局部复制篡改检测方法。
背景技术
近年来,图像已经成为人们生活中必不可少的一部分,尤其是在新闻、医学、军事、情报等领域中,图像作为信息的载体和数字证据发挥着不可或缺的重要作用。随着计算机网络技术和便携式数码设备的高速发展以及应用范围的不断扩大,给人们的生活和工作带来效率的同时,也使得图像的篡改变得更为容易。图像处理软件的飞速增长更是催生了大量没有明显痕迹的篡改图像,篡改图像出现的频率越来越高,人们不依靠技术手段,仅仅凭借肉眼去辨认的可能性愈发降低。大多数篡改图像往往会被篡改者用于不良的用途,对整个社会造成极大影响,特别是同源局部复制的篡改,近几年来泛滥出现于各种场合之中,因此急需可靠的篡改检测方法来确定候选图像的可信度,图像篡改检测技术便由此应运而生。
同源局部复制篡改是一种十分常见的图像篡改技术,篡改者往往将图像中的一部分内容进行复制后再粘贴到原图像上的另一部分,用以达到掩盖特定对象的目的。为了使这种篡改更难被检测,篡改者通常会对复制部分进行一系列诸如压缩、旋转、缩放的操作后再进行粘贴。目前,针对同源局部复制篡改的检测方法可以分为两大类:基于块的检测方法和基于特征点的检测方法。其中基于块的检测方法具有定位精度高的优势,且对经过噪声、压缩等处理过的图像检测效果好,但该类方法检测的时间复杂度极高;而基于特征点的检测方法通过提取待检测图像的特征点,并使用相应的局部特征描述子进行特征点之间的匹配,使得该方法的检测时间复杂度得到了大幅度的降低,但该类方法对于平滑的篡改区域检测效果并不理想。因此,现有的大多数篡改检测方法仍然存在检测精度提高,而时间复杂度也随之提高,或是时间复杂度降低,但检测精度也随之降低的问题。
发明内容
本发明是为了解决现有的复制粘贴篡改检测技术所存在的上述技术问题,提供一种基于超像素多特征匹配的同源局部复制篡改检测方法。
本发明的技术解决方案是:一种基于超像素多特征匹配的同源局部复制篡改检测方法,其特征在于按照如下步骤进行:
约定:I指待检测图像;I1指经过高斯平滑滤波预处理的图像;BEMD是指二维经验模式分解算法;SLIC是指超像素分割算法;矩阵superMatrix用于存放所有超像素分块的特征信息;Rg2NN是指逆序广义2近邻算法;BBF为优化查找算法;RANSAC是指去除错误匹配的参数估计方法;ZNCC是立体匹配算法;
a. 初始设置
读取待检测图像I,对I进行高斯平滑滤波预处理得到图像I1;
b. 自适应初始化超像素分块个数
b.1 对图像I进行BEMD分解,得到4个固有模式函数和残差
b.2 按照下式计算残差能量和固有模式函数能量,以及固有模式函数能量占总能量比例,设置近似表示为纹理程度比例
b.3 求取超像素个数,其中,表示分块相关系数;
c. 图像超像素分块
利用SLIC算法,按照步骤b对图像I1进行分割,标记每个超像素分块;
d. 超像素分块特征计算
d.1 选择一超像素分块,采用查色表颜色量化技术对超像素分块做量化,计算查色表中颜色与超像素块颜色之间的欧拉距离
具有最小欧拉距离的颜色即为的量化颜色,统计超像素每一分块的量化颜色,并将最多的前3种主要颜色作为该分块的颜色特征;
d.2 然后使用纹理矩分析方法,求得均值、方差、扭曲度、峰度和熵5个纹理特征如下:
其中为不同的图像灰度级;
d.3 每一超像素分块的特征由3个颜色特征与5个纹理特征组成,重复d.1~d.2,将所有超像素分块的特征信息放入矩阵superMatrix中;
e. Rg2NN快速匹配
e.1 使用BBF算法计算每个分块的近邻特征块;
e.2 计算每一超像素分块和其他分块的欧氏距离并按照从小到大的顺序排列,得到
e.3 依次逆序计算特征描述子欧氏距离之间的比值,若,且,则说明待检测超像素分块和距其个分块均特征相似,取值为0.7;
f. 后处理
f.1 使用SIFT算法进行特征点提取,若一对匹配块中一个块的特征点数为,另一个块的特征点数为,则得到点对
f.2 对每一对匹配块上的匹配点对进行RANSAC操作;
f.3 将所有具有相同仿射变化的匹配点对放到同一矩阵中,再次利用RANSAC方法来消除错误匹配;
f.4 最后使用ZNCC算法及形态学方法对匹配的区域进行标记。
本发明首先对待检测图像进行高斯平滑滤波预处理,使用BEMD算法求取图像纹理程度比例,自适应初始化超像素分块个数;其次,利用基于SLIC超像素分割算法获取图像超像素分块,使用查色表颜色量化技术和纹理矩分析方法得到超像素分块的特征信息;然后,利用Rg2NN算法和BBF算法进行超像素分块的特征匹配;最后,利用SIFT特征点提取算法、RANSAC方法、ZNCC算法、形态学方法等进行后处理。实验结果表明,本发明的方法由于使用SLIC算法、颜色特征、纹理特征以及Rg2NN快速匹配算法,并结合了SIFT算法解决了图像不规则分块的后处理问题,能够在保证分块检测精确度的同时大幅度降低时间复杂度。
与现有技术相比,本发明具有以下有益效果:
第一,使用基于SLIC超像素的图分割算法对待检测图像进行分块,具有运行速度、生成超像素的紧凑度、轮廓保持方面的优势;
第二,结合图像的颜色特征与纹理特征,很好的刻画描述了超像素块的特征信息,在抵抗攻击方面更具有鲁棒性;
第三,使用Rg2NN和BBF快速匹配算法,提高了特征匹配速度以及篡改检测精度;
第四,结合SIFT特征点提取算法与RANSAC方法,有效地解决了不规则分块的后处理问题。
附图说明
图1为本发明实施例篡改区域标记图。
图2为本发明实施例中间具体步骤结果图。
图3为本发明实施例FUA库的部分篡改检测结果图。
图4为本发明实施例GRIP库的部分篡改检测结果图。
图5为本发明实施例Ard库的部分篡改检测结果图。
图6为本发明实施例的流程图。
具体实施方式
本发明的方法如图6所示共包括五个阶段:自适应初始化超像素分块个数、图像超像素分块、超像素分块特征计算、Rg2NN快速匹配和后处理。
约定:I指待检测图像;I1指经过高斯平滑滤波预处理的图像;BEMD是指二维经验模式分解算法;SLIC是指超像素分割算法;矩阵superMatrix用于存放所有超像素分块的特征信息;Rg2NN是指逆序广义2近邻算法;BBF为优化查找算法;RANSAC是指去除错误匹配的参数估计方法;ZNCC是立体匹配算法;
a. 初始设置
读取待检测图像I,对I进行高斯平滑滤波预处理得到图像I1;
b. 自适应初始化超像素分块个数
b.1 对图像I进行BEMD分解,得到4个固有模式函数和残差
b.2 按照下式计算残差能量和固有模式函数能量,以及固有模式函数能量占总能量比例,设置近似表示为纹理程度比例
b.3 求取超像素个数,其中,表示分块相关系数;
c. 图像超像素分块
利用SLIC算法,按照步骤b对图像I1进行分割,标记每个超像素分块;
d. 超像素分块特征计算
d.1 选择一超像素分块,采用查色表颜色量化技术对超像素分块做量化,计算查色表中颜色与超像素块颜色之间的欧拉距离
具有最小欧拉距离的颜色即为的量化颜色,统计超像素每一分块的量化颜色,并将最多的前3种主要颜色作为该分块的颜色特征;
d.2 然后使用纹理矩分析方法,求得均值、方差、扭曲度、峰度和熵5个纹理特征如下:
其中为不同的图像灰度级;
d.3 每一超像素分块的特征由3个颜色特征与5个纹理特征组成,重复d.1~d.2,将所有超像素分块的特征信息放入矩阵superMatrix中;
e. Rg2NN快速匹配
e.1 使用BBF算法计算每个分块的近邻特征块;
e.2 计算每一超像素分块和其他分块的欧氏距离并按照从小到大的顺序排列,得到
e.3 依次逆序计算特征描述子欧氏距离之间的比值,若,且,则说明待检测超像素分块和距其个分块均特征相似,取值为0.7;
f. 后处理
f.1 使用SIFT算法进行特征点提取,若一对匹配块中一个块的特征点数为,另一个块的特征点数为,则得到点对
f.2 对每一对匹配块上的匹配点对进行RANSAC操作;
f.3 将所有具有相同仿射变化的匹配点对放到同一矩阵中,再次利用RANSAC方法来消除错误匹配;
f.4 最后使用ZNCC算法及形态学方法对匹配的区域进行标记。
实验测试和参数设置:
实验是在Windows10×64位环境下,使用Matlab7.12.0执行的,实验所涉及到的三个图像库分别为FAU、GRIP和Ard,均已公开,可以在网上自行搜索下载。
图1为本发明实施例篡改区域标记图。
图2为本发明实施例中间具体步骤结果图。
图3为本发明实施例FUA库的部分篡改检测结果图。
图4为本发明实施例GRIP库的部分篡改检测结果图。
图5为本发明实施例Ard库的部分篡改检测结果图。
图6为本发明实施例流程图。
实验结果表明,本发明的方法由于使用SLIC算法、颜色特征、纹理特征以及Rg2NN快速匹配算法,并结合了SIFT算法解决了图像不规则分块的后处理问题,能够在保证分块检测精确度的同时大幅度降低时间复杂度。

Claims (1)

1.一种基于超像素多特征匹配的同源局部复制篡改检测方法,其特征在于按如下步骤进行:
约定:I指待检测图像;I1指经过高斯平滑滤波预处理的图像;BEMD是指二维经验模式分解算法;SLIC是指超像素分割算法;矩阵superMatrix用于存放所有超像素分块的特征信息;Rg2NN是指逆序广义2近邻算法;BBF为优化查找算法;RANSAC是指去除错误匹配的参数估计方法;ZNCC是立体匹配算法;
a.初始设置
读取待检测图像I,对I进行高斯平滑滤波预处理得到图像I1;
b. 自适应初始化超像素分块个数
b.1 对图像I进行BEMD分解,得到4个固有模式函数和残差
b.2 按照下式计算残差能量和固有模式函数能量以及固有模式函数能量占总能量比例,设置近似表示为纹理程度比例
b.3 求取超像素个数,其中,表示分块相关系数;
c. 图像超像素分块
利用SLIC算法,按照步骤b对图像I1进行分割,标记每个超像素分块;
d. 超像素分块特征计算
d.1 选择一超像素分块,采用查色表颜色量化技术对超像素分块做量化,计算查色表中颜色与超像素块颜色之间的欧拉距离
具有最小欧拉距离的颜色即为的量化颜色,统计超像素每一分块的量化颜色,并将最多的前3种主要颜色作为该分块的颜色特征;
d.2 然后使用纹理矩分析方法,求得均值、方差、扭曲度、峰度和熵5个纹理特征如下:
其中为不同的图像灰度级;
d.3 每一超像素分块的特征由3个颜色特征与5个纹理特征组成,重复d.1~d.2,将所有超像素分块的特征信息放入矩阵superMatrix中;
e. Rg2NN快速匹配
e.1 使用BBF算法计算每个分块的近邻特征块;
e.2 计算每一超像素分块和其他分块的欧氏距离并按照从小到大的顺序排列,得到
e.3 依次逆序计算特征描述子欧氏距离之间的比值,若,且,则说明待检测超像素分块和距其个分块均特征相似,取值为0.7;
f. 后处理
f.1 使用SIFT算法进行特征点提取,若一对匹配块中一个块的特征点数为,另一个块的特征点数为,则得到点对
f.2 对每一对匹配块上的匹配点对进行RANSAC操作;
f.3 将所有具有相同仿射变化的匹配点对放到同一矩阵中,再次利用RANSAC方法来消除错误匹配;
f.4 最后使用ZNCC算法及形态学方法对匹配的区域进行标记。
CN201711365742.4A 2017-12-18 2017-12-18 基于超像素多特征匹配的同源局部复制检测方法 Active CN108109141B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711365742.4A CN108109141B (zh) 2017-12-18 2017-12-18 基于超像素多特征匹配的同源局部复制检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711365742.4A CN108109141B (zh) 2017-12-18 2017-12-18 基于超像素多特征匹配的同源局部复制检测方法

Publications (2)

Publication Number Publication Date
CN108109141A true CN108109141A (zh) 2018-06-01
CN108109141B CN108109141B (zh) 2021-11-19

Family

ID=62209871

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711365742.4A Active CN108109141B (zh) 2017-12-18 2017-12-18 基于超像素多特征匹配的同源局部复制检测方法

Country Status (1)

Country Link
CN (1) CN108109141B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112991408A (zh) * 2021-04-19 2021-06-18 湖南大学 大场景高分遥感图像自适应区域多特征配准方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1900970A (zh) * 2006-07-20 2007-01-24 中山大学 一种鲁棒的图像区域复制篡改检测方法
CN101707711A (zh) * 2009-11-03 2010-05-12 上海大学 一种基于压缩域的视频序列Copy-Move篡改检测方法
CN102521821A (zh) * 2011-10-24 2012-06-27 南开大学 数字图像中一种自动识别和定位篡改区域的方法
CN102956026A (zh) * 2011-08-26 2013-03-06 鲍丽燕 基于不变矩的Copy-Move型篡改图像盲检测方法
CN106683074A (zh) * 2016-11-03 2017-05-17 中国科学院信息工程研究所 一种基于雾霾特性的图像篡改检测方法
KR101755980B1 (ko) * 2016-01-29 2017-07-10 세종대학교산학협력단 스케일 스페이스 표현에 기반한 유사 영역 검출 방법 및 그 장치

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1900970A (zh) * 2006-07-20 2007-01-24 中山大学 一种鲁棒的图像区域复制篡改检测方法
CN101707711A (zh) * 2009-11-03 2010-05-12 上海大学 一种基于压缩域的视频序列Copy-Move篡改检测方法
CN102956026A (zh) * 2011-08-26 2013-03-06 鲍丽燕 基于不变矩的Copy-Move型篡改图像盲检测方法
CN102521821A (zh) * 2011-10-24 2012-06-27 南开大学 数字图像中一种自动识别和定位篡改区域的方法
KR101755980B1 (ko) * 2016-01-29 2017-07-10 세종대학교산학협력단 스케일 스페이스 표현에 기반한 유사 영역 검출 방법 및 그 장치
CN106683074A (zh) * 2016-11-03 2017-05-17 中国科学院信息工程研究所 一种基于雾霾特性的图像篡改检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WANG, XIANG-YANG等: "A new keypoint-based copy-move forgery detection for small smooth regions", 《MULTIMEDIA TOOLS & APPLICATIONS》 *
李岩等: "一种基于逆序广义2近邻的图像多重复制粘贴篡改检测算法", 《电子与信息学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112991408A (zh) * 2021-04-19 2021-06-18 湖南大学 大场景高分遥感图像自适应区域多特征配准方法及系统

Also Published As

Publication number Publication date
CN108109141B (zh) 2021-11-19

Similar Documents

Publication Publication Date Title
CN105608671B (zh) 一种基于surf算法的图像拼接方法
CN103345758B (zh) 基于dct统计特征的jpeg图像区域复制篡改盲检测方法
CN108335331B (zh) 一种钢卷双目视觉定位方法及设备
CN108562589A (zh) 一种对磁路材料表面缺陷进行检测的方法
CN102541954B (zh) 一种商标检索方法及系统
CN111445459B (zh) 一种基于深度孪生网络的图像缺陷检测方法及系统
CN108898132B (zh) 一种基于形状上下文描述的太赫兹图像危险品识别方法
CN104778701A (zh) 一种基于rgb-d传感器的局部图像描述方法
Su et al. A novel forgery detection algorithm for video foreground removal
CN110047063B (zh) 一种物料掉落的检测方法、装置、设备及存储介质
CN109816676A (zh) 一种拼接图像篡改检测方法
CN107240086B (zh) 一种基于积分图算法的织物瑕疵检测方法
CN110738216A (zh) 基于改进surf算法的药品识别方法
CN108122225A (zh) 基于自适应特征点的数字图像篡改检测方法
Muhammad et al. A non-intrusive method for copy-move forgery detection
CN105740872A (zh) 图像特征提取方法及其装置
CN106228557A (zh) 基于二次jpeg压缩的图像篡改区域检测方法
CN110321869A (zh) 基于多尺度融合网络的人员检测和提取方法
CN114187247A (zh) 一种基于图像配准的安瓿瓶印刷字缺陷检测方法
CN109145914A (zh) 一种基于超复数傅里叶变换和均值偏移的显著性目标检测方法
CN108109141A (zh) 基于超像素多特征匹配的同源局部复制检测方法
CN113052234A (zh) 一种基于图像特征和深度学习技术的玉石分类方法
CN109544614B (zh) 一种基于图像低频信息相似度的匹配图像对识别的方法
CN108960285B (zh) 一种分类模型生成方法、舌体图像分类方法及装置
CN108010076B (zh) 一种面向密集工业棒材图像检测的端面外观建模方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant