clear all
addpath(genpath(pwd))
dirpath = uigetdir('./','Enter directory with the 2dseq file');
dseqpath=[dirpath '/2dseq'];
visupath=[dirpath '/visu_pars'];

%% load separate reconstructed image:
obj=ImageDataObject(dseqpath, visupath);
obj=ImageDataObject(dseqpath, visupath);
DataSize = size(obj.data);

%uloženie pôvo
for l = 1:DataSize(end)
        for j = 1:DataSize(end-1)
            data1{1,1}{j,l} = obj.data(:,:,1,1,j,l);
end
end
Z=data1{1,1};

%uloženie TI časov do premennej
for i=1:l
    TI1(i)=obj.FrameGroupParams(1,i).VisuAcqInversionTime;
    TI = [TI1;TI1]; 

save('Z.mat','Z','TI','l')