vertices = skygrid.geodesic(4 * np.pi * u.sr / n_vertices_target,
                                class_='III')
n_vertices = len(vertices)

ax = plt.axes(projection='astro globe', center='0d 25d')
plt.suptitle('Class III')
ax.set_title(f'{n_vertices} vertices (goal was {n_vertices_target})')
ax.plot_coord(vertices, '.')
ax.grid()