sift sift特征提取

sift特征提取，包括三部分，特征点提取，特征点描述，特征点匹配

function descriptors = do_descriptor（octave oframes  sigma0 Ssminvarargin）

for k=1:2:length（varargin）
	switch lower（varargin{k}）
      case ‘magnif‘
        magnif = varargin{k+1} ;      
      case ‘numspatialbins‘
        NBP = varargin{k+1} ;       
      case  ‘numorientbins‘
        NBO = varargin{k+1} ;   
      otherwise
        error（[‘Unknown parameter ‘ varargin{k} ‘.‘]） ;
     end
end 
                                   
num_spacialBins = NBP;
num_orientBins = NBO;
key_num = size（oframes 2）;
% 计算图像梯度 
[M N s_num] = size（octave）; 
descriptors = [];
magnitudes = zeros（M N s_num）;
angles = zeros（M N s_num）;

for si = 1: s_num
    img = octave（::si）;
    dx_filter = [-0.5 0 0.5];
    dy_filter = dx_filter‘;
    gradient_x = imfilter（img dx_filter）;
    gradient_y = i

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

     文件     506173  2011-11-23 12:30  sift\1.jpg

     文件     263222  2011-11-23 12:57  sift\2.bmp

     文件       4156  2010-05-20 13:02  sift\do_descriptor.m

     文件        323  2010-05-20 12:30  sift\do_diffofg.m

     文件       3278  2010-05-20 12:51  sift\do_extrefine.m

     文件       2365  2010-05-20 12:30  sift\do_gaussian.m

     文件       1181  2010-05-20 12:41  sift\do_localmax.m

     文件       1753  2010-05-20 12:57  sift\do_orientation.m

     文件       2477  2011-12-07 22:09  sift\do_sift.m

     文件        101  2010-05-21 10:53  sift\smooth.m

     文件     842240  2010-05-23 13:36  sift\片段.shs

     目录          0  2011-12-07 22:03  sift

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

              1627269                    12