%?Example?3?Affine?registration??
%?Load?images??
clc;
clear?all;??
close?all;??
%%%發(fā)現(xiàn)只有匹配點(diǎn)對(duì)遍布在圖像中，匹配效果才會(huì)變好??

%
I11=（imread（‘11.jpg‘））;%紅外圖像??
I22=（imread（‘22.jpg‘））;%可見(jiàn)光圖像??
I10=rgb2gray（I11）;??
I20=rgb2gray（I22）;??
I1=imadjust（I10）;??
I2=imadjust（I20）;??


I33=（imread（‘111.png‘））;%紅外圖像???
I30=rgb2gray（I33）;??
%I31=imadjust（I30）;??

%?I1=im2double（edge_detection_1（I1））;??
%?I2=im2double（edge_detection_1（I2））;??
I1?=?im2double（edge（I1‘canny‘））;???
I2?=?im2double（edge（I2‘canny‘））;??
I30?=?im2double（edge（I30‘canny‘））;
figureimshow（I1）;??
figureimshow（I2）;??
figureimshow（I30）;?

%
%?
%?I11?=?imread（‘111-111.jpg‘）;
%?I22?=?imread（‘222-222.jpg‘）;
%?%?I10=rgb2gray（I11）;??
%?%?I20=rgb2gray（I22）;
%?
%?I1=imadjust（I11）;??
%?I2=imadjust（I22）;
%?
%?I1?=?im2double（I1）;???
%?I2?=?im2double（I2）;


??
%?Get?the?Key?Points??
Options.upright=true;??
Options.tresh=0.0001;??
Ipts1=OpenSurf（I1Options）;??%1*268??
Ipts2=OpenSurf（I2Options）;??


%?Put?the?landmark?descriptors?in?a?matrix??
D1?=?reshape（[Ipts1.descriptor]64[]）;??%64*268??
D2?=?reshape（[Ipts2.descriptor]64[]）;??
??
%?Find?the?best?matches??
err=zeros（1length（Ipts1））;?????
cor1=1:length（Ipts1）;??
cor2=zeros（1length（Ipts1））;??
for?i=1:length（Ipts1）??
????distance=sum（（D2-repmat（D1（:i）[1?length（Ipts2）]））.^21）;%用D2的每一列分別去見(jiàn)D1的每一列??
????[err（i）cor2（i）]=min（distance）;??
?????
end??
??
%?Sort?matches?on?vector?distance適量距離排序??
[err?ind]=sort（err）;??
cor1=cor1（ind）;??
cor2=cor2（ind）;??
??
%?Make?vectors?with?the?coordinates?of?the?best?matches??
Pos1=[[Ipts1（cor1）.y]‘[Ipts1（cor1）.x]‘];??
Pos2=[[Ipts2（cor2）.y]‘[Ipts2（cor2）.x]‘];??
Pos1_1=Pos1;??
Pos2_1=Pos2;??
Pos1（:3）=1;?Pos2（:3）=1;Pos1_1（:3）=1;?Pos2_1（:3）=1;??
Pos1=Pos1（1:20:）;??
Pos2=Pos2（1:20:）;??
Pos11=zeros（size（Pos1））;??
Pos22=zeros（size（Pos2））;??

%?Calculate?affine?matrix??
for?i=1:20??
?????k（i）=Pos2（i2）/（Pos2（i1）-Pos1（i1））;??
end??
tic;??
%%??
%自適應(yīng)確定k的范圍??
positive_num=find（k>0）;??
negative_num=find（k<0）;??
if?length（positive_num）>length（negative_num）??
????flag=1?;?
????num_err=negative_num;??
????num_correct=positive_num;??
else??
????flag=0;??
????num_err=positive_num;??
????num_correct=negative_num;??
end??
k（num_err）=[];??
k_1=round（k）;??
k_2=unique（k_1）;%求出k_1有哪些數(shù)，剔除重復(fù)數(shù)據(jù)??
temp=zeros（length（k_2）2）;%第一列存數(shù)，第二列存該數(shù)的個(gè)數(shù)%后邊還需要對(duì)個(gè)數(shù)排序??
for?i=1:length（k_2）??
????temp_num=0;??
????temp（i1）=k_2（i）;??
???for?j=1:length（k）??
???????if?temp（i1）==k_1（j）??
??????????temp_num=temp_num+1;???
???????end??
???end??
???temp（i2）=temp_num;??
end??
%temp??
??
[temp1temp1_location]=sort（temp（:2）‘descend‘）??
xu=temp1（2）;??
temp=temp（temp1_location）??
k_1=zeros（size（12））;??
k_1（1）=temp（1）;??
k_1（2）=temp（2）;??
k_1=sort（k_1‘descend‘）;??
k_11=k_1（11）+0.5;??
k_12=k_1（11）-0.5;??
k_21=k_1（12）-0.5??
k_22=k_1（12）+0.5??
count=0;??
??
for?i=1:20??
?????k=Pos2（i2）/（Pos2（i1）-Pos1（i1））;???%%%%可不可以把數(shù)據(jù)歸一化??
?????if

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-03-12?18:40??1\
?????文件????????5854??2017-12-21?14:05??1\11.jpg
?????文件????????6409??2017-12-21?14:06??1\22.jpg
?????目錄???????????0??2018-03-12?18:40??1\OpenSURF_version1c\
?????文件????????3778??2010-09-07?00:27??1\OpenSURF_version1c\OpenSurf.m
?????目錄???????????0??2018-03-12?18:40??1\OpenSURF_version1c\SubFunctions\
?????文件?????????603??2010-07-27?04:56??1\OpenSURF_version1c\SubFunctions\FastHessian_BuildDerivative.m
?????文件?????????738??2010-07-30?04:32??1\OpenSURF_version1c\SubFunctions\FastHessian_Responselayer.m
?????文件????????1699??2010-07-30?00:10??1\OpenSURF_version1c\SubFunctions\FastHessian_buildResponselayer.m
?????文件????????1877??2010-07-21?00:37??1\OpenSURF_version1c\SubFunctions\FastHessian_buildResponseMap.m
?????文件????????2353??2010-07-30?23:43??1\OpenSURF_version1c\SubFunctions\FastHessian_getIpoints.m
?????文件?????????455??2010-07-27?04:56??1\OpenSURF_version1c\SubFunctions\FastHessian_getLaplacian.m
?????文件?????????450??2010-07-27?04:56??1\OpenSURF_version1c\SubFunctions\FastHessian_getResponse.m
?????文件????????2478??2010-07-30?23:02??1\OpenSURF_version1c\SubFunctions\FastHessian_interpolateExtremum.m
?????文件????????1680??2010-07-30?22:56??1\OpenSURF_version1c\SubFunctions\FastHessian_isExtremum.m
?????文件?????????791??2010-07-29?22:58??1\OpenSURF_version1c\SubFunctions\IntegralImage_BoxIntegral.m
?????文件?????????494??2010-07-30?04:32??1\OpenSURF_version1c\SubFunctions\IntegralImage_HaarX.m
?????文件?????????556??2010-07-30?04:37??1\OpenSURF_version1c\SubFunctions\IntegralImage_HaarY.m
?????文件?????????796??2010-07-29?22:42??1\OpenSURF_version1c\SubFunctions\IntegralImage_IntegralImage.m
?????文件????????1089??2010-07-30?23:39??1\OpenSURF_version1c\SubFunctions\PaintSURF.m
?????文件????????1554??2010-07-31?00:39??1\OpenSURF_version1c\SubFunctions\SurfDescriptor_DecribeInterestPoints.m
?????文件????????3636??2010-07-31?00:40??1\OpenSURF_version1c\SubFunctions\SurfDescriptor_GetDescriptor.m
?????文件????????2921??2010-07-31?00:07??1\OpenSURF_version1c\SubFunctions\SurfDescriptor_GetOrientation.m
?????目錄???????????0??2018-03-12?18:40??1\OpenSURF_version1c\TestImages\
?????文件??????239755??2010-09-07?00:22??1\OpenSURF_version1c\TestImages\lena1.png
?????文件??????211716??2010-09-07?00:22??1\OpenSURF_version1c\TestImages\lena2.png
?????文件????????5296??2010-07-20?05:38??1\OpenSURF_version1c\TestImages\test.png
?????文件??????123228??2010-07-27?01:06??1\OpenSURF_version1c\TestImages\testc1.png
?????文件??????122063??2010-07-27?01:06??1\OpenSURF_version1c\TestImages\testc2.png
?????目錄???????????0??2018-03-12?18:40??1\OpenSURF_version1c\WarpFunctions\
?????文件????????9721??2010-09-07?01:04??1\OpenSURF_version1c\WarpFunctions\affine_warp.m
............此處省略5個(gè)文件信息