SIFT复制粘贴篡改检测

matlab代码，用SIFT特征点检测图像中复制粘贴篡改，需安装vl_feat库

clc;
clear;
tic;
color=imread（‘40.png‘）;%1\
result_name=‘F:\2016 lvjing\picture\000008.jpg‘;%1\
if size（size（color）2）==3
    im=rgb2gray（color）;
else
    im=color;
end
im=single（im）;

[fd] = vl_sift（im） ;
pos=f（1:2:）‘;
scale=f（3:）‘;
descr=double（d‘）;
count=size（descr1）; %特征点个数
S=[];
D=[];
color=double（color）;      

for i=1:size（descr1） %所有特征点循环
    k=1;       
    distance=zeros（1（size（descr1）））;
    for j=1:size（descr1）
        distance（k）=norm（descr（i:）-descr（j:））;
        k=k+1;
    end
        [distanceposition]=sort（distance）;
        m=distance（2）/distance（3）;
        if （m<0.3） 
            S=[S;pos（i:）];
            D=[D;pos（position（2）:）];                
        end
end
correctPoints = ransac_stitch（S（:1）S（:2）D（:1）D（:2））;   
figureimshow（color/256）;
hold on;
for i = 1: size（S1）
% Rif find（correctPoints（:）==i）
    line（[S（i1）D（i1）] ...
         [S（i2）D（i2）] ‘Color‘ ‘g‘‘LineWidth‘2）;
% end
  
end
%plot（pos（:1）‘pos（:2）‘‘*‘）;
plot（S（:1）‘S（:2）‘‘*‘）;
plot（D（:1）‘D（:2）‘‘*‘）;
hold off;


saveas（gcfresult_name）%

match=size（S1）;
toc;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      152076  2017-06-13 16:26  SIFT源码\0014.jpg
     文件      153882  2017-06-13 16:34  SIFT源码\0016.jpg
     文件    12738707  2011-05-08 17:03  SIFT源码\40.png
     文件        1112  2017-12-11 17:23  SIFT源码\detectdistance.m
     文件        1015  2017-03-19 15:00  SIFT源码\ransac_stitch.m
     目录           0  2017-12-11 18:05  SIFT源码\