manipulation de données binaires (optimisation)

**tubaas** · 02/07/2012, 10h21

Bonjour
je manipule des données binaires et je suis parvenu à un code fonctionnel mais pas efficace.
Voici mon problème : dans un fichier binaire, j'ai des blocs de 2048 données utiles codées sur 32bits. mais dont seuls les 3 premiers octets doivent être interprétés. Chaque bloc est précédé de 4 octets inintéressants. Voici un extrait du fichier en hexadécimal :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
08 04 00 00 70 F4 15 01 00 00 00 00 00 00 00 00  %inintéressant
DA FF FF 00 DA FF FF 00 E1 FF FF 00 DB FF FF 00 %données : pour les 4 premiers octets, la valeur est en fait DA FF FF codée en little endian soit 0xFFFFDA et ainsi de suite
D8 FF FF 00 D9 FF FF 00 D3 FF FF 00 D8 FF FF 00 
E1 FF FF 00 DD FF FF 00 DA FF FF 00 D9 FF FF 00 
DC FF FF 00 DC FF FF 00 DA FF FF 00 E1 FF FF 00 
E0 FF FF 00 E4 FF FF 00 E6 FF FF 00 D8 FF FF 00 
DB FF FF 00 E0 FF FF 00 DD FF FF 00 DD FF FF 00...

Si ja'i 2 voies (variable nbvoies dans le code ci-dessous), je dois ensuite entrelacer 24bits du premier bloc avec 24 du second 2048 fois (longueur des blocs) puis passer aux blocs 3 et 4 et ainsi de suite. Les 24bits représentant une valeur interprétable, je leur applique éventuellement un facteur multiplicatif.
Voici mon code fonctionnel :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
    fid=fopen(filepath,'r');
    % lecture et remise en ordre données audio en vue d'écriture fichier
    % wav
    wavdata=zeros(1,nbentrees*nbvoies*2048);
    for i=0:1:nbvoies*nbentrees-1
        fread(fid,16);
        start_ind=fix(i/2)*nbvoies*2048+mod(i,nbvoies)+1;
        for j=1:2048
            wavdata(1,start_ind+nbvoies*(j-1))=fread(fid,1,'bit24');
            fread(fid,1,'bit8');
        end
    end
    % GAIN
    wavdata=wavdata*2^pow;
    fclose(fid);
 
    %écriture des données dans un nouveau fichier
    fid2=fopen(filepath2,'w');
    fwrite(fid2,wavdata,'bit24');
    fclose(fid2);

Je me demande s'il n'y aurait pas une manière plus efficace, par exemple en récupérant chacun des octets dans une matrice, en réorganisant cette matrice et en interprétant les valeurs effectives codées en 24 bits seulement ensuite mais je ne sais pas trop comment m'y prendre.
J'espère avoir été suffisamment clair.
Merci !

**Jerome Briot** · 02/07/2012, 20h57

Essaie ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
C = {'08' '04' '00' '00' '70' 'F4' '15' '01' '00' '00' '00' '00' '00' '00' '00' '00'
'DA' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'E1' 'FF' 'FF' '00' 'DB' 'FF' 'FF' '00'
'D8' 'FF' 'FF' '00' 'D9' 'FF' 'FF' '00' 'D3' 'FF' 'FF' '00' 'D8' 'FF' 'FF' '00' 
'E1' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'D9' 'FF' 'FF' '00' 
'DC' 'FF' 'FF' '00' 'DC' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'E1' 'FF' 'FF' '00' 
'E0' 'FF' 'FF' '00' 'E4' 'FF' 'FF' '00' 'E6' 'FF' 'FF' '00' 'D8' 'FF' 'FF' '00' 
'DB' 'FF' 'FF' '00' 'E0' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00'
'08' '04' '00' '00' '70' 'F4' '15' '01' '00' '00' '00' '00' '00' '00' '00' '00'
'DA' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'E1' 'FF' 'FF' '00' 'DB' 'FF' 'FF' '00'
'D8' 'FF' 'FF' '00' 'D9' 'FF' 'FF' '00' 'D3' 'FF' 'FF' '00' 'D8' 'FF' 'FF' '00' 
'E1' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'D9' 'FF' 'FF' '00' 
'DC' 'FF' 'FF' '00' 'DC' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'E1' 'FF' 'FF' '00' 
'E0' 'FF' 'FF' '00' 'E4' 'FF' 'FF' '00' 'E6' 'FF' 'FF' '00' 'D8' 'FF' 'FF' '00' 
'DB' 'FF' 'FF' '00' 'E0' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00'
'08' '04' '00' '00' '70' 'F4' '15' '01' '00' '00' '00' '00' '00' '00' '00' '00'
'DA' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'E1' 'FF' 'FF' '00' 'DB' 'FF' 'FF' '00'
'D8' 'FF' 'FF' '00' 'D9' 'FF' 'FF' '00' 'D3' 'FF' 'FF' '00' 'D8' 'FF' 'FF' '00' 
'E1' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'D9' 'FF' 'FF' '00' 
'DC' 'FF' 'FF' '00' 'DC' 'FF' 'FF' '00' 'DA' 'FF' 'FF' '00' 'E1' 'FF' 'FF' '00' 
'E0' 'FF' 'FF' '00' 'E4' 'FF' 'FF' '00' 'E6' 'FF' 'FF' '00' 'D8' 'FF' 'FF' '00' 
'DB' 'FF' 'FF' '00' 'E0' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00' 'DD' 'FF' 'FF' '00'};
 
X = cellfun(@hex2dec,C);
 
fid = fopen('data.bin','w');
    fwrite(fid,X.','uint8');   
fclose(fid);
 
n = 1;
fid = fopen('data.bin','r');
while ~feof(fid)
    s = fseek(fid,16,'cof');
    if s==-1
        break
    end
    wavdata{n} = fread(fid,[4,6],'ubit24',8).';
    n = n+1;
end
fclose(fid);

**tubaas** · 03/07/2012, 09h21

Bonjour,
merci pour la réponse

. En fait c'est l'argument SKIP de la fonction FREAD que je n'avais pas compris. Par contre la prochaine fois évitez de me faire une frayeur en remplaçant mon 'bit24' par un 'ubit24'

(c'est un peu ce qui posait problème car en fait le quatrième octet est toujours nul donc si j'avais été en non signé, j'aurais pu l'interpréter comme un 'ubit32').
Pour les curieux voici le nouveau code (en rouge les lignes qui ont changé)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
    fid=fopen(filepath,'r');
    % lecture et remise en ordre données audio en vue d'écriture fichier
    % wav
    wavdata=zeros(1,nbentrees*nbvoies*2048);
    for i=0:1:nbvoies*nbentrees-1
        s = fseek(fid,16,'cof');
            if s==-1
            break
        end
        start_ind=fix(i/2)*nbvoies*2048+mod(i,nbvoies)+1;
        wavdata(1,start_ind:nbvoies:start_ind+2047*nbvoies) = fread(fid,2048,'bit24',8).';
    end
    % GAIN
    wavdata=wavdata*2^pow;
    fclose(fid);

Connaissant le nombre d’occurrences dans mon fichier, j'ai pensé qu'il était préférable de conserver ma boucle for à la boucle while qui ajoute un test.Pour être tout à fait complet, en faisant la moyenne du temps d'exécution de la boucle sur 15 exécutions je passe de 14,5s

avec la première méthode à 0,5s

avec la seconde.

PS : une petite question comme je reprends MATLAB après quelques temps d'interruption, je ne me rappelle plus la différence entre

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

a'

et

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

a.'

?

Encore merci.

**Jean Dumoncel** · 03/07/2012, 09h26

Bonjour,

voir la

:
Quelle est la différence entre ' et .' ?

**tubaas** · 03/07/2012, 09h29

Merci...

**Jerome Briot** · 03/07/2012, 09h42

Dans ce cas, la ligne avec FREAD se simplifie en :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wavdata(1,start_ind:nbvoies:start_ind+2047*nbvoies) = fread(fid,[1,2048],'bit24',8);

manipulation de données binaires (optimisation)

MATLAB

Discussions similaires

Partager

Partager