%?Example?3?Affine?registration
%?Load?images
I1=im2double（imread（‘TestImages/lena1.png‘））;
I2=im2double（imread（‘TestImages/lena2.png‘））;

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

%?Put?the?landmark?descriptors?in?a?matrix
D1?=?reshape（[Ipts1.descriptor]64[]）;
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）;
????[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=Pos1（1:30:）;
Pos2=Pos2（1:30:）;

%?Show?both?images
I?=?zeros（[size（I11）?size（I12）*2?size（I13）]）;
I（:1:size（I12）:）=I1;?I（:size（I12）+1:size（I12）+size（I22）:）=I2;
figure?imshow（I）;?hold?on;

%?Show?the?best?matches
plot（[Pos1（:2）?Pos2（:2）+size（I12）]‘[Pos1（:1）?Pos2（:1）]‘‘-‘）;
plot（[Pos1（:2）?Pos2（:2）+size（I12）]‘[Pos1（:1）?Pos2（:1）]‘‘o‘）;

%?Calculate?affine?matrix
Pos1（:3）=1;?Pos2（:3）=1;
M=Pos1‘/Pos2‘;

%?Add?subfunctions?to?Matlab?Search?path
functionname=‘OpenSurf.m‘;
functiondir=which（functionname）;
functiondir=functiondir（1:end-length（functionname））;
addpath（[functiondir?‘/WarpFunctions‘]）

%?Warp?the?image
I1_warped=affine_warp（I1M‘bicubic‘）;

%?Show?the?result
figure
subplot（131）?imshow（I1）;title（‘Figure?1‘）;
subplot（132）?imshow（I2）;title（‘Figure?2‘）;
subplot（133）?imshow（I1_warped）;title（‘Warped?Figure?1‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1776??2012-05-29?22:53??surf\example3.m

?????文件???????1335??2010-09-06?07:21??surf\license.txt

?????文件???????3778??2010-09-06?12:27??surf\OpenSurf.m

?????文件?????641272??2010-11-08?17:14??surf\opensurf.png

?????文件????????603??2010-07-26?16:56??surf\SubFunctions\FastHessian_BuildDerivative.m

?????文件???????1699??2010-07-29?12:10??surf\SubFunctions\FastHessian_buildResponselayer.m

?????文件???????1877??2010-07-20?12:37??surf\SubFunctions\FastHessian_buildResponseMap.m

?????文件???????2353??2010-07-30?11:43??surf\SubFunctions\FastHessian_getIpoints.m

?????文件????????455??2010-07-26?16:56??surf\SubFunctions\FastHessian_getLaplacian.m

?????文件????????450??2010-07-26?16:56??surf\SubFunctions\FastHessian_getResponse.m

?????文件???????2478??2010-07-30?11:02??surf\SubFunctions\FastHessian_interpolateExtremum.m

?????文件???????1680??2010-07-30?10:56??surf\SubFunctions\FastHessian_isExtremum.m

?????文件????????738??2010-07-29?16:32??surf\SubFunctions\FastHessian_Responselayer.m

?????文件????????791??2010-07-29?10:58??surf\SubFunctions\IntegralImage_BoxIntegral.m

?????文件????????494??2010-07-29?16:32??surf\SubFunctions\IntegralImage_HaarX.m

?????文件????????556??2010-07-29?16:37??surf\SubFunctions\IntegralImage_HaarY.m

?????文件????????796??2010-07-29?10:42??surf\SubFunctions\IntegralImage_IntegralImage.m

?????文件???????1089??2010-07-30?11:39??surf\SubFunctions\PaintSURF.m

?????文件???????1554??2010-07-30?12:39??surf\SubFunctions\SurfDescriptor_DecribeInterestPoints.m

?????文件???????3636??2010-07-30?12:40??surf\SubFunctions\SurfDescriptor_GetDescriptor.m

?????文件???????2921??2010-07-30?12:07??surf\SubFunctions\SurfDescriptor_GetOrientation.m

?????文件????2902447??2012-03-30?09:38??surf\surf.pdf

?????文件?????163573??2011-08-27?17:20??surf\TestImages\1.tif

?????文件?????156898??2011-09-06?09:32??surf\TestImages\1wallis.tif

?????文件?????158580??2011-09-06?09:32??surf\TestImages\1wallis_referrence.tif

?????文件?????139826??2011-09-02?16:57??surf\TestImages\1_contrast?enhancement.tif

?????文件?????163573??2011-08-27?17:21??surf\TestImages\1_referrence.tif

?????文件?????152606??2011-09-06?09:36??surf\TestImages\2wallis.tif

?????文件?????141124??2011-09-06?09:42??surf\TestImages\2wallis_referrence.tif

?????文件?????239755??2010-09-06?12:22??surf\TestImages\lena1.png

............此處省略14個文件信息