def ConsistencyStrat2():
# Walk through all the lower directories
for root, dirs, files in os.walk(CSV_OUTPUT_DIR):
# We only care about directories that have files
if(len(files) != 0):
files.sort()
# The directories with files have the .csv files
# Get the pressure data from them and store in list
pressureList = []
testerName = ""
deviceName = ""
for fileName in files:
seed = int(fileName.split(":")[0])
#Only seeds above 1000 are extra long paths likely to create
#bit strings of atleast 128 bits
if( seed > 1000 ):
with open(join(root, fileName), "rb") as csvfile:
respReader = csv.reader(csvfile)
dataStarted = False
for row in respReader:
if(dataStarted):
pressureList.append(float(row[2]))
elif(row[0] != "ChallengeX" and row[0] != "X"):
testerName = row[2]
deviceName = row[3]
if(row[0] == "X" and row[1] == "Y" and row[2] == "PRESSURE"):
dataStarted = True
#Calculate moving average n = 5 for use in arbiter
movingAvg5 = simpleMovingAverage(pressureList, 5)
#Build the bit string and convert it to a bytearray() type for writing
#to a binary file
bitString = ""
for i in range(0, len(pressureList)):
bitString += str(arbiter(pressureList[i], movingAvg5[i]))
byteArr = convBitStrToByteArr(bitString)
#Copy the directory structure of the OutputCSVs folder
deviceName = basename(dirname(root))
testerName = basename(root)
#If the directory doesn't exist make it
att_path = os.path.join(GENERATED_OUTPUT_DIR, "ConsistencyStrat2", deviceName, testerName, fileName + ".bin")
if not os.path.exists(os.path.dirname(att_path)):
os.makedirs(os.path.dirname(att_path))
#Write binary file
with open(att_path, "wb") as outputFile:
outputFile.write(byteArr)
pressureList = []
testerName = ""
deviceName = ""
def Strat3():
# Walk through all the lower directories
for root, dirs, files in os.walk(CSV_OUTPUT_DIR):
# We only care about directories that have files
if(len(files) != 0):
files.sort()
# The directories with files have the .csv files
# Get the pressure data from them and store in list
pressureList = []
for fileName in files:
with open(join(root, fileName), "rb") as csvfile:
respReader = csv.reader(csvfile)
dataStarted = False
for row in respReader:
if(dataStarted):
pressureList.append(float(row[2]))
if(row[0] == "X" and row[1] == "Y" and row[2] == "PRESSURE"):
dataStarted = True
#Calculate moving average n = 5 for use in arbiter
movingAvg10 = simpleMovingAverage(pressureList, 10)
#Build the bit string and convert it to a bytearray() type for writing
#to a binary file
bitString = ""
for i in range(0, len(pressureList)):
bitString += str(arbiter(pressureList[i], movingAvg10[i]))
byteArr = convBitStrToByteArr(bitString)
#Copy the directory structure of the OutputCSVs folder
deviceName = basename(dirname(root))
testerName = basename(root)
#If the directory doesn't exist make it
att_path = os.path.join(GENERATED_OUTPUT_DIR, "Strat3", deviceName, testerName, "responseBinary")
if not os.path.exists(os.path.dirname(att_path)):
os.makedirs(os.path.dirname(att_path))
#Write binary file
with open(att_path, "wb") as outputFile:
outputFile.write(byteArr)
def NormalizedStrat4():
# Walk through all the lower directories
for root, dirs, files in os.walk(CSV_OUTPUT_DIR):
# We only care about directories that have files
if(len(files) != 0):
files.sort()
# The directories with files have the .csv files
# Get the pressure data from them and store in list
challengeList = []
pressureList = []
respList = []
testerName = ""
deviceName = ""
for fileName in files:
seed = int(fileName.split(":")[0])
#Only seeds above 2000 are the single line challenges
if( seed >= 2000 ):
print fileName
with open(join(root, fileName), "rb") as csvfile:
respReader = csv.reader(csvfile)
responseStarted = False
for row in respReader:
if(responseStarted):
respX = float(row[0])
respY = float(row[1])
respList.append((respX, respY)) #Append the points as x,y tuples
pressureList.append(float(row[2]))
elif(row[0] != "ChallengeX" and row[0] != "X"):
challX = float(row[0])
challY = float(row[1])
challengeList.append((challX, challY))
testerName = row[2]
deviceName = row[3]
if(row[0] == "X" and row[1] == "Y" and row[2] == "PRESSURE"):
responseStarted = True
challengeLength = 0
isX = False
for i in range(0, len(challengeList)-1):
yLength = fabs(challengeList[i][1] - challengeList[i+1][1])
xLength = fabs(challengeList[i][0] - challengeList[i+1][0])
if(xLength > yLength):
challengeLength = xLength
isX = True
else:
challengeLength = yLength
isX = False
normalizedDist = challengeLength/32.0
normalizedList = []
currPt = challengeList[0]
for i in range(0, 32):
normalizedList.append(currPt)
if(isX):
currPt = (currPt[0]+normalizedDist, currPt[1])
else:
currPt = (currPt[0], currPt[1]+normalizedDist)
normalizedPressList = []
print len(respList)
print len(pressureList)
for normPt in normalizedList:
if(isX):
closestLeftPt = respList[0]
closestRightPt = respList[-1]
closestLeftPtInd = 0
closestRightPtInd = -1
if(respList[0][0] > normPt[0]):
closestLeftPt = respList[-1]
closestRightPt = respList[0]
closestLeftPtInd = -1
closestRightPtInd = 0
for i in range(0, len(respList)):
if(respList[i][0] <= normPt[0]): #if it's to the left
if((normPt[0] - respList[i][0]) < (normPt[0] - closestLeftPt[0])):
closestLeftPt = respList[i]
closestLeftPtInd = i
if(respList[i][0] > normPt[0]): #if it's to the left
if((respList[i][0]-normPt[0]) < (closestRightPt[0]-normPt[0])):
closestRightPt = respList[i]
closestRightPtInd = i
if(closestRightPt == closestLeftPt):
closestLeftPtInd = closestLeftPtInd - 1
closestLeftPt = respList[closestLeftPtInd]
deltaX = closestRightPt[0]-closestLeftPt[0]
deltaY = closestRightPt[1]-closestLeftPt[1]
angle = fabs(atan(deltaY/deltaX))
deltaX1 = normPt[0]-closestLeftPt[0]
h1 = deltaX1/cos(angle)
h = deltaX/cos(angle)
h2 = h-h1
interpPress = (h1/h)*pressureList[closestLeftPtInd]+(h2/h)*pressureList[closestRightPtInd]
normalizedPressList.append(interpPress)
else:
closestLowerPt = respList[0]
closestUpperPt = respList[-1]
closestLowerPtInd = 0
closestUpperPtInd = -1
if(respList[0][1] > normPt[1]):
closestLowerPt = respList[-1]
closestUpperPt = respList[0]
closestLowerPtInd = -1
closestUpperPtInd = 0
for i in range(0, len(respList)):
if(respList[i][1] <= normPt[1]): #if it's below
if((normPt[1] - respList[i][1]) < (normPt[1] - closestLowerPt[1])):
closestLowerPt = respList[i]
closestLowerPtInd = i
if(respList[i][1] > normPt[1]): #if it's above
if((respList[i][1]-normPt[1]) < (closestUpperPt[1]-normPt[1])):
closestUpperPt = respList[i]
closestUpperPtInd = i
if(closestUpperPt == closestLowerPt):
closestLowerPtInd = closestLowerPtInd - 1
closestLowerPt = respList[closestLowerPtInd]
deltaX = closestUpperPt[0]-closestLowerPt[0]
deltaY = closestUpperPt[1]-closestLowerPt[1]
angle = fabs(atan(deltaX/deltaY))
deltaY1 = normPt[1]-closestLowerPt[1]
h1 = deltaY1/cos(angle)
h = deltaY/cos(angle)
h2 = h-h1
interpPress = (h1/h)*pressureList[closestLowerPtInd]+(h2/h)*pressureList[closestUpperPtInd]
normalizedPressList.append(interpPress)
#Calculate moving average n = 5 for use in arbiter
movingAvg5 = simpleMovingAverage(normalizedPressList, 40)
#Build the bit string and convert it to a bytearray() type for writing
#to a binary file
bitString = ""
for i in range(0, len(normalizedPressList)):
bitString += str(arbiter(normalizedPressList[i], movingAvg5[i]))
byteArr = convBitStrToByteArr(bitString)
#Copy the directory structure of the OutputCSVs folder
deviceName = basename(dirname(root))
testerName = basename(root)
#If the directory doesn't exist make it
att_path = os.path.join(GENERATED_OUTPUT_DIR, "NormalizedStrat4", deviceName, testerName, fileName + ".bin")
if not os.path.exists(os.path.dirname(att_path)):
os.makedirs(os.path.dirname(att_path))
#Write binary file
with open(att_path, "wb") as outputFile:
outputFile.write(byteArr)
pressureList = []
respList = []
challengeList = []
testerName = ""
deviceName = ""
def NormalizedStrat3():
# Walk through all the lower directories
for root, dirs, files in os.walk(CSV_OUTPUT_DIR):
# We only care about directories that have files
if(len(files) != 0):
files.sort()
# The directories with files have the .csv files
# Get the pressure data from them and store in list
challengeList = []
pressureList = []
respList = []
testerName = ""
deviceName = ""
for fileName in files:
seed = int(fileName.split(":")[0])
#Only seeds above 1000 are extra long paths likely to create
#bit strings of atleast 128 bits
if( seed > 1000 ):
print fileName
with open(join(root, fileName), "rb") as csvfile:
respReader = csv.reader(csvfile)
responseStarted = False
for row in respReader:
if(responseStarted):
respX = float(row[0])
respY = float(row[1])
respList.append((respX, respY)) #Append the points as x,y tuples
pressureList.append(float(row[2]))
elif(row[0] != "ChallengeX" and row[0] != "X"):
challX = float(row[0])
challY = float(row[1])
challengeList.append((challX, challY))
testerName = row[2]
deviceName = row[3]
if(row[0] == "X" and row[1] == "Y" and row[2] == "PRESSURE"):
responseStarted = True
challengeLength = 0
for i in range(0, len(challengeList)-1):
challengeLength += distBetweenPts(challengeList[i], challengeList[i+1])
normalizedDist = challengeLength/32.0
normalizedList = []
currPt = challengeList[0]
for i in range(0, len(challengeList)-1):
currPt = challengeList[i]
nextPt = challengeList[i+1]
while(distBetweenPts(currPt, nextPt) > normalizedDist):
normalizedList.append(currPt)
currPt = calcNextPt(currPt, nextPt, normalizedDist)
normalizedList.append(challengeList[-1])
normalizedPressList = []
print len(respList)
print len(pressureList)
for normPt in normalizedList:
closestPtIndex = 0
secondClosestPtIndex = 0
closestPt = respList[0]
secondClosestPt = respList[0]
closestPtDist = distBetweenPts(normPt, closestPt)
for i in range(0, len(respList)):
if(distBetweenPts(normPt, respList[i]) <= closestPtDist):
closestPtDist = distBetweenPts(normPt, respList[i])
secondClosestPt = closestPt
secondClosestPtIndex = closestPtIndex
closestPt = respList[i]
closestPtIndex = i
secondClosestPtDist = distBetweenPts(normPt, secondClosestPt)
distBetweenClosestPts = secondClosestPtDist + closestPtDist
normPressure = (closestPtDist/distBetweenClosestPts)*pressureList[closestPtIndex] + (secondClosestPtDist/distBetweenClosestPts)*pressureList[secondClosestPtIndex] #weighted average weighted towards the closer point
normalizedPressList.append(normPressure)
#Calculate moving average n = 5 for use in arbiter
movingAvg5 = simpleMovingAverage(normalizedPressList, 5)
#Build the bit string and convert it to a bytearray() type for writing
#to a binary file
bitString = ""
for i in range(0, len(normalizedPressList)):
bitString += str(arbiter(normalizedPressList[i], movingAvg5[i]))
byteArr = convBitStrToByteArr(bitString)
#Copy the directory structure of the OutputCSVs folder
deviceName = basename(dirname(root))
testerName = basename(root)
#If the directory doesn't exist make it
att_path = os.path.join(GENERATED_OUTPUT_DIR, "NormalizedStrat3", deviceName, testerName, fileName + ".bin")
if not os.path.exists(os.path.dirname(att_path)):
os.makedirs(os.path.dirname(att_path))
#Write binary file
with open(att_path, "wb") as outputFile:
outputFile.write(byteArr)
pressureList = []
respList = []
challengeList = []
testerName = ""
deviceName = ""
def NormalizedStrat2():
# Walk through all the lower directories
for root, dirs, files in os.walk(CSV_OUTPUT_DIR):
# We only care about directories that have files
if(len(files) != 0):
files.sort()
# The directories with files have the .csv files
# Get the pressure data from them and store in list
challengeList = []
pressureList = []
respList = []
testerName = ""
deviceName = ""
for fileName in files:
seed = int(fileName.split(":")[0])
#Only seeds above 1000 are extra long paths likely to create
#bit strings of atleast 128 bits
if( seed > 1000 ):
print fileName
with open(join(root, fileName), "rb") as csvfile:
respReader = csv.reader(csvfile)
responseStarted = False
for row in respReader:
if(responseStarted):
respX = float(row[0])
respY = float(row[1])
respList.append((respX, respY)) #Append the points as x,y tuples
pressureList.append(float(row[2]))
elif(row[0] != "ChallengeX" and row[0] != "X"):
challX = float(row[0])
challY = float(row[1])
challengeList.append((challX, challY))
testerName = row[2]
deviceName = row[3]
if(row[0] == "X" and row[1] == "Y" and row[2] == "PRESSURE"):
responseStarted = True
challengeLength = 0
for i in range(0, len(challengeList)-1):
challengeLength += distBetweenPts(challengeList[i], challengeList[i+1])
normalizedDist = challengeLength/32.0
normalizedList = []
currPt = challengeList[0]
for i in range(0, len(challengeList)-1):
currPt = challengeList[i]
nextPt = challengeList[i+1]
while(distBetweenPts(currPt, nextPt) > normalizedDist):
normalizedList.append(currPt)
currPt = calcNextPt(currPt, nextPt, normalizedDist)
normalizedList.append(challengeList[-1])
normalizedPressList = []
print len(respList)
print len(pressureList)
lastPtIndex = 0
for normPt in normalizedList:
lastIndexChanged = False
closestPt = respList[lastPtIndex]
closestPtDist = distBetweenPts(normPt, closestPt)
if(len(respList[lastPtIndex:lastPtIndex+100]) < 100):
lastPtIndex = len(respList) - 100 #I always want to ensure we atleast search through 20 points
for i in range(0, len(respList)):#lastPtIndex+100): #The assumption is the next closest point won't be 20 points away from the last one. This reduces search time and prevents the situation of two intersecting pressure lines mismatching
if(distBetweenPts(normPt, respList[i]) <= closestPtDist):
closestPtDist = distBetweenPts(normPt, respList[i])
closestPt = respList[i]
if(i != lastPtIndex):
lastPtIndex = i
lastIndexChanged = True
print "%s, lastPtInd = %d, normPt = %s, closestPt = %s, dist = %f" % (lastIndexChanged, lastPtIndex, normPt, closestPt, closestPtDist)
normalizedPressList.append(pressureList[lastPtIndex])
if not lastIndexChanged:
lastPtIndex = lastPtIndex+1
#Calculate moving average n = 5 for use in arbiter
movingAvg5 = simpleMovingAverage(normalizedPressList, 5)
#Build the bit string and convert it to a bytearray() type for writing
#to a binary file
bitString = ""
for i in range(0, len(normalizedPressList)):
bitString += str(arbiter(normalizedPressList[i], movingAvg5[i]))
byteArr = convBitStrToByteArr(bitString)
#Copy the directory structure of the OutputCSVs folder
deviceName = basename(dirname(root))
testerName = basename(root)
#If the directory doesn't exist make it
att_path = os.path.join(GENERATED_OUTPUT_DIR, "NormalizedStrat2", deviceName, testerName, fileName + ".bin")
if not os.path.exists(os.path.dirname(att_path)):
os.makedirs(os.path.dirname(att_path))
#Write binary file
with open(att_path, "wb") as outputFile:
outputFile.write(byteArr)
pressureList = []
respList = []
challengeList = []
testerName = ""
deviceName = ""