Discussion :

[yabbiyou]

Bonjour,

Afin de filtrer des images j ai besoin de creer des masques en formes d anneaux.

Pour cela je cree deux cercle pleins de differents diametre que je soustrait l un a l autre. Pour pouvoir multiplier mon mask par la matrice a analyser je dois discretiser l interieur de mes cercles sur les entiers:
voici le code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
for x=-round(r/2):round(r/2)
    for y=-round(r/2):round(r/2)
        u=x^2+y^2;
        if u<(r/2)^2
            g=g+1;
            pts(g,1)=x+xo;
            pts(g,2)=y+yo;
            mat(x+xo,y+yo)=1;
        else bibi=0;
        end
    end
end

Le code marche bien, mais le gros soucis c est que les images sur lesquelles je travail sont de tres grandes tailles (4000*5000), et que la creation des mask prend un temps tres (trop!) important.
Est ce que vous avez une idee (ou des idees!) pour ameliorer l algo pour acceler les calculs.

Merci d avance
Pierre

[Jerome Briot]

=> Création d'une image 3D qui contient une sphère (Niveaux gris)

A noter que dans la solution que je donnais dans cette discussion, le FIND est inutile :

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
% Diamètre du disque et taille de la matrice
d=30;
 
% Création de la matrice 2D remplie de 0
X=zeros(d);
 
% Tous les couples de coordonnées des pixels de X
[col,lig] = meshgrid(1:d); 
 
% Translation du centre de la matrice en (0,0)
col(:) = col(:)-d/2-0.5; 
lig(:) = lig(:)-d/2-0.5;
 
% Recherche des indices des pixels contenu dans le disque
% idx=sqrt(col(:).^2+lig(:).^2)<=d/2
idx = (col(:).^2+lig(:).^2)<=(d^2/4);
 
% Mise à 1 de la valeur de ces pixels
X(idx) = 1;
 
% Affichage
figure
imagesc(X)
axis image

[yabbiyou]

Merci pour la reponse, (rapide et efficace...comme toujours)
J ai change un petit detail pour permettre la creation d un petit cercle dans une grande matrice, si ca interesse quelqu un...

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
function [ mat]=points_cercle(r,pts_mat)
 
% Diamètre du disque et taille de la matrice
d=2*r;
 
% Création de la matrice 2D remplie de 0
X=zeros(pts_mat);
 
% Tous les couples de coordonnées des pixels de X
[col,lig] = meshgrid(1:pts_mat); 
 
% Translation du centre de la matrice en (0,0)
col(:) = col(:)-pts_mat/2-0.5; 
lig(:) = lig(:)-pts_mat/2-0.5;
 
% Recherche des indices des pixels contenu dans le disque
% idx=sqrt(col(:).^2+lig(:).^2)<=d/2
idx = (col(:).^2+lig(:).^2)<=(d^2/4);
 
% Mise à 1 de la valeur de ces pixels
X(idx) = 1;
mat=X;
% Affichage
% figure
% imshow(mat2gray(mat))
% axis image