vtcwriteheader.m

VTC file is the file format of the Brainvoyager. The FMRI matchine scan human s brain and get the vt

function [cols, rows, slices] = vtcWriteHeader(fid, NoOfVolumes, protocol, TR, VTC_resolution, ver, XStart, XEnd,...
                                                                                YStart, YEnd, ZStart, ZEnd, HemDelay, HemDelta, HemTau,...
                                                                                SegmentSize, SegmentOffset, fmr);
% syntax: vtcWriteHeader(fid, NoOfVolumes, protocol, TR, VTC_resolution, ver, XStart, XEnd,...
%                                                                                 YStart, YEnd, ZStart, ZEnd, HemDelay, HemDelta, HemTau,...
%                                                                                 SegmentSize, SegmentOffset, fmr);
% default values: protocol = '<none>'; TR = 3000; VTC_resolution = 3; ver =
% 2; XStart = 57; XEnd = 231; YStart = 52; YEnd = 172; ZStart = 59; ZEnd
% =197; HemDelay = 1; HemDelta = 2.5000; HemTau = 1.2500; SegmentSize = 10;
% SegmentOffset = 0; fmr = 'dummy.fmr';
% 
% must specify: fid and NoOfVolumes
% 
% This function was written by:
% Hagar Gelbard
% Rafi Malach's Lab
% Weizmann Institute of Science
% Rehovot, Israel
% [e-mail]

if nargin < 18, fmr = ['dummy.fmr' 0]; end                                                                                                                                                        
if nargin < 17, SegmentOffset = 0; end                                                                            
if nargin < 16, SegmentSize = 10; end
if nargin < 15, HemTau = 1.2500; end
if nargin < 14, HemDelta = 2.5000; end
if nargin < 13, HemDelay = 1; end
if nargin < 12, ZEnd =197; end
if nargin < 11, ZStart = 59; end
if nargin < 10, YEnd = 172; end
if nargin < 9, YStart = 52; end
if nargin < 8, XEnd = 231; end
if nargin < 7, XStart = 57; end
if nargin < 6, ver = 2; end
if nargin < 5, VTC_resolution = 3; end
if nargin < 4, TR = 3000; end
if nargin < 3, protocol = ['<none>' 0]; end

fwrite(fid, ver, 'short');
fwrite(fid, char(fmr), 'char');
fwrite(fid, char(protocol), 'char');
fwrite(fid, NoOfVolumes, 'short');
fwrite(fid, VTC_resolution, 'short');
fwrite(fid, XStart, 'short');
fwrite(fid, XEnd, 'short');
fwrite(fid, YStart, 'short');
fwrite(fid, YEnd, 'short');
fwrite(fid, ZStart, 'short');
fwrite(fid, ZEnd, 'short');
fwrite(fid, HemDelay, 'short');
fwrite(fid, TR, 'float');
fwrite(fid, HemDelta, 'float');
fwrite(fid, HemTau, 'float');
fwrite(fid, SegmentSize, 'short');
fwrite(fid, SegmentOffset, 'short');

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