use astropy for writing and fitsio for reading