读取siftgeo格式文件的matlab程序

% This function reads a siftgeo binary file %读取siftgeo格式的二进制文件

%

% Usage: [v, meta] = siftgeo_read (filename, maxdes)

%   filename    the input filename

%   maxdes      maximum number of descriptors to be loaded 要加载的描述子最大数量

%   (default=unlimited)%默认无限制

%

% Returned values %返回值

%   v           the sift descriptors (1 descriptor per line) %每行一个sift描述子

%   meta        meta data for each descriptor, i.e., per line: %每个描述子的元数据

%               x, y, scale, angle, mi11, mi12, mi21, mi22, cornerness

function [v, meta] = siftgeo_read (filename, maxdes)

if nargin < 2 %nargin:number of arguments input %nargout:number of arguments ouput

  maxdes = 100000000; %1亿

end

% open the file and count the number of descriptors

fid = fopen (filename, 'r');

fseek (fid, 0, 1); %fseek(fid, 0, 'eof'); 跳到文件末尾 %'bof' or -1   Beginning of file;  'cof' or  0   Current position in file;  'eof' or  1   End of file

n = ftell (fid) / (9 * 4 + 1 * 4 + 128); %n是descriptors的数量 %return the current position (number of bytes from the file beginning)

fseek (fid, 0, -1); %fseek(fid, 0, 'bof') 跳到文件开头

if n > maxdes

  n = maxdes;

end;

% first read the meta information associated with the descriptor

meta = zeros (n, 9, 'single'); %n*9矩阵

v = zeros (n, 128, 'single'); %n*128矩阵

% read the elements

for i = 1:n %n是要读区的descriptors数量

  meta(i,:) = fread (fid, 9, 'float'); %(float)*9 元数据信息

  d = fread (fid, 1, 'int'); %(int)*1

  v(i,:) = fread (fid, d, 'uint8=>single'); %(uint8=>single)*d 描述子

end

fclose (fid);