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