tableau de valeurs sans doublons à arrondis

Version imprimable

Bonjour à tous,

Je suis en train d'écrire un code lisant un fichier stl (fichier CAO maillé surfacique) et interprétant son contenu.
Le contexte : je veux créer un tableau de points (une liste de coordonnées xyz de points = tableau "n x 3" avec n= nbre de points) mais certains points apparaissent plusieurs fois dans le tableau. J'enlève donc les doublons dans la liste en utilisant le code ci après :
Code:

1 2 3 4 5 6 tableau_des_points = [[-300.0, 0.0, 700.0], [-300.0, 100.0, 700.0]], [[-300.0, 0.0, 700.0], [-325.88189999999997, 0.0, 796.59259999999995]], [[-428.7901, 0.0, 665.49080000000004], [-454.67200000000003, 0.0, 762.08339999999998]], [[-300.0, 100.0, 700.0], [-300.0, 1.421085e-014, 700.0]], [[-300.0, 66.666659999999993, 700.0], [-325.88189999999997, 66.666659999999993, 796.59259999999995]], [[-325.88189999999997, 0.0, 796.59259999999995], [-351.7638, 0.0, 893.18520000000001]], [[-325.88189999999997, 0.0, 796.59259999999995], [-300.0, 1.421085e-014, 700.0]], [[-351.5573, 0.0, 720.69449999999995], [-377.43920000000003, 0.0, 817.28700000000003]], [[-377.64569999999998, 0.0, 989.77779999999996], [-351.7638, 0.0, 893.18520000000001]], [[-308.62729999999999, 22.22222, 732.19749999999999], [-334.50920000000002, 22.22222, 828.79010000000005]] tableau_des_points_sans_doublon = [] for k in tableau_des_points: if k not in tableau_des_points_sans_doublon: tableau_des_points_sans_doublon.append(k)
Le problème est que certaines valeurs des coordonnées sont arrondies. Exemple : un même point peut être sous la forme (100.0, 0.0, 10.0) et sous la forme (99.9999999, 1.4e-12, 10.000001) car la conversion STL (réalisée par CATIA V5) gère mal les arrondis.
Quelqu'un pourrait-il m'indiquer une astuce pour que mon code reconnaisse ces deux écritures différentes du même point comme identiques et ne les ajoute donc pas au tableau sans doublons. L'idéal serait qu'une variable me permette de définir à combien de chiffres après la virgule on peut considérer que les points sont identiques.
Merci d'avance,

Fred

Une méthode possible

Voici une méthode à tester et à adapter à ton cas :

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
#!/usr/bin/env python
#coding=utf-8
 
tableauDeNbres = [300, 300.000, (300+10**(-14)),1,2,3,0,10**(-14)]
tableauNettoye = []
 
for unNbre in tableauDeNbres:
    if not unNbre in tableauNettoye:
        tableauNettoye.append(unNbre)
 
print '1ere Methode :'
print '\t',
print tableauNettoye
 
# Ici on demande à Python d'arrondir les nombres
# en gardant 3 chiffres après la virgule.
tableauDeNbres = [float("%.3f" %float(x)) for x in tableauDeNbres]
tableauNettoye = []
 
for unNbre in tableauDeNbres:
    if not unNbre in tableauNettoye:
        tableauNettoye.append(unNbre)
 
print '2eme Methode :'
print '\t',
print tableauNettoye

On obtient :

Code:

1
2
3
4
1ere Methode :
	[300, 1, 2, 3, 0, 1e-14]
2eme Methode :
	[300.0, 1.0, 2.0, 3.0, 0.0]

Merci Rambc,

En retravaillant un peu, j'ai trouvé une solution que je pense assez semblable à la tienne (mais un peu moins propre je pense). La voici :
Code:

1 2 3 4 5 for k in tableau_des_points: if [eval("%.5f" %k[0]), eval("%.5f" %k[1]), eval("%.5f" %k[2])] not in tableau_des_points_sans_doublon: tableau_des_points_sans_doublon.append([eval("%.5f" %k[0]), eval("%.5f" %k[1]), eval("%.5f" %k[2])]) print 'tableau_des_points_sans_doublon', tableau_des_points_sans_doublon
Peux-tu me donner ton avis sur cette méthode stp? De toutes façons, je vais prendre la tienne...
D'autre part, quelle est la différence entre

Code:

if not k in tableau

et

Code:

if k not in tableau

? Y a-t-il une différence?

Merci d'avance,

Fred

Aucune différence dans la sémantique.
Les deux testent que l'objet "k" ne soit pas dans le tableau.
D'ailleurs, les deux deviennent le même bytecode:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
>>> def test1():
...     return k not in tableau
...
>>> def test2():
...     return not k in tableau
...
>>> from dis import dis
>>> dis(test1)
  2           0 LOAD_GLOBAL              0 (k)
              3 LOAD_GLOBAL              1 (tableau)
              6 COMPARE_OP               7 (not in)
              9 RETURN_VALUE
>>> dis(test2)
  2           0 LOAD_GLOBAL              0 (k)
              3 LOAD_GLOBAL              1 (tableau)
              6 COMPARE_OP               7 (not in)
              9 RETURN_VALUE

14/10/2009, 17h17
fraid49

Ok, c'est noté, merci Antoine!
14/10/2009, 17h43
rambc

Citation:

Envoyé par fraid49

En retravaillant un peu, j'ai trouvé une solution que je pense assez semblable à la tienne (mais un peu moins propre je pense).

Peux-tu me donner ton avis sur cette méthode stp?

float("%.3f" %(1**(-14))) et eval("%.3f" %(1**(-14))) renvoient la même chose. Ceci étant dit, la 1ère instruction utilise juste une conversion, tandis que la seconde lance une évaluation de calcul, ce qui demande plus de travail en général. A choisir, je prendrais moi aussi la 1ère méthode.
15/10/2009, 09h19
fraid49

C'est noté pour le temps de travail (de toutes façons, ta méthode est beaucoup plus propre que la mienne).
Merci,

Fred