matlab实现sift提取特征点

SIFT算法的相关实现 自己用matlab写的sift特征点检测与匹配程序，运行能通过

function [databasehandle] = add_descriptors_to_database（ im pos scale orient desc existing_database ）

% [database] = add_descriptors_to_database（ im pos orient scale desc existing_database ）
%
% Add an image and its descriptors to a database or create a new database if
% one is not specified.
% 
% Input:
% im - the image.
% pos orient scale desc - keypoint information from the SIFT function.
% existing_database - database to add the descriptors to.  If omitted a 
%   new database is created.
%
% Output:
% database - the database structure.
% handle - the integer handle of the descriptor block in the database.
%
% Thomas F. El-Maraghi
% May 2004

if exist（‘existing_database‘‘var‘）
   database = existing_database;
else
   database = empty_descriptor_database;
end

database.num_im = database.num_im + 1;
database.im{database.num_im} = im;
database.index = [database.index; database.num_im*ones（size（pos1）1）];
if nargout == 2
   handle = database.num_im;
end

database.pos = [database.pos; pos];
database.scale = [database.scale; scale];
database.orient = [database.orient; orient];
database.desc = [database.desc; desc];

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

     文件       1148  2009-01-06 21:23  SIFT\add_descriptors_to_database.m

     文件       2546  2008-12-25 10:34  SIFT\display_keypoints.m

     文件        328  2008-12-25 10:34  SIFT\empty_descriptor_database.m

     文件       1131  2009-01-06 21:24  SIFT\find_nearest_neighbours.m

     文件       3578  2009-01-06 21:25  SIFT\fit_robust_affine_transform.m

     文件        642  2009-01-06 21:22  SIFT\gaussian_filter.m

     文件        442  2008-12-25 10:34  SIFT\get_linear_index.m

     文件        588  2008-12-25 10:34  SIFT\get_multiple_indices.m

     文件       7690  2009-01-06 21:26  SIFT\hough.m

     文件       2505  2009-01-06 21:26  SIFT\imWarpAffine.m

     文件       1573  2009-01-07 10:09  SIFT\pgmread.m

     文件       1189  2009-01-06 21:27  SIFT\read_keypoint_file.m

     文件       1662  2009-01-06 21:28  SIFT\read_lowe_keypoints.m

     文件        492  2008-12-25 10:34  SIFT\resizeImageFig.m

     文件      30545  2009-03-05 20:47  SIFT\SIFT.m

     文件       2703  2008-12-25 10:34  SIFT\SIFT_from_cache.m

     文件      12343  2009-03-05 21:32  SIFT\tutorial.m

     文件        856  2008-12-25 10:34  SIFT\write_keypoint_file.m

     目录          0  2009-05-12 00:04  SIFT

----------- ---------  ---------- -----  ----

                71961                    19