Discussion :

[piter_pam]

Bonjour a tous

Quand un indice est un vecteur colonne, on peut directement récupérer d'un tableau les valeurs qui correspondent a cet indice.

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
>> tableau = rand(8,1)
 
tableau =
 
    0.9501
    0.2311
    0.6068
    0.4860
    0.8913
    0.7621
    0.4565
    0.0185
 
>> indice= [5;8]; tableau(indice, 1)
 
ans =
 
    0.8913
    0.0185

Si indice est un vecteur ligne, on doit transposer le résultat pour récupérer les valeurs sous forme de vecteur ligne.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
>> indice= [5 8];
tableau(indice, 1)
 
ans =
 
    0.8913
    0.0185
>> tableau(indice, 1)'
 
ans =
 
    0.8913    0.0185

Mais si l'indice est de taille n-dimensionnel (typiquement bidim), comment récupérer les valeurs sous la même forme que l'indice ?


Voici mon problème
J'ai une variable bkg: un tableau de taille (60,1,3). Il s'agit d'une sorte de colormap, la 3eme dimension code les valeurs R G et B.

J'ai un tableau d'indice A, de taille (300, 800), qui contient "numero" de la couleur a utiliser pour chaque pixel

Et un tableau im=zeros(600,800,3), a afficher avec la commande image. (la partie du bas sera noire)

Code :

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

1
2
3
4
5
>> size(bkg(round(max(min(A,length(bkg)),1)),1,:))
 
ans =
 
      240000           1           3

donc je ne peux pas affecter directement

Code :

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

>>im(1:300,:,:) = bkg(round(max(min(A,length(bkg)),1)),1,:)

Bien sûr, je peux le faire avec des boucles for, mais on doit pouvoir faire plus propre.

Zavez vous une idée ? Une commande qui permet de faire ça ?

M'ci

[Jerome Briot]

Je ne vois pas où si situe le problème avec l'orientation des indices dans ton premier message :

Code :

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

1
2
3
4
5
6
7
8
9
tableau = rand(8,1);
 
indice= [5;8];
 
tableau(indice)
 
indice= [5 8];
 
tableau(indice)

donne le même résultat :

Code :

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

1
2
3
4
5
6
7
8
9
10
ans =
 
    0.5309
    0.8200
 
 
ans =
 
    0.5309
    0.8200

Ensuite, pour ton problème, voici une solution :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
r = size(A,1);
c = size(A,2);
n = numel(A);
 
img = zeros(r,c,3);
 
img(1:n) = bkg(A,1,1);
img(n+1:2*n) = bkg(A,1,2);
img(2*n+1:end) = bkg(A,1,3);
 
img(600,800,3) = 0;

ce qui se résume en une seule ligne à l’aide de la fonction ind2rgb :

Code :

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

1
2
3
img = ind2rgb(A,bkg);
 
img(600,800,3) = 0;

[piter_pam]

bonjour Dut

Je ne connaissais pas cette façon d'indexer les matrices. Ça se rapproche du langage de bas niveau, avec un pointeur, une valeurs de décollage et c'est au programmeur de savoir comment sont arrangées ses données en mémoire. C'est intéressant.

Les premiers exemples étaient là pour introduire le sujet. Soit une matrice quelconque, (ligne, colonne, m*n, m*n*o, ....); comment remplacer les éléments de cette matrice (les indexes) par les valeurs contenus dans une autre matrice (en fct de l'index). Pour les deux premiers cas, on a une façon assez simple.

La première solution que tu me donnes

Code :

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

1
2
3
img(1:n) = bkg(A,1,1);
img(n+1:2*n) = bkg(A,1,2);
img(2*n+1:end) = bkg(A,1,3);

fonctionne très bien si la taille d'img = la taille d'A
Dans mon cas, img est 2 fois plus grande que A. Donc les index ne tombent pas en face.
(Je suis passé par une variable intermédiaire toto de la taille d'A puis scr(1:end/2,:,: ) = toto; mais ça gaspille de la mémoire pour rien).

Par contre,ça ne fonctionne que si tu connais le nombre et la taille des dimensions de A, bkg et img. (Dans mon cas c'est fixé donc pas de problème, mais si j'avais besoin de créer une fonction plus générique, ça ne marcherai pas.

Pareil pour la 2nd solution, ça marche très bien dans mon cas (merci à toi) mais ce n'est pas extrapolable a n dimensions de tailles inconnues.

(En plus j'ai été obligé de modifier ma colormap pour supprimer la 2eme dimension(qui est singleton))


Pour la beauté de l'art, y a-t-il une autre façon de faire ça ?

M'ci

[le fab]

Par contre,ça ne fonctionne que si tu connais le nombre et la taille des dimensions de A, bkg et img. (Dans mon cas c'est fixé donc pas de problème, mais si j'avais besoin de créer une fonction plus générique, ça ne marcherai pas.

Code :

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

doc size

[Jerome Briot]

Code :

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

1
2
3
img(1:n) = bkg(A,1,1);
img(n+1:2*n) = bkg(A,1,2);
img(2*n+1:end) = bkg(A,1,3);

fonctionne très bien si la taille d'img = la taille d'A
Dans mon cas, img est 2 fois plus grande que A. Donc les index ne tombent pas en face.
(Je suis passé par une variable intermédiaire toto de la taille d'A puis scr(1:end/2,:,: ) = toto; mais ça gaspille de la mémoire pour rien).

Je ne comprend pas bien. Tu veux dire que la variable img est définie avant d'exécuter ce morceau de code ?

Dans ce cas, soit tu la détruis à l'aide de clear, puis tu appliques la méthode que j'ai donnée précédemment, soit tu utilises cette autre méthode :

Code :

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

1
2
3
4
r = size(A,1);
c = size(A,2);
 
img(1:r,:,:) = reshape(bkg(A,1,:),r,c,[]);

[piter_pam]

reshape, c'est exactement ce qu'il me manquait

Merci les gars