@@ -0,0 +1,817 @@
+#!/usr/bin/env python3%0A# -*- coding: utf-8 -*-%0A%0A# Documentation: http://docs.astropy.org/en/stable/io/fits/index.html%0A%0Aimport argparse%0Afrom astropy.io import fits%0A%0A%0A# PARSE OPTIONS ###############################################################%0A%0Aparser = argparse.ArgumentParser(description=%22An astropy snippet%22)%0A%0Aparser.add_argument(%22filearg%22, nargs=1, metavar=%22FILE%22, help=%22the FITS file to process%22)%0A%0Aargs = parser.parse_args()%0A%0Afile_path = args.filearg%5B0%5D%0A%0A%0A# READ DATA ###################################################################%0A%0A# Open the FITS file%0Ahdu_list = fits.open(file_path)%0A%0A# Print the content of the FITS file (HDU headers)%0Ahdu_list.info()%0A%0Afor hdu in hdu_list:%0A    data = hdu.data   # %22hdu.data%22 is a Numpy Array%0A    print(data)%0A    print(data.shape)%0A%0A# Close the FITS file%0Ahdu_list.close()%0A%0A

Add a snippet (Python/Astropy).