Discussion :

[callofduty]

Bonsoir, je dispose d'un fichier texte contenant 3 colonnes (time, altitude, valeurs de Ln(Pr2))mal séparées. Afin d'éviter plusieurs répétitions au niveau de temps et d'altitude, je voudrais l'organiser comme le montre cet exemple:

T      A     Ln(Pr2)                                                   
0.2    2.4   -1.4                                                          
0.2    2.5    2.0                                                           
.                                                       
.                                                                 
0.2    2.9   0.09                                                        
0.3    2.4    -8.04
0.3    2.5    -0.087
...     2.9     3.0

vers

A/T  0.2       0.3
2.4 -1.4       -8.04
2.5  2.0        2.5
.
.
2.9  0.09   3.0

Le résultat doit être affecter dans un fichier texte. Voila, j'espere que j'ai bien illustré mon idée.

[Jerome Briot]

A tester :

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('data.txt','r');
X = textscan(fid,'%f%f%f','headerlines',1,'collectoutput',true);
fclose(fid);
 
uT = unique(X{1}(:,1));
uA = unique(X{1}(:,2));
 
M = reshape(X{1}(:,3), numel(uA),[]);
 
fid = fopen('new_data.txt','wt');
 
fprintf(fid, 'A/T ');
fprintf(fid, '%.2f ', uT);
fprintf(fid, '\n');
 
fmt = [repmat('%.2f ', 1, numel(uT)) '%.2f\n'];
 
fprintf(fid, fmt, [uA(:) M].');
 
fclose(fid);

[callofduty]

Merci Dut, lors de la phase d'exécution j'ai eu ce message d'erreur:

Error using reshape
Product of known dimensions, 999, not divisible into total number of
elements, 1132089.

[Jerome Briot]

Donne nous le fichier texte mis dans une archive zip

[callofduty]

Je pense pas que sera possible. La taille d'archive est de 6.6 Mo a part si je vous envoie une portion de ce fichier.

[Jerome Briot]

Les valeurs dans les deux premières colonnes doivent former une grille rectangulaire et régulière

Il ne doit pas y avoir de valeurs "manquantes"

Par exemple :

T      A                                                  
0.2    2.4                                                          
0.2    2.5                                                         
0.2    2.6
0.2    2.7
0.2    2.8                                                 
0.2    2.9                                                       
0.3    2.4                                                         
0.3    2.5                                                          
0.3    2.6
0.3    2.7 
0.3    2.8                                                 
0.3    2.9 

[Jerome Briot]

Les valeurs ne sont pas exactement les mêmes dans le morceau de fichier que tu nous donnes :

-0.0497341 14.985 -12.068

-0.0329838 14.955 -10.8144

-0.0162163 14.955 -10.9264

0.00055 15 -13.505

Elles ne forment donc pas une grille régulière

[callofduty]

J'ai compris ce que vous voulez me dire, les données sont enregistrées comme le montre le fichier texte le montre.

[Jerome Briot]

Il faut donc interpoler les valeurs sur une grille régulière.

Voici le principe :

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
fid = fopen('data1.txt','rt');
fgetl(fid);
M = fscanf(fid,'%f',[3, inf]).';
fclose(fid);
 
t = linspace(min(M(:,1)), max(M(:,1)), 4);
a = linspace(min(M(:,2)), max(M(:,2)), 100);
 
[x,y] = meshgrid(t,a);
 
Mi = griddata(M(:,1),M(:,2),M(:,3),x,y,'cubic');
 
fid = fopen('new_data.txt','wt');
 
fprintf(fid, 'A/T ');
fprintf(fid, '%.2f ', t);
fprintf(fid, '\n');
 
fmt = [repmat('%.2f ', 1, numel(t)) '%.2f\n'];
 
fprintf(fid, fmt, [a(:) Mi].');
 
fclose(fid);
 
edit new_data.txt
 
figure
plot3(M(:,1),M(:,2),M(:,3),'rd');
hold on
surf(x,y,Mi,'edgecolor','none');

A toi de l'adapter à ton problème.