Discussion :

[matt67]

Bonjour,

Ayant fait l'acquisition SCAN d'un fantôme de tête contenant une grille régulière 3D de billes de 3 mm espacées d'1 cm, je souhaite extraire la grille de points correspondant aux centres de gravité des billes.

J'ai du mal à trouver la méthode à adapter...
Ci-dessous le code actuel et ci-joint l'image sur laquelle je souhaite détecter les centres de gravité :

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
%% Récupération des informations de l'image DICOM traitée %%
 
    Ima1info=dicominfo(name);
    Ima1= dicomread(Ima1info);       
    figure (1)
    imshow(Ima1,[],'InitialMagnification','fit');
 
    %% Fitrage adaptatif du bruit blanc %%
 
    Ima2 = wiener2(Ima1,[5 5]);
    figure (2)
    imshow(Ima2,[],'InitialMagnification','fit');
 
    %% Définition de l'élément structurant et opération bottom hat
 
    se = strel('disk',10,0);  %ES forme disque de rayon lissé
    Ima3 = imbothat(Ima2,se);
    figure (3)
    imshow(Ima3,[],'InitialMagnification','fit');
 
    %% Hessian filter %%
    I=double(Ima3);
    Ivessel=FrangiFilter2D(I);
    figure (4)
    subplot(1,2,1), imshow(Ima1,[]);
    subplot(1,2,2), imshow(Ivessel,[0 0.25]);

Je récupère actuellement l'image DICOM, que je traite pour bruit blanc puis bottom hat et filtre Hessian. On obtient une image où les centres de gravité sont bien visibles (points) mais il y a encore beaucoup trop d'éléments parasitaires autour...

Merci de toute aide que vous pourrez m'apporter !!!