base_url = "https://www.mpi-hd.mpg.de/hfm/HESS/hgps/data/"
for filename in hgps_filenames:
url = base_url + filename
path = hgps_data_path / filename
if path.exists():
print("Already downloaded: {}".format(path))
else:
print("Downloading {} to {}".format(url, path))
urlretrieve(url, str(path))
hgps_data_path.mkdir(parents=True, exist_ok=True)
hgps_data_download()
print("\n\nFiles at {} :\n".format(hgps_data_path.absolute()))
for path in hgps_data_path.iterdir():
print(path)
# ## Catalog with Astropy
#
# ### FITS file content
#
# Let's start by just opening up `hgps_catalog_v1.fits.gz` and looking at the content.
#
# Note that ``astropy.io.fits.open`` doesn't work with `Path` objects yet,
# so you have to call `str(path)` and pass a string.
# In[ ]:
path = hgps_data_path / "hgps_catalog_v1.fits.gz"
hdu_list = fits.open(str(path))
