rapport de stage janus : radiographie par les rayons ... - HAL - IN2P3
rapport de stage janus : radiographie par les rayons ... - HAL - IN2P3
rapport de stage janus : radiographie par les rayons ... - HAL - IN2P3
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
CHAPITRE 4 : TRAVAIL SUR LA RECONSTITUTION DE L’IMAGE :<br />
a) Voxelisation :<br />
Comme on l’a vu précé<strong>de</strong>mment la reconstitution <strong>de</strong> l’image nécessite<br />
<strong>de</strong> <strong>par</strong>tager le volume à <strong>radiographie</strong>r en voxels.<br />
Les voxels sont <strong>de</strong> petits cubes d’arête Δx= Δy=Δz qui ont le même rôle<br />
que <strong>les</strong> pixels.<br />
Pour repérer <strong>les</strong> voxels dans l’espace on leur attribue 3 entiers :(I, J,K).<br />
Chaque triplet définit une zone <strong>de</strong> l’espace délimité <strong>par</strong> le voxel lui<br />
correspondant. Ici <strong>par</strong> exemple le centre du volume est représenté <strong>par</strong> (0,<br />
0,0). Si on se déplace du point central (0, 0,0) situé au centre du voxel (0,<br />
0,0) d’une distance égale à Δx dans <strong>les</strong> trois directions du repère cartésien<br />
on se trouvera au centre du voxel (1, 1,1).<br />
Comme on veut connaître dans quel voxel se trouve le point d’approche<br />
le plus près on exprime (I, J, K) en fonction <strong>de</strong> (x, y, z) cordonnées du point<br />
d’approche le plus près qui seront calculées dans la suite.<br />
On trouve facilement la relation liant I et x <strong>par</strong> un raisonnement sur<br />
l’axe (0,x). On déduit J en fonction <strong>de</strong> y et K en fonction <strong>de</strong> z <strong>par</strong> le même<br />
raisonnement sur <strong>les</strong> <strong>de</strong>ux autres axes du repère cartésien. Voilà ces trois<br />
relations :<br />
⎛ Δx<br />
⎞<br />
⎜ x +<br />
I sgn ( x) E 2 ⎟<br />
= ⎜ ⎟<br />
⎜ Δx<br />
⎟<br />
⎝ ⎠<br />
⎛ Δy<br />
⎞<br />
⎜ y +<br />
J sgn ( y) E 2 ⎟<br />
= ⎜ ⎟<br />
⎜ Δy<br />
⎟<br />
⎝ ⎠<br />
⎛ Δz<br />
⎞<br />
⎜ z +<br />
K sgn ( z) E 2 ⎟<br />
= ⎜ ⎟<br />
⎜ Δz<br />
⎟<br />
⎝ ⎠<br />
Où E est la <strong>par</strong>tie entière.<br />
36