problème de comparaison de vectors

Bonjour à tous. je vous expose mon problème qui me tracasse en ce moment:

En fait j'ai à la base 3 csv qui stockent des coordonnées de points en x et y et leur niveau de gris correspondant. en gros ca donne ca

cordx cordy niveau gris
50 150 3810

Je lis mes 3 csv et stocke les valeurs dans un vector par csv (donc un vector 3d qui stocke en colonne 0 la valeur en x, en colonne 1 la valeur en y et en colonne 2 la valeur de niveau de gris.

Ce que je voudrais arriver à faire, c'est trouver les points communs aux 3 csv et les stocker dans un vector résultat avec les 3 niveaux de gris correspondant. En fait j'ai 3 coupes d'images et je dois récupérer les points communs au 3 pour avoir les niveaux de gris dans les 3 dimensions.

Je rajoute à mon explication que les vectors n'ont pas la même taille (ils peuvent l'avoir mais ca veut dire que l'utilisateur qui a tracé la coupe sur l'image a fait 3 fois le même tracé donc il est très fort :mrgreen:). Voilà si quelqu'un a une idée je suis preneur car là je sèche totalement.

Merci d'avance à tous ceux qui pourront m'aider ;)

Bonjour,

La première chose à faire, c'est de rendre les images comparables, c'est à dire de les mettres toutes à la même taille.

Une fois que toutes les images sont à la même taille, il devient facile de les parcourir de manière uniforme pour les comparer.

Code:

---

1 2 3 4 5 6 7 8 9

for(int i; i < width; ++i){ for(int j = 0; j < height; ++j){ // cas vector 2d if((vvImage1[x][y] == vvImage2[x][y]) && (vvImage1[x][y] == vvImage3[x][y])) { // Stoker la valeur dans la nouvelle image } } }

---

Citation:

---

En fait j'ai 3 coupes d'images et je dois récupérer les points communs au 3 pour avoir les niveaux de gris dans les 3 dimensions.

---

Un tour sur le forum traitement d'images serai à mon avis la première chose à faire avant de passer à l'implémentation.

Je ne sais pas si ça peut t'aider, mais en standard tu as des fonctions comme std::set_intersection (aussi std::set_union, std::set_difference ) qui te permettent de récupérer les points communs (ou les différences ) entre deux containers triés ( à ne pas oublier ) mais de tailles différentes.

Tu peux faire quelque chose comme ça:

Code:

---

1 2 3 4 5 6

std::vector< TaClasse > v1, v2, results;   std::set_difference( v1.begin(), v1.end(), v2.begin(), v2.end(), std::back_inserter( results ) );

---

Il faut que ta classe puisse être comparable ( avec un operateur < ou autre fonction à donner comme dernier paramètre à set_difference )

Merci de ta réponse. En fait ce matin j'ai eu une idée, au lieu d'avoir un vector < vector <int> > j'ai stocké mes points sous forme de CPoint (classe de la mfc qui en gros récupère la valeur en x et y et en fait un point en mémoire).

Du coup je me retrouve avec cette structure:

vector

CPoint(x;y)
CPoint(x;y)
CPoint(x;y)
CPoint(x;y)
...
...
...

Je pense que la comparaison sera plus simple là? A partir de ces poins je peux faire un multimap et comparer les clés? Le problème c'est que les points identiques sont pas forcément au même endroit dans le tableau d'où la difficulté de la comparaison en fait.

Citation:

---

Je pense que la comparaison sera plus simple là? A partir de ces poins je peux faire un multimap et comparer les clés? Le problème c'est que les points identiques sont pas forcément au même endroit dans le tableau d'où la difficulté de la comparaison en fait.

---

tu peux toujours utiliser set_differencet, set_ etc... pour avoir tes résultats.

ok donc la fonction intersection va m'être utile, maintenant le soucis est la surcharge de < qui m'est demandée. Comment surcharger une classe des MFC? je peux directement le faire comme ca où je veux? Je suis pas sûr :mouarf: quelles merdes ces mfc... :mrgreen:

Citation:

---

Envoyé par joesatriani

ok donc la fonction intersection va m'être utile, maintenant le soucis est la surcharge de < qui m'est demandée. Comment surcharger une classe des MFC? je peux directement le faire comme ca où je veux? Je suis pas sûr :mouarf: quelles merdes ces mfc... :mrgreen:

---

Peut-être qu'elle existe déjà?

Si tu peux le mettre dans un container trié, il y a des chances qu'elle existe non?

Sinon tu peux créer ta fonction de comparaison où tu veux oui ( et pas forcément un operator< en utilisant les versions des fonctions/containers qui s'appellent avec une fonction de comparaison en plus )

En fait mes points sont triés car les csv sont triés et vu que je remplis mes vectors à partir de ces csv j'ai donc des vectors triés. Le soucis c'est qu'à la compilation il me dit:

c:\Program Files\Microsoft Visual Studio .NET 2003\Vc7\include\algorithm(2384): error C2678: '<' binaire : aucun opérateur trouvé qui accepte un opérande de partie gauche de type 'std::allocator<_Ty>::value_type' (ou il n'existe pas de conversion acceptable)
with
[
_Ty=CPoint
]
donc j'en déduis qu'il veut la surcharge de '<' de la classe CPoint nan?

Citation:

---

Envoyé par joesatriani

En fait mes points sont triés car les csv sont triés et vu que je remplis mes vectors à partir de ces csv j'ai donc des vectors triés. Le soucis c'est qu'à la compilation il me dit:

c:\Program Files\Microsoft Visual Studio .NET 2003\Vc7\include\algorithm(2384): error C2678: '<' binaire : aucun opérateur trouvé qui accepte un opérande de partie gauche de type 'std::allocator<_Ty>::value_type' (ou il n'existe pas de conversion acceptable)
with
[
_Ty=CPoint
]
donc j'en déduis qu'il veut la surcharge de '<' de la classe CPoint nan?

---

oui tu peux créer un:

Code:

---

1 2 3 4

bool operator<( const CPoint & p1, const CPoint & p2 ) { }

---

à un endroit approprié.

ouaip j'y avais pensé par contre l'autre soucis c'est que dans mon appel de fonction:

set_intersection(X.begin(),X.end(), Y.begin(), Y.end(), results);

je n'ai pas vraiment de points à donner en paramètre à donner à l'opérateur :roll:

bon la surcharge a marché mais j'ai encore une erreur ...lol:

c:\Program Files\Microsoft Visual Studio .NET 2003\Vc7\include\algorithm(2389): error C2676: '++' : 'std::vector<_Ty>' binaire ne définit pas cet opérateur ou une conversion vers un type acceptable pour l'opérateur prédéfini
with
[
_Ty=CPoint
]

j'avoue que celle-ci me laisse pantois 8O j'avoue ne pas comprendre pourquoi il veut un opérateur ++ sur mes CPoint... lol


EDIT: oups je crois que c'est parce que j'ai oublié ca std::back_inserter

je teste ca en revenant au bureau staprem et si ca plante encore je reviendrai par ici ;D en tout cas merci beaucoup pour votre aide ;)

Bonjour,

si tu as des soucis avec les set_xxx, tu peux jeter un coup d'oeil ici ;)

yes merci bien à tous vous m'enlevez une sacrée épine du pied ca va me permettre de finir mon stage tranquilou en rendant quelque chose qui marche :mrgreen: je teste tout ca à 14h et je vous tiens au jus :yaisse2:

A priori ca marche donc je finis mes tests et si j'ai un problème je sais où venir. merci beaucoup ;)

Effectivement ca a bien marché merci beaucoup. J'ai un bug marant à vous faire partager. J'ai stocké dans un vector <vector <int> > des points avec leur niveau de gris

X Y nvg1 nvg2 nvg3
.
.
.
.
.
.
nbpoints

donc dans mon tableau tout va bien si je demande d'afficher dans une messagebox la valeur de mes niveaux de gris dans une boucle ca s'affiche nickel. Par contre quand je demande d'écrire cette même valeur dans un fichier, j'ai toujours la meme valeur de niveau de gris qui s'écrit (la première en l'occurrence). J'avouerai que là, je ne comprends pas d'où vient le problème.

Code:

---

1 2 3 4 5 6 7 8 9

for (int i = 0; i< ecx.size(); i++) { nvg = ecx[i][2];   // si je demande d'afficher les valeurs, j'ai bien les bonnes valeurs ici.   ficres<<nvg;// alors qu'ici dans le fichier j'ai toujours la même qui s'écrit.   }

---

si vous avez une idée je suis preneur car là j'avoue que je comprends pas du tout ce qui se passe.

EDIT: j'ai encore plus drôle:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

for (int i = 0; i< ecx.size(); i++) { int test = ecx[i][2]; int test1 = ecx[i][3]; int test2 = ecx[i][4]; if (i == ecx.size()-2) { char c[10]; itoa(test,c,10); CString MFCString; MFCString = c; AfxMessageBox(MFCString); itoa(test2,c,10); MFCString = c; AfxMessageBox(MFCString); itoa(test3,c,10); MFCString = c; AfxMessageBox(MFCString); }   ficres<<ecx[i][2]<<";"<<ecx[i][3]<<";"<<ecx[i][4]<<";"<<"\n";   }

---

et bein devinez quoi? avec ce code là, ca m'écrit les bonnes valeurs dans le fichier.... je pige rien du tout sur ce coup là lol :yaisse2:

bon finalement je vais utiliser une autre méthode. Je vous explique le topo:

J'ai 3 map <CPoint,int>. Dans celles ci, j'ai mes points en clés et à chaque point j'attribue un niveau de gris. Jusque là ca va. Ce que je voudrais faire, c'est avoir les points en commun des 3 map. Avec set_intersection, a priori c'est impossible donc si quelqu'un avait une idée, je suis preneur ;) parce que là je bloque dessus depuis 2 jours et je dors pas beaucoup à me creuser la tête sur ce bug :(

un bug? non. Impossible? c'est toi qui le dis :)

Citation:

---

Envoyé par joesatriani

Ce que je voudrais faire, c'est avoir les points en commun des 3 map.

---

set_union() ?

le problème si je fais set_union ca va me retourner l'union des map et non l'intersection ?

En fait j'ai trouvé une idée qui me parait bien

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

for (it = Xm.begin(); it!=Xm.end(); it++) { pair<CPoint,int> t; if (Ym.find((*it).first) != NULL) { if (Zm.find((*it).first) != NULL) { t.first = (*it).first; t.second = 0; res.insert(t); } } }

---

malheureusement ca me donne rien de bien bon , en gros ca me retourne ce que j'ai dans Xm ce qui n'est pas logique du tout aux vues du code.

oui oops pardon, pas réveillé :oops:

Pourquoi set_intersection ne te conviens pas?

bein en fait avec des maps je vois pas comment utiliser set_intersection avec des maps, avec les vectors j'y arrive mais pas avec les map :s. Je pige pas pourquoi mon code me retourne pas ce que je veux d'ailleurs pour le coup....

Sinon, ton code est erroné, car:

Code:

---

(Ym.find((*it).first) != NULL)

---

sera toujours vrai.
idem pour:

Code:

---

(Zm.find((*it).first) != NULL)

---

les fonctions find() de la STL renvoient toujours un itérateur, donc elles ne renvoient jamais NULL. Si l'élément recherché n'est pas trouvé, cette fonction retourne l'itérateur spécial end().

Il faut donc remplacer tes tests par:

Code:

---

if ( Ym.find((*it).first) != Ym.end() )

---

et

Code:

---

if ( Zm.find((*it).first) != Zm.end() )

---

Alors là je dois dire merci, nan nan pire que ca MERCIIIIIIIIIIIIIII. Comment débloquer une situation qui paraissait désespérée depuis lundi. Plein de bisous :king:

Problème résolu à 100 %, merci beaucoup ca va me faire rattraper mon retard.

bon ok pourtant... enfin bref

Ce n'est effectivement pas évident d'utiliser les algorithmes set_xxx avec des maps. Je n'y suis pas encore arrivé en tout cas.

En effet, le back_inserter utilise la fonction membre push_back du conteneur de destination, et map n'a pas de push_back. J'ai essayé avec un simple insert_iterator, mais je ne comprend pas trop ce qu'il se passe, ça ne fonctionne pas. Dès que j'aurais un peu de temps, j'y regarderai de plus près.

En attendant, pour ceux que ça intéresse, voici le code que j'utilise pour essayer de trouver une solution (code qui est censé solutionner le problème de ce topic):

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 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 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71

#include <iostream> #include <map> #include <algorithm> #include <iterator>   using namespace std;   struct MyPoint { MyPoint( unsigned int px = 0, unsigned int py = 0 ) : x(px), y(py) {} unsigned int x; unsigned int y; };   ostream& operator << ( ostream & os, const MyPoint & pt ){ os << "x: " << pt.x << " , y: " << pt.y; return os;}   bool operator < ( const MyPoint & left, const MyPoint & right ) { if ( left.x < right.x ) return true; if ( left.x == right.x ) return ( left.y < right.y ); return false; }   bool operator == ( const MyPoint & right, const MyPoint & left ) { return ( right.x == left.x && right.y == left.y ); }   typedef pair<MyPoint, unsigned int> PPair;   bool operator < ( const PPair & left, const PPair & right ) { return ( left.first < right.first ); }   bool operator == ( const PPair & left, const PPair & right ) { return ( left.first == right.first ); }   ostream& operator << ( ostream & os, const PPair & ppair ){ os << ppair.first << " , value: " << ppair.second << endl; return os;}   typedef map<MyPoint, unsigned int> PMap; typedef PMap::iterator PMapIt;   void PrintPMap( const PMap & pm, ostream & os = cout ) { copy( pm.begin(), pm.end(), ostream_iterator<PPair>(os)); }   int main(int argc, char** argv){ PMap m1, m2, m3; m1.insert( PPair( MyPoint( 10, 20 ), 128 ) ); m1.insert( PPair( MyPoint( 15, 10 ), 212 ) ); m1.insert( PPair( MyPoint( 5, 30 ), 124 ) ); m1.insert( PPair( MyPoint( 10, 40 ), 65 ) ); cout << "map1:" << endl; PrintPMap( m1 );   m2.insert( PPair( MyPoint( 15, 20 ), 89 ) ); m2.insert( PPair( MyPoint( 15, 10 ), 54 ) ); m2.insert( PPair( MyPoint( 5, 30 ), 212 ) ); m2.insert( PPair( MyPoint( 12, 30 ), 165 ) ); cout << endl << "map2:" << endl; PrintPMap( m2 );   set_intersection( m1.begin(), m1.end(), m2.begin(), m2.end(), std::insert_iterator<PMap>( m3, m3.begin() ) ); //ici ça ne fonctionne pas, aucun élément n'est inséré dans m3 cout << endl << "map3:" << endl; PrintPMap( m3 );   cout << "end" << endl; cin.get(); return 0; }

---

Evidemment, dans ce cas précis, personnellement, j'aurais intégré le valeur du niveau de gris dans la classe MyPoint, mais j'aimerais bien parvenir à faire fonctionner ces fonctions de set_xxx avec des maps :aie:

edit: une solution pourrait consister à créer notre propre inserter. Mais j'aimerais trouver plus simple.

insert_iterator est sensé inséré à un endroit spécifié, notion qui n'est pas supportée par les maps...

arf...

vois-tu une solution alors?

Je ne comprends pas très bien.

un std::inserter utilise la méthode insert qui pour un set ou un map peut prendre un deuxième paramètre qui est une position dans la map.

Ce paramètre n'est qu'une indication, il ne précise pas que l'on va insérer à un endroit précis, c'est un "hint".

Ces algorithmes sont fait pour fonctioner avec des sets, donc il n'y a pas de raison que ça ne marche pas.

Perso j'aurais ptet fait une struct/classe avec un CPoint et un niveau de gris, une manière de trier tout ça et envoyer tout ça dans un set.

Hmm... je n'y comprend plus rien :aie:

Je n'ai pas la norme sous les yeux, mais dans l'implémentation de <iterator> que j'ai sous la main (gcc 4.1.2), l'operateur d'affectation du insert_iterator est:

Code:

---

1 2 3 4 5 6

insert_iterator& operator=(const typename _Container::const_reference __value) { iter = container->insert(iter, __value); ++iter; return *this; }

---

Le premier paramètre (iter) est un iterateur, et __value la valeur à insérer.

D'apres cgi.com, la map implémente une fonction membre

Code:

---

iterator insert(iterator pos, const value_type& x)

---

Donc ça devrait marcher non? :koi:
Mais alors pourquoi le code que j'ai posté ci-dessus ne fonctionne pas?


edit: nikko34, je n'avais pas vu ton post :oops:, mais en fait nous disons la même chose.

edit2:

Citation:

---

Envoyé par nikko34

Perso j'aurais ptet fait une struct/classe avec un CPoint et un niveau de gris, une manière de trier tout ça et envoyer tout ça dans un set.

---

Je suis d'accord, c'est en substance ce que j'avais dit quelques posts auparavant. Mais on ne connait pas suffisamment le contexte pour être sûr.

Salut.
Y as un truc que je ne pige pas.
Tu as trois images et tu fais du traitement dessus. C'est bien Çà?

Pourquoi ne pas traiter ton problème par des images???

Tu parcoure tes trois images en même temps et tu fait ce que tu as à faire.

Ce que j'ai du mal à comprendre est que parle de point 2D et ensuite de 3D :koi: