% Plot the results of XRMC_NoCrimeData_comp.m. 
%
% Arguments:
% N    resolution of the previous computation is N x N
%
% Samuli Siltanen March 2011

function XRMC_NoCrimeData_plot(N)

% Load precomputed results at resolution N from file
loadcommand = ['load data/XRMC_NoCrime', num2str(N), ' N m mnc mncn measang target err_sup err_squ'];
eval(loadcommand)


% Take a look at the measurement WITH inverse crime
figure(2)
clf
imagesc(m)
colormap gray
axis square
axis off

% Take a look at the measurement WITHOUT inverse crime
figure(3)
clf
imagesc(mnc)
colormap gray
axis square
axis off

% Take a look at the difference between the two measurements
figure(4)
clf
imagesc(abs(m-mnc))
colormap gray
axis square
axis off
title(['Sup norm error ', num2str(round(err_sup*100)), '%; square norm error ', ...
    num2str(round(err_squ*100)), '%'])

% Take a look at the noisy measurement WITHOUT inverse crime
figure(5)
clf
imagesc(mncn)
colormap gray
axis square
axis off
title('noisy measurement WITHOUT inverse crime')

% Show actual phantom (target)
figure(6)
clf
imagesc(target)
colormap gray
axis square
axis off
title(['Shepp-Logan phantom at size ', num2str(N), 'x', num2str(N)])