Project

General

Profile

Data format » History » Version 4

Anonymous, 03/05/2020 03:25 PM

1 1 Anonymous
h1. Data format
2
3 3 Anonymous
La fonction de lecture des données est la suivantes, pour un fichier:
4
5 1 Anonymous
<pre><code class="python">
6
def readFileADLINKCh0( filename, samplesPerFile, timeSerie_A ):
7
    fd = open(filename,'rb')
8
    
9
    dum = np.fromfile(fd, dtype = np.uint16)
10
11
    timeSerie_A[:] = dum[:]
12
    
13
    fd.close()
14
</code></pre>
15 3 Anonymous
16
On fait une boucle pour lire tous les fichiers et créer un grand vecteur:
17 1 Anonymous
18
<pre><code class="python">
19 4 Anonymous
timeSerie_A = np.zeros( samplesPerFile )
20
adc_A = np.zeros( (nbFiles, samplesPerFile) )
21 1 Anonymous
lastBuffer = (nbFiles-1) * buffersPerFile + firstBuffer
22
for k in range( nbFiles ):
23
    
24
    # read the data
25
    nb = str( int( k * buffersPerFile + firstBuffer) )
26
    filename = data_dir + "/record" + nb + ".bin"
27
    readFileADLINKCh0( filename, samplesPerFile, timeSerie_A )
28
    
29
    print("k = {}, block = {} / {}".format(k, nb, lastBuffer))
30
    
31
    adc_A[ k, : ] = (timeSerie_A + offset) * scalingFactor
32
</code></pre>
33 2 Anonymous
34
Les données étant lues, on remet en forme le vecteur pour ne garder que le signal correspondant aux rampes montantes:
35
36
<pre><code class="python">
37
A_reshaped = adc_A.reshape(nbFiles * rampsPerFile, samplesPerRamp)
38
</code></pre>