SURF算法代码 matlab

可以实现SURF特征点的提取和匹配运算，经过测试可以运行。供给大家分享学习。

% Example 2 Corresponding points
% Load images
  I1=imread（‘TestImages/testc1.png‘）;
  I2=imread（‘TestImages/testc2.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）;
% 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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1211  2014-06-04 16:50  SURF\example2.m

     文件       1776  2014-06-04 16:50  SURF\example3.m

     文件       1335  2014-06-04 16:50  SURF\license.txt

     文件       3778  2014-06-04 16:50  SURF\OpenSurf.m

     文件        603  2014-06-04 16:50  SURF\SubFunctions\FastHessian_BuildDerivative.m

     文件       1699  2014-06-04 16:50  SURF\SubFunctions\FastHessian_buildResponselayer.m

     文件       1877  2014-06-04 16:50  SURF\SubFunctions\FastHessian_buildResponseMap.m

     文件       2353  2014-06-04 16:50  SURF\SubFunctions\FastHessian_getIpoints.m

     文件        455  2014-06-04 16:50  SURF\SubFunctions\FastHessian_getLaplacian.m

     文件        450  2014-06-04 16:50  SURF\SubFunctions\FastHessian_getResponse.m

     文件       2478  2014-06-04 16:50  SURF\SubFunctions\FastHessian_interpolateExtremum.m

     文件       1680  2014-06-04 16:50  SURF\SubFunctions\FastHessian_isExtremum.m

     文件        738  2014-06-04 16:50  SURF\SubFunctions\FastHessian_Responselayer.m

     文件        791  2014-06-04 16:50  SURF\SubFunctions\IntegralImage_BoxIntegral.m

     文件        494  2014-06-04 16:50  SURF\SubFunctions\IntegralImage_HaarX.m

     文件        556  2014-06-04 16:50  SURF\SubFunctions\IntegralImage_HaarY.m

     文件        796  2014-06-04 16:50  SURF\SubFunctions\IntegralImage_IntegralImage.m

     文件       1089  2014-06-04 16:50  SURF\SubFunctions\PaintSURF.m

     文件       1554  2014-06-04 16:50  SURF\SubFunctions\SurfDescriptor_DecribeInterestPoints.m

     文件       3636  2014-06-04 16:50  SURF\SubFunctions\SurfDescriptor_GetDescriptor.m

     文件       2921  2014-06-04 16:50  SURF\SubFunctions\SurfDescriptor_GetOrientation.m

     文件     574557  2014-06-04 16:50  SURF\TestImages\lena1.png

     文件     578897  2014-06-04 16:50  SURF\TestImages\lena2.png

     文件       5296  2014-06-04 16:50  SURF\TestImages\test.png

     文件     123228  2014-06-04 16:50  SURF\TestImages\testc1.png

     文件     122063  2014-06-04 16:50  SURF\TestImages\testc2.png

     文件       9721  2014-06-04 16:50  SURF\WarpFunctions\affine_warp.m

     目录          0  2014-06-04 16:50  SURF\SubFunctions

     目录          0  2014-06-04 16:50  SURF\TestImages

     目录          0  2014-06-04 16:50  SURF\WarpFunctions

............此处省略4个文件信息