Discussion :

[mapotam]

Bonjour,

j ai un fichier .geom (http://people.sc.fsu.edu/~burkardt/data/off/off.html) contenant les vertices et faces d un objet 3D.
Il a la forme suivante:

OFF
8815 17579 0
16.7722 -200.773 -791.821
16.9744 -200.869 -791.787
17.2209 -200.972 -791.768
17.4559 -201.05 -791.785
...
%la ligne suivante est la 8818eme (8815+2+1)
3 0 47 1
3 1 47 2
3 2 47 48
3 2 48 3
3 3 48 49
3 3 49 4
3 4 49 50
3 4 50 5
3 5 50 51
3 5 51 6
3 6 51 52
3 6 52 7
3 7 52 53
3 7 53 8
3 8 53 54
....
3 8 53 54
% La derniere ligne est la 26396eme.

Je souhaite recuperer les coordonnees des vertices et des faces.

pour cela, j ai fait la fonction suivante:

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
40
41
42
43
44
45
46
47
48
49
50
51
 
function [vertex,faces] = ReadGeomFile(file);
 
%%
% Purpose:  given a file.geom, extract a vector containing the vertex, and
%           a vector containing the faces
%
%   file        file.geom
%   vertex      vector of vertex
%                   x y z 
%                   for vertex 0, 1, ..., vertex_count-1
%   faces       vector of faces
%                   (n) v1 v2 ... vn, 
%                   (the number of vertices) and the vertex indices for each face.
%
% HK 2010/01/04
 
%% Processing
 
% Progress bar
h_waitbar = waitbar(0,'Please wait...'); set(h_waitbar,'windowstyle','modal');
 
f = fopen(file);
 
% Initialization: get the number of vertex, faces and edges
[nb_vertex,nb_faces,nb_edges] = textread(file, ' %f %f %f ', 1, 'headerlines', 1);
 
% Get the vertex vector
 vertex=zeros(nb_vertex,3);
 for i=1:nb_vertex
     [x,y,z] = textread(file, '%f %f %f ', 1, 'headerlines', i+1);
     vertex(i,1)=x;
     vertex(i,2)=y;
     vertex(i,3)=z;
 end
 
% Get the faces vector
faces=zeros(nb_faces,3);
for i=1:nb_faces
    %Update waitbar
    for j=1:300,
        waitbar((i+j/300)/nb_faces,h_waitbar)
    end
    set(h_waitbar,'windowstyle','modal');
    [n, v1,v2,v3] = textread(file, '%f %f %f %f ', 1, 'headerlines', (i+1)+nb_vertex);
 
    % Assign v+1 because first vertex vector indice is 1 instead of 0
    faces(i,1)=v1+1;
    faces(i,2)=v2+1;
    faces(i,3)=v3+1;
end

Ce code marche tres bien..... pour de petit fichier !
Mais pour un fichier de cette dimension, cela prend un temps fou.

Comme ce code doit être intégré dans une fonction dans une GUI, il faut que ce soit rapide.

Quelqu'un aurait-il une idée d'optimisation?

Merci d'avance!

[Jerome Briot]

Tu devrais pouvoir simplement faire ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
fid = fopen('data.geom','rt');
 
fgetl(fid);
 
N = fscanf(fid,'%d%*d%*d');
 
vertices = fscanf(fid,'%f',[3,N]).';
 
faces = fscanf(fid,'%d',[4,inf]).';
 
fclose(fid);

Enfin c'est l'idée...

Par contre le vrai format OFF est un peu plus compliqué (gestion de la couleur et de la transparence me semble-t-il).

[mapotam]

Merci beaucoup Dut!
Ca marche tres bien, et instantanement
J utilise toujours 'textread', mais c est vrai que je devrais penser a 'fscanf' plus souvent.

Et en effet, le .format OFF peut etre plus complexe, celui que j ai est basic.

[mapotam]

J ai encore un petit probleme.
Il y a une erreur dans ma matrcie 'vertices'.
Le premier element de la premiere colonne est toujours faux.
Par exemple, si dans mon fichier, j ai

23.0443 -240.517 -820.986

le code me renvoie

0.0443 -240.517 -820.986

J ai essaye avec plusieurs fichiers differents. J ai bien compris chaque ligne du code, mais je ne vois pas d ou provient l erreur.

[mapotam]

Ok, alors avec

Code :

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

1
2
% Get the number of vertices, faces and edges
N = fscanf(fid,'%d%d%d',[1,3])

au lieu de

Code :

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

1
2
% Get the number of vertices, faces and edges
N = fscanf(fid,'%d%d%d')

il n y a plus l erreur.