IH/iihmsp13/exp/plot_attaques.py

   1 import matplotlib.pyplot as mtpl
   2
   3
   4
   5
   6 nbexp=20
   7
   8 l_rifqi_spat=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.6138,0.8568,0.9132,0.9392,0.9444,0.9456,0.9512,0.9486,0.9526,0.9528,0.9536,0.9548,0.953,0.9536,0.9522,0.9528,0.9562,0.9556,0.9526,0.9532,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001,0.1906,0.7642,0.9164,0.9564,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0.1378,0.9564,0.1076,0.2862,0.2826,0.4354,0.7418,0.5414,0.602,0.6134,0.6736,0,0.0008,0.0002,0.0004,0.0002,0.0062,0.0004,0.0006,0.0004,0.0012,0,0.9288,0.6318,0.901,0.6288,0.893,0.6852,0.8784,0.6322,0.6448,0,0.77,0.637,0.8036,0.5488,0.8926,0.6954,0.9082,0.7478,0.9278,0.0072,0.9538,0.9544,0.9542,0.9522,0.9234,0.954,0.9524,0.9536,0.9544,0,0.9552,0.9552,0.9546,0.9552,0.7168,0.9552,0.9552,0.9546,0.9552,0.0952,0.9564,0.6782,0.1348,0.0122,0.0002,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
   9
  10
  11 l = l_rifqi_spat
  12
  13
  14 atq_rot =l[0:20]
  15 atq_rot_x=[3*(k+1) for k in range(nbexp)]
  16
  17 atq_comp =l[20:40]
  18 atq_comp_x=[100-3*(k+1) for k in range(nbexp)]
  19 #print atq_comp,atq_comp_x
  20
  21 atq_jp2 =l[40:60]
  22 atq_jp2_x=[100-5*(k+1) for k in range(nbexp)]
  23
  24 atq_dec=l[60:80]
  25 atq_dec_x=[(0.05*(k+1))**2 for k in range(nbexp)]
  26
  27 atq_flou=l[80:100]
  28 atq_flou=[0.1414,0.156,0.1396,0.15,0.0218,0.149,0.1674,0.2098,0.1034,0.2412,0.12,0.133,0.5606,0.5694,0.2058,0.5356,0.8128,0.8092,0.863,0.9442,0.8892,0.902,0.8986,0.9152,0,0.018,0.0216,0.0242,0.0252,0.0584,0.0338,0.0432,0.037,0.0434,0,0.0416,0.0372,0.0446,0.0436,0.0668,0.0472,0.0306,0.0374,0.0454,0,0.0462,0.04,0.0448,0.0396,0.919]
  29
  30 atq_flou_x=[0.04*(k+1) for k in range(50)]
  31
  32 atq_contrast=l[100:120]
  33
  34 atq_contrast=[0.949,0.9836,0.9414,0.9748,0.0048,0.9848,0.92,0.9208,0.9092,0.964,0.9294,0.9464,0.949,0.984,0.0426,0.9664,0.8988,0.9616,0.9108,0.9336,0.8754,0.9176,0.7738,0.7584,0,0.7854,0.906,0.9402,0.9458,0.9546,0.9628,0.968,0.942,0.953,0.057,0.9916,0.9622,0.9732,0.967,0.9888,0.9586,0.9824,0.968,0.9872,0.1848,0.9926,0.9766,0.993,0.974,0.9892]
  35 atq_contrast_x=[+0.04*(k+1) for k in range(50)]
  36 lm=[]
  37 for j in range(len(atq_contrast)):
  38     if atq_contrast[j] < 0.4:
  39         lm+=[atq_contrast_x[j]]
  40 print lm
  41
  42
  43
  44 atq_redim=l[120:140]
  45 atq_redim_x=[0.5+0.05*(k+1) for k in range(nbexp)]
  46
  47 atq_lsb=l[140:160]
  48 atq_lsb_x=[1-float(k)/(nbexp) for k in range(nbexp)]
  49
  50 atq_ci2_again=[0,0,0,0,0,0,0.01,0.01,0.022,0.048,0.083,0.167,0.26,0.362,0.491,0.618,0.701,0.773,0.841,0.871,0.902,0.926,0.938,0.948,0.955,0.959,0.969,0.97,0.972,0.978,0.98,0.979,0.982,0.981,0.983,0.987,0.988,0.989,0.99,0.99]
  51 atq_ci2_again_x=[0.1*30*(k+1)/40 for k in range(40)]
  52
  53
  54
  55
  56 fig = mtpl.figure()
  57
  58 ax_rot = fig.add_subplot(111)
  59 ax_rot.plot(atq_rot_x,atq_rot, '-')
  60 #mtpl.legend(('dwt(neg)','dwt(fl)', 'dct(neg)', 'dct(fl)'),'lower right', shadow=True)
  61 ax_rot.set_xlabel('Rotation Angle (degree)')
  62 ax_rot.set_ylabel('Fermi-Dirac measure betw. Marks')
  63 mtpl.savefig("atq-rot.pdf")
  64
  65
  66 fig = mtpl.figure()
  67
  68 ax_comp = fig.add_subplot(111)
  69 ax_comp.plot(atq_comp_x,atq_comp, '-')
  70 #mtpl.legend(('dwt(neg)','dwt(fl)', 'dct(neg)', 'dct(fl)'),'center left', shadow=True)
  71 ax_comp.set_xlabel('JPEG Compression Ratio')
  72 ax_comp.set_ylabel('Fermi-Dirac measure betw. Marks')
  73 mtpl.savefig("atq-jpg.pdf")
  74
  75
  76 fig = mtpl.figure()
  77
  78
  79 ax_jp2 = fig.add_subplot(111)
  80 ax_jp2.plot(atq_jp2_x,atq_jp2, '-')
  81 ax_jp2.set_xlabel('JPEG 2000 Compression Ratio')
  82 ax_jp2.set_ylabel('Fermi-Dirac measure betw. Marks')
  83 mtpl.savefig("atq-jp2.pdf")
  84
  85
  86
  87 fig = mtpl.figure()
  88 ax_dec = fig.add_subplot(111)
  89 ax_dec.plot(atq_dec_x,atq_dec, '-')
  90 ax_dec.set_xlabel('Percentage of Cropping')
  91 ax_dec.set_ylabel('Fermi-Dirac measure betw. Marks')
  92 mtpl.savefig("atq-dec.pdf")
  93
  94
  95
  96 fig = mtpl.figure()
  97 ax_flou = fig.add_subplot(111)
  98 ax_flou.plot(atq_flou_x,atq_flou, '-')
  99 ax_flou.set_xlabel('Noise Ratio')
 100 ax_flou.set_ylabel('Fermi-Dirac measure betw. Marks')
 101 mtpl.savefig("atq-flou.pdf")
 102
 103
 104 fig = mtpl.figure()
 105 ax_contrast = fig.add_subplot(111)
 106 ax_contrast.plot(atq_contrast_x,atq_contrast, '-')
 107 ax_contrast.set_xlabel('Contrast Ratio')
 108 ax_contrast.set_ylabel('Fermi-Dirac measure betw. Marks')
 109 mtpl.savefig("atq-contrast.pdf")
 110
 111
 112 fig = mtpl.figure()
 113 ax_redim = fig.add_subplot(111)
 114 ax_redim.plot(atq_redim_x,atq_redim, '-')
 115 ax_redim.set_xlabel('Redim Ratio')
 116 ax_redim.set_ylabel('Fermi-Dirac measure betw. Marks')
 117 mtpl.savefig("atq-redim.pdf")
 118
 119
 120 fig = mtpl.figure()
 121 ax_lsb = fig.add_subplot(111)
 122 ax_lsb.plot(atq_lsb_x,atq_lsb, '-')
 123 ax_lsb.set_xlabel('LSB Ratio')
 124 ax_lsb.set_ylabel('Fermi-Dirac measure betw. Marks')
 125 mtpl.savefig("atq-lsb.pdf")
 126
 127
 128 fig = mtpl.figure()
 129 ax_ci2 = fig.add_subplot(111)
 130 ax_ci2.plot(atq_ci2_again_x,atq_ci2_again, '-')
 131 ax_ci2.set_xlabel('number of iteration per pixel')
 132 ax_ci2.set_ylabel('Fermi-Dirac measure betw. Marks')
 133 mtpl.savefig("atq-ci2.pdf")
 134
 135
 136
 137
 138
 139
 140
 141 #mtpl.show()
 142
 143
 144 l_leven=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0026,0.01,0.0174,0.2658,0.2868,0.2906,0.2932,0.2938,0.2944,0.2946,0.293,0.2936,0.2948,0.294,0.294,0.295,0.294,0.2946,0.2932,0.294,0.294,0.2944,0.2946,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.009,0.1032,0.2562,0.2878,0.294,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0012,0.0234,0.091,0.2268,0.294,0.1994,0.2452,0.2412,0.2572,0.2834,0.2658,0.269,0.2702,0.2756,0,0.012,0.0144,0.0158,0.0162,0.0324,0.0182,0.023,0.02,0.0282,0.0052,0.2846,0.2736,0.2824,0.2612,0.2792,0.259,0.272,0.2428,0.2184,0,0.2274,0.2544,0.2814,0.279,0.2778,0.281,0.2906,0.285,0.2912,0.1012,0.2942,0.2944,0.2948,0.2946,0.2934,0.2942,0.2944,0.2946,0.2944,0,0.2942,0.2942,0.2942,0.2938,0.2836,0.2942,0.2942,0.2936,0.2938,0.2148,0.294,0.2824,0.2268,0.1478,0.0814,0.04,0.0166,0.0064,0,0,0,0,0,0,0,0,0,0,0,0]
 145
 146
 147 l_rifqi=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.6138,0.8568,0.9132,0.9392,0.9444,0.9456,0.9512,0.9486,0.9526,0.9528,0.9536,0.9548,0.953,0.9536,0.9522,0.9528,0.9562,0.9556,0.9526,0.9532,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001,0.1906,0.7642,0.9164,0.9564,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0.1378,0.9564,0.1076,0.2862,0.2826,0.4354,0.7418,0.5414,0.602,0.6134,0.6736,0,0.0008,0.0002,0.0004,0.0002,0.0062,0.0004,0.0006,0.0004,0.0012,0,0.9288,0.6318,0.901,0.6288,0.893,0.6852,0.8784,0.6322,0.6448,0,0.77,0.637,0.8036,0.5488,0.8926,0.6954,0.9082,0.7478,0.9278,0.0072,0.9538,0.9544,0.9542,0.9522,0.9234,0.954,0.9524,0.9536,0.9544,0,0.9552,0.9552,0.9546,0.9552,0.7168,0.9552,0.9552,0.9546,0.9552,0.0952,0.9564,0.6782,0.1348,0.0122,0.0002,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]