[regionprops] Problème avec une boucle ou une notation

Version imprimable

16/11/2007, 13h28
lilyla

[regionprops] Problème avec une boucle ou une notation
Bonjour,

J'essaye de faire tourner les lignes suivantes:
Code:

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,

Citation:
Envoyé par lilyla
Bonjour,

J'essaye de faire tourner les lignes suivantes:
Code:

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:

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

16/11/2007, 16h09
lilyla

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

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

Si reg est calculé à partir de REGIONPROPS :

Code:

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

19/11/2007, 09h58
lilyla

Il me renvoie une erreur :

Citation:

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 ;)
19/11/2007, 10h48
Jerome Briot

A la place de ceci :

Code:

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

Je ferais simplement :

Code:

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

Et le code que j'ai donné précédemment devrais fonctionner
19/11/2007, 10h52
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 !
19/11/2007, 10h53
Jerome Briot

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

Au final :
Code:

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:

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:

Enfin!

Code:

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,

Citation:

Envoyé par lilyla

Code:

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:

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

21/11/2007, 10h01
lilyla

Je corrige cela et merci encore pour ton aide