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

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