surf算法对图像的配准实时的取证

surf算法对图像的配准实时的取证，经本人验证，该程序好用无误，可放心使用

% Example 2 Corresponding points
% Load images
clc
clear all
close all
  I1=imread（‘TestImages/abcd.jpg‘）;
  I2=imread（‘TestImages/abc.jpg‘）;
% Get the Key Points
  Options.upright=true;
  Options.tresh=0.0001;
  tic
  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）*2 size（I13）]）;
  I（:1:size（I12）:）=I1; I（:size（I12）+1:size（I12）+size（I22）:）=I2;
  figure imshow（I/255）; 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）
  end
  fprintf（‘global in %f \n‘toc）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-05-19 15:29  surf\
     文件        1804  2013-05-19 15:54  surf\K-mean.asv
     文件        1771  2013-05-18 11:41  surf\K-mean.m
     文件     1362910  2013-05-16 16:08  surf\OpenSURF_version1c.zip
     文件        3778  2010-09-07 00:27  surf\OpenSurf.m
     目录           0  2013-05-16 16:09  surf\SubFunctions\
     文件         603  2010-07-27 04:56  surf\SubFunctions\FastHessian_BuildDerivative.m
     文件         738  2010-07-30 04:32  surf\SubFunctions\FastHessian_Responselayer.m
     文件        1699  2010-07-30 00:10  surf\SubFunctions\FastHessian_buildResponselayer.m
     文件        1877  2010-07-21 00:37  surf\SubFunctions\FastHessian_buildResponseMap.m
     文件        2353  2010-07-30 23:43  surf\SubFunctions\FastHessian_getIpoints.m
     文件         455  2010-07-27 04:56  surf\SubFunctions\FastHessian_getLaplacian.m
     文件         450  2010-07-27 04:56  surf\SubFunctions\FastHessian_getResponse.m
     文件        2478  2010-07-30 23:02  surf\SubFunctions\FastHessian_interpolateExtremum.m
     文件        1680  2010-07-30 22:56  surf\SubFunctions\FastHessian_isExtremum.m
     文件         791  2010-07-29 22:58  surf\SubFunctions\IntegralImage_BoxIntegral.m
     文件         494  2010-07-30 04:32  surf\SubFunctions\IntegralImage_HaarX.m
     文件         556  2010-07-30 04:37  surf\SubFunctions\IntegralImage_HaarY.m
     文件         796  2010-07-29 22:42  surf\SubFunctions\IntegralImage_IntegralImage.m
     文件        1089  2010-07-30 23:39  surf\SubFunctions\PaintSURF.m
     文件        1554  2010-07-31 00:39  surf\SubFunctions\SurfDescriptor_DecribeInterestPoints.m
     文件        3636  2010-07-31 00:40  surf\SubFunctions\SurfDescriptor_GetDescriptor.m
     文件        2921  2010-07-31 00:07  surf\SubFunctions\SurfDescriptor_GetOrientation.m
     目录           0  2013-05-18 09:46  surf\TestImages\
     文件      409708  2013-05-14 15:18  surf\TestImages\a.jpg
     文件      410556  2013-05-16 15:26  surf\TestImages\a1.jpg
     文件      422268  2013-05-15 20:08  surf\TestImages\aa.jpg
     文件      400870  2013-05-14 18:42  surf\TestImages\abc.jpg
     文件      400870  2013-05-14 18:42  surf\TestImages\abcd.jpg
     文件      319686  2013-05-09 18:40  surf\TestImages\cc.jpg
     文件     1488221  2013-05-15 11:01  surf\TestImages\d.jpg
............此处省略21个文件信息