Discussion :

[grohino]

Bonjour à tous,
Je travaille sur un projet d'imagerie, plus précisément de reconstruction d'objet 3D (méthode de décalage de phase, projection de lumière structurée).

mon probleme est le suivant:

je cherche à convertir une matrice de dimension (3,m,n) en une matrice de dimension (m,n), où les colonnes ont été remplacées par des chaines de caracteres.
les colonnes de la matrice 3D originale sont des permutations du triplet {1,2,3}.

mon but est d'attribuer à chacune de ces permutations
un numero spécifique{1,...,6},et ainsi d'obtenir une sorte de 'carte' de la matrice 3D d'origine.

Comme je travaille avec des images d'au moins 600x800, et que l'analyse se fait quasiment en temps réel, je ne dois pas avoir de boucles dans mon algorithme.


merci de voter aide, à bientot

[duf42]

Bonjour,

Je suis pas sur d'avoir bien compris ton problème, mais voici une petite idée:

Code :

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

1
2
3
4
 
M1 = arrayfun(@(x) num2str(x), M(1,:,:), 'UniformOutput', false);
M2 = arrayfun(@(x) num2str(x), M(2,:,:), 'UniformOutput', false);
M3 = arrayfun(@(x) num2str(x), M(3,:,:), 'UniformOutput', false);

Sachant que arrayfun te permet d'éviter la boucle pour la conversion des num en string.

Voilà, j'espère que ca pourra t'aider...

Bon courage,
Duf

[mr_samurai]

Salut,

Une autre petite idée :

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
 
p = 1000;
q = 1000;
 
M = zeros(3,p,q);
for i=1:p
    for j=1:q
        M(:,i,j) = randperm(3);
    end
end
 
tic 
 
M1 = reshape(M(1,:,:)+ 5*M(2,:,:)+7*M(3,:,:),p,q);
 
M1(M1==32)= 1; % cas [1 2 3]
M1(M1==30)= 2; % cas [1 3 2]
M1(M1==22)= 3; % cas [3 1 2]
M1(M1==20)= 4; % cas [3 2 1]
M1(M1==28)= 5; % cas [2 1 3]
M1(M1==24)= 6; % cas [2 3 1]
 
toc
% En 0.14 secondes sur ma machine pour 1000*1000

++ bonne chance

[grohino]

merci Mr samurai
par contre, je suis pas bien sûr d'avoir compris cette ligne:

Code :

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

1
2
3
 
 
M1 = reshape(M(1,:,:)+ 5*M(2,:,:)+7*M(3,:,:),p,q);

avec la suite du code, je suppose que les coefficients 1,5 et 7 servent à lever l'indétermination sur la somme des termes de la permutation?(sinon on aurait toujours 6 comme résultat?)

[grohino]

merci Mr samurai
par contre, je suis pas bien sûr d'avoir compris cette ligne:

Code :

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

1
2
3
 
 
M1 = reshape(M(1,:,:)+ 5*M(2,:,:)+7*M(3,:,:),p,q);

avec la suite du code, je suppose que les coefficients 1,5 et 7 servent à lever l'indétermination sur la somme des termes de la permutation?(sinon on aurait toujours 6 comme résultat?)

voilà ce qui arrive quand on répond trop vite à un message...j'ai oublié une chose: dans quelle matrice restockes-tu les informations extraites?
ma question n'est peut etre pas futée, mais j'ai le cerf-volant en ce moment...

[grohino]

ok...

c'était vraiment une question stupide, mais tant qu'on l'a pas posée...

les élés sont stockés dans M1, j'avais un doute la dessus, j'ai toujours un peu peur de me planter en laissant un meme nom pour des modifications de matrice.

merci beaucoup à vous!

[mr_samurai]

Salut,

Si on avais mis uniquement :

Code :

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

1
2
 
M1 = M(1,:,:)+ 5*M(2,:,:)+7*M(3,:,:);

M1 serait une matrice 3D de taille [1,p,q].
Alors que :

Code :

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

M1 = reshape(M(1,:,:)+ 5*M(2,:,:)+7*M(3,:,:),p,q);

produit une matrice 2D de taille [p,q].

Plus clairement , shape c'est forme . Et la fonction RESHAPE réarrange les données dans une autre forme (en p lignes et q colonnes dans notre cas).

++ bonne chance