construire un voxel

**Akramou** · 12/02/2012, 22h42

Salut,
je travaille sur la reconstruction des images 2D( des projections) en 3D, dans un premier lieu de mon travail je doit avoir une images volumique initiale avec des valeurs quelconques.
ma question est comment déclarer un voxel en c++,(càd construire un voxel et après le remplir, quelque chose comme "float volume[depth][rows][cols] )
merci

**pyros** · 13/02/2012, 09h30

Salut,

Je ne comprends pas trop ta question. Un voxel n'est rien d'autre qu'une valeur (float, int, double, ... selon ce que tu manipule), au même titre qu'un pixel d'une image (je considère que tes images sont en niveau de gris, ce qui est généralement le cas pour de la reconstruction 3D).
Pour moi, tu aurais quelque chose dans le genre :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
class Volume
{
  Volume(int width, int height, int depth)
  : m_values(width * height * depth)
  {
  }
 
  private:
    std::vector<float> m_values; // float, int, ... selon le type que tu utilise
}

**DonQuiche** · 13/02/2012, 13h56

Une telle approche serait extrêmement gourmande en mémoire et CPU, on peut mieux faire. En effet, à moins d'avoir besoin de visualiser des plans de coupe, on n'a besoin, pour une ligne donnée, que des données de deux voxels : celui d'entrée et celui de sortie, avec leurs coordonnées. Autrement dit, un tableau bidimensionnel de structures type { int color1, int z1, int color2, int z2 }.

Problème : certains objets complexes peuvent avoir plusieurs segments par ligne (cas d'une cruche, on peut avoir une ligne passant par le corps, puis ressortant dans le vide et repénétrant plus loin dans l'anse). Soit on garde l'approche précédente en subdivisant les objets complexes en objets simples (la cruche est divisée en deux morceaux), soit on stocke plusieurs structures comme celle donnée plus haut par ligne.

Enfin, même si l'on voulait des plans de coupe, on peut au moins faire des économies en virant les voxels vides. Par exemple en décrivant une ligne via une structure du type : { int z1, int z2, int[] voxels }, ou via une séquence de telles structures.

EDIT : J'ai oublié le pb des surfaces planes. Grosso modo on reprend la dernière solution en virant les voxels invisibles.