Discussion :

[safojava]

Bonsoir,
J'ai besoin d'utiliser l'ACP sur un ensemble d'images après faire la projection des composants (en matlab)
J'ai extrait les nouvelles composantes mais comment les projeter par exemple juste la 1ère composante et 2ème
voici le code pour le moment

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
moy = mean(matFace);% matface matrice d'images
[nligne ncol] =size(matFace);
vunit=linspace(1,1,nligne)';
Moy=matFace-vunit*moy;
ect=std(matFace);
ecar=diag(ect)^-1 ;
%%% matrice centrée réduite
Xcr=Moy*ecar;
%%%%matrice de covariance
Cov=(Xcr'*Xcr)/nligne;
[vec val] = eig(Cov) ;%vecteur et valeur propre
col=diag(val);
U=[col vec];
U=flipud(U);
U=U(:,2:ncol+1);
 
CP=Xcr*U;
for i=1:ncol
    ele=CP(:,i);
    elt=reshape(ele,30,20);
    figure
    imshow(elt),
end

[Franck Dernoncourt]

Pourquoi n'utilises-tu pas la fonction princomp() ? http://www.mathworks.com/help/toolbo.../princomp.html

[Jean Dumoncel]

Pourquoi n'utilises-tu pas la fonction princomp() ? http://www.mathworks.com/help/toolbo.../princomp.html

Tout le monde ne possède pas la toolbox statistic et l'ACP n'est pas très compliquée à coder avec les outils de base de Matlab.


@sofajava : je ne comprends pas bien ce que tu cherches à faire avec ces lignes :

Code :

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

1
2
3
U=[col vec];
U=flipud(U);
U=U(:,2:ncol+1);

Il suffit de réordonner dans l'ordre décroissant les valeurs propres :

Code :

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

1
2
3
4
5
[val,idx] = sort(diag(val),'descend');
vec = vec(:,idx);
 
% Projection :
CP=Xcr*vec;

[safojava]

Bonsoir,
merci pour vos réponses:
@Frank: j'ai déjà la toolbox mais il faut la programmer
@magelan: oui c'est la même chose, juste ajouter vecteur propre et après les ordonnées j'ai pas le problème ici mais les problèmes que j'ai sont:
en premier lieu comment donner un pourcentage pour juste travailler avec les premiers vecteurs (grands) est-ce qu'il y a une formule ou non ?
Le 2ème problème si j'ai une nouvelle image comment la projeter avec les autres images est-ce que je dois refaire le même processus ou il faut un changement quelque part?
Merci d'avance.

[Jean Dumoncel]

Sur la matrice des valeurs propres val correctement ordonnées, pour avoir le pourcentage cumulé :

Code :

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

cumsum(diag(val)) / sum(diag(val))

le 2 eme pb si j'ai une image nouveau comment la projeté avec les autres images est ce que je refaire la mm processus ou il faut un changement kelle que part et mercii d'avance

Tout dépend ce que tu veux faire, si tu veux projeter cette nouvelle image il faut multiplier par les vecteurs propres trouvés précédemment, si tu veux utiliser cette image dans le calcul de tes composantes principales, il faut refaire le calcul avec cette nouvelle donnée.

[safojava]

Merci mais comment par exemple juste utiliser les 5 premiers val d'après le pourcentage ?
Pour l'image nouveau il faut calculer la composant principale...
C'est une question bête, mais comment puis-je dessiner les composantes 1 et 2?
Je fais

Code :

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

1
2
3
s1=CP(:,1);
s2=CP(:,2);
figure;plot(s1,'*',s2,'+')

Mais ça ne marche pas.
Je voudrais représenter la 1ère composante avec * et la 2ème avec +.

[Jean Dumoncel]

Merci mais comment par exemple juste utiliser les 5 premiers val d'après le pourcentage ?

Il te suffit de conserver uniquement les 5 premières composantes de ton tableau contenant les projections. Avec ta notation :
s1a5=CP(:,1:5);

C'est une question bête, mais comment puis-je dessiner les composantes 1 et 2?
Je fais

Code :

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

1
2
3
s1=CP(:,1);
s2=CP(:,2);
figure;plot(s1,'*',s2,'+')

Mais ça ne marche pas.
Je voudrais représenter la 1ère composante avec * et la 2ème avec +.

Tu veux les représenter sur un axe 1D? Une solution rapide est de mettre les ordonnées à 0 :

Code :

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

plot(s1,zeros(size(s1)),'*',s2,zeros(size(s2)),'+')

[safojava]

Merci c'est très gentil.
Je veux essayer avec la première solution, mais pour dessiner les composants en 2D pour voir la dispersion des variables.

[Jean Dumoncel]

Je n'ai pas compris ce que tu a écrits...

Habituelement, on trace la deuxième composante en fonction de la première (ou la troisième en fonction de la premier ou la troisième en fonction de la deuxième, etc ou 3 composante sur un graphe 3D) avec plot3:

Code :

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

plot(s1,s2,'+')

[safojava]

okk moi je cherche par exemple de faire la différence entre par exemple le composant 1 et 2 de faire + et * s'il est possible

[Jean Dumoncel]

Désole je ne vois pas de quoi tu veux parler.

[safojava]

d'accord, pas grave, peut être c'est mon problème j'ai pas bien compris la représentation des composants 1 et 2 , on représente la composante par rapport a l'autre ou on présente les 2 composants?

[Jean Dumoncel]

Essaie de faire un effort sur la mise en forme de tes réponses, le langage sms, ça ne donne pas envie de te lire et ça n'aide pas du tout à te comprendre...

Tout est possible mais ce qui se fait généralement, c'est de représenter en 2D avec la 1ere en abscisse, et la 2eme composante en ordonnée, et d'attribuer des marqueurs à tes différentes classes d'images par exemple (si le but est bien de trouver les composantes les plus discriminantes par rapport à différentes classes d'image).

[safojava]

ok merci pour tout