surf的matlab代码

SURF角点检测算法是对SIFT的一种改进，主要体现在速度上，效率更高。它和SIFT的主要区别是图像多尺度空间的构建方法不同。作为尺度不变特征变换算法（Sift算法）的加速版，Surf算法在适中的条件下完成两幅图像中物体的匹配基本实现了实时处理，其快速的基础实际上只有一个——积分图像haar求导

% 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个文件信息