hdulist[1] = table print("Axis:") for name in hdulist[1].columns.names: if name.islower() and not name == 'labs': ax = np.unique(hdulist[1].data.field(name)) print(name, len(ax), min(ax), max(ax)) teffs = hdulist[1].data.field('teff') loggs = hdulist[1].data.field('logg') ebvs = hdulist[1].data.field('ebv') #for logg in np.unique(loggs): #keep = loggs==logg #plt.figure() #plt.title(logg) #plt.plot(teffs[keep],ebvs[keep],'ko',ms=2) keep = hdulist[1].data.field('teff') > 0 logger.info('Removing {}/{} false entries'.format(sum(-keep), len(keep))) #print len(ff[1].data),sum(keep) hdulist[1].data = hdulist[1].data[keep] hdulist.close() #} if __name__ == "__main__": logger = loggers.get_basic_logger() imethod, iargs, ikwargs = argkwargparser.parse() output = globals()[imethod](*iargs, **ikwargs)
from ivs.aux import argkwargparser from ivs.inout import ascii #-- if no arguments are given, we just do a test run if not sys.argv[1:]: doctest.testmod() pl.show() sys.exit() #-- if arguments are given, we assume the user wants to run one of the # functions with arguments given in the command line # EXAMPLES: # $:> python freqanalyse.py find_frequency infile=test.dat full_output=True # $:> python freqanalyse.py time_frequency infile=test.dat full_output=True else: method,args,kwargs = argkwargparser.parse() print "Running method %s with arguments %s and keyword arguments %s"%(method,args,kwargs) if '--help' in args or 'help' in args or 'help' in kwargs: sys.exit() full_output = kwargs.get('full_output',False) times,signal = ascii.read2array(kwargs.pop('infile')).T[:2] out = globals()[method](times,signal, **kwargs) #-- when find_frequency is called if method=='find_frequency' and full_output: print pl.mlab.rec2txt(out[0],precision=8) pl.figure() pl.subplot(211) pl.plot(out[1][0],out[1][1],'k-') pl.subplot(212) pl.plot(times,signal,'ko',ms=2)
return (order[match_success],match_success),(order[-match_success],-match_success) def photometry2str(master): """ String representation of master record array @param master: master record array containing photometry @type master: numpy record array """ master = master[np.argsort(master['photband'])] txt = '#%19s %12s %12s %12s %10s %12s %12s %11s %s\n'%('PHOTBAND','MEAS','E_MEAS','UNIT','CWAVE','CMEAS','E_CMEAS','UNIT','SOURCE') txt+= '#=========================================================================================================================\n' for i,j,k,l,m,n,o,p in zip(master['photband'],master['meas'],master['e_meas'],master['unit'],master['cwave'],master['cmeas'],master['e_cmeas'],master['source']): txt+='%20s %12g %12g %12s %10.0f %12g %12g erg/s/cm2/AA %s\n'%(i,j,k,l,m,n,o,p) return txt if __name__=="__main__": logger = loggers.get_basic_logger() method,args,kwargs = argkwargparser.parse() print_help = '--help' in args or '-h' in args if print_help: help(globals()[method]) else: master = globals()[method](*args,**kwargs) print(photometry2str(master)) #xmatch_vizier('II/169/main','B/pastel/pastel')
table.header.update(key,hdulist[1].header[key]) hdulist[1] = table print "Axis:" for name in hdulist[1].columns.names: if name.islower() and not name=='labs': ax = np.unique(hdulist[1].data.field(name)) print name,len(ax),min(ax),max(ax) teffs = hdulist[1].data.field('teff') loggs = hdulist[1].data.field('logg') ebvs = hdulist[1].data.field('ebv') #for logg in np.unique(loggs): #keep = loggs==logg #plt.figure() #plt.title(logg) #plt.plot(teffs[keep],ebvs[keep],'ko',ms=2) keep = hdulist[1].data.field('teff')>0 logger.info('Removing {}/{} false entries'.format(sum(-keep),len(keep))) #print len(ff[1].data),sum(keep) hdulist[1].data = hdulist[1].data[keep] hdulist.close() #} if __name__=="__main__": logger = loggers.get_basic_logger() imethod,iargs,ikwargs = argkwargparser.parse() output = globals()[imethod](*iargs,**ikwargs)