Rotation d'un element structurant

Version imprimable

30/07/2010, 16h41
ImagingAllthe

Rotation d'un element structurant

Bonjour,
je souhaiterai faire roter un ES de type strel
C'est possible via rot90, mais ca limite les angles de rotation.
Je vais le coder, mais peut etre y as t il une facon direct de le faire ...

ImagingAll
31/07/2010, 09h03
Jean Dumoncel

Bonjour,

quel element structurant précisément?
02/08/2010, 10h06
ImagingAllthe

un element type rectangle ...

Code:

b1 = strel('line',50,2);

Pour faire une ouverture par des segments de droite d'épaisseur.

ImagingAll
02/08/2010, 11h04
kmaniche

Bonjour;

Tu peux par exemple utiliser imrotate :

Code:

>> doc imrotate

Citation:

Envoyé par kmaniche

Tu peux par exemple utiliser imrotate :

Un exemple :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
function se = strel2(len,wid,deg)
 
if nargin==0
    len = 12;
    wid = 4;
    deg = 45;
end
 
nhood = ones(wid,len);
 
nhood = imrotate(nhood,deg);
 
nhood = logical(nhood);
 
se = strel('arbitrary', nhood);

Merci Dut,
Ca donne ca avec une tite fonction, mais c'est un peu lourd en calcul...

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
function [IMGout] = RadialOpenrect_SA(IMG,len,wid,step,affichages)
% INPUT 
% [IMGout] = RadialOpenrect_SA(IMG,len,wid,step,affichages)
% affichages: show image result if ==1 
 
 IMGout=im2uint8(zeros(size(IMG)));
 
for x=0:step:180
      deg = x;
      nhood = ones(wid,len);
      nhood = imrotate(nhood,deg);
      nhood = logical(nhood);
      se = strel('arbitrary', nhood);
      radialopen=imopen(IMG,se);
      RadialopeningIMG=max(im2uint8(IMGout),im2uint8(radialopen));
end
 
if (affichages==1)   
    figure, imshow(RadialopeningIMG);
else
end
 
 
end

Bonne journée,

ImagingAll

Dans un premier temps, la boucle devrait pouvoir se simplifier ainsi :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
nhood = ones(wid,len);
nhood = logical(nhood);
 
n = round(180/step);
deg = linspace(0,180,n);
 
for k = 1:n
      nhoodk = imrotate(nhood,deg(n),'nearest');
      se = strel('arbitrary', nhoodk);
      radialopen = imopen(IMG,se);
      RadialopeningIMG = max(im2uint8(IMGout),im2uint8(radialopen));
end

oui c'est un peu plus rapide et sans cette coquille:
deg(n)->deg(k)
Ca donnne ca:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
nhood = ones(wid,len);
nhood = logical(nhood);
 
n = round(180/step);
deg = linspace(0,180,n);
 
for k = 1:n
      nhoodk = imrotate(nhood,deg(k),'nearest');
      se = strel('arbitrary', nhoodk);
      radialopen = imopen(IMG,se);
      RadialopeningIMG = max(im2uint8(IMGout),im2uint8(radialopen));
end
 
if (affichages==1)   
    figure, imshow(RadialopeningIMG);
else
end

ImagingAll

03/08/2010, 11h28
Jerome Briot

L'autre solution consiste à pré-calculer les éléments structurants (si possible) et à les stocker sous forme d'un tableau de cellules dans un fichier mat qui serait appelé en début de fonction (comme une sorte d'abaque).

Cela économiserait les fonction IMROTATE et STREL dans la boucle FOR-END
05/08/2010, 17h25
ImagingAllthe

Avec les images en niveau de gris ... c'est super long.
Alors je nettoie.

ImagingAll