import pyfar as pf
coords = pf.samplings.sph_lebedev(sh_order=10)
result = coords.get_nearest_sph(0, 0, 1, 45, show=True)