1""" Create example sphere points """
2
3import numpy as np
4
5
6def _make_pts():
7    """ Make points around sphere quadrants """
8    thetas = np.arange(1, 4) * np.pi/4
9    phis = np.arange(8) * np.pi/4
10    north_pole = (0, 0, 1)
11    south_pole = (0, 0, -1)
12    points = [north_pole, south_pole]
13    for theta in thetas:
14        for phi in phis:
15            x = np.sin(theta) * np.cos(phi)
16            y = np.sin(theta) * np.sin(phi)
17            z = np.cos(theta)
18            points.append((x, y, z))
19    return np.array(points)
20
21
22sphere_points = _make_pts()
23