gaussian_smooth_line.py

def smooth_line(x, y, sigma=5, interp_points=300):
    from scipy.ndimage import gaussian_filter1d

    x_interp = np.linspace(x.min(), x.max(), interp_points)

    # The gaussian filter assumes spacing 1, so we have to adjust sigma
    interp_diff = x_interp[1] - x_interp[0]
    effective_sigma = sigma / interp_diff
    
    y_interp = np.interp(x_interp, x, y)
    y_smooth = gaussian_filter1d(y_interp, sigma=effective_sigma)

    return x_interp, y_smooth