Markus Feldmann 2009-05-26 extraktpeaks2
von Anonymous- SNIPPET_TEXT:
-
- def extrahierePeaks(self,starttime,endtime,timescale,nitems):
- f_output = open('../../../peakswerteundzeit.dat','wb')
- f_p_1 = open('../../../peakdetektor1.hex','rb')
- f_p_2 = open('../../../peakdetektor2.hex','rb')
- f_p_3 = open('../../../peakdetektor3.hex','rb')
- f_pv_1 = open('../../../peakband1.hex','rb')
- f_pv_2 = open('../../../peakband2.hex','rb')
- f_pv_3 = open('../../../peakband3.hex','rb')
- a_p_1 = numpy.fromfile(f_p_1,dtype='b',count=-1)
- a_p_2 = numpy.fromfile(f_p_2,dtype='b',count=-1)
- a_p_3 = numpy.fromfile(f_p_3,dtype='b',count=-1)
- a_pv_1 = numpy.fromfile(f_pv_1,dtype='f',count=-1)
- a_pv_2 = numpy.fromfile(f_pv_2,dtype='f',count=-1)
- a_pv_3 = numpy.fromfile(f_pv_3,dtype='f',count=-1)
- i = 0
- a_time = []
- a_a = []
- a_b = []
- a_c = []
- a_ar = []
- a_ai = []
- a_br = []
- a_bi = []
- a_cr = []
- a_ci = []
- while i <= nitems-1:
- if (a_p_1[i] or a_p_2[i] or a_p_3[i]):
- a_time += [starttime+timescale*i]
- a_a.append(a_p_1[i])
- # print "Mit i ist ",i," ergibt sich für a_p_1[i] ",a_p_1[i],\
- # "und für a_a ",a_a[:]
- a_b.append(a_p_2[i])
- a_c.append(a_p_3[i])
- a_ar.append(a_pv_1[i*2])
- a_ai.append(a_pv_1[i*2+1])
- a_br.append(a_pv_2[i*2])
- a_bi.append(a_pv_2[i*2+1])
- a_cr.append(a_pv_3[i*2])
- a_ci.append(a_pv_3[i*2+1])
- i += 1
- cPickle.dump(a_time,f_output)
- cPickle.dump(a_a,f_output)
- cPickle.dump(a_b,f_output)
- cPickle.dump(a_c,f_output)
- cPickle.dump(a_ar,f_output)
- cPickle.dump(a_ai,f_output)
- cPickle.dump(a_br,f_output)
- cPickle.dump(a_bi,f_output)
- cPickle.dump(a_cr,f_output)
- cPickle.dump(a_ci,f_output)
- f_output.flush()
- f_output.close()
Quellcode
Hier kannst du den Code kopieren und ihn in deinen bevorzugten Editor einfügen. PASTEBIN_DOWNLOAD_SNIPPET_EXPLAIN