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)
