Discussion :

[lilyla]

Bonjour,

Voici mon code

Code :

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

1
2
3
4
5
6
7
8
affiche=zeros(20,4); affiche(:,1)=1:20;affiche(:,2)=1;
for k=1:N
    reg=regionprops(L{k},'Centroid'); 
    affiche(:,2*k+1:2*k+2)=cat(1, reg.Centroid); % Merci à Dut en passant
end
for k=1:N
    j1=2*k+1;    x(:,j1)=affiche(:,j1);
end

Affiche est une matrice de 20 lignes et 2*k+2 colonnes (ici k=3):

Code :

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

1
2
3
4
5
affiche =
 
    1.0000    1.0000   50.2185  328.7130   49.5573  329.0815
    ...         ...       ...        ...        ...      ... 
   20.0000    1.0000  617.0956   34.4563  617.8324   34.2253

Je souhaiterais récupérer les colonnes d'indices impaires en vecteur xk sachant que je n'utilise pas les 2 1ères colonnes. Les colonnes d'indices paires correspondant à yk mais là le problème est le même!

Donc, lil me renvoie une matrice pour x que je veux vecteur ...?:

Code :

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

1
2
3
4
x =
         0         0   50.2185         0   49.5573
        ...        ...    ...         ...     ...
         0         0  617.0956         0  617.8324

En fait, je souhaiterais récupérer les N vecteurs x et les N vecteurs y (qui correspondent aux colonnes k>3 de ma matrice 'affiche' dans le but de calculer une distance car x et y sont des coordonnées!

Merci de vos conseils

[Jerome Briot]

une petite lecture préalable : Introduction à la gestion des matrices

[lilyla]

Je l'ai lu et plus d'une fois, j'ai même penser à utiliser la fonction reshape mais je ne pense pas que pour mon problème elle m'aide grandement.

Peut-être que je me trompe : je ne suis pas une experte, en fait je fais comme avec Matlab : du pas à pas pour avancer dans mon programme (et mes idées aussi d'ailleurs)

[Jerome Briot]

Essaie déjà ceci :

Code :

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

1
2
3
4
5
6
affiche = [1.0000    1.0000   50.2185  328.7130   49.5573  329.0815 
   20.0000    1.0000  617.0956   34.4563  617.8324   34.2253]
 
x=affiche(:,3:2:end)
 
y=affiche(:,4:2:end)

[lilyla]

j'essaye de mettre ça dans mon programme alors!

Merci et je reviens ici si je n'y arrive pas !

Par contre, je ne pensais pas qu'on pouvait faire i:pas:end comme ça dans la matrice, mais un test ne m'aurait rien couté merci beaucoup de ton aide!

Et en plus je ne fais pas de boucle ça c'est cool mais en fait j'étais persuadée que j'avais besoin d'une boucle !

[FR119492]

Salut !

Pourquoi ne fais-tu pas deux matrices distinctes, l'une avec les colonnes impaires et l'autre avec les colonnes paires, vu que, sur la base de ton exemple, elles semblent contenir des données de nature différente?

Jean-Marc Blanc

[lilyla]

En fait avec les lignes données par Dut:

Code :

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

1
2
x=affiche(:,3:2:end)
y=affiche(:,4:2:end)

, mon problème est résolu

Par contre, je souhaiterais savoir comment faire pour ne pas avoir ce problème de 1ères lignes que j'ai eu car j'ai à nouveau le problème :

Code :

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

1
2
3
4
5
6
xs=zeros(NbObjets{1},1);
ys=zeros(NbObjets{1},1);
for k=1:N-1
    xs=[xs x(:,k+1)-x(:,k)];
    ys=[ys y(:,k+1)-y(:,k)];
end

et pour xs il me renvoie la matrice que je veux sauf que j'ai en plus une ligne de zéros (qui correspond à mon initialisation mais sans ça, il ne peut fonctionner)... Donc si quelqu'un peut m'expliquer comment m'y prendre autrement, je suis preneuse

[Jerome Briot]

Tu dois initialiser xs et ys comme des vecteurs vides :

Code :

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

1
2
3
4
5
6
xs=[];
ys=[];
for k=1:N-1
    xs=[xs x(:,k+1)-x(:,k)];
    ys=[ys y(:,k+1)-y(:,k)];
end

Ou alors faire la préallocation de mémoire pour xs et ys et utiliser l'indexage dans la boucle :

Code :

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

1
2
3
4
5
6
7
8
9
xs=zeros(<Taille finale de xs>);
% Si j'ai tout suivi  : xs=zeros(size(x,1),N-1);
ys=zeros(<Taille finale de ys>);
% Si j'ai tout suivi  : ys=zeros(size(y,1),N-1);
 
for k=1:N-1
    xs(:,k)=x(:,k+1)-x(:,k);
    ys(:,k)=y(:,k+1)-y(:,k);
end

C'est l'idée

[lilyla]

MERCI!

[Jerome Briot]

Donc maintenant que tu as compris, tu peux remplacer directement la boucle FOR-END par ceci :

Code :

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

1
2
3
 
xs=diff(x,[],2)
ys=diff(y,[],2)