Eigen et (très) grosse matrice

Bonjour,

C'est encore moi et j'ai encore un problème avec Eigen :)

Je travail sur une appli qui manipule des images. A partir de mon image je créer une matrice dont les dimensions sont N*N avec N étant le nombre de pixels de mon image. Du coup pour une image 100*200 pixels, ma matrice doit avoir pour dimension 20000*20000. Lors de l'allocation j'ai alors :

Code:

---

Exception non gérée à 0x7c812afb dans LocalAdjustment.exe*: Exception Microsoft C++*: std::bad_alloc à l'emplacement mémoire 0x0012f53c..

---

J'ai testé avec une image plus petite 50*50 pixels et là pas de soucis même si il y a un petit temps d'attente lors de l'allocation.

Bref du coup je ne sais pas comment faire, mon but final étant de résoudre un système linéaire AX = B avec A étant la matrice N*N et X et B des vecteurs colonne de dimension N.

Sinon je gère mes images avec OpenCV, peut être y a-t-il des fonctions intégrées à OpenCV?

Merci pour ceux qui on eu le courage de me lire jusque là!

Citation:

---

Envoyé par jexxo

Bonjour,

C'est encore moi et j'ai encore un problème avec Eigen :)

Je travail sur une appli qui manipule des images. A partir de mon image je créer une matrice dont les dimensions sont N*N avec N étant le nombre de pixels de mon image. Du coup pour une image 100*200 pixels, ma matrice doit avoir pour dimension 20000*20000. Lors de l'allocation j'ai alors :

Code:

---

Exception non gérée à 0x7c812afb dans LocalAdjustment.exe*: Exception Microsoft C++*: std::bad_alloc à l'emplacement mémoire 0x0012f53c..

---

---

20000*20000 = 400 000 000 d'elements. Pas etonnant. Si ta matrice est creuse, essaie de trouver une interface qui te permet de n'en donner que les elements non nuls, sinon, passe en 64 bits.

Pour information (matrices creuses d'Eigen) :
http://eigen.tuxfamily.org/dox/TutorialSparse.html

Sinon, pour compléter ce que dit Jean-Marc, essaie de voir si ton gros système peut pas se ramener à la résolution de plusieurs "petits" systèmes linéaires indépendants, peut-être ?

Ok merci je vais essayer les matrices creuses étant donné que sur chaque ligne au maximum 5 valeurs sont non nulles (le pixel considéré et ses 4 voisins).

A vrai dire, c'est pour toi la seule alternative jouable... Bonne nouvelle, ça va aller a priori plus vite :D

Ah bah voilà tout est bien qui finit bien 8-)