讀取siftgeo格式檔案的matlab程式_siftgeo

% 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);