Discussion :

[Pixcoder]

Bonjour,

J'ai une image volumique en niveau de gris sous forme d'un tableau tridimensionnel d'unsigned char ptr donc chaque élément correspond à un niveau de gris (blanc=niveau maximal, noir=niveau nul).

Je désire trier mes voxels selon leur niveau de gris. J'ai pensé à une std::map dont les clés sont les adresses des voxels dans ptr et dont les valeurs sont les niveaux de gris. Voici un exemple :

Code :

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

1
2
3
4
5
6
7
 
std::map<unsigned char*, unsigned char, classcomp> mamap;
//remplissage de la map
for(int i=0; i<3; i++)
 for(int j=0; j<2; j++)
  for(int k=0; k<4; k++)
    mamap[&ptr[i][j][k]] = ptr[i][j][k];

où classcomp est défini comme suit :

Code :

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

1
2
3
4
5
struct classcomp {
  bool operator()( unsigned char* s1, unsigned char* s2 ) {
        return (*s1)<(*s2);
      }
};

Après ceci, le premier constat est que la taille de mamap vaut 1!!! C'est pas bon! La taille devrait être la taille du volume (=nbre total de voxels).

Où pensez-vous que se situe mon erreur?

Merci de m'aider, je ne comprends pas
Pix

[JulienDuSud]

Insérer les voxels selon leur adresse mémoire est à priori très douteux. Mais bon, si c'est ce que tu veux faire...

Une adresse n'est pas un unsigned char, c'est un unsigned int.

Donc, pour ce faire, tu dois faire comme ça:

Code :

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

1
2
3
4
5
6
std::map<unsigned int, unsigned char, classcomp> mamap;
//remplissage de la map
for(int i=0; i<3; i++)
 for(int j=0; j<2; j++)
  for(int k=0; k<4; k++)
    mamap[(unsigned int)&ptr[i][j][k]] = ptr[i][j][k];

[Pixcoder]

Salut,

Merci pour ta proposition, mais je pense que ca ne peut pas marcher :
étant donné que je veux trier les voxels par niveau de gris croissant rien qu'en les insérant dans la map (qui trie en interne par ordre croissant selon l'opérateur de comparaison, lors de l'insertion), c'est l'opérateur de comparaison qui est important et en mettant "unsigned int" comme type d'adresse, je ne peux plus accéder à la valeur pointée! Et donc plus comparer les niveaux de gris!

Tu comprends ce que j'ai voulu faire? Mais bon, ça ne marche pas...

[JulienDuSud]

Tu peux ré-effectuer un cast avec l'adresse pour avoir la valeur, mais c'est vraiment trop barbare (déjà que je ne mettrais jamais le code que j'ai moi même posté plus haut dans un de mes code sources...)

Tu ferais mieux alors d'avoir un vecteur, simplement. Tu ajoutes tes voxels un par un avec un push_back, et ensuite tu fais un sort avec ton foncteur.

[Invité]

Faudrait pas faire dans l'autre sens? Si tes niveaux de gris sont les valeurs de tes pointeurs, ce doit être eux les clefs de ta map (la map est triée par clefs, pas par valeurs)

Donc, faudrait un truc comme

Code :

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

1
2
 
mamap[ptr[i][j][k]]=&(ptr[i][j][k]);

(et peut être mettre les coordonnées i,j,k au lieu des adresses comme valeurs de la map...)

non?



Francois

[Pixcoder]

J'ai essayé ta proposition :

Code :

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

1
2
3
4
5
6
7
8
9
 
std::vector<unsigned char *> monvector;
//remplissage du vecteur
for(int i=0; i<3; i++)
 for(int j=0; j<2; j++)
  for(int k=0; k<4; k++)
    monvector.push_back(&ptr[i][j][k]);
 
std::sort(monvector.begin(), monvector.end(), uc_cpr);

avec uc_cpr défini comme suit :

Code :

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

1
2
3
4
 
bool uc_cpr( unsigned char * s1, unsigned char * s2 ) {
  return (*s1)<=(*s2);
}

et le programme plante, je ne comprends pas bien..

[Pixcoder]

fcharton : le souci est que pour un niveau de gris donné, plusieurs voxels peuvent l'avoir... Donc en inversant, je perds des données!

[JulienDuSud]

J'ai essayé ta proposition :

Code :

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

1
2
3
4
5
6
7
8
9
 
std::vector<unsigned char *> monvector;
//remplissage du vecteur
for(int i=0; i<3; i++)
 for(int j=0; j<2; j++)
  for(int k=0; k<4; k++)
    monvector.push_back(&ptr[i][j][k]);
 
std::sort(monvector.begin(), monvector.end(), uc_cpr);

avec uc_cpr défini comme suit :

Code :

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

1
2
3
4
 
bool uc_cpr( unsigned char * s1, unsigned char * s2 ) {
  return (*s1)<=(*s2);
}

et le programme plante, je ne comprends pas bien..

Normal, uc_cpr doit renvoyer vrai si le 1er argument est strictement inférieur, et non inférieur ou égal !

Ce code fonctionne chez moi:

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
#include <vector>
#include <iostream>
#include <algorithm>
 
using namespace std;
 
bool uc_cpr( unsigned char * s1, unsigned char * s2 ) {
  return (*s1)<(*s2);
}
 
 
int main()
{
	vector< unsigned char * > v;
 
	unsigned char ptr[10][10][10];
 
	for( int i = 0; i < 10; ++i )
		for( int j = 0; j < 10; ++j )
			for( int k = 0; k < 10; ++k )
				ptr[i][j][k] = i+(i*j+k)+(i+j*k);
 
	for( int i = 0; i < 10; ++i )
		for( int j = 0; j < 10; ++j )
			for( int k = 0; k < 10; ++k )
				v.push_back( &ptr[i][j][k] );
 
	cout << "pre order:" << endl;
	for( unsigned int i = 0; i < v.size(); ++i )
		cout << '0' + *v[i] << endl;
 
	sort( v.begin(), v.end(), uc_cpr );
 
	cout << "post order:" << endl;
	for( unsigned int i = 0; i < v.size(); ++i )
		cout << '0' + *v[i] << endl;
}