from scipy.spatial import ConvexHull
points = np.random.rand(30, 2)   # 30 random points in 2-D
hull = ConvexHull(points)

# The convex hull is represented as a set of N-1 dimensional simplices,
# which in 2-D means line segments. The storage scheme is exactly the
# same as for the simplices in the Delaunay triangulation discussed
# above.

# We can illustrate the above result:

import matplotlib.pyplot as plt
plt.plot(points[:,0], points[:,1], 'o')
for simplex in hull.simplices:
    plt.plot(points[simplex,0], points[simplex,1], 'k-')
plt.show()