Discussion :

[matt67]

Bonjour,

J'utilise un code qui me permet en connaissant la matrice de rotation (3x3) et celle de translation (vecteur 3x1) de recaler des points sur des autres. Le problème est que matrice de points Pinterb est de taille 3x70000 (points en 3D) et que lorsque j'utilise la transformation pour recaler chacun de mes points dans une boucle for, le programme rame et plante finalement... Voici le code que j'utilise ci-dessous :

Code :

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

1
2
3
for i=1:size(Pinterpb,2)          %Calcul des points recalés                
   Q(:,i)=inv(R)*(Pinterpb(:,i)-q(5:7));
end

Pinterpb est ma matrice de points 3x70000 à recaler, R la matrice de rotation et q le vecteur translation : c'est une transformation inverse en fait que j'applique.

Merci de votre aide s'il existe un moyen pour optimiser et ne pas faire planter l'ordi !

[Caro-Line]

1ère chose à faire : as-tu penser à pré-allouer ta matrice :

Code :

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

1
2
 
Q=zeros(size(Pinterb));

[matt67]

Merci, impeccable , j'avais oublié de le faire pour celle-ci...

[Jerome Briot]

2ème chose :

Code :

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

1
2
3
4
5
6
 
invR=inv(R);
Q=zeros(size(Pinterb));
for i=1:size(Pinterpb,2)          %Calcul des points recalés                
   Q(:,i)=invR*(Pinterpb(:,i)-q(5:7));
end

Ou mieux, l'utilisation de l'opérateur \ (backslash)