debianforum.de

Markus Feldmann 2009-05-26 extraktpeaks3

SNIPPET_TEXT:

1. def extrahierePeaks(self,starttime,endtime,timescale,nitems):
2.         f_output = open('../../../peakswerteundzeit.dat','wb')
3.         f_p_1 = open('../../../peakdetektor1.hex','rb')
4.         f_p_2 = open('../../../peakdetektor2.hex','rb')
5.         f_p_3 = open('../../../peakdetektor3.hex','rb')
6.         f_pv_1 = open('../../../peakband1.hex','rb')
7.         f_pv_2 = open('../../../peakband2.hex','rb')
8.         f_pv_3 = open('../../../peakband3.hex','rb')
9.         a_p_1 = numpy.fromfile(f_p_1,dtype='b')
10.         a_p_2 = numpy.fromfile(f_p_2,dtype='b')
11.         a_p_3 = numpy.fromfile(f_p_3,dtype='b')
12.         a_pv_1 = numpy.fromfile(f_pv_1,dtype=numpy.complex64)
13.         a_pv_2 = numpy.fromfile(f_pv_2,dtype=numpy.complex64)
14.         a_pv_3 = numpy.fromfile(f_pv_3,dtype=numpy.complex64)
15.         criterion = (a_p_1==1) | (a_p_2==1) | (a_p_3==1)
16.         i = 0
17.         a_time = []
18. #         print a_p_1[0:4]
19. #         print a_p_2[0:4]
20. #         print a_p_3[0:4]
21. #         print a_pv_1[0:4]
22. #         print a_pv_2[0:4]
23. #         print a_pv_3[0:4]
24.         #print "Timescale*1e6 ist ",timescale*10000
25.         while i <= nitems-1:
26.             if (a_p_1[i] or a_p_2[i] or a_p_3[i]):
27.                 a_time += [timescale*i*10000]
28.             i += 1
29. #         print a_time[0:4]
30.         cPickle.dump((a_time[criterion],
31.                       a_p_1[criterion],
32.                       a_p_2[criterion],
33.                       a_p_3[criterion],
34.                       a_pv_1[criterion],
35.                       a_pv_2[criterion],
36.                       a_pv_3[criterion]),f_output)
37.         f_output.close()

Quellcode

Hier kannst du den Code kopieren und ihn in deinen bevorzugten Editor einfügen. PASTEBIN_DOWNLOAD_SNIPPET_EXPLAIN