pip install astropy python from astropy.io import fits data = fits.getdata('data.fits') python shape = data.shape python mean_value = np.mean(data) python from astropy import coordinates as coord from astropy import units as u coord_obj = coord.SkyCoord(ra=10.68458*u.degree, dec=41.26917*u.degree, frame='icrs') galactic_coord = coord_obj.transform_to(coord.Galactic)