annotate scripts/timeline_view_movement.py @ 1814:a6d3b7b64943

Bug fixes in scripts: display plays in last timeline segment of each fragment (post last move); call sub-processes of generate_report.py in same folder (pass on argument); indentation in LaTeX file.
author Brecht De Man <b.deman@qmul.ac.uk>
date Wed, 16 Sep 2015 13:31:40 +0100
parents e4d060951e47
children adbaff4bf3f0 235594325b84
rev   line source
b@1807 1 #!/usr/bin/python
b@1807 2
b@1807 3 import xml.etree.ElementTree as ET
b@1807 4 import os # list files in directory
b@1807 5 import sys # command line arguments
b@1807 6 import matplotlib.pyplot as plt # plots
b@1807 7 import matplotlib.patches as patches # rectangles
b@1807 8
b@1807 9
b@1807 10 # COMMAND LINE ARGUMENTS
b@1807 11
b@1807 12 assert len(sys.argv)<3, "timeline_view_movement takes at most 1 command line argument\n"+\
b@1807 13 "Use: python timeline_view_movement.py [XML_files_location]"
b@1807 14
b@1807 15 # XML results files location
b@1807 16 if len(sys.argv) == 1:
b@1807 17 folder_name = "../saves" # Looks in 'saves/' folder from 'scripts/' folder
b@1807 18 print "Use: python timeline_view_movement.py [XML_files_location]"
b@1807 19 print "Using default path: " + folder_name
b@1807 20 elif len(sys.argv) == 2:
b@1807 21 folder_name = sys.argv[1] # First command line argument is folder
b@1807 22
b@1807 23 # check if folder_name exists
b@1807 24 if not os.path.exists(folder_name):
b@1807 25 #the file is not there
b@1807 26 print "Folder '"+folder_name+"' does not exist."
b@1807 27 sys.exit() # terminate script execution
b@1807 28 elif not os.access(os.path.dirname(folder_name), os.W_OK):
b@1807 29 #the file does exist but write privileges are not given
b@1807 30 print "No write privileges in folder '"+folder_name+"'."
b@1807 31
b@1807 32
b@1807 33 # CONFIGURATION
b@1807 34
b@1807 35 # Folder where to store timelines
b@1807 36 timeline_folder = folder_name + '/timelines_movement/' # Stores in 'saves/timelines_movement/' by default
b@1807 37
b@1807 38 # Font settings
b@1807 39 font = {'weight' : 'bold',
b@1807 40 'size' : 16}
b@1807 41 plt.rc('font', **font)
b@1807 42
b@1807 43 # Colormap for to cycle through
b@1807 44 colormap = ['b', 'g', 'c', 'm', 'y', 'k']
b@1807 45
b@1807 46 # figure size
b@1807 47 fig_width = 25
b@1807 48 fig_height = 10
b@1807 49
b@1807 50
b@1807 51 # CODE
b@1807 52
b@1807 53 # create timeline_folder if not yet created
b@1807 54 if not os.path.exists(timeline_folder):
b@1807 55 os.makedirs(timeline_folder)
b@1807 56
b@1807 57 # get every XML file in folder
b@1807 58 for file in os.listdir(folder_name):
b@1807 59 if file.endswith(".xml"):
b@1807 60 tree = ET.parse(folder_name + '/' + file)
b@1807 61 root = tree.getroot()
b@1807 62 subject_id = file[:-4] # drop '.xml'
b@1807 63
b@1807 64 previous_audioholder_time = 0 # time spent before current audioholder
b@1807 65 time_offset = 0 # test starts at zero
b@1807 66
b@1807 67 # ONE TIMELINE PER PAGE - make new plot per page
b@1807 68
b@1807 69 # get list of all page names
b@1807 70 for audioholder in root.findall("./audioholder"): # iterate over pages
b@1807 71 page_name = audioholder.get('id') # get page name
b@1809 72 plot_empty = True # check if any data is plotted
b@1807 73
b@1807 74 if page_name is None: # ignore 'empty' audio_holders
b@1807 75 print "Skipping empty audioholder name from "+subject_id+"."
b@1807 76 break
b@1807 77
b@1807 78 # subtract total audioholder length from subsequent audioholder event times
b@1807 79 audioholder_time_temp = audioholder.find("./metric/metricresult/[@id='testTime']")
b@1807 80 if audioholder_time_temp is not None:
b@1807 81 audioholder_time = float(audioholder_time_temp.text)
b@1807 82 else:
b@1807 83 print "Skipping audioholder without total time specified from "+subject_id+"."
b@1807 84 break
b@1807 85
b@1807 86 # get audioelements
b@1807 87 audioelements = audioholder.findall("./audioelement")
b@1807 88
b@1807 89 # sort alphabetically
b@1807 90 data = []
b@1807 91 for elem in audioelements: # from http://effbot.org/zone/element-sort.htm
b@1807 92 key = elem.get("id")
b@1807 93 data.append((key, elem))
b@1807 94 data.sort()
b@1807 95
b@1807 96 N_audioelements = len(audioelements) # number of audio elements for this page
b@1807 97 increment = 0 # increased for every new audioelement
b@1807 98
b@1807 99 # get axes handle
b@1807 100 fig = plt.figure(figsize=(fig_width, fig_height))
b@1807 101 ax = fig.add_subplot(111)
b@1807 102
b@1807 103 # for page [page_name], print comments related to fragment [id]
b@1807 104 #for tuple in data:
b@1807 105 # audioelement = tuple[1]
b@1807 106 for tuple in data:
b@1807 107 audioelement = tuple[1]
b@1807 108 if audioelement is not None: # Check it exists
b@1807 109 audio_id = str(audioelement.get('id'))
b@1807 110
b@1807 111 # break if no initial position or move events registered
b@1807 112 initial_position_temp = audioelement.find("./metric/metricresult/[@name='elementInitialPosition']")
b@1807 113 if initial_position_temp is None:
b@1807 114 print "Skipping "+page_name+" from "+subject_id+": does not have initial positions specified."
b@1807 115 break
b@1807 116
b@1808 117 # get move events, initial and eventual position
b@1807 118 initial_position = float(initial_position_temp.text)
b@1807 119 move_events = audioelement.findall("./metric/metricresult/[@name='elementTrackerFull']/timepos")
b@1807 120 final_position = float(audioelement.find("./value").text)
b@1807 121
b@1808 122 # get listen events
b@1808 123 start_times_global = []
b@1808 124 stop_times_global = []
b@1808 125 listen_events = audioelement.findall("./metric/metricresult/[@name='elementListenTracker']/event")
b@1808 126 for event in listen_events:
b@1808 127 # get testtime: start and stop
b@1808 128 start_times_global.append(float(event.find('testtime').get('start'))-time_offset)
b@1808 129 stop_times_global.append(float(event.find('testtime').get('stop'))-time_offset)
b@1808 130
b@1807 131 # display fragment name at start
b@1807 132 plt.text(0,initial_position+0.02,audio_id,color=colormap[increment%len(colormap)]) #,rotation=45
b@1807 133
b@1807 134 # previous position and time
b@1807 135 previous_position = initial_position
b@1807 136 previous_time = 0
b@1807 137
b@1808 138 # assume not playing at start
b@1808 139 currently_playing = False # keep track of whether fragment is playing during move event
b@1814 140
b@1807 141 # draw all segments except final one
b@1807 142 for event in move_events:
b@1809 143 # mark this plot as not empty
b@1809 144 plot_empty = False
b@1809 145
b@1808 146 # get time and final position of move event
b@1807 147 new_time = float(event.find("./time").text)-time_offset
b@1807 148 new_position = float(event.find("./position").text)
b@1808 149
b@1808 150 # get play/stop events since last move until current move event
b@1808 151 stop_times = []
b@1808 152 start_times = []
b@1808 153 # is there a play and/or stop event between previous_time and new_time?
b@1808 154 for time in start_times_global:
b@1808 155 if time>previous_time and time<new_time:
b@1808 156 start_times.append(time)
b@1808 157 for time in stop_times_global:
b@1808 158 if time>previous_time and time<new_time:
b@1808 159 stop_times.append(time)
b@1808 160 # if no play/stop events between move events, find out whether playing
b@1808 161
b@1808 162 segment_start = previous_time # first segment starts at previous move event
b@1808 163
b@1808 164 # draw segments (horizontal line)
b@1808 165 while len(start_times)+len(stop_times)>0: # while still play/stop events left
b@1808 166 if len(stop_times)<1: # upcoming event is 'play'
b@1808 167 # draw non-playing segment from segment_start to 'play'
b@1808 168 currently_playing = False
b@1808 169 segment_stop = start_times.pop(0) # remove and return first item
b@1808 170 elif len(start_times)<1: # upcoming event is 'stop'
b@1808 171 # draw playing segment (red) from segment_start to 'stop'
b@1808 172 currently_playing = True
b@1808 173 segment_stop = stop_times.pop(0) # remove and return first item
b@1808 174 elif start_times[0]<stop_times[0]: # upcoming event is 'play'
b@1808 175 # draw non-playing segment from segment_start to 'play'
b@1808 176 currently_playing = False
b@1808 177 segment_stop = start_times.pop(0) # remove and return first item
b@1808 178 else: # stop_times[0]<start_times[0]: upcoming event is 'stop'
b@1808 179 # draw playing segment (red) from segment_start to 'stop'
b@1808 180 currently_playing = True
b@1808 181 segment_stop = stop_times.pop(0) # remove and return first item
b@1808 182
b@1808 183 # draw segment
b@1808 184 plt.plot([segment_start, segment_stop], # x-values
b@1808 185 [previous_position, previous_position], # y-values
b@1808 186 color='r' if currently_playing else colormap[increment%len(colormap)],
b@1808 187 linewidth=3
b@1808 188 )
b@1808 189 segment_start = segment_stop # move on to next segment
b@1808 190 currently_playing = not currently_playing # toggle to draw final segment correctly
b@1808 191
b@1808 192 # draw final segment (horizontal line) from last 'segment_start' to current move event time
b@1808 193 plt.plot([segment_start, new_time], # x-values
b@1807 194 [previous_position, previous_position], # y-values
b@1808 195 # color depends on playing during move event or not:
b@1808 196 color='r' if currently_playing else colormap[increment%len(colormap)],
b@1807 197 linewidth=3
b@1807 198 )
b@1808 199
b@1807 200 # vertical line from previous to current position
b@1807 201 plt.plot([new_time, new_time], # x-values
b@1807 202 [previous_position, new_position], # y-values
b@1808 203 # color depends on playing during move event or not:
b@1808 204 color='r' if currently_playing else colormap[increment%len(colormap)],
b@1807 205 linewidth=3
b@1807 206 )
b@1807 207
b@1807 208 # update previous_position value
b@1807 209 previous_position = new_position
b@1807 210 previous_time = new_time
b@1807 211
b@1814 212
b@1814 213
b@1814 214 # draw final horizontal segment (or only segment if audioelement not moved)
b@1807 215 # horizontal line from previous time to end of audioholder
b@1814 216
b@1814 217 # get play/stop events since last move until current move event
b@1814 218 stop_times = []
b@1814 219 start_times = []
b@1814 220 # is there a play and/or stop event between previous_time and new_time?
b@1814 221 for time in start_times_global:
b@1814 222 if time>previous_time and time<audioholder_time-time_offset:
b@1814 223 start_times.append(time)
b@1814 224 for time in stop_times_global:
b@1814 225 if time>previous_time and time<audioholder_time-time_offset:
b@1814 226 stop_times.append(time)
b@1814 227 # if no play/stop events between move events, find out whether playing
b@1814 228
b@1814 229 segment_start = previous_time # first segment starts at previous move event
b@1814 230
b@1814 231 # draw segments (horizontal line)
b@1814 232 while len(start_times)+len(stop_times)>0: # while still play/stop events left
b@1814 233 # mark this plot as not empty
b@1814 234 plot_empty = False
b@1814 235 if len(stop_times)<1: # upcoming event is 'play'
b@1814 236 # draw non-playing segment from segment_start to 'play'
b@1814 237 currently_playing = False
b@1814 238 segment_stop = start_times.pop(0) # remove and return first item
b@1814 239 elif len(start_times)<1: # upcoming event is 'stop'
b@1814 240 # draw playing segment (red) from segment_start to 'stop'
b@1814 241 currently_playing = True
b@1814 242 segment_stop = stop_times.pop(0) # remove and return first item
b@1814 243 elif start_times[0]<stop_times[0]: # upcoming event is 'play'
b@1814 244 # draw non-playing segment from segment_start to 'play'
b@1814 245 currently_playing = False
b@1814 246 segment_stop = start_times.pop(0) # remove and return first item
b@1814 247 else: # stop_times[0]<start_times[0]: upcoming event is 'stop'
b@1814 248 # draw playing segment (red) from segment_start to 'stop'
b@1814 249 currently_playing = True
b@1814 250 segment_stop = stop_times.pop(0) # remove and return first item
b@1814 251
b@1814 252 # draw segment
b@1814 253 plt.plot([segment_start, segment_stop], # x-values
b@1814 254 [previous_position, previous_position], # y-values
b@1814 255 color='r' if currently_playing else colormap[increment%len(colormap)],
b@1814 256 linewidth=3
b@1814 257 )
b@1814 258 segment_start = segment_stop # move on to next segment
b@1814 259 currently_playing = not currently_playing # toggle to draw final segment correctly
b@1814 260
b@1814 261 # draw final segment (horizontal line) from last 'segment_start' to current move event time
b@1814 262 plt.plot([segment_start, audioholder_time-time_offset], # x-values
b@1807 263 [previous_position, previous_position], # y-values
b@1814 264 # color depends on playing during move event or not:
b@1814 265 color='r' if currently_playing else colormap[increment%len(colormap)],
b@1807 266 linewidth=3
b@1807 267 )
b@1807 268
b@1814 269 # plt.plot([previous_time, audioholder_time-time_offset], # x-values
b@1814 270 # [previous_position, previous_position], # y-values
b@1814 271 # color=colormap[increment%len(colormap)],
b@1814 272 # linewidth=3
b@1814 273 # )
b@1814 274
b@1807 275 # display fragment name at end
b@1807 276 plt.text(audioholder_time-time_offset,previous_position,\
b@1807 277 audio_id,color=colormap[increment%len(colormap)]) #,rotation=45
b@1807 278
b@1807 279 increment+=1 # to next audioelement
b@1807 280
b@1807 281 last_audioholder_duration = audioholder_time-time_offset
b@1807 282 time_offset = audioholder_time
b@1807 283
b@1809 284 if not plot_empty: # if plot is not empty, show or store
b@1809 285 # set plot parameters
b@1809 286 plt.title('Timeline ' + file + ": "+page_name)
b@1809 287 plt.xlabel('Time [seconds]')
b@1809 288 plt.xlim(0, last_audioholder_duration)
b@1809 289 plt.ylabel('Rating') # default
b@1809 290 plt.ylim(0, 1) # rating between 0 and 1
b@1807 291
b@1809 292 #y-ticks: labels on rating axis
b@1809 293 label_positions = []
b@1809 294 label_text = []
b@1809 295 scale_tags = root.findall("./BrowserEvalProjectDocument/audioHolder/interface/scale")
b@1809 296 scale_title = root.find("./BrowserEvalProjectDocument/audioHolder/interface/title")
b@1809 297 for tag in scale_tags:
b@1809 298 label_positions.append(float(tag.get('position'))/100) # on a scale from 0 to 100
b@1809 299 label_text.append(tag.text)
b@1809 300 if len(label_positions) > 0: # if any labels available
b@1809 301 plt.yticks(label_positions, label_text) # show rating axis labels
b@1809 302 # set label Y-axis
b@1809 303 if scale_title is not None:
b@1809 304 plt.ylabel(scale_title.text)
b@1807 305
b@1809 306 #plt.show() # uncomment to show plot; comment when just saving
b@1809 307 #exit()
b@1807 308
b@1809 309 plt.savefig(timeline_folder+subject_id+"-"+page_name+".pdf", bbox_inches='tight')
b@1809 310 plt.close()
b@1807 311