new_table in astropy pyfits is deprecated so now I use BinTableHDU.from_columns to do the same thing