-
Notifications
You must be signed in to change notification settings - Fork 7
/
velocityTest.py
139 lines (109 loc) · 4.05 KB
/
velocityTest.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
###############################################################################
# Event Detection Algorithm Suite
# Copyright (C) 2012 Gian Perrone (http://github.com/gian)
#
# Permission to use, copy, modify, and distribute this software and its
# documentation for any purpose and without fee is hereby granted,
# provided that the above copyright notice appear in all copies and that
# both the copyright notice and this permission notice and warranty
# disclaimer appear in supporting documentation, and that the name of
# the above copyright holders, or their entities, not be used in
# advertising or publicity pertaining to distribution of the software
# without specific, written prior permission.
#
# The above copyright holders disclaim all warranties with regard to
# this software, including all implied warranties of merchantability and
# fitness. In no event shall the above copyright holders be liable for
# any special, indirect or consequential damages or any damages
# whatsoever resulting from loss of use, data or profits, whether in an
# action of contract, negligence or other tortious action, arising out
# of or in connection with the use or performance of this software.
###############################################################################
from detect.sample import Sample
from detect.sample import FileSampleStream
from detect.dispersion import *
from detect.velocity import *
from detect.hmm import *
from detect.prefix import *
from detect.aoi import *
from detect.movingaverage import *
from detect.srr import *
from detect.intersamplevelocity import *
from detect.sgfilter import *
from detect.blinkfilter import *
print "============= I-VT MovingAverage test ==============="
stream = FileSampleStream('testData/UH27_img_vy_labelled_MN.txt')
fstream = IntersampleVelocity(MovingAverageFilter(BlinkFilter(stream),9))
idt = Velocity(fstream, 0.004)
fixations = []
for i in idt:
print i
fixations.append(i)
verifStream = BlinkFilter(FileSampleStream('testData/UH27_img_vy_labelled_MN.txt'))
taggedEvents = []
cfix = []
eventType = []
for i in verifStream:
if i.eventType == 1:
cfix.append(i)
eventType.append(1)
else:
eventType.append(0)
if len(cfix) == 0:
continue
print ("Fixation of length: " + str(len(cfix)) + " starting at sample " + str(cfix[0].index))
p = idt.centroid(cfix)
f = EFixation(p, len(cfix), cfix[0], cfix[-1])
taggedEvents.append(f)
cfix = []
matchedSamples = 0
errorSamples = 0
for f in fixations:
s = f.start.index
for i in range(s,s+f.length):
if eventType[i] == 1:
matchedSamples = matchedSamples + 1
else:
errorSamples = errorSamples + 1
mPct = matchedSamples / float(len(eventType))
ePct = errorSamples / float(len(eventType))
print "Matched Samples: " + str(matchedSamples) + " (" + str(mPct * 100) + "%)"
print "Error Samples: " + str(errorSamples) + " (" + str(ePct * 100) + "%)"
print "============= I-VT SG test ==============="
stream = FileSampleStream('testData/UH27_img_vy_labelled_MN.txt')
fstream = SGFilter(IntersampleVelocity(BlinkFilter(stream)), 21, 2)
idt = Velocity(fstream, 0.005)
fixations = []
for i in idt:
print i
fixations.append(i)
verifStream = BlinkFilter(FileSampleStream('testData/UH27_img_vy_labelled_MN.txt'))
taggedEvents = []
cfix = []
eventType = []
for i in verifStream:
if i.eventType == 1:
cfix.append(i)
eventType.append(1)
else:
eventType.append(0)
if len(cfix) == 0:
continue
print ("Fixation of length: " + str(len(cfix)) + " starting at sample " + str(cfix[0].index))
p = idt.centroid(cfix)
f = EFixation(p, len(cfix), cfix[0], cfix[-1])
taggedEvents.append(f)
cfix = []
matchedSamples = 0
errorSamples = 0
for f in fixations:
s = f.start.index
for i in range(s,s+f.length):
if eventType[i] == 1:
matchedSamples = matchedSamples + 1
else:
errorSamples = errorSamples + 1
mPct = matchedSamples / float(len(eventType))
ePct = errorSamples / float(len(eventType))
print "Matched Samples: " + str(matchedSamples) + " (" + str(mPct * 100) + "%)"
print "Error Samples: " + str(errorSamples) + " (" + str(ePct * 100) + "%)"