init;

seti.rBall = 0.05;
seti.qBall = 0.8;

guideSetGrid;
% Evaluate predefined contrast referenceBall2D in the grid
q2D = referenceBall2D(seti.gridROI(1,:), seti.gridROI(2,:), seti);
x = seti.gridROI(1,:); y = seti.gridROI(2,:);
figure(1); imagesc(x,y,real(reshape(q2D,[seti.nROI seti.nROI]))); axis xy; colorbar;
% set(gca,'FontSize',20); axis square; print(1,'-depsc','guideContrast2D.eps');

guideSetGrid3D;
q3D = referenceBall3D(seti.gridROI(1,:), seti.gridROI(2,:), seti.gridROI(3,:), seti);
seti.faceColor3D = 'blue';
figure(2); contourPlotROI(q3D, seti, 'real');
% set(gca,'FontSize',20); axis square; print(2,'-depsc','guideContrast3D.eps');