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源碼\