% 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, cornernessfunction [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 descriptorsfid = 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 filen = 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 descriptormeta = zeros (n, 9, 'single'); %n*9矩阵v = zeros (n, 128, 'single'); %n*128矩阵% read the elementsfor 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 描述子endfclose (fid);