Calculate Volume Profile, Point of Control (POC) & Value Area High/Low in Python

import pandas as pd
import numpy as np

def calculate_volume_profile(df, bins=50):
    df = df.copy()
    price_min, price_max = df['Low'].min(), df['High'].max()
    bin_edges = np.linspace(price_min, price_max, bins + 1)
    bin_centers = (bin_edges[:-1] + bin_edges[1:]) / 2
    
    binned = pd.cut(df['Close'], bins=bin_edges, labels=bin_centers, include_lowest=True)
    vol_profile = df.groupby(binned, observed=False)['Volume'].sum()
    poc_price = float(vol_profile.idxmax())
    return {'POC_Price': poc_price, 'Volume_Profile_Series': vol_profile}