vtcmovetoafterheaderposition.m

来自「VTC file is the file format of the Brain」· M 代码 · 共 64 行

function [NoOfVolumes, cols, rows, slices, protocol, TR] = vtcMoveToAfterHeaderPosition(fid);
% This function was written by:
% Hagar Gelbard
% Rafi Malach's Lab
% Weizmann Institute of Sceince
% Rehovot, Israel
% [e-mail]
% 
% syntax: [NoOfVolumes, cols, rows, slices, protocol, TR] = vtcMoveToAfterHeaderPosition(fid);
% 
% this func moves the given vtc fid to the after header position
% (i.e., where the data starts)
% it returns NoOfVolumes, cols, rows & slices, which are frequently used
% parameters
% 
% This function was written by:
% Hagar Gelbard
% Rafi Malach's Lab
% Weizmann Institute of Science
% Rehovot, Israel
% [e-mail]

% READ HEADER

% read version number
ver = fread(fid,1,'short');
	
% read fmr name
fmr = [];
a = 'a';
while a ~= 0
  a = fread(fid,1,'char');
  fmr = [fmr a];
end

% read protocol name
protocol = [];
a = 'a';
while a ~= 0
  a = fread(fid,1,'char');
  protocol = [protocol a];
end

% read important parameters
NoOfVolumes = fread(fid,1,'short');
VTC_resolution  = fread(fid,1,'short');
XStart = fread(fid,1,'short');
XEnd = fread(fid,1,'short');
YStart = fread(fid,1,'short');
YEnd = fread(fid,1,'short');
ZStart = fread(fid,1,'short');
ZEnd = fread(fid,1,'short');
HemDelay = fread(fid,1,'short');
TR = fread(fid,1,'float');
HemDelta = fread(fid,1,'float');
HemTau = fread(fid,1,'float');
SegmentSize = fread(fid,1,'short');
SegmentOffset = fread(fid,1,'short');

% end of header

cols = (XEnd-XStart)/VTC_resolution;
rows = (YEnd-YStart)/VTC_resolution;
slices = (ZEnd-ZStart)/VTC_resolution;