1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73
|
#!/usr/bin/python
# Programme tracé de courbe OHZ
############### IMPORTS MODULES ###############################
import matplotlib.pyplot as plt
import numpy as np
import os
import scipy as sc
import Tkinter # Gestion import de fichiers
import tkFileDialog #
from xlwt import Workbook
from os.path import basename, splitext
from matplotlib.pyplot import figure,show
############### CHARGE LES FICHIERS ##########################
os.chdir=('/home/benjamin/Documents/doc_python/')
root= Tkinter.Tk()
path = tkFileDialog.askopenfilename(parent=root,title='Please select a file') # Selection du fichier
Data=np.loadtxt(filename,skiprows=7) # Lecture des donnees
filename= splitext(basename(path))[0] # Extraction du nom du fichier
num_voie=filename[1:3] # determination de la voie qui emet
############### REPRESENTATION ###############################
fech=5e6 # Hz
Npp=2048
t_final=Npp/fech
time=np.arange(0,t_final,1/fech)
choix_voie_visu=raw_input("Entrez la valeur de la voie a visualiser : ")
os.system("pause")
fig = figure(1) # Trace de la voie qui emet et des voies qui recoient
subplot(211)
plt.plot(time,Data[10:,num_voie]) ,xlabel('Temps') ,ylabel('Amplitude')
title('Signaux lors emission sur la voie '+ num_voie)
def onpick(event):
thisline = event.artist
xdata = thisline.get_xdata()
ydata = thisline.get_ydata()
ind = event.ind
print 'onpick points:', zip(xdata[ind],ydata[ind])
fig.canvas.mpl_connect('pick_event',onpick)
plt.show()
subplot(212)
#os.system("pause")
plt.plot(time,Data[10:,choix_voie_visu],'r')
time_S_min=1.2*time_P # Détermination de l'intervalle de temps arrivée onde S
time_S_max=1.6*time_P
figure(2)
plt.plot(time,Data[10:,choix_voie_visu],'r')
hold(True)
axvline(time_S_min,color='r') # Definition de l'intervalle
axvline(time_S_max,color='r')
############### SAUVEGARDE ###############################
resultat_pointeP = open("resultat_pointeP.txt","w") # Ouvre le fichier des resutats
resultat_pointeP.write("haha ca fonctionne ! et oui !") # Ecriture dans le fichier
resultat_pointeP.close() # Ferme le fichier
# os.chmod('path','mode') # Change les droits d'acces |
Partager