vtcaddmeanpoint.m

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

function newNoOfVolumes = vtcAddMeanPoints(path, vtc, noOfPoints, location);

% syntax: vtcAddMeanPoints(path, vtc, noOfPoints, location);
% default values: noOfPoints =1, location = 'end'
%
% this func adds 'noOfPoints' points with the timecourse mean
% in the 'location' ('beginning'/'end' ) of each timecourse
% 
% This function was written by:
% Hagar Gelbard
% Rafi Malach's Lab
% Weizmann Institute of Science
% Rehovot, Israel
% [e-mail]

if nargin < 4, location = 'end'; end
if nargin < 3, noOfPoints = 1; end

% open the file for reading  
fid = fopen([path vtc '.vtc'],'r');
% open a new file for writing the cropped file
suffix = ['_' num2str(noOfPoints) 'meanPoints_' location];
nfid = fopen([path vtc suffix '.vtc'],'w');

% READ VTC HEADER
[NoOfVolumes, cols, rows, slices, protocol,TR] = vtcMoveToAfterHeaderPosition(fid);
% WRITE THE HEADER OF THE NEW VTC
newNoOfVolumes = NoOfVolumes + noOfPoints;
vtcWriteHeader(nfid, newNoOfVolumes,protocol,TR);

voxelsNo = rows * cols * slices;

% timecourse manipulation

for voxel = 1:voxelsNo
   % read time course of current voxel
    fseek(fid, 0, 'cof');
    TC = fread(fid, NoOfVolumes, 'uint16');
    if strcmp(location,'beginning')
        newTC = [ones(noOfPoints,1)*mean(TC); TC];
    else
        newTC = [TC; ones(noOfPoints,1)*mean(TC)];
    end
    fwrite(nfid, newTC, 'uint16');
end    
fclose('all');