Discussion :

[lilyla]

Bonjour,

J'essaye de faire tourner les lignes suivantes:

Code :

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

1
2
3
4
5
6
for k=1:N
   le=length(reg{1}); 
   for i=1:le
      centre{k}=[reg{k}(i).Centroid(1)] % vecteur à revoir
   end
end

ce n'est qu'une partie du code...

Il s'agit donc de la ligne en gras qui ne revoie pas un vecteur comme on le souhaiterait (dans le but de l'intégrer par la suite dans une matrice).

Merci de votre aide
Cordialement,

[LordPeterPan2]

Bonjour,

J'essaye de faire tourner les lignes suivantes:

Code :

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

1
2
3
4
5
6
7
8
9
10
      for k=1:N
le=length(reg{1}); 
  A = (1:le)';         B = ones(le,1);        C= zeros(le,1);         D= zeros(le,1);
    
for i=1:le
    centre{k}=[reg{k}(i).Centroid(1)] % vecteur à revoir
end
 
%     matrice = [A B C D]
end

ce n'est qu'une partie du code...

Il s'agit donc de la ligne en gras qui ne revoie pas un vecteur comme on le souhaiterait (dans le but de l'intégrer par la suite dans une matrice).

Merci de votre aide
Cordialement,

Ta variable reg{k}(i).Centroid(1) ne peut en aucun cas etre un vecteur (au mieux un cellule contenant un vecteur).

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
 
>> Q=cell(4,1);
>> Q{1}=[1 2 3];
>> Q{2}=rand(2,4);
>> Q(1)
 
ans = 
 
    [1x3 double]
 
>> Q{1}
 
ans =
 
     1     2     3

[lilyla]

oui en testant il me renvoie que ça contient des 'cells'

Je n'arrive pas à voir ce que je dois faire...
Merci encore

[Jerome Briot]

Si reg est calculé à partir de REGIONPROPS :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
for k=1:N
   le=length(reg{1}); 
   A = (1:le)';         
   B = ones(le,1);       
   C= zeros(le,1);        
   D= zeros(le,1);
 
      centre{k}=[reg.Centroid];
% ou 
% centre{k}=cat(1, reg.Centroid);
 
 
%     matrice = [A B C D]
end

[lilyla]

Il me renvoie une erreur :

Attempt to reference field of non-structure array.
Error in ==> pour_tester at 56
centre{k}=[reg.Centroid];

et j'ai la même erreur si j'utilise l'autre ligne ...

J'utilise effectivement la fonction regionprops.

En fait, j'ai besoin de créer une matrice avec :
- en 1ère colonne : la numérotationj des objets (k)
- en 2nd colonne : une ligne de 1 (pour l'instant)
- en 3ème colonne: les coordonnées en -x- du centre de l'objet n°k
- en 4ème la même chose en -y-...

Et je rencontre un problème pour les 2 dernières colonnes que j'essaye de résoudre depuis plusieurs jours tout de même ... Et malgré vos pistes (merci beaucoup déjà!) , rien n'y fait!

Je vous remercie encore de l'aide que vous pourrez m'apporter j'espère et bonne semaine à tous

Avec le programme complet, peux-être que cela sera plus clair... ainsi que mes explications

[Jerome Briot]

A la place de ceci :

Code :

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

reg{1}=regionprops(L{k},'All');

Je ferais simplement :

Code :

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

reg=regionprops(L{k},'All');

Et le code que j'ai donné précédemment devrais fonctionner

[lilyla]

En changeant ce que tu me dis:

centre = [1x40 double]

Je n'arrive pas à comprendre pourquoi je n'obtiens pas mon vecteur!
Merci de ton aide !

[Jerome Briot]

La variable "centre" est un tableau de cellules => : Les tableaux de cellules (cell array)

[lilyla]

Au final :

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
24
25
26
27
28
29
30
31
32
    N=input('quel est le nombre total d''images? ');
for k=1:N
    fname=sprintf('levures + 50mS par m + 3Vpp +44MHz_t00%d(c1).jpg',k);
    I{k}=imread(fname);
end

se=strel('disk',10);
for k=1:N
    background{k}=imclose(I{k},se); 
    I2{k}=imsubtract(background{k},I{k}); 
    level=graythresh(I2{k});
    bw{k}=im2bw(I2{k},level); 
    bw{k}=bwareaopen(bw{k},200);   
    bw{k}=imclearborder(bw{k}); 
    fill{k}=imfill(bw{k},'holes'); 
        
    b{k}=bwboundaries(fill{k});
    L{k}=bwlabel(fill{k}); 
    reg=regionprops(L{k},'All');
    le=length(reg); 
    
    num = (1:le)'; 
    test = ones(le,1);

    for i=1:le
         centre{k}=cat(1, reg.Centroid);
    end   
            
    X=centre{k}(:,1);
    Y=centre{k}(:,2);
    matrice = [num test X Y]
end

Bon alors tout ceci fonctionne pour une image (k=1)
Après je rencontre un problème de dimension, ce qui est normal mais je ne vois pas comment le résoudre!

en fait c'est la ligne :

Code :

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

reg=regionprops(L{k},'All');

qui pose problème apparemment.

Je voudrais arriver à créer une matrice 'centre' qui correspondrait aux coordonnées x,y de chaque image comme suit:
objet 1: x(k) y(k) x(k+1) y(k+1)
objet i: x(k) y(k) x(k+1) y(k+1)
objet le: x(k) y(k) x(k+1) y(k+1)


J'espère que vous voyez à peu près ce que je souhaite faire...

Et voici sinon les images sur lesquelles je travaille:

[lilyla]

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
24
25
26
27
28
29
30
31
32
33
34
35
36
N=input('quel est le nombre total d''images? ');
tail_disq=10; 
seuil=800; 
petite_aire=200;
 
for k=1:N
    fname=sprintf('levures + 50mS par m + 3Vpp +44MHz_t00%d(c1).jpg',k);
    I{k}=imread(fname);
end
 
se=strel('disk',tail_disq);
for k=1:N
    background{k}=imclose(I{k},se); 
    I2{k}=imsubtract(background{k},I{k}); 
    level=graythresh(I2{k});
    bw{k}=im2bw(I2{k},level); 
    bw{k}=bwareaopen(bw{k},petite_aire);   
    bw{k}=imclearborder(bw{k}); 
    fill{k}=imfill(bw{k},'holes'); 
 
    [etiquette{k},NbObjets{k}]=bwlabel(fill{k},8);
    b{k}=bwboundaries(fill{k});
    L{k}=bwlabel(fill{k}); 
end
 
affiche=[(1:NbObjets{k})' ones(NbObjets{k},1)]; 
for k=1:N
    reg=regionprops(L{k},'All');
    num = (1:NbObjets{k})'; 
    test = ones(NbObjets{k},1);
 
    centre{k}=cat(1, reg.Centroid);
    X{k}=centre{k}(:,1);
    Y{k}=centre{k}(:,2);
    affiche=[affiche X{k} Y{k}]
end

Et un grand merci aux personnes qui ont contribué à la résolution de ce problème qui n'est peut-être pas le dernier!

Bien cordialement,

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
affiche=[(1:NbObjets{k})' ones(NbObjets{k},1)]; 
for k=1:N
    reg=regionprops(L{k},'All');
    num = (1:NbObjets{k})'; 
    test = ones(NbObjets{k},1);
 
    centre{k}=cat(1, reg.Centroid);
    X{k}=centre{k}(:,1);
    Y{k}=centre{k}(:,2);
    affiche=[affiche X{k} Y{k}]
end

Cette partie s'écrit "simplement" :

Code :

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

1
2
3
4
5
6
7
8
9
10
affiche=zeros(NbObjets{1},6);
affiche(:,1)=1:NbObjets{1};
affiche(:,2)=1;
 
for k=1:N
 
    reg=regionprops(L{k},'Centroid');
    affiche(:,2*k+1:2*k+2)=cat(1, reg.Centroid);
 
end

[lilyla]

Je corrige cela et merci encore pour ton aide