|
|||||||
|
|
|
![]() |
|
|
Strumenti |
|
|
#1 |
|
Junior Member
Iscritto dal: Jan 2011
Messaggi: 17
|
[Python] doppi valori sullo stesso asse x
Salve a tutti.
Ho sviluppato il seguente codice per plottare la figura in allegato: Codice:
rc('font',**{'family':'serif','serif':['Helvetica']})
rc('ps',usedistiller='xpdf')
rc('text', usetex=True)
#------------------------------------------------------
tmin=56203
tmax=56250
data=pyfits.open('http://heasarc.gsfc.nasa.gov/docs/swift/results/transients/weak/GX304-1.orbit.lc.fits')
time = data[1].data.field(0)/86400. + data[1].header['MJDREFF'] + data[1].header['MJDREFI']
rate = data[1].data.field(1)
error = data[1].data.field(2)
data.close()
cond = ((time > tmin) & (time < tmax)) | ((time + 132.5 > tmin) & (time + 132.5 < tmax)) | ((time + 265 > tmin) & (time + 265 < tmax) )
time=time[cond]
rate=rate[cond]
error=error[cond]
#..plot the current light curve
errorbar(time, rate, error, fmt='r.', capsize=0)
gca().xaxis.set_major_formatter(FormatStrFormatter('%5.1f'))
#..overplot a previous outburst
errorbar(time + 132.5, rate, error, fmt='b.', capsize=0)
errorbar(time + 265, rate, error, fmt='k.', capsize=0)
axis([tmin,tmax,-0.00,0.45])
xlabel('Time, MJD')
savefig("sync.eps",orientation='portrait',papertype='a4',format='eps')
os.system('gv sync.eps')
La figura rappresenta delle curve di luce (conteggi verso il tempo) di una sorgente periodica. I picchi rossi sono l'emissione piu' recente. Quelli blu quella passata, e quello nero quella prima ancora. Il tempo sull'asse x e' il tempo relativo alla curva di luce piu' recente. Quello che voglio fare io e' riportare sull'asse x anche i tempi degli altri due picchi, magari coi rispettivi colori, SOTTO i numeri che gia' stanno sull'asse x. In pratica, ad esempio, sotto il valore 56210 deve starci 56110 (diverso colore), e lo stesso per gli altri ticks. Qualche suggerimento please? |
|
|
|
|
| Strumenti | |
|
|
Tutti gli orari sono GMT +1. Ora sono le: 11:55.



















