python from astropy import units as u from astropy.coordinates import SkyCoord from astropy.io import fits from astropy.table import Table c = SkyCoord(ra=10.68458*u.degree, dec=41.26917*u.degree, frame='icrs') galactic = c.transform_to('galactic') data = fits.getdata('example.fits') table = Table(data) sorted_table = table[np.argsort(table['column_name'])]