vtc2vmp.m

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

function vtc2vmp(path, vtc, output);
% syntax: vtc2vmp(path, vtc, output);
% this func gets a vtc with 1 value at each voxel (1 point TC), and make a
% VMP from it (actually write this value (normalized between 0-1) in a
% correlation map format)
% 
% This function was written by:
% Hagar Gelbard
% Rafi Malach's Lab
% Weizmann Institute of Science
% Rehovot, Israel
% [e-mail]

datestr(now)

% open the vtc file for reading  
fid = fopen([path vtc '.vtc'],'r');
[NoOfVolumes, cols, rows, slices] = vtcMoveToAfterHeaderPosition(fid);

voxelsNo = cols*rows*slices;
vmpVec = zeros(1,voxelsNo);

for k = 1:voxelsNo
       fseek(fid, 0, 'cof');
        vmpVec(k) = fread(fid, 1, 'uint16');
end

vmpVec = vmpVec / max(vmpVec);
writeVMP(path, output, vmpVec, 3, 10 ,0.1,1);

fclose('all');
datestr(now)