clear;
close all;
clc;

load('diffusion_data');

%% Visualize values of a single voxel
scatter3(bvecs(1,:), bvecs(2,:), bvecs(3,:), 50, (1-t), 'filled');


%% DTI Reconstruction
for i=1:numel(bvals),
    DTIdata(i).VoxelData = single(data(:,:,:,i));
    DTIdata(i).Gradient = bvecs(:,i);
    DTIdata(i).Bvalue = bvals(i);
end

DTens = dti_recon(DTIdata, 223);






