from scipy.spatial import ConvexHull # ... generate 'points', a 2D numpy array hull = ConvexHull(points) # to plot the hull vertices = hull.vertices plt.plot(points[vertices,0], points[vertices,1], 'r', lw = 2)