MakeExampleGridsPlot.m

clear all;

baseFolder = pwd;

doSave = false;
saveFolder = pwd;
if doSave
    dfws = get(0,'DefaultFigureWindowStyle');
    set(0,'DefaultFigureWindowStyle','normal');
end

Nsim = 5;
R = 4.4/Nsim;
FOVinplane = 0.22;

sliceSpacing = FOVinplane/2/Nsim;

[~, ~, ~, ~, ~, ~, ~, fh] = ...
    ComputeMostHexagonalSampling(R, FOVinplane, sliceSpacing, Nsim, [], [], [], true, 0, [], 0.001);

if doSave
    set(fh{1},'Position',[29 174 1479 588])
    export_fig(fullfile(saveFolder, 'ExampleGrids'),'-png','-transparent','-a1','-r300',fh{1});

    set(fh{2},'Position',[845.8  418.6  372.8  244.8])
    hLegend = findobj(gcf, 'Type', 'Legend');
    hLegend.Position = [0.2517    0.2370    0.5359    0.2708];
    export_fig(fullfile(saveFolder, 'ExampleGridsPackingDensity'),'-png','-transparent','-a1','-r300',fh{2});
    set(0,'DefaultFigureWindowStyle',dfws);
end