surf算法matlab源码

赵小川书的代码，直接可以运行，实现特征匹配

% Example 3 Affine registration
% Load images
clear all;

I1=im2double（imread（‘TestImages/lena1.png‘））;
I2=im2double（imread（‘TestImages/lena2.png‘））;
% I1=im2double（imread（‘sky\A1.png‘））;
% I2=im2double（imread（‘sky\A2.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（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-03-19 20:41  OpenSURF_version1c\
     文件        3778  2010-09-07 00:27  OpenSURF_version1c\OpenSurf.m
     目录           0  2013-03-19 19:20  OpenSURF_version1c\SubFunctions\
     文件         603  2010-07-27 04:56  OpenSURF_version1c\SubFunctions\FastHessian_BuildDerivative.m
     文件         738  2010-07-30 04:32  OpenSURF_version1c\SubFunctions\FastHessian_Responselayer.m
     文件        1699  2010-07-30 00:10  OpenSURF_version1c\SubFunctions\FastHessian_buildResponselayer.m
     文件        1877  2010-07-21 00:37  OpenSURF_version1c\SubFunctions\FastHessian_buildResponseMap.m
     文件        2353  2010-07-30 23:43  OpenSURF_version1c\SubFunctions\FastHessian_getIpoints.m
     文件         455  2010-07-27 04:56  OpenSURF_version1c\SubFunctions\FastHessian_getLaplacian.m
     文件         450  2010-07-27 04:56  OpenSURF_version1c\SubFunctions\FastHessian_getResponse.m
     文件        2478  2010-07-30 23:02  OpenSURF_version1c\SubFunctions\FastHessian_interpolateExtremum.m
     文件        1680  2010-07-30 22:56  OpenSURF_version1c\SubFunctions\FastHessian_isExtremum.m
     文件         791  2010-07-29 22:58  OpenSURF_version1c\SubFunctions\IntegralImage_BoxIntegral.m
     文件         494  2010-07-30 04:32  OpenSURF_version1c\SubFunctions\IntegralImage_HaarX.m
     文件         556  2010-07-30 04:37  OpenSURF_version1c\SubFunctions\IntegralImage_HaarY.m
     文件         796  2010-07-29 22:42  OpenSURF_version1c\SubFunctions\IntegralImage_IntegralImage.m
     文件        1089  2010-07-30 23:39  OpenSURF_version1c\SubFunctions\PaintSURF.m
     文件        1554  2010-07-31 00:39  OpenSURF_version1c\SubFunctions\SurfDescriptor_DecribeInterestPoints.m
     文件        3636  2010-07-31 00:40  OpenSURF_version1c\SubFunctions\SurfDescriptor_GetDescriptor.m
     文件        2921  2010-07-31 00:07  OpenSURF_version1c\SubFunctions\SurfDescriptor_GetOrientation.m
     目录           0  2013-03-19 20:24  OpenSURF_version1c\TestImages\
     文件       75288  2011-04-17 06:11  OpenSURF_version1c\TestImages\043.png
     文件       81918  2011-04-17 06:11  OpenSURF_version1c\TestImages\044.png
     文件       41472  2012-10-06 17:08  OpenSURF_version1c\TestImages\Thumbs.db
     文件       15904  2006-01-19 11:35  OpenSURF_version1c\TestImages\fish1.png
     文件       15428  2006-01-19 11:37  OpenSURF_version1c\TestImages\fish2.png
     文件      239755  2010-09-07 00:22  OpenSURF_version1c\TestImages\lena1.png
     文件      211716  2010-09-07 00:22  OpenSURF_version1c\TestImages\lena2.png
     文件       19555  2013-03-19 20:24  OpenSURF_version1c\TestImages\lena20.jpg
     文件      113463  2013-03-19 20:24  OpenSURF_version1c\TestImages\lena21.jpg
     文件        5296  2010-07-20 05:38  OpenSURF_version1c\TestImages\test.png
............此处省略12个文件信息