unew = np.arange(0, 1.01, 0.01)
out = interpolate.splev(unew, tck)
plt.figure()
plt.plot(x, y, out[0], out[1])
plt.show()
Approximating data with a multi segment cubic bezier curve and a distance as well as a curvature contraint