goFigure.m

load('syntheticDataset_01.mat');

K = 3;

tmp_xy(:,:,1) = squeeze(max(syntheticDataset(1).phantomKymo(:,:,:,1),[],2));
tmp_xy(:,:,2) = squeeze(max(syntheticDataset(11).phantomKymo(:,:,:,1),[],2));
tmp_xy(:,:,3) = squeeze(max(syntheticDataset(21).phantomKymo(:,:,:,1),[],2));

figure('Position',[300,200,900,300]);
layout1 = tiledlayout(1,3,'TileSpacing','compact','Padding','tight');
layout2 = tiledlayout(layout1,1,3,'TileSpacing','tight','Padding','tight');
layout2.Layout.Tile = 1;
layout2.Layout.TileSpan = [1 1];

axVel = nexttile(layout1,2,[1 1]);
plotTitle = 'CSF_{vel}^{3}';
drawSyntheticSSF(syntheticDatasetVel_d,K,syntheticDataset,plotTitle);
view(axVel,-40,20)

axRval = nexttile(layout1,3,[1 1]);
plotTitle = 'CSF_{rnd}^{3}';
drawSyntheticSSF(syntheticDatasetRVal_d,K,syntheticDataset,plotTitle);
view(axRval,-40,20)

for i = 1:3
    ax = nexttile(layout2);
    imagesc(tmp_xy(:,:,i));
    ylim([1 160])
    xlim([1 24])
    yticks([30 80 130]);
    xticks([6 12 18]);
    clim([0,7]);

    if i == 1
        ax.YLabel.String = 'xy \Downarrow x';
    end

    if i == 2
        ax.XLabel.String = 'time';
    end

end

cm = colormap();    
cm(1,:) = [1,1,1];
colormap(cm);
cb = colorbar;
cb.Label.String = 'Velocity SSF';

% outfile = 'phantom.png';
% exportgraphics(gcf,outfile,'Resolution',300);