annotate scripts/timeline_view.py @ 1491:5fd7e8f4313f

Add interfaces section to WAC paper
author Dave <djmoffat@users.noreply.github.com>
date Mon, 12 Oct 2015 12:22:42 +0100
parents 0b095f66de65
children 1b6fa37d46a4 235594325b84
rev   line source
djmoffat@718 1 #!/usr/bin/python
djmoffat@718 2
djmoffat@718 3 import xml.etree.ElementTree as ET
djmoffat@718 4 import os # list files in directory
djmoffat@718 5 import sys # command line arguments
djmoffat@718 6 import matplotlib.pyplot as plt # plots
djmoffat@718 7 import matplotlib.patches as patches # rectangles
djmoffat@718 8
djmoffat@718 9 # COMMAND LINE ARGUMENTS
djmoffat@718 10
djmoffat@718 11 assert len(sys.argv)<3, "timeline_view takes at most 1 command line argument\n"+\
djmoffat@718 12 "Use: python timeline_view.py [XML_files_location]"
djmoffat@718 13
djmoffat@718 14 # XML results files location
djmoffat@718 15 if len(sys.argv) == 1:
djmoffat@718 16 folder_name = "../saves" # Looks in 'saves/' folder from 'scripts/' folder
djmoffat@718 17 print "Use: python timeline_view.py [XML_files_location]"
djmoffat@718 18 print "Using default path: " + folder_name
djmoffat@718 19 elif len(sys.argv) == 2:
djmoffat@718 20 folder_name = sys.argv[1] # First command line argument is folder
djmoffat@718 21
djmoffat@718 22 # check if folder_name exists
djmoffat@718 23 if not os.path.exists(folder_name):
djmoffat@718 24 #the file is not there
djmoffat@718 25 print "Folder '"+folder_name+"' does not exist."
djmoffat@718 26 sys.exit() # terminate script execution
djmoffat@718 27 elif not os.access(os.path.dirname(folder_name), os.W_OK):
djmoffat@718 28 #the file does exist but write privileges are not given
djmoffat@718 29 print "No write privileges in folder '"+folder_name+"'."
djmoffat@718 30
djmoffat@718 31
djmoffat@718 32 # CONFIGURATION
djmoffat@718 33
djmoffat@718 34 # Folder where to store timelines
djmoffat@718 35 timeline_folder = folder_name + '/timelines/' # Stores in 'saves/timelines/'
djmoffat@718 36
djmoffat@718 37 # Font settings
djmoffat@718 38 font = {'weight' : 'bold',
djmoffat@718 39 'size' : 16}
djmoffat@718 40 plt.rc('font', **font)
djmoffat@718 41
djmoffat@718 42 # Colormap for to cycle through
djmoffat@718 43 colormap = ['b', 'r', 'g', 'c', 'm', 'y', 'k']
djmoffat@718 44
djmoffat@718 45 # if enabled, x-axis shows time per audioholder, not total test time
djmoffat@718 46 show_audioholder_time = True
djmoffat@718 47
djmoffat@718 48 # bar height (<1 to avoid overlapping)
djmoffat@718 49 bar_height = 0.6
djmoffat@718 50
djmoffat@718 51 # figure size
djmoffat@718 52 fig_width = 25
djmoffat@718 53 fig_height = 5
djmoffat@718 54
djmoffat@718 55
djmoffat@718 56 # CODE
djmoffat@718 57
djmoffat@718 58 # create timeline_folder if not yet created
djmoffat@718 59 if not os.path.exists(timeline_folder):
djmoffat@718 60 os.makedirs(timeline_folder)
djmoffat@718 61
djmoffat@718 62 # get every XML file in folder
djmoffat@718 63 for file in os.listdir(folder_name):
djmoffat@718 64 if file.endswith(".xml"):
djmoffat@718 65 tree = ET.parse(folder_name + '/' + file)
djmoffat@718 66 root = tree.getroot()
djmoffat@718 67 subject_id = file[:-4] # drop '.xml'
djmoffat@718 68
djmoffat@718 69 time_offset = 0 # test starts at zero
djmoffat@718 70
djmoffat@718 71 # ONE TIMELINE PER PAGE - make new plot per page
djmoffat@718 72
djmoffat@718 73 # get list of all page names
djmoffat@718 74 for audioholder in root.findall("./audioholder"): # iterate over pages
djmoffat@718 75 page_name = audioholder.get('id') # get page name
djmoffat@718 76 plot_empty = True # check if any data is plotted
djmoffat@718 77
djmoffat@718 78 if page_name is None: # ignore 'empty' audio_holders
djmoffat@718 79 break
djmoffat@718 80
djmoffat@718 81 # SORT AUDIO ELEMENTS ALPHABETICALLY
djmoffat@718 82 audioelements = audioholder.findall("./audioelement")
djmoffat@718 83
djmoffat@718 84 data = []
djmoffat@718 85 for elem in audioelements: # from http://effbot.org/zone/element-sort.htm
djmoffat@718 86 key = elem.get("id")
djmoffat@718 87 data.append((key, elem))
djmoffat@718 88 data.sort()
djmoffat@718 89
djmoffat@718 90 N_audioelements = len(audioelements) # number of audio elements for this page
djmoffat@718 91 increment = 0 # increased for every new audioelement
djmoffat@718 92 audioelements_names = [] # store names of audioelements
djmoffat@718 93
djmoffat@718 94 # get axes handle
djmoffat@718 95 fig = plt.figure(figsize=(fig_width, fig_height))
djmoffat@718 96 ax = fig.add_subplot(111) #, aspect='equal'
djmoffat@718 97
djmoffat@718 98 # for page [page_name], print comments related to fragment [id]
djmoffat@718 99 for tuple in data:
djmoffat@718 100 audioelement = tuple[1]
djmoffat@718 101 if audioelement is not None: # Check it exists
djmoffat@718 102 audio_id = str(audioelement.get('id'))
djmoffat@718 103 audioelements_names.append(audio_id)
djmoffat@718 104
djmoffat@718 105 # for this audioelement, loop over all listen events
djmoffat@718 106 listen_events = audioelement.findall("./metric/metricresult/[@name='elementListenTracker']/event")
djmoffat@718 107 for event in listen_events:
djmoffat@718 108 # mark this plot as not empty
djmoffat@718 109 plot_empty = False
djmoffat@718 110
djmoffat@718 111 # get testtime: start and stop
djmoffat@718 112 start_time = float(event.find('testtime').get('start'))-time_offset
djmoffat@718 113 stop_time = float(event.find('testtime').get('stop'))-time_offset
djmoffat@718 114 # event lines:
djmoffat@718 115 ax.plot([start_time, start_time], # x-values
djmoffat@718 116 [0, N_audioelements+1], # y-values
djmoffat@718 117 color='k'
djmoffat@718 118 )
djmoffat@718 119 ax.plot([stop_time, stop_time], # x-values
djmoffat@718 120 [0, N_audioelements+1], # y-values
djmoffat@718 121 color='k'
djmoffat@718 122 )
djmoffat@718 123 # plot time:
djmoffat@718 124 ax.add_patch(
djmoffat@718 125 patches.Rectangle(
djmoffat@718 126 (start_time, N_audioelements-increment-bar_height/2), # (x, y)
djmoffat@718 127 stop_time - start_time, # width
djmoffat@718 128 bar_height, # height
djmoffat@718 129 color=colormap[increment%len(colormap)] # colour
djmoffat@718 130 )
djmoffat@718 131 )
djmoffat@718 132
djmoffat@718 133 increment+=1 # to next audioelement
djmoffat@718 134
djmoffat@718 135 # subtract total audioholder length from subsequent audioholder event times
djmoffat@718 136 audioholder_time = audioholder.find("./metric/metricresult/[@id='testTime']")
djmoffat@718 137 if audioholder_time is not None and show_audioholder_time:
djmoffat@718 138 time_offset = float(audioholder_time.text)
djmoffat@718 139
djmoffat@718 140 if not plot_empty:
djmoffat@718 141 # set plot parameters
djmoffat@718 142 plt.title('Timeline ' + file + ": "+page_name)
djmoffat@718 143 plt.xlabel('Time [seconds]')
djmoffat@718 144 plt.ylabel('Fragment')
djmoffat@718 145 plt.ylim(0, N_audioelements+1)
djmoffat@718 146
djmoffat@718 147 #y-ticks: fragment IDs, top to bottom
djmoffat@718 148 plt.yticks(range(N_audioelements, 0, -1), audioelements_names) # show fragment names
djmoffat@718 149
djmoffat@718 150
djmoffat@718 151 #plt.show() # uncomment to show plot; comment when just saving
djmoffat@718 152 #exit()
djmoffat@718 153
djmoffat@718 154 plt.savefig(timeline_folder+subject_id+"-"+page_name+".pdf", bbox_inches='tight')
djmoffat@718 155 plt.close()
djmoffat@718 156
djmoffat@718 157 #TODO: if 'nonsensical' or unknown: dashed line until next event
djmoffat@718 158 #TODO: Vertical lines for fragment looping point
djmoffat@718 159