%% fig6.m: script to display fig6:

% variables in fig6.mat:
% imMap15 - FDG-guided MAPEM reconstruction with nMax = 15
% imNlm35 - FDG-NLM filtered MLEM reconstruction with nMax = 35
% imNlm125 - FDG-NLM filtered MLEM reconstruction with nMax = 125

load('fig6.mat')

sliceNum = 64;
sliceNumSag = 120;

figure
subplot(2,3,3),petdisplay(imMap15(:,:,sliceNum));
pettitle('FDG-guided MAPEM, nMax = 15')

subplot(2,3,6),petdisplay(squeeze(imMap15(:,sliceNumSag,:))');
pettitle('FDG-guided MAPEM, nMax = 15')

subplot(2,3,1),petdisplay(imNlm35(:,:,sliceNum));
pettitle('MLEM + FDG-NLM, nMax = 35 ')

subplot(2,3,4),petdisplay(squeeze(imNlm35(:,sliceNumSag,:))');
pettitle('MLEM + FDG-NLM, nMax = 35 ')

subplot(2,3,2),petdisplay(imNlm125(:,:,sliceNum));
pettitle('MLEM + FDG-NLM, nMax = 125 ')

subplot(2,3,5),petdisplay(squeeze(imNlm125(:,sliceNumSag,:))');
pettitle('MLEM + FDG-NLM, nMax = 125 ')

setColourAxes(gcf,[0,0.12])
