Big Sur 2021-09-25 thru 26

import pandas as pd
from gpx_converter import Converter

from haversine import haversine, Unit

files = [
    'big-sur-2021-09-25-thru-26/salmon-creek-station-to-north-buckeye.gpx',
    'big-sur-2021-09-25-thru-26/north-buckeye-camp-to-villa-creek-camp.gpx',
    'big-sur-2021-09-25-thru-26/track-92621-100917am.gpx'
]

def transform_gpx_data(filename):
    df = (Converter(input_file=filename)
          .gpx_to_dataframe())
    df = df[df['time'].dt.day == df['time'].dt.day.values[0]].reset_index(drop=True)
    df['time'] = df['time'].apply(lambda x: x.tz_convert('US/Pacific'))
    df['seconds_delta'] = (((df['time'].shift(-1)-df['time'])
                            .fillna(pd.Timedelta(seconds=0))
                            .astype(int)/1000000000)
                           .astype(int))
    df['human_date'] = df['time'].dt.strftime('%Y-%m-%d')
    df['human_time'] = df['time'].dt.strftime('%I:%M:%S %p')
    
    df['altitude_feet'] = round(df['altitude'] * 3.280839895).astype('int')
    
    for i in range(df.shape[0]-1):
        start = df.at[i,   'latitude'], df.at[i,   'longitude']
        end =   df.at[i+1, 'latitude'], df.at[i+1, 'longitude']
        distance = round(haversine(start,
                                   end,
                                   unit=Unit.FEET),1)
        df.at[i, 'distance_feet'] = distance

        altitude_change = df.at[i+1, 'altitude_feet'] - df.at[i, 'altitude_feet']
        df.at[i, 'altitude_change'] = altitude_change
    
    df['speed_mph'] = ((df['distance_feet'] / df['seconds_delta']) * (3600/5280)).round(1)
    
    df = df[['time', 'human_date', 'human_time', 'seconds_delta',
             'latitude', 'longitude', 'altitude', 'altitude_feet',
             'distance_feet', 'speed_mph', 
             'altitude_feet','altitude_change']].copy()
    return df

for file in files:
    print(file)
    d = transform_gpx_data(file)
    string = '''{}, {} - {}
Start: ({}, {}),   End: ({}, {})
{} GPS datapoints
{} duration
{:3.2f} miles @ {:3.1f} avg MPH
{:3.0f}/{:3.0f} feet total/net elevation change'''.format(
        d['human_date'].min(), 
        d.iloc[0]['human_time'], 
        d.iloc[d.shape[0]-1]['human_time'], 
        d.iloc[0]['latitude'], 
        d.iloc[0]['longitude'],
        d.iloc[d.shape[0]-1]['latitude'], 
        d.iloc[d.shape[0]-1]['longitude'],
        d.shape[0],
        str(d['time'].max()-d['time'].min())[7:12],
        d['distance_feet'].sum()/5280, 
        d['speed_mph'].mean(),
        d['altitude_change'].abs().sum(), 
        d['altitude_change'].sum())
    print(string + '\n')