traveltime_traveltime_per_pc4

traveltime_traveltime_per_pc4

pc = gpd.read_file(os.path.join('PC_4-shp','PC4.shp'))
df = pd.read_csv(OUTPUT)
gdf = gpd.GeoDataFrame(df, 
                       geometry=gpd.points_from_xy(df.EndLon, df.EndLat))
                            [['geometry', 'DurationMin']]

def get_avg_traveltime(pc4):
    v = gdf[gdf.within(pc4['geometry'] )].mean(skipna=True)
    if math.isnan(v):
        return np.NaN
    else:
        return int(v[0])
    
pc['Traveltime'] = pc.apply(lambda x: get_avg_traveltime(x), axis=1)