Cercle pour image processing

**totoc1001** · 14/12/2006, 03h27

Bonsoir,

J'aurais besoin d'un petit conseil;
Je dois "nettoyer" une image en y appliquant un cercle et en effaçant tout ce qui y est a l'extérieur.
Mais j'ai un peu de mal avec les cercles, quelqu'un pourrait-il me dire si il existe une fonction déjà toute faite pour ça, ou comment faire ?
Je ne me souviens plus comment appliquer un cercle.
Mon algo serait le suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
Pour tous les points X de mon image
     Pour tous les points Y de mon image
          Si on est a l'intérieur du cercle
               Copier la valeur de l'ancienne image dans la nouvelle
          Sinon (Si on est à l'extérieur du cercle)
               La valeur de la nouvelle image sera 0
          fin si
     fin pour Y
fin pour X

Mon problème c'est le SI on est a l'intérieur du cercle.
Je connais le rayon du cercle, et aussi son centre.

Merci beaucoup
Thomas

**Kangourou** · 14/12/2006, 14h16

tu regardes si la distance euclidienne de ton point au centre du cercle est inferieure au rayon.
distance euclidienne + sqrt(dx*dx+dy*dy), avec dx et dy les differences de coordonnnees en x et y.

je crois pas qu'il y ait de fct toute faite, tu peut regarder 'strel' qui te genere un element structurant, mais tu auras surement un peu d'emballage a faire.

A+

**totoc1001** · 14/12/2006, 14h55

Je vais essayer de faire ca, merci.
Si j'ai des soucis je reviens ;-)

Merci encore
Thomas

**totoc1001** · 15/12/2006, 02h22

Si ca peut aider quelqu'un, voici ma solution, elle est extrèmement simple et fonctionne parfaitement:

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
 
function TestCircle
    Radius=500;
    A=500;
    B=500;
    for y=1:1000
        for x=1:1000
            CircleEquation=sqrt(((x-A)^2)+((y-B)^2));
            if(CircleEquation>Radius)
                NewImage(y,x)=0;
            else
                NewImage(y,x)=255;
            end
        end
    end
    figure(1)
    colormap('gray');
    imshow(NewImage);
 
end

Celle-ci permet déterminer les points a l'intérieur du cercle de ceux a l'extérieur.

Bye
THomas

**tug83** · 15/12/2006, 11h02

merci pour ton programme, juste une remarque si je peux me permettre, l'executer tel quel met 22 secondes sur ma machine, en préallouant la mémoire j'arrive à aller 100 fois plus vite , environ 0.2 secondes

juste avant ta boucle for j'ecris:

Code :

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

NewImage=zeros(1000);

**Kangourou** · 15/12/2006, 18h18

on peut aussi vectoriser le code :

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
 
 
% nb de points dans chaque direction
Nx=1000;
Ny=1000;
A=500; B=500;R=500;
 
% vecteur de coordonnees dans caque directions
lx = 1:Nx;
ly = 1:Ny;
 
% matrice des coordonnees de chaque pixel, centree sur le cercle
x = repmat(reshape(lx, [1 Nx]), [Ny 1]) - A;
y = repmat(reshape(ly, [Ny 1]), [1 Nx]) - B;
 
% on cree une image binaire comme resultat d'une operation de comparaison
img = (x*x+y*y)<R*R;

Par contre attention, c'est plus gourmand en memoire...

**totoc1001** · 15/12/2006, 22h10

Envoyé par tug83

merci pour ton programme, juste une remarque si je peux me permettre, l'executer tel quel met 22 secondes sur ma machine, en préallouant la mémoire j'arrive à aller 100 fois plus vite , environ 0.2 secondes

juste avant ta boucle for j'ecris:

Code :

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

NewImage=zeros(1000);

Merci pour le conseil, je ne savais pas le fait de ne pas définir la taille de la matrice pouvait tellement ralentir le temps d'éxécution.
Je te remercie, je vais peut etre pouvoir gagner pas mal de temps sur mes scripts, je n'initialise que très rarement la taille des tableaux.
Simplement parce que comme ca je suis sur d'avoir la matrice la plus petite possible.

Merci pour le tuyau ;-)

**tug83** · 16/12/2006, 15h11

tu pourrais regarder la note technique suivante pour comprendre l'importance de la préallocation:
http://www.mathworks.fr/support/solu...lution=1-18150

**salseropom** · 18/12/2006, 10h44

Envoyé par totoc1001

Si ca peut aider quelqu'un, voici ma solution, elle est extrèmement simple et fonctionne parfaitement:

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
 
function TestCircle
    Radius=500;
    A=500;
    B=500;
    for y=1:1000
        for x=1:1000
            CircleEquation=sqrt(((x-A)^2)+((y-B)^2));
            if(CircleEquation>Radius)
                NewImage(y,x)=0;
            else
                NewImage(y,x)=255;
            end
        end
    end
    figure(1)
    colormap('gray');
    imshow(NewImage);
 
end