for lat in range(-90, 100, 10):
    for lon in range(-180, 190, 10):
        point = np.radians([lon, lat])
round_trip = vectorToGeogr(sphericalToVector(point))
assert np.allclose(point, round_trip)
Matplotlib custom projection: How to transform points