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