[voronoin] Diagramme de Voronoi dans un cube

Bonjour tous,

j'aimerai utiliser la fonction "voronoin" de matlab :
http://www.mathworks.fr/fr/help/matl.../voronoin.html

pour faire le maillage d'un cube comme vous pouvez voir en pièce jointe.

Le soucis est que j'ai lu la doc et je n'ai rien compris, je trouve que les exemples sont pourri...

Pourriez vous s'il vous plait me montrer comment faire ce genre
de choses avec matlab ?

je vous remercie pour votre aide

Je viens de trouver deux discussions très intéressantes :
http://www.developpez.net/forums/d44...noi-3d-resolu/
et
http://www.developpez.net/forums/d12...hlight=voronoi

Du coup, voici ce que j'ai fais:
1) j'ai donné les coordonnées des sommets de mon cube dans X
2) j'ai appelé la fonction de voronoi pour qu'elle me fasse un maillage dans ce domaine (avec l'option Qbb car j'ai une dimension=3)
3) j'ai mis l'option 'QJ' mais je ne sais pas vraiment pourquoi... le message d'erreur optenu me semblait moins grave que si je mets rien... :aie:
4) ensuite à partir des données renvoyées par 'voronoin' j'essai de faire une enveloppe convexe grace à 'convhulln' et tracer cette enveloppe grace a patch.

Code:

---

1 2 3 4 5 6

clear all;close all;clc; X=[0 0 0;1 0 0;1 1 0;0 1 0;0 0 1;1 0 1;1 1 1;0 1 1]; [V,C] = voronoin(X,{'Qbb','QJ'}); K = convhulln(V,{'Qt','Pp'});   patch('faces', K, 'Vertices',V)

---

le soucis est que je ne comprends pas vraiment la logique des deux fonctions "voronoin" et "convhulln" du coup je ne sais pas trop quoi leur donner en entrée et avec quelle options....

j'ai l'impression qu'il serait mieux correct de donner "C" en 'input' de la fonction convhulln mais d'après la doc il n'est pas spécifié que cette fonction prend des cell en entrée :cry:

pourriez vous me débloquer pour ce petit bout de code svp ?

merci

Je ne suis pas un spécialiste mais en suivant le lien donné dans le deuxième message de cette discussion : voronoi cells in bounded surface or volume

Et en faisant simplement :

Code:

---

1 2 3 4

seeds = 2*(rand(100,3)-0.5); Options.pbound = unitbox(3, 1); V = mpt_voronoi(seeds, Options); plot(V)

---

J'obtiens un résultat assez semblable.

Il faudra sans doute lire la documentation de la toolbox pour améliorer le résultat.

merci beaucoup DUT!!!!!!!!!!!!
c'est exactement ce que je cherche à faire :ccool:

par contre ton code utilise une toolbox special "MPT" et je n'arrive pas à l'installer....

d'après ce que j'ai compris il faut aller sur le site:
http://control.ee.ethz.ch/~mpt

ensuite prendre le dossier "MPT", le mettre sur un endroit de notre disque et suivre les instructions de cette pages :
http://control.ee.ethz.ch/~mpt/docs/install.php

le soucis est que même en procédant comme ceci matlab ne reconnait pas les fonctions de cette toolbox...

pourrais tu m'en dire plus sur l'installation de MPT, comment as tu fais de ton côté ?

merci :ccool::ccool::ccool:

Bonjour,

Je n'ai eu aucun problème en suivant ces mêmes directives... :?
Qu'as-tu fait précisément ?

Voir aussi la FAQ Comment utiliser une Toolbox ?

je viens de refaire les directives et à présent ça fonctionne.... :aie:
je pense que ça devait être un problème avec 'savepath' je n'avais pas cette commande auparavant, où peut être que mon 'Path' etait pas top (espace dans nom fichier...)

bref, maintenant c'est réglé et ça fonctionne, merci beaucoup tous les 2 :ccool: