Discussion :

[AlphaJ]

Bonjour à tous,
Je cherche une méthode simple afin de construire une image 3D (grey scale) qui contient une sphère ?
la fonction sphere me donne des X, Y, Z mais sous forme matrice 2D pour chacun ! comment je fais pour obtenir les coordonnées en pixels x, y, z ?

Merci beaucoup en avance

[Jerome Briot]

Qu'entends-tu par "image 3D" ?

[AlphaJ]

L'image 3D est tout simplement une matrice 3D ceci est une reconstruction 3D d une série d'images consecutives, pour simplifier les chose je suppose une matrice de 30*30*30.

[Jerome Briot]

En 2D (plus simple pour comprendre) :

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=find(sqrt(col(:).^2+lig(:).^2)<=d/2);
 
% Mise à 1 de la valeur de ces pixels
X(idx)=1;
 
% Affichage
figure
imagesc(X)
axis image

Et en 3D :

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

[Kangourou]

Salut,

utilise la fonction meshgrid pour generer les coordonnées x, y, et z de chaque pixel,
puis tu peux calcuelr pour chaque pixel si il appartient ou non a la sphere :

Code :

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

1
2
3
 
[x y z] = meshgrid(1:30, 1:30, 1:30);
img = (x-15).^2+(y-15).^2+(z-15).^2 < 10^2;

pour une sphere de rayon 10 centree autoru du pixel (15,15,15).

A+

[AlphaJ]

Merci beaucoup, après vos réponses je peux bien utiliser le mesh