MATLAB SIFT 代码

matlab 实现的 sift 变换 的代码，包含整个过程的详细步骤。

function dss = diffss（ss）
% DIFFSS  Difference of scale space
%   DSS=DIFFSS（SS） returns a scale space DSS obtained by subtracting
%   consecutive levels of the scale space SS.
%
%   In SIFT this function is used to compute the difference of
%   Gaussian scale space from the Gaussian scale space of an image.
%
%   See also GAUSSIANSS（） PDF:SIFT.USER.SS.

% AUTORIGHTS
% Copyright （c） 2006 The Regents of the University of California.
% All Rights Reserved.
% 
% Created by Andrea Vedaldi
% UCLA Vision Lab - Department of Computer Science
% 
% Permission to use copy modify and distribute this software and its
% documentation for educational research and non-profit purposes
% without fee and without a written agreement is hereby granted
% provided that the above copyright notice this paragraph and the
% following three paragraphs appear in all copies.
% 
% This software program and documentation are copyrighted by The Regents
% of the University of California. The software program and
% documentation are supplied “as is“ without any accompanying services
% from The Regents. The Regents does not warrant that the operation of
% the program will be uninterrupted or error-free. The end-user
% understands that the program was developed for research purposes and
% is advised not to rely exclusively on the program for any reason.
% 
% This software embodies a method for which the following patent has
% been issued: “Method and apparatus for identifying scale invariant
% features in an image and use of same for locating an object in an
% image“ David G. Lowe US Patent 6711293 （March 23
% 2004）. Provisional application filed March 8 1999. Asignee: The
% University of British Columbia.
% 
% IN NO EVENT SHALL THE UNIVERSITY OF CALIFORNIA BE LIABLE TO ANY PARTY
% FOR DIRECT INDIRECT SPECIAL INCIDENTAL OR CONSEQUENTIAL DAMAGES
% INCLUDING LOST PROFITS ARISING OUT OF THE USE OF THIS SOFTWARE AND
% ITS DOCUMENTATION EVEN IF THE UNIVERSITY OF CALIFORNIA HAS BEEN
% ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. THE UNIVERSITY OF
% CALIFORNIA SPECIFICALLY DISCLAIMS ANY WARRANTIES INCLUDING BUT NOT
% LIMITED TO THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
% A PARTICULAR PURPOSE. THE SOFTWARE PROVIDED HEREUNDER IS ON AN “AS IS“
% BASIS AND THE UNIVERSITY OF CALIFORNIA HAS NO OBLIGATIONS TO PROVIDE
% MAINTENANCE SUPPORT UPDATES ENHANCEMENTS OR MODIFICATIONS.

dss.smin = ss.smin ;
dss.smax = ss.smax-1 ;
dss.omin = ss.omin ;
dss.O = ss.O ;
dss.S = ss.S ;
dss.sigma0 = ss.sigma0 ;

for o=1:dss.O
  % Can be done at once but it turns out to be faster
  % this way
  [MNS] = size（ss.octave{o}） ;
  dss.octave{o} = zeros（MNS-1） ;
  for s=1:S-1
    dss.octave{o}（::s） = ...
        ss.octave{o}（::s+1） - ss.octave{o}（::s） ;
  end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2006-08-04 10:29  data\
     文件      464954  2006-04-05 04:45  data\img3.jpg
     文件     1420810  2006-04-05 04:45  data\img3.sift
     文件      467808  2006-04-05 04:45  data\img5.jpg
     文件     1585635  2006-04-05 04:45  data\img5.sift
     文件      275111  2006-05-01 19:51  data\landscape-a.jpg
     文件      304723  2006-05-01 19:51  data\landscape-b.jpg
     文件        7620  2006-04-05 04:45  data\test.png
     文件       34162  2006-04-05 04:45  data\test.sift
     文件        2797  2006-05-12 21:43  diffss.m
     目录           0  2006-08-04 12:59  doc\
     文件         246  2006-04-11 07:37  doc\bibliography.bib
     目录           0  2006-04-27 23:03  doc\figures\
     文件         550  2006-04-27 23:02  doc\figures\sift-descriptor-AI.bb
     文件      475918  2006-04-27 23:02  doc\figures\sift-descriptor-AI.pdf
     文件        3275  2006-04-27 23:00  doc\figures\sift-descriptor-SAVED.tex
     文件           8  2006-04-27 23:03  doc\figures\sift-descriptor.aux
     文件       13430  2006-04-27 23:03  doc\figures\sift-descriptor.log
     文件      443880  2006-04-27 23:03  doc\figures\sift-descriptor.pdf
     文件        3275  2006-04-27 23:03  doc\figures\sift-descriptor.tex
     文件          13  2006-04-27 23:03  doc\figures\sift-descriptor.wrm
     文件       48557  2001-07-05 02:24  doc\figures\warmread.sty
     文件       21807  2006-05-13 07:16  doc\sift.tex
     文件        6174  2006-04-11 07:39  doc\visionlab.sty
     文件        7920  2006-05-12 21:44  gaussianss.m
     文件        2359  2006-05-12 21:45  imreadbw.m
     文件        2994  2006-05-01 22:36  imsmooth.c
     文件        1990  2006-04-05 08:42  LICENSE
     文件        2117  2006-08-04 10:31  Makefile
     文件        2178  2006-05-01 22:36  mexutils.c
     文件        4360  2006-05-18 17:36  plotmatches.m
............此处省略24个文件信息