def getDirectionality(track, **kwargs): FIELD_VECTOR = kwargs['gradientVector'] xDiffs = gen.getDifferenceArray(track.x) yDiffs = gen.getDifferenceArray(track.y) directedness = [] #calculate instantenous directionality between frames for dx, dy in zip(xDiffs, yDiffs): vector = [dx, dy] normVector = gen.normalizeVector(vector) if normVector != 0: normFieldVector = gen.normalizeVector(FIELD_VECTOR) directionality = gen.dotProduct(vector, FIELD_VECTOR) / normVector * normFieldVector directedness.append(directionality) else: directedness.append(0) return gen.nanMean(directedness)
