%% fig1.m: script to display fig1:

% variables in fig1.mat:
% gtFdg - high resolution ground truth FDG image volume
% gtMet - high resolution ground truth methionine image volume
% gtMu  - high resolution ground truth attenuation map

load('fig1.mat')
%%
figure
subplot(1,3,1),petdisplay(gtFdg(:,:,166));
subplot(1,3,2),petdisplay(gtMet(:,:,166));
subplot(1,3,3),petdisplay(gtMu(:,:,166));
colormap(gca,gray(256))