% Data from: http://dx.doi.org/10.1109/TIP.2012.2235847
% Just showing how to plot this...

l = {'FISTA [36]','l_0-SPAR [37]','IDD-BM3D [42]','ASDS-Reg [21]','NCSR'};
d1 = [28.37 27.10 29.21 28.70 29.68
      29.04 29.86 31.20 30.80 31.08
      26.82 26.97 28.56 28.08 28.62
      31.99 32.98 34.44 34.03 34.31
      29.11 29.34 31.06 31.22 31.95
      28.33 28.72 29.70 29.92 29.96
      25.75 26.42 27.98 27.86 28.10
      27.75 28.11 29.48 29.72 30.28
      28.43 28.66 29.62 29.48 29.66
      26.49 26.30 29.38 28.59 29.98];
d2 = [0.9119 0.8879 0.9287 0.9053 0.9271
      0.8858 0.9094 0.9304 0.9236 0.9294
      0.8627 0.8689 0.9007 0.8950 0.9026
      0.9017 0.9225 0.9369 0.9337 0.9415
      0.9002 0.9262 0.9364 0.9306 0.9411
      0.8798 0.9063 0.9197 0.9256 0.9254
      0.8375 0.8691 0.9014 0.9088 0.9117
      0.8775 0.8951 0.9167 0.9208 0.9293
      0.8813 0.9066 0.9200 0.9203 0.9220
      0.8958 0.8776 0.9295 0.9075 0.9341];
d3 = [30.36 30.73 29.83 30.84
      29.36 31.68 30.27 31.49
      26.81 28.17 27.29 28.34
      31.50 34.08 31.87 33.63
      31.23 32.89 32.93 33.39
      29.47 31.45 30.36 31.26
      25.03 27.19 27.05 27.91
      29.65 31.66 31.91 32.27
      29.42 29.99 28.95 30.16
      29.36 31.40 30.62 31.57];
d4 = [0.9452 0.9442 0.9126 0.9381
      0.9024 0.9426 0.9064 0.9371
      0.8845 0.9136 0.8637 0.9078
      0.8968 0.9359 0.8978 0.9333
      0.9290 0.9561 0.9576 0.9587
      0.9011 0.9430 0.9058 0.9389
      0.8415 0.8986 0.8881 0.9088
      0.9256 0.9496 0.9491 0.9551
      0.9057 0.9373 0.9039 0.9331
      0.9393 0.9512 0.9304 0.9508];

set(gcf,'position',[390,878,980,400])
set(gcf,'defaultaxesfontsize',14)

subplot(2,2,1)
boxplot(d1)
set(gca,'xticklabel',l)
view([90 90])
box on
title('9\times{}9 Uniform Blur, \sigma_n = \surd{}2, PSNR')

subplot(2,2,2)
boxplot(d2)
set(gca,'xticklabel',l)
view([90 90])
box on
title('9\times{}9 Uniform Blur, \sigma_n = \surd{}2, FSIM')
set(gca,'ylim',[0.82,0.96])

subplot(2,2,3)
boxplot(d3)
set(gca,'xticklabel',l([1,3:5]))
view([90 90])
box on
title('Gaussian Blur, \sigma_n = \surd{}2, PSNR')

subplot(2,2,4)
boxplot(d4)
set(gca,'xticklabel',l([1,3:5]))
view([90 90])
box on
title('Gaussian Blur, \sigma_n = \surd{}2, FSIM')
set(gca,'ylim',[0.82,0.96])

print -dpng -r0 PSNR_big_table_graph.png

%%% ALT 2:

set(gcf,'position',[390,878,980,400])
set(gcf,'defaultaxesfontsize',14)

subplot(2,2,1)
plot([0,5],[0,0],'k:')
hold on
boxplot(repmat(d1(:,5),1,4)-d1(:,1:4))
set(gca,'xticklabel',l(1:4))
view([90 90])
box on
title('9\times{}9 Uniform Blur, \sigma_n = \surd{}2, PSNR')
set(gca,'ylim',[-3,4])

subplot(2,2,2)
plot([0,5],[0,0],'k:')
hold on
boxplot(repmat(d2(:,5),1,4)-d2(:,1:4))
set(gca,'xticklabel',l(1:4))
view([90 90])
box on
title('9\times{}9 Uniform Blur, \sigma_n = \surd{}2, FSIM')

subplot(2,2,3)
plot([0,4],[0,0],'k:')
hold on
boxplot(repmat(d1(:,4),1,3)-d3(:,1:3))
set(gca,'xticklabel',l([1,3:4]))
view([90 90])
box on
title('Gaussian Blur, \sigma_n = \surd{}2, PSNR')
set(gca,'ylim',[-3,4])

subplot(2,2,4)
plot([0,4],[0,0],'k:')
hold on
boxplot(repmat(d4(:,4),1,3)-d4(:,1:3))
set(gca,'xticklabel',l([1,3:4]))
view([90 90])
box on
title('Gaussian Blur, \sigma_n = \surd{}2, FSIM')

print -dpng -r0 PSNR_big_table_graph2.png
