from ligo.skymap.tool.ligo_skymap_plot_volume import main
from astropy.utils.data import download_file
filename = download_file('https://gwosc.org/s/events/GW170104/P1500227/bayestar.fits.gz', cache=True)
main([filename, '--radecdist', '135', '37', '500'])