Matlab实现SIFT三幅图像拼接

使用Matlab实现寻找出SIFT特征点，然后进行匹配，最终完成三幅图像拼接。

function n2 = dist2（x c）
%DIST2	Calculates squared distance between two sets of points.
%
%	Description
%	D = DIST2（X C） takes two matrices of vectors and calculates the
%	squared Euclidean distance between them.  Both matrices must be of
%	the same column dimension.  If X has M rows and N columns and C has
%	L rows and N columns then the result has M rows and L columns.  The
%	I Jth entry is the  squared distance from the Ith row of X to the
%	Jth row of C.
%
%	See also
%	GMMACTIV KMEANS RBFFWD
%

%	Copyright （c） Christopher M Bishop Ian T Nabney （1996 1997）

[ndata dimx] = size（x）;
[ncentres dimc] = size（c）;
if dimx ~= dimc
	error（‘Data dimension does not match dimension of centres‘）
end

n2 = （ones（ncentres 1） * sum（（x.^2）‘ 1））‘ + ...
  		ones（ndata 1） * sum（（c.^2）‘1） - ...
  		2.*（x*（c‘））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2010-07-25 00:31  mosaic\
     文件        6148  2010-07-25 00:30  mosaic\.DS_Store
     目录           0  2010-07-25 00:32  __MACOSX\
     目录           0  2010-07-25 00:32  __MACOSX\mosaic\
     文件          82  2010-07-25 00:30  __MACOSX\mosaic\._.DS_Store
     文件         846  2008-10-26 12:21  mosaic\dist2.m
     文件        1274  2010-07-24 23:48  mosaic\ext_desc.m
     文件        1559  2010-07-24 20:28  mosaic\harris.m
     文件        1069  2008-10-26 12:30  mosaic\hnormalise.m
     文件        2493  2008-10-26 12:30  mosaic\homography2d.m
     文件        2380  2008-10-26 12:29  mosaic\iscolinear.m
     文件       69455  2007-02-26 22:39  mosaic\keble_a.jpg
     文件       69819  2007-02-26 22:39  mosaic\keble_b.jpg
     文件       68821  2007-02-26 22:39  mosaic\keble_c.jpg
     文件        4314  2010-07-24 23:55  mosaic\mosaic.m
     文件      372769  2010-07-24 23:41  mosaic\mosaic.pdf
     文件          82  2010-07-24 23:41  __MACOSX\mosaic\._mosaic.pdf
     文件        2243  2008-10-26 12:26  mosaic\normalise2dpts.m
     目录           0  2009-11-08 20:13  mosaic\paris\
     文件      181033  2009-11-05 13:09  mosaic\paris\paris_a.jpg
     目录           0  2010-07-25 00:32  __MACOSX\mosaic\paris\
     文件         201  2009-11-05 13:09  __MACOSX\mosaic\paris\._paris_a.jpg
     文件      180911  2009-11-05 13:09  mosaic\paris\paris_b.jpg
     文件         201  2009-11-05 13:09  __MACOSX\mosaic\paris\._paris_b.jpg
     文件      177405  2009-11-05 13:08  mosaic\paris\paris_c.jpg
     文件         201  2009-11-05 13:08  __MACOSX\mosaic\paris\._paris_c.jpg
     文件         201  2009-11-08 20:13  __MACOSX\mosaic\._paris
     文件        2282  2008-10-26 12:29  mosaic\randomsample.m
     文件        9953  2008-10-26 12:29  mosaic\ransac.m
     文件        5050  2008-10-26 12:19  mosaic\ransacfithomography.m
     文件        4538  2008-11-14 14:45  mosaic\vgg_gui_H.m
............此处省略1个文件信息