Traitement d'image par ACP

Tout ce qui concerne les différents logiciels mathématiques (Matlab, Maxima, Scilab, R, Maple, Mathematica, ...) et les calculatrices scientifiques.
Règles du forum
Merci d'éviter le style SMS dans vos messages et de penser à utiliser la fonction Recherche avant de poster un message. Pour joindre des fichiers à vos messages, consulter ce sujet.
> Penser à utiliser les balises Code pour poster du code.
> Préciser le nom du logiciel concerné dans le titre de votre message.

Traitement d'image par ACP

Messagepar mariona » Samedi 09 Février 2008, 13:33

Bonjour,

je doit effectuer une analyse en composantes principales pour des images. J'ai écrit le script et souhaiterai savoir si ce dernier est correcte:

Code: Tout sélectionner
clear all
%%% on charge les données %%%
database=load('essai_realdefau_1_kHz.dat');


xmin=250
xmax=500
ymin=250
ymax=500



%%%%%%%%%%%reduction de ma matrice de données%%%%%%%%%%%%%%%%%
for j=xmin:xmax
    for k=ymin:ymax
        database_reduit(k-ymin+1,j-xmin+1)=database(k,j);

        [nb_ex,nb_feat] = size(database_reduit);
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%Calcul de la matrice de donnée centrée et normée%%%%%%%%%%%%%%
for j=1:nb_feat
    rowsum = sum(database_reduit(:, j))/nb_ex;
    ect=std(database_reduit(:, j));
    for k=1:nb_ex
        database1=(database_reduit(k,j)-rowsum)/(ect)^0,5;
       database2_reduit(k,j)=database1;
    end
end
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%% on calcule la matrice de covariance %%%
S = cov(database2_reduit);

%%% on calcule les vecteurs propres et les valeurs propres %%%
[vec_p, L] = eig(S);

%%% on trie les vecteurs propres selon les valeurs propres %%%
[val_p, ind] = sort(diag(L));
val_p = flipud(val_p);
ind = flipud(ind);
vec_p2 = vec_p(:,ind);
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%Calcul de la matrice d'estimation%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%(Matrice vecteurs propores*Matrice données centrée%%%%%%%%%%%%%%%%%
S_chapo=vec_p2'*database2_reduit
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%coeff=corr2(S_chapo,database2_reduit)


Puis, je souhaiterai représenter une "nouvelle" image, en ne tenant plus compte du premier vecteur propre (seulement de tous les autres). Je ne sais pas comment faire?

Merci pour votre aide!
mariona
Utilisateur
 
Messages: 2
Inscription: Samedi 09 Février 2008, 12:26
Statut actuel: Post-bac | Ecole d'ingénieur

Publicité

Re: Traitement d'image par ACP

Messagepar keane » Vendredi 25 Avril 2014, 22:21

bonsoir mariona je suis étudiant en master 2 telecom je suis entrin de faire mon projet de fin d'etude est justement votre message il ma interpeller par ce que c le titre de mon projet et j un très grand problème c que je maîtrise pas beaucoup le matlab et je pêne a trouver le code de mon programme (réduction de donne d'une image en utilisant l'acp) alors es-que si possible de m aider est de me propose un code détaille ça sera vraiment très gentille de votre part (*@yahoo.fr) merci !! :D :D
Dernière édition par MB le Dimanche 27 Avril 2014, 16:21, édité 1 fois.
Raison: Suppression de l'adresse mail.
keane
Utilisateur
 
Messages: 1
Inscription: Vendredi 25 Avril 2014, 21:53
Statut actuel: Post-bac | Master

Re: Traitement d'image par ACP

Messagepar rebouxo » Samedi 26 Avril 2014, 13:07

Ça, cela s'appelle du déterrage de sujet.

Olivier
A line is a point that went for a walk. Paul Klee
Par solidarité, pas de MP
rebouxo
Modérateur
 
Messages: 6703
Inscription: Mercredi 15 Février 2006, 13:18
Localisation: le havre
Statut actuel: Actif et salarié | Enseignant

Re: Traitement d'image par ACP

Messagepar MB » Dimanche 27 Avril 2014, 16:22

Éviter également d'indiquer votre adresse mail dans le message.
MB (Pas d'aide en Message Privé)
Merci d'utiliser $\LaTeX$ (voir ici) et d'éviter le style SMS pour la lisibilité des messages.
MB
Administrateur
 
Messages: 6744
Inscription: Samedi 28 Mai 2005, 13:23
Localisation: Créteil
Statut actuel: Actif et salarié | Enseignant


Retourner vers Logiciels mathématiques

 


  • Articles en relation
    Réponses
    Vus
    Dernier message

Qui est en ligne

Utilisateurs parcourant ce forum: Aucun utilisateur enregistré et 2 invités