04.07.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!