annotate scripts/timeline_view_movement.py @ 22:1f375b7d75fd tip

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