Discussion :

[Newenda]

Bonjour,

Voici mon petit problème. J'ai un fichier .raw (ouverture avec Excel, WordPad etc..) avec une header de 32 lignes avant 4 colonnes de données mais également avec un "header de fin"...

Ce qui fait que quand j'essaye de l'ouvrir de la sorte :

Code :

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

[DATA{1:4}] = textread('rawdata0001.raw','%f%f%f%f','delimiter','; \b\t','headerlines',32);

il m'envoit le message d'erreur suivant :

??? Error using ==> dataread
Trouble reading floating point number from file (row 511, field
1) ==> FFT Successful\n

Error in ==> textread at 176
[varargout{1:nlhs}]=dataread('file',varargin{:});

où FFT Successful représente bien la 1ère ligne du Header de fin.

Le but est simple : faire abstraction du "header de fin"

Merci pour votre aide (et puis bonne année aussi!)

[Jerome Briot]

Je ferais comme ceci :

Code :

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

1
2
3
4
X = textread('rawdata0001.raw','%s','delimiter','\n','headerlines',32);
 
idx = find(strncmp(X,'FFT',3));
X(idx:end) = [];

Il "suffit" ensuite de traiter le tableau de cellules X

[Newenda]

Merci ! Bon j'ai essayé ensuite de traiter X mais c'est un cauchemars !

>> X{1}()

ans =

0 5.578701e-04 0.000000e+00 4.376011e-05 2.406050e+02

>> X{45}

ans =

44 5.291265e-04 -6.214043e-06 3.739349e-05 2.377500e+02

tous les valeurs d'une même colonne de X non pas forcément le même nombre de chiffre significatif, et quand les valeurs sont négatives ça ajouter un caractère, les espaces entre les chaque colonne n'est pas non plus forcément égal à 1.

Bref impossible de séparer les valeurs par indexage numérique. Il faudrait un indexage logique du type: séparer les nombres dont les caractères ne se touchent pas mais là je coince. Je n'arrive même pas à défnir les espaces :

>> X{1}(2)

ans =

Merci d'avance !

[tubaas]

salut
si dans une ligne le nombre de décimal est le même je te propose ceci à adapter

Code :

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

1
2
3
a='gqyri rqgio hoqsf';
a=reshape([a ' ']',find(isspace(a),1),(length(a)+1)/find(isspace(a),1))'
a=a(:,1:end-1)

mais il doit être plus efficace de le faire directement à la lecture en utilisant TEXTSCAN avec comme délimiteur de fin de ligne ('endofline') '\n' et d'utiliser les espace comme délimiteurs ('delimiter').
tu auras alors un cell array à deux dimensions

[Jerome Briot]

Il serait préférable de commencer par une méthode de plus bas niveau :

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
fid = fopen('rawdata0001.raw','rt');
 
for n=1:32
    fgetl(fid);
end
 
n = 1;
while ~feof(fid)
 
    str = fgetl(fid);
 
    if strncmp(str,'FFT',3)
        break
    end
 
    % Extraction des données à partir de str ici
    data(n,:) = sscanf(str,'%f',[1,4]);
    n = n+1;
 
end
 
fclose(fid);

Enfin c'est l'idée de base... essaies avec tes propres données.
Tu pourras ensuite passer à une version plus "optimisée"

Si vraiment tu n'y arrives pas, joint un fichier raw dans une archive zip ou rar à ton prochain message

[Newenda]

Merci Tubaas, j'imagine quand tu parles de nombre de décimal tu veux parler également du nombre d'espace entre les valeurs: malhereusement celui ci n'est pas identique à chaque ligne, c'est à dire que si je rajoute un espace (en rouge) :

Code :

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

1
2
3
a='gqyri_ rqgio hoqsf';
a=reshape([a ' ']',find(isspace(a),1),(length(a)+1)/find(isspace(a),1))'
a=a(:,1:end-1)

plus rien ne marche...

Merci Dut mais problème également :

??? Subscripted assignment dimension mismatch.

Error in ==> RheoOSC_dec2010_auto at 41
data(n, = sscanf(str,'%f',[1,4]);

Je vous joint un fichier à ouvrir... merci...désolé..

[tubaas]

ce que je voulais te montrer c'est qu'une fois que tu as lu le fichier tu peux toujours t'en tirer. à toi de construire ton algorithme.
par exemple au lieu d'utiliser

Code :

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

find(isspace(a),1)

tu peux utiliser

Code :

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

find(isspace(a))

et récupèrer ce qu'il y a entre les espaces.
Cependant le mieux est de le faire dès la lecture du fichier mais je te laisse voir avec Dut pour ça.

[Jerome Briot]

Merci Dut mais problème également :

Il y a une ligne remplie de 4 tabulations avant l'entete final.

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
fid = fopen('rawdata0001.raw','rt');
 
for n=1:31
    fgetl(fid);
end
 
n = 1;
while ~feof(fid)
 
    str = fgetl(fid);
 
    if ~isstrprop(str(1),'digit')
        break
    end
 
    data(n,:) = sscanf(str,'%f',[1,4]);
    n = n+1;
 
end
 
fclose(fid);

Par contre, si tu connais le nombre d'échantillons à l'avance (ici 256) il y a encore plus simple...

Et sinon il n'y a pas 4 mais 5 colonnes de données... si tu ne veux garder que les 4 dernières :

Code :

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

data(n,:) = sscanf(str,'%*d%f%f%f%f',[1,4]);

[Newenda]

Merci Dut, ça marche parfaitement... Merci Tubaas aussi, on en apprends tous les jours !

[Newenda]

Hum finalement petit problème. En voulant ouvrir un autre fichier identique :

rawdata0002.raw, il y a un message d'erreur :

??? Attempted to access str(1); index out of bounds because
numel(str)=0.

Comprends pas du tout le pourquoi du comment...

MErci d'avance

PS: en pièce jointe le fichier rawdata0002

[Jerome Briot]

Dans le premier fichier, les lignes "vides" contenaient des tabulations pas dans le second...

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
fid = fopen('rawdata0001.raw','rt');
 
for n=1:31
    fgetl(fid);
end
 
n = 1;
while ~feof(fid)
 
    str = fgetl(fid);
 
    if isempty(str) || ~isstrprop(str(1),'digit')
        break
    end
 
    data(n,:) = sscanf(str,'%*d%f%f%f%f',[1,4]);
    n = n+1;
 
end
 
fclose(fid);

[Newenda]

Ha oui ! super !

Mais le rawdata0002 passe toujours pas, cette fois à cause du dernier str !
Alors que dans rawdata0001 le dernier str était vide, dans rawdata0002 il vaut: str = ''.
Et donc problème pour créer data !

??? Subscripted assignment dimension mismatch.

Error in ==> RheoOSC_dec2010_auto2 at 33
data(n,: ) = sscanf(str,'%f',[1,5]);

Peut on éviter d'enregistrer le dernier str? ou alors créer une autre condition pour que str = '' soit virer?

[Jerome Briot]

J'ai édité mon dernier code avant que tu ne répondes...

[Newenda]

Merci ! maintenant c'est bon, pour tous les fichiers !