% plot fsc

t = importdata('fsc_postprocess.txt');
f = figure;
plot((t(:,2)),t(:,4),'linewidth',2)
xlabel(['Resolution [1/' char(197) ']'],'fontsize',12)
ylabel('Fourier shell correlation','fontsize',12)
yline(0.143)
text(0.23,0.18,['4.5' char(197)])
set(gca, 'YTick', unique([0.143, get(gca, 'YTick')]));
grid on
print(f,'fsc_plot.png','-dpng','-r600');