Distance_calculator.py

#################
# Name: Jason Grant
# ID: 12430732
# Description: IOT Assignment 2: Calculates the Summary Items (Distance, Time, Avg Speed)
#################

from math import radians, cos, sin, asin, sqrt
import glob
import os 
from datetime import datetime, date
import BlynkLib
import os.path
import subprocess


BLYNK_AUTH = '_1YBrbat_TJksBX_p4ni9jz5gr3q62so'
# initialize Blynk
blynk = BlynkLib.Blynk(BLYNK_AUTH)

now = datetime.now() # current date and time
date_time = now.strftime("%m_%d_%Y_%H_%M_%S")
##dateTimeObj = datetime.strptime(date_time, "%m_%d_%Y_%H_%M_%S").time()

##Calculate the distance (in various units) between two points on Earth using their latitude and longitude.
def haversine(lon1, lat1, lon2, lat2):
    lon1 = radians(lon1)
    lat1 = radians(lat1)
    lon2 = radians(lon2)
    lat2 = radians(lat2)

    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
    distance = 2 * asin(sqrt(a)) * 6371 #6371 is the radius of the Earth
    return distance

##https://stackoverflow.com/questions/903853/how-do-you-extract-a-column-from-a-multi-dimensional-array
def column(matrix, i):
    return [row[i] for row in matrix]

def calctime(time):
    date_Fudge = date.today()
    datetime1 = datetime.combine(date_Fudge, time[0])
    datetime2 = datetime.combine(date_Fudge, time[-1])
    time_elapsed = datetime2 - datetime1
    #timeRes = time[-1] - time[0]
    print(time_elapsed)
    return time_elapsed


def calcDistanceDelta_Sum(data):
    latCol = column(data,0)
    lonCol = column(data,1)
    x=0
    distCalcList = []
    while x<(len(latCol)-1):
        distanceCalc = haversine(lonCol[x], latCol[x], lonCol[x+1], latCol[x+1])
        ###print(distanceCalc)
        distCalcList.append(distanceCalc)
        x+=1

    distanceTotal = sum(distCalcList)
    ##print(distCalcList)
    ##print(distanceTotal)
    return distanceTotal, distCalcList

def avgSpeed(data, timeCalcRes):
    totalDist, intervDist = calcDistanceDelta_Sum(data)
    timeCalcRes_secs = timeCalcRes.total_seconds()
    ##print(timeCalcRes_secs)
    avgSpeed = totalDist / timeCalcRes_secs
    return avgSpeed

def summaryToHTML(time, dist, avgSpeed, date_time ):
    index= open("index.html", 'r').read().format(htmlTime = time,htmldist=dist, htmlAvgSpeed=avgSpeed, htmlDate_Time=date_time)
    print(index)
    
    datareplace = open("index.html", 'wt')
    datareplace.write(index)
    datareplace.close

##https://stackoverflow.com/questions/39327032/how-to-get-the-latest-file-in-a-folder##
list_of_files = glob.glob('/home/pi/Assignment_2/ResultsFolder/*.txt') # * means all if need specific format then *.csv 
latest_file = max(list_of_files, key=os.path.getctime)

ins = open(latest_file, 'r')
data = []
timedata = []
for line in ins:
    number_strings = line.split()                                   # Split the line on runs of whitespace
    numbers = [float(n) for n in number_strings]                    # Convert to floats
    data.append(numbers)                                            # Add the "row" to your list.
    ##print(data) # [[1, 3, 4], [5, 5, 6]]
ins.close   
for x in data:
    strngTime = str(x[:][2])
    dateTimeObj = datetime.strptime(strngTime, '%H%M%S.0').time()
    timedata.append(dateTimeObj)
    

timeCalcRes = calctime(timedata)
timeCalcRes_secs = timeCalcRes.total_seconds()
print('time calc res :' + str(timeCalcRes_secs))
totalDist, distDelta = calcDistanceDelta_Sum(data)
print('Total Dist : ' + str(totalDist))

avgSpeedRes = avgSpeed(data, timeCalcRes)
print('AVG Speed :  '+ str(avgSpeedRes))

file2 = open(f'ResultsFolder/GPS_Stats_1.txt', 'w')
file2.write(f"Time: {timeCalcRes_secs:.5f} Total Distance: {totalDist:.5f} Average Speed:{avgSpeedRes:.5f} \n")
file2.close
file2.seek(0)
blynk.virtual_write(1, totalDist)
blynk.virtual_write(3, timeCalcRes_secs)
blynk.virtual_write(4, avgSpeedRes)
summaryToHTML(timeCalcRes_secs, totalDist, avgSpeedRes, date_time)
process = subprocess.run(f"python3 fileopenandRead_test.py {date_time}", shell = True)