voisinage d'objets sur Matlab

[invite105e858b]

Bonjour à tous

j'ai un problème depuis une semaine et j'ai pas arrivé à trouver une solution et je compte énormément sur votre aide.

J'ai besoin d'un code en matlab qui me permet de détecter le voisinage de chaque objets étiqueté dans une image binaire. En effet je voulais enregistrer chaque voisin trouvé de l'objet d’intérêt dans une structure car je voudrais après calculer le barycentre des voisins de chaque objet.
Par exemple pour l'objet d'étiquette 3 ses voisins seront les objets d'étiquettes 5 et 7.

j'ai commencé avec ce code:

%%%%%%%%%%%%%%%%##################Extraction des zones d'interet dans l'image binaire
 
[nb,C]=i4mat_components (w, h,I);%%Décomposition en composantes annexes: la fonction retourne C: la matrice des labels et nb : le nombre des objets dans l'image

%%%La boucle suivante permet de determiner la matrice des vecteurs
%%%caractéristiques de chaque objets.
 for i=1:nb
     
     %[r,c]=find(C==i);
%A=zeros(r,c);
     % A=[r,c];
     A=zeros(size(C));
     A(C==i)=1;
    % figure,imagesc(A)
     [M]=feature_vec(A); 
    % Ap=double (A);


%ici commence le code qui m'interesse

 labelled(:,i )=i;
  % assume labelled array is called "labelled"
label = i;  % label of object to find neighbours of
object = labelled == label;
se = ones(3);   % 8-connectivity for neighbours - could be changed
neighbours = imdilate(object, se) & ~object;
neighbourLabels = unique(labelled(neighbours));
disp(neighbourLabels);

 
B(:,i )= abs(M);


 end;

mais ce code ne m'affiche pas le bon résultat
j'espère que vous pouvez m'aidez et merci d'avance
-----

[invite2d7144a7]

Bonjour,

demande d'aide : en voilà un titre explicite ...

[invite105e858b]

salut whoami,

en principe on est là pour ça, aider et demander l'aide non ?

[invite2d7144a7]

Bonjour,

salut whoami,

en principe on est là pour ça, aider et demander l'aide non ?

Oui, mais on peut exiger des titres en rapport avec la demande.

En conséquence à ta réaction, je vais m'abstenir de t'aider.

[A voir en vidéo sur Futura]

[invite105e858b]

Bonjour,

c'est la première fois que je partage un truc sur un forum et je sais pas qu'un tel titre peux vous déranger. En fait j'ai juste traduit mon besoin.
je te remercie en tout cas

[Kissagogo27]

Bonjour, préciser par exemple sur quoi vous avez besoin d'aide par exemple

Besoin d'aide sur Logiciel Matlab un modérateur peut modifier le titre en conséquence

[invite105e858b]

Bonjour,

merci pour ton conseil je vais le prendre en considération prochainement

[JPL]

c'est la première fois que je partage un truc sur un forum et je sais pas qu'un tel titre peux vous déranger.

Quand tu t'es inscrit le formulaire t'a demandé de lire la charte du forum et de cocher une case pour certifier que tu l'avais lue. Si tu l'avais lue tu aurais vu :

Les titres des messages doivent être explicites.

Tu aurais vu aussi pourquoi j'ai écrit ce message en vert.

[Antoane]

Bonjour,

Pour moi, ce n'est pas clair...
Qu'est ce qu'une image binaire ? Un exemple ?
qui sont w, h et I ?
Il serait bon de donner des infos sur les fonction que tu utilises lorsqu'elles ne sont pas natives de matlab (i4mat_components, feature_vec, etc).

Présenter un code bien indenté (sélectionne ton code dans l'éditeur Matlab et presse CTRL+I ), sans lignes inutiles (commentées) mais avec des commentaires aiderait.

[invite105e858b]

Bonjour,

je vous remercie pour votre réponse

une image binaire est une image pour laquelle chaque pixel ne peut avoir pour valeur que 0 ou 1.
J'ai pas expliqué ces fonctions : i4mat_components et feature_vec car j'en ai pas besoin de les utiliser pour la partie de détection de voisinage de chaque objet. Mais comme tu me l'avais demandé voici tout le code :

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Algorithme de la fonction de détection des objets rares
%%%%image est l'input de la fonction DetectionAlgo
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear all

I=imread('D:\PFE\matlab_rech\images\imgSynthetic.tif');
sz=size(I);
%% nombre de lignes
w=sz(1,1);
% nombre de colonnes
h=sz(1,2);
% taille de la fenêtre
tailleF= 30
%%fixer la taille de la fenètre
nr = floor(w/tailleF);
nc = floor(h/tailleF);

%%%Calcul des inavaiants de Hu pour chaque fenètre de calcul suceptible de
%%%contenir un objet.
for r = 1:nr
    for c = 1:nc
        subimg = I((r-1)*tailleF+1:tailleF*r,(c-1)*tailleF+1:tailleF*c,:);
        
    end
end
%%%%%%%%%%%%%%%%##################Extraction des zones d'interet dans l'image binaire

[nb,C]=i4mat_components (w, h,I);%%Décomposition en composantes annexes: la fonction retourne C: la matrice des labels et nb : le nombre des objets dans l'image


%%%La boucle suivante permet de determiner la matrice des vecteurs
%%%caractéristiques de chaque objets.

for i=1:nb
    
    A=zeros(size(C));
    A(C==i)=1;
    
    [M]=feature_vec(A); %calculer les invariants de Hu pour chaque objet A de l'image
    
    labelled(:,i )=i;

    % assume labelled array is called "labelled"
    label = 1;  % label of object to find neighbours of
    object = labelled == label;
    se = ones(3);   % 8-connectivity for neighbours - could be changed
    neighbours = imdilate(object, se) & ~object;
    neighbourLabels = unique(labelled(neighbours));
    disp(neighbourLabels);
    
    
    B(:,i )= abs(M);
    
    
end;
figure,imagesc(C),colorbar


%%%%%%%%%%%%%%%%##################Calcul du Vecteur Caractéristique de chaque objet de
%%%%%%%%%%%%%%%%l'image

[invite9dc7b526]

Je crois qu'il faudrait surtout que tu expliques ce que sont un "objet" et un "voisinage". Je pense deviner qu'un objet est une composante connexe (mais même ça demande à être précisé) de l'ensemble des pixels de valeur 1 de l'image, mais qu'est-ce qu'un voisinage?

[invite105e858b]

Comme vous l'avez dis, un objet est une composante connexe où tout les pixels ont la valeur 1.

Pour le voisinage c'est tout autre objet qui peut être assez proche de l'objet d’intérêt car après cette étape je dois calculer la distance euclidienne de chaque objet de l'image avec que ses voisins pour dégager enfin la dispersion géométrique de chaque objet dans l'image.

[invite9dc7b526]

Donc tu te donnes une distance, disons D0, telle que deux objets distants de moins de D0 soient considérés comme voisins (?) la distance entre deux objets étant le minimum des distances entre pixels des deux objets (?) ou la moyenne (?)

[invite105e858b]

plutôt je crois qu'il faute travailler avec une fenêtre glissante de taille (3*3) ou (5*5) peut être comme j'ai déjà commencer à faire dans cette partie :

     labelled(:,i )=i;
    % assume labelled array is called "labelled"
    label = 1;  % label of object to find neighbours of
    object = labelled == label;
    se = ones(3);   % 8-connectivity for neighbours - could be changed
    neighbours = imdilate(object, se) & ~object;
    neighbourLabels = unique(labelled(neighbours));
    disp(neighbourLabels);

mais le problème que ce code ne fonctionne que sur le dernier objet ( de numéro nb)

[invite9dc7b526]

A mon humble avis, avant de commencer à programmer, il faudrait définir clairement ce qui est recherché.