matlab surf图像配准 图像检索

% Example 3 Affine registration
% Load images
I1=im2double（imread（‘TestImages/lena1.png‘））;
I2=im2double（imread（‘TestImages/lena2.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）;

% Make vectors with the coordinates of the best matches
Pos1=[[Ipts1（cor1）.y]‘[Ipts1（cor1）.x]‘];
Pos2=[[Ipts2（cor2）.y]‘[Ipts2（cor2）.x]‘];
Pos1=Pos1（1:30:）;
Pos2=Pos2（1:30:）;

% 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）; hold on;

% Show the best matches
plot（[Pos1（:2） Pos2（:2）+size（I12）]‘[Pos1（:1） Pos2（:1）]‘‘-‘）;
plot（[Pos1（:2） Pos2（:2）+size（I12）]‘[Pos1（:1） Pos2（:1）]‘‘o‘）;

% Calculate affine matrix
Pos1（:3）=1; Pos2（:3）=1;
M=Pos1‘/Pos2‘;

% Add subfunctions to Matlab Search path
functionname=‘OpenSurf.m‘;
functiondir=which（functionname）;
functiondir=functiondir（1:end-length（functionname））;
addpath（[functiondir ‘/WarpFunctions‘]）

% Warp the image
I1_warped=affine_warp（I1M‘bicubic‘）;

% Show the result
figure
subplot（131） imshow（I1）;title（‘Figure 1‘）;
subplot（132） imshow（I2）;title（‘Figure 2‘）;
subplot（133） imshow（I1_warped）;title（‘Warped Figure 1‘）;

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

     文件       1776  2012-05-29 22:53  surfexample3.m

     文件       1335  2010-09-06 07:21  surflicense.txt

     文件       3778  2010-09-06 12:27  surfOpenSurf.m

     文件     641272  2010-11-08 17:14  surfopensurf.png

     文件        603  2010-07-26 16:56  surfSubFunctionsFastHessian_BuildDerivative.m

     文件       1699  2010-07-29 12:10  surfSubFunctionsFastHessian_buildResponselayer.m

     文件       1877  2010-07-20 12:37  surfSubFunctionsFastHessian_buildResponseMap.m

     文件       2353  2010-07-30 11:43  surfSubFunctionsFastHessian_getIpoints.m

     文件        455  2010-07-26 16:56  surfSubFunctionsFastHessian_getLaplacian.m

     文件        450  2010-07-26 16:56  surfSubFunctionsFastHessian_getResponse.m

     文件       2478  2010-07-30 11:02  surfSubFunctionsFastHessian_interpolateExtremum.m

     文件       1680  2010-07-30 10:56  surfSubFunctionsFastHessian_isExtremum.m

     文件        738  2010-07-29 16:32  surfSubFunctionsFastHessian_Responselayer.m

     文件        791  2010-07-29 10:58  surfSubFunctionsIntegralImage_BoxIntegral.m

     文件        494  2010-07-29 16:32  surfSubFunctionsIntegralImage_HaarX.m

     文件        556  2010-07-29 16:37  surfSubFunctionsIntegralImage_HaarY.m

     文件        796  2010-07-29 10:42  surfSubFunctionsIntegralImage_IntegralImage.m

     文件       1089  2010-07-30 11:39  surfSubFunctionsPaintSURF.m

     文件       1554  2010-07-30 12:39  surfSubFunctionsSurfDescriptor_DecribeInterestPoints.m

     文件       3636  2010-07-30 12:40  surfSubFunctionsSurfDescriptor_GetDescriptor.m

     文件       2921  2010-07-30 12:07  surfSubFunctionsSurfDescriptor_GetOrientation.m

     文件    2902447  2012-03-30 09:38  surfsurf.pdf

     文件     163573  2011-08-27 17:20  surfTestImages1.tif

     文件     156898  2011-09-06 09:32  surfTestImages1wallis.tif

     文件     158580  2011-09-06 09:32  surfTestImages1wallis_referrence.tif

     文件     139826  2011-09-02 16:57  surfTestImages1_contrast enhancement.tif

     文件     163573  2011-08-27 17:21  surfTestImages1_referrence.tif

     文件     152606  2011-09-06 09:36  surfTestImages2wallis.tif

     文件     141124  2011-09-06 09:42  surfTestImages2wallis_referrence.tif

     文件     239755  2010-09-06 12:22  surfTestImageslena1.png

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