manipulation de données binaires (optimisation)

Version imprimable

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:

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:

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 !

Essaie ceci :

Code:

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);

Bonjour,
merci pour la réponse :ccool:. 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:

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 :calim2: avec la première méthode à 0,5s 8-) 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:

a'

et

Code:

a.'

?

Encore merci.

03/07/2012, 09h26
Jean Dumoncel

Bonjour,

voir la :faq: :
Quelle est la différence entre ' et .' ?
03/07/2012, 09h29
tubaas

Merci...
03/07/2012, 09h42
Jerome Briot

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

Code:

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