SIMART : un programme MATLAB de simulation numérique des artefacts de susceptibilité magnétique en IRM.

Ce programme Matlab met en oeuvre l'algorithme 1 présenté à la section 7.2.2 de la thèse. Nous considérons ici le cas où l'implant métallique est une bille métallique de rayon r₀ et de susceptibilité magnétique X_m. Nous disposons dans le cas de la bille uniformément aimanté d'une formule donnant le champ magnétique induit ce qui facilite la mise en oeuvre d'un programme de simulation.

La modélisation mathématique du processus de création des artefacts de susceptibilité magnétique ayant permis d'obtenir l'algorithme est décrit dans la thèse. Cette modélisation est présentée de manière plus concise dans cet article.

Dans le fichier artef.m, les variables suivantes sont modifiables:

global r0; r0 = 1e-2;     % rayon de la bille (en mètre)
global Xm; Xm = 3e-4;     % susceptibilité magnétique du matériau
global  C;  C = [0 0 0];  % origine du plan de coupe
global B0; B0 = 1;        % intensité du champ B_0 en tesla
global ec; ec = 3e-3;     % épaisseur de coupe en mètre
global  g;  g = 1e-2;     % intensité des gradients en Tesla par mètre
global ns; ns = [0 1 0];  % direction du gradient de sélection de coupe
global nr; nr = [0 0 1];  % direction du gradient de lecture
global np; np = [1 0 0];  % direction du gradient de phase
global Lr; Lr = 6e-2;     % longueur de la coupe dans le direction nr du gradient 
                            de lecture en mètre
global Lp; Lp = 6e-2;     % longueur de la coupe dans le direction np du gradient 
                            de phase en mètre

Pour avoir la signification précise de ces variables, consulter la modélisation présenté dans cet article.

Le programme

Il est constitué des fichiers suivants:

artef.m : le programme principal.
voxel.m
fndrc.m
polnom.m
F.m
calBdB.m

Ces fichiers sont à copier dans un même répertoire. Lancez une session Matlab puis la commande artef.

Exemple de simulation

Voici un exemple d'image d'un artefact engendré par une bille métallique obtenue par simulation numérique et une image IRM expérimentale.

Figure 1: Image brute obtenue par simulation (à gauche) et image IRM (à droite).

retour à la page d'accueil.

Dernière modification : 4 avril 1999
(c) Copyright 1999 - Stéphane Balac - Dept. de Mathematiques - Université de Bretagne Occidentale