Modified histogram shell plot script, added python dynamics plot script

webrtc/modules/remote_bitrate_estimator/test/plot_dynamics.py

+#!/usr/bin/env python
+# Copyright (c) 2015 The WebRTC project authors. All Rights Reserved.
+#
+# Use of this source code is governed by a BSD-style license
+# that can be found in the LICENSE file in the root of the source
+# tree. An additional intellectual property rights grant can be found
+# in the file PATENTS.  All contributing project authors may
+# be found in the AUTHORS file in the root of the source tree.
+
+# This script is used to plot simulation dynamics.
+# Able to plot each flow separately. Other plot boxes can be added,
+# currently one for Throughput, one for Latency and one for Packet Loss.
+
+import matplotlib
+import matplotlib.pyplot as plt
+import math
+import numpy
+import re
+import sys
+import pprint
+
+
+class Variable:
+  def __init__(self, variable):
+    self._ID = variable[0]
+    self._xlabel = variable[1]
+    self._ylabel = variable[2]
+    self._subplot = variable[3]
+    self._y_max = variable[4]
+    self._samples = dict()
+
+  def getID(self):
+    return self._ID
+
+  def getXLabel(self):
+    return self._xlabel
+
+  def getYLabel(self):
+    return self._ylabel
+
+  def getSubplot(self):
+    return self._subplot
+
+  def getYMax(self):
+    return self._y_max
+
+  def getNumberOfFlows(self):
+    return len(self._samples)
+
+
+  def addSample(self, line):
+    groups = re.search(r'/\d_(\w+)#\d@(\w*)', line)
+    var_name = groups.group(1)
+    alg_name = groups.group(2)
+
+    if alg_name not in self._samples.keys():
+      self._samples[alg_name] = {}
+
+    if var_name not in self._samples[alg_name].keys():
+      self._samples[alg_name][var_name] = []
+
+    sample = re.search(r'(\d+\.\d+)\t([-]?\d+\.\d+)', line)
+    s = (sample.group(1),sample.group(2))
+    self._samples[alg_name][var_name].append(s)
+
+def plotVar(v, ax, show_legend, show_x_label):
+  if show_x_label:
+    ax.set_xlabel(v.getXLabel(), fontsize='large')
+  ax.set_ylabel(v.getYLabel(), fontsize='large')
+
+  for alg in v._samples.keys():
+    print "alg %s: %d" % (alg, len(v._samples[alg]))

[stefan-webrtc, 2015/07/21 12:16:37]

Can probably be moved

[magalhaesc, 2015/07/21 17:05:03]

Done.

+    i = 1
+    for series in v._samples[alg].keys():
+      x = [sample[0] for sample in v._samples[alg][series]]
+      y = [sample[1] for sample in v._samples[alg][series]]
+      x = numpy.array(x)
+      y = numpy.array(y)
+      line = plt.plot(x, y, label=alg, linewidth=4.0)
+      colormap = {'Available1':'#AAAAAA',
+                  'Available2':'#AAAAAA',
+                  'GCC1':'#80D000',
+                  'GCC2':'#008000',
+                  'GCC3':'#00F000',
+                  'GCC4':'#00B000',
+                  'GCC5':'#70B020',
+                  'NADA1':'#0000AA',
+                  'NADA2':'#A0A0FF',
+                  'NADA3':'#0000FF',
+                  'NADA4':'#C0A0FF',
+                  'NADA5':'#9060B0',
+                  'TCP1':'#AAAAAA',
+                  'TCP2':'#AAAAAA',
+                  'TCP3':'#AAAAAA',
+                  'TCP4':'#AAAAAA',
+                  'TCP5':'#AAAAAA',
+                  'TCP6':'#AAAAAA',
+                  'TCP7':'#AAAAAA',
+                  'TCP8':'#AAAAAA',
+                  'TCP9':'#AAAAAA',
+                  'TCP10':'#AAAAAA',}
+
+      plt.setp(line, color=colormap[alg + str(i)])
+      if alg.startswith('Available'):
+        plt.setp(line, linestyle='--')
+      plt.grid(True)
+
+      x1, x2, y1, y2 = plt.axis()
+      if v.getYMax() >= 0:
+        y2 = v.getYMax()
+      plt.axis((0, x2, 0, y2))
+      i += 1
+
+    if show_legend:
+      legend = plt.legend(loc='upper right', shadow=True, fontsize='large', ncol=len(v._samples))
+
+if __name__ == '__main__':
+
+  variables = [
+          #id, xlabel, ylabel,  subplot
+          #('Send_estimate', "Time (s)", "Throughput (kbps)", 1, 4000),
+          ('Throughput_kbps', "Time (s)", "Throughput (kbps)", 1, 4000),
+          ('Delay_ms', "Time (s)", "One-way Delay(ms)", 2, 500),
+          ('Packet_Loss', "Time (s)", "Packet Loss Ratio", 3, 1.0),
+          #('threshold', "Time (s)", "Threshold (ms)",  3),
+          #('offset', "Time (s)", "Offset (ms)", 3),

[stefan-webrtc, 2015/07/21 12:16:37]

This should be cleaned up

[magalhaesc, 2015/07/21 17:05:03]

Done.

+          ]
+
+  var = []
+
+  # Create objects.
+  for variable in variables:
+    var.append(Variable(variable))
+
+  # Add samples to the objects.
+  for line in sys.stdin:
+    if line.startswith("[ RUN      ]"):
+        test_name = re.search('\.(\w+)', line).group(1)
+    if line.startswith("PLOT"):
+      for v in var:
+        if v.getID() in line:
+          v.addSample(line)
+
+  matplotlib.rcParams.update({'font.size': 20})
+
+  # Plot variables.
+  fig = plt.figure()
+  #fig.suptitle(test_name)

[stefan-webrtc, 2015/07/21 12:16:37]

Remove?

[magalhaesc, 2015/07/21 17:05:02]

Done.

+
+  print "Plotting \"%s\"" % (test_name)
+
+  # Offest and threshold on the same plot.
+  n = var[-1].getSubplot()
+  i = 0
+  for v in var:
+    ax = fig.add_subplot(n, 1, v.getSubplot())
+    plotVar(v, ax, i == 0, i == n - 1)
+    i += 1
+
+  #fig.savefig(test_name+".jpg")
+  plt.show()
+