Contents
init; seti.contrast = 'corner2D'; seti.rotation = 20; seti = setGeomSim(seti); figure(1); imagesc(real(seti.G(seti.qROIexact))); axis xy; colorbar; % set(gca,'FontSize',20); axis square; print(1,'-depsc','guideForward1.eps'); [FFqMeas,FFqROI,seti] = forward(seti,seti.qROIexact); figure(2); imagesc(real(seti.G(FFqROI(:,7)))); axis xy; colorbar; % 7-th transmitter % set(gca,'FontSize',20); axis square; print(2,'-depsc','guideForward2.eps'); figure(3); imagesc(real(FFqMeas)); axis xy; colorbar; % set(gca,'FontSize',20); axis square; print(3,'-depsc','guideForward3.eps');
add noise...
seti.delta = 0.05; seti.whichNoise = 'normal'; seti.seed = 10; [seti, FmeasDelta] = addNoise(seti, FFqMeas); figure(4); imagesc(real(FFqMeas)); axis xy; colorbar; set(gca,'FontSize',20); axis square; print(4,'-depsc','guideForward4.eps');
Alternative
init; % Set a grid in 2D seti.qBall = 0.8; seti.rBall = 0.05; guideSetGrid; % Contrast qROI = referenceBall2D(seti.gridROI(1,:), seti.gridROI(2,:), seti); % Forward [FFqMeas,FFqROI,seti] = forward(seti,reshape(qROI,[seti.nROI^seti.dim 1]));