%?Example?2?Corresponding?points
%?Load?images
I1?=?imread（‘D:\image\SPOT\SPOT_low_3.tif‘）;
I2?=?imread（‘D:\image\SPOT\SPOT_low_4.tif‘）;

imshow（I1）;
%?Get?the?Key?Points
??Options.upright=true;
??Options.tresh=0.0001;
%?Options.verbose?=?true;
??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）;
%?Show?both?images
??I?=?zeros（[size（I11）?size（I12）+size（I22）?size（I13）]）;
??I（:1:size（I12）:）=I1;?I（1:size（I21）size（I12）+1:size（I12）+size（I22）:）=I2;
??figureimshow（uint8（I））;?hold?on;
%?Show?the?best?matches
??for?i=1:30
??????c=rand（13）;
??????plot（[Ipts1（cor1（i））.x?Ipts2（cor2（i））.x+size（I12）][Ipts1（cor1（i））.y?Ipts2（cor2（i））.y]‘-‘‘Color‘c）
??????plot（[Ipts1（cor1（i））.x?Ipts2（cor2（i））.x+size（I12）][Ipts1（cor1（i））.y?Ipts2（cor2（i））.y]‘o‘‘Color‘c）
??????k1（i?:）?=?[Ipts1（cor1（i））.x?Ipts1（cor1（i））.y];
??????k2（i?:）?=?[Ipts2（cor2（i））.x?Ipts2（cor2（i））.y];
??end

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

?????文件??????19062??2016-03-28?14:22??新建文件夾\1.jpg

?????文件??????39760??2016-04-05?14:21??新建文件夾\2.jpg

?????文件???????1401??2014-10-23?14:27??新建文件夾\example2.m

?????文件???????1547??2016-04-05?16:45??新建文件夾\example3.m

?????文件??????19782??2014-09-16?09:16??新建文件夾\hs_err_pid1348.log

?????文件??????19276??2014-11-06?14:30??新建文件夾\hs_err_pid6024.log

?????文件??????19947??2014-11-13?16:12??新建文件夾\hs_err_pid7852.log

?????文件???????1335??2010-09-06?19:21??新建文件夾\license.txt

?????文件????????655??2014-10-14?18:34??新建文件夾\myexample.m

?????文件???????3778??2010-09-07?00:27??新建文件夾\OpenSurf.m

?????文件????????603??2010-07-27?04:56??新建文件夾\SubFunctions\FastHessian_BuildDerivative.m

?????文件???????1765??2014-09-22?14:18??新建文件夾\SubFunctions\FastHessian_buildResponselayer.m

?????文件???????1892??2014-09-17?16:20??新建文件夾\SubFunctions\FastHessian_buildResponseMap.m

?????文件???????2431??2014-10-14?10:02??新建文件夾\SubFunctions\FastHessian_getIpoints.m

?????文件????????455??2010-07-27?04:56??新建文件夾\SubFunctions\FastHessian_getLaplacian.m

?????文件????????450??2010-07-27?04:56??新建文件夾\SubFunctions\FastHessian_getResponse.m

?????文件???????2493??2014-09-22?15:30??新建文件夾\SubFunctions\FastHessian_interpolateExtremum.m

?????文件???????1701??2014-09-16?16:59??新建文件夾\SubFunctions\FastHessian_isExtremum.m

?????文件????????738??2010-07-30?04:32??新建文件夾\SubFunctions\FastHessian_Responselayer.m

?????文件?????198316??2014-09-11?15:36??新建文件夾\SubFunctions\img_1.tif

?????文件?????201156??2014-09-11?15:36??新建文件夾\SubFunctions\img_2.tif

?????文件????????820??2014-09-19?14:27??新建文件夾\SubFunctions\IntegralImage_BoxIntegral.m

?????文件????????494??2010-07-30?04:32??新建文件夾\SubFunctions\IntegralImage_HaarX.m

?????文件????????556??2010-07-30?04:37??新建文件夾\SubFunctions\IntegralImage_HaarY.m

?????文件????????796??2010-07-29?22:42??新建文件夾\SubFunctions\IntegralImage_IntegralImage.m

?????文件???????1089??2010-07-30?23:39??新建文件夾\SubFunctions\PaintSURF.m

?????文件???????1554??2014-10-14?09:11??新建文件夾\SubFunctions\SurfDescriptor_DecribeInterestPoints.m

?????文件???????3636??2014-09-23?19:45??新建文件夾\SubFunctions\SurfDescriptor_GetDescriptor.m

?????文件???????2921??2010-07-31?00:07??新建文件夾\SubFunctions\SurfDescriptor_GetOrientation.m

?????文件?????239755??2010-09-07?00:22??新建文件夾\TestImages\lena1.png

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