def write_absid_file(self, outfil=None): from astropy.table import Column from astropy.table.table import Table wrest = self.lines.keys() wrest.sort() if outfil is None: outfil = self.absid_file # Columns cols = [Column(np.array(wrest), name='WREST')] clm_nms = self.lines[wrest[0]].analy.keys() for clm_nm in clm_nms: clist = [self.lines[iwrest].analy[clm_nm] for iwrest in wrest] cols.append( Column(np.array(clist), name=clm_nm) ) cols.append( Column(np.ones(len(cols[0]))*self.zabs, name='ZABS') ) table = Table(cols) prihdr = fits.Header() prihdr['COMMENT'] = "Above are the data sources" prihdu = fits.PrimaryHDU(header=prihdr) table_hdu = fits.BinTableHDU.from_columns(np.array(table.filled())) thdulist = fits.HDUList([prihdu, table_hdu]) thdulist.writeto(outfil,clobber=True) print('Wrote AbsID file: {:s}'.format(outfil))
def write_absid_file(self, outfil=None): from astropy.table import Column from astropy.table.table import Table wrest = self.lines.keys() wrest.sort() if outfil is None: outfil = self.absid_file # Columns cols = [Column(np.array(wrest), name='WREST')] clm_nms = self.lines[wrest[0]].analy.keys() for clm_nm in clm_nms: clist = [self.lines[iwrest].analy[clm_nm] for iwrest in wrest] cols.append( Column(np.array(clist), name=clm_nm) ) cols.append( Column(np.ones(len(cols[0]))*self.zabs, name='ZABS') ) table = Table(cols) prihdr = fits.Header() prihdr['COMMENT'] = "Above are the data sources" prihdu = fits.PrimaryHDU(header=prihdr) table_hdu = fits.BinTableHDU.from_columns(np.array(table.filled())) thdulist = fits.HDUList([prihdu, table_hdu]) thdulist.writeto(outfil,clobber=True) print('Wrote AbsID file: {:s}'.format(outfil))
def write_absid_file(self, outfil=None): '''TO BE DEPRECATED (probably) Writes portions of the LLS lines to a FITS table. ''' from astropy.table import Column from astropy.table.table import Table # #wrest = self.lines.keys() #wrest.sort() wrest = [line.wrest.value for line in self.lines] if outfil is None: outfil = self.absid_file # Columns cols = [Column(np.array(wrest), name='WREST')] clm_nms = self.lines[0].analy.keys() for clm_nm in clm_nms: if clm_nm == 'spec': continue clist = [line.analy[clm_nm] for line in self.lines] cols.append( Column(np.array(clist), name=clm_nm) ) cols.append( Column(np.ones(len(cols[0]))*self.zabs, name='ZABS') ) #from PyQt4 import QtCore #QtCore.pyqtRemoveInputHook() #xdb.set_trace() #QtCore.pyqtRestoreInputHook() table = Table(cols) prihdr = fits.Header() prihdr['COMMENT'] = "Above are the data sources" prihdu = fits.PrimaryHDU(header=prihdr) table_hdu = fits.BinTableHDU.from_columns(np.array(table.filled())) thdulist = fits.HDUList([prihdu, table_hdu]) thdulist.writeto(outfil,clobber=True) print('Wrote AbsID file: {:s}'.format(outfil))
def write_absid_file(self, outfil=None): '''TO BE DEPRECATED (probably) Writes portions of the LLS lines to a FITS table. ''' from astropy.table import Column from astropy.table.table import Table # #wrest = self.lines.keys() #wrest.sort() wrest = [line.wrest.value for line in self.lines] if outfil is None: outfil = self.absid_file # Columns cols = [Column(np.array(wrest), name='WREST')] clm_nms = self.lines[0].analy.keys() for clm_nm in clm_nms: if clm_nm == 'spec': continue clist = [line.analy[clm_nm] for line in self.lines] cols.append(Column(np.array(clist), name=clm_nm)) cols.append(Column(np.ones(len(cols[0])) * self.zabs, name='ZABS')) #from PyQt4 import QtCore #QtCore.pyqtRemoveInputHook() #xdb.set_trace() #QtCore.pyqtRestoreInputHook() table = Table(cols) prihdr = fits.Header() prihdr['COMMENT'] = "Above are the data sources" prihdu = fits.PrimaryHDU(header=prihdr) table_hdu = fits.BinTableHDU.from_columns(np.array(table.filled())) thdulist = fits.HDUList([prihdu, table_hdu]) thdulist.writeto(outfil, clobber=True) print('Wrote AbsID file: {:s}'.format(outfil))