Discussion :

[lilyla]

Bonjour,

Voici mon code utilisé pour visualiser l'axe majeur sur mes objets détectés dans les images afin de me rendre compte de la fiabilité du paramètre "orientation". Mes objets étant très petits je souhaite les visualiser un à un AVEC leur axe majeur correspondant. Ce que, donc, je n'arrive pas à faire dans ce cas précis! La partie à revoir est en gras dans mon programme.

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=1; tail_disq=10; seuil=100; petite_aire=25; se = strel('disk',tail_disq); 
for k=1:N 
    fname = sprintf('20060516 rot 50khz 8vpp_t%03d(c1).jpg',k); 
    I{k} = imread(fname); 
    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');  [L{k},NbObjets{k}] = bwlabel(fill{k},8); 
    b{k} = bwboundaries(fill{k});  rego{k} = regionprops(L{k},'All'); 
end

for k=1:N    
    for i=1:NbObjets{k}
        cellule{k}(i).x = rego{k}(i).Centroid(1); cellule{k}(i).y = rego{k}(i).Centroid(2);
        cellule{k}(i).orientation = rego{k}(i).Orientation; 
        cellule{k}(i).maj = rego{k}(i).MajorAxisLength; cellule{k}(i).min = rego{k}(i).MinorAxisLength;
        x1{k}(i) = cellule{k}(i).x + cellule{k}(i).maj*cos(-cellule{k}(i).orientation*pi/180)/2;
        y1{k}(i) = cellule{k}(i).y + cellule{k}(i).maj*sin(-cellule{k}(i).orientation*pi/180)/2;
        x2{k}(i) = cellule{k}(i).x - cellule{k}(i).maj*cos(-cellule{k}(i).orientation*pi/180)/2;
        y2{k}(i) = cellule{k}(i).y - cellule{k}(i).maj*sin(-cellule{k}(i).orientation*pi/180)/2;
    end
end

for k=1 % pour visualiser les objets un par un dans une image donnée
    for i=1:NbObjets{k}
        bx{i}=zeros(1,4);
        bx{i}=rego{k}(i).BoundingBox;
        ob{i}=I{k}(bx{i}(2)-.5:bx{i}(2)+bx{i}(4)-.5,bx{i}(1)-.5:bx{i}(1)+bx{i}(3)-.5,:);
        figure(i)
        imagesc(ob{i}), title(['cellule n° ' num2str(i)])
        colormap('gray')
        figure(i)
        hold on
        plot([x1{k}(i) x2{k}(i)],[y1{k}(i) y2{k}(i)])
    end
end

Quelqu'un aurait une idée une suggestion voire la solution à ce problème?

Merci d'avance de l'aide apportée qui me fait avancer à grand pas!

[lilyla]

Voici une image sur laquelle j'aimerais pouvoir obtenir N nouvelles images (N étant le nombre d'objet détecté dans mon image) qui seraient chacun un zoom d'un objet avec le tracé de l'axe majeur correspondant que je récupère avec les coordonnées de 2 points (x1, x2, y1 et y2).

Faut peut-être aussi revoir la taille de la fenêtre prise en compte dans ob{i} mais je ne sais pas trop...

Merci !

[Jerome Briot]

Utilise l'information BoundingBox donnée par REGIONPROPS

Quelque chose comme :

Code :

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

1
2
3
4
5
R = reg.BoundingBox;
 
subimg = bw(R(1):R(1)+R(3),R(2):R(2)+R(4));
 
imagesc(subimg);

C'est l'idée...

[lilyla]

Si je lance ces quelques lignes pour une image et un objet j'ai déjà l'erreur suivante:

Code :

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

1
2
3
4
5
6
Warning: Integer operands are required for colon operator when used as index.
> In visualisation at 39
??? Index exceeds matrix dimensions.
 
Error in ==> visualisation at 39
        subimg = bw(R(1):R(1)+R(3),R(2):R(2)+R(4));

Mais, j'ai déjà dans mon code ob{i} qui est à quelque chose près ce que tu me dis non? Et ce que je voudrais c'est rajouter à chacune de ces images (qui sont donc des zooms sur objets) LA droite correspondant à l'axe majeur de l'objet sur lequel je zoome...

Merci de votre aide !