def _rnd_interpolate(xy, num_points, ordered=False):
tck, u = splprep([
xy[:, 0],
xy[:, 1]],
s=0
)
unew = random(num_points)
if sort:
unew = sort(unew)
out = splev(unew, tck)
return column_stack(out)