$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
1. Reconstruire des surfaces corticales 3D
Cette première partie du protocole utilise l'oléoduc FreeSurfer standard comme décrit dans le Wiki ( http://surfer.nmr.mgh.harvard.edu/fswiki ). Notez que les commandes détaillées ici décrire un moyen de réaliser les reconstructions de surface corticale, mais les commandes équivalentes peuvent également être utilisés.
- Importer les premières IRM DICOM dans FreeSurfer et vérifier la qualité de l'image (par exemple, que l'orientation est correcte, le contraste suffisant et les images n'est pas déplacé). Ce procédé utilise les commandes suivantes (remplacez le texte entre <...> (inclus) avec des valeurs appropriées à une instance spécifique, et «#» désigne les commentaires):
mksubjdirs # Créer l'architecture du dossier utilisé par FreeSurfer
cd / Mri # aller dans le dossier IRM de votre sujet
mri_convert cm tkmedit 001.mgz # de visualiser le volume converti
- Créer des modèles en trois dimensions mailles corticale 8,9. Afin de faire face au problème de la enterrés sillons, FreeSurfer première crée un volume de matière blanche unitaire, qui est utilisé comme point de départ pour la première gris-blanc de surface. Cette surface est ensuite optimisé en fonction du gradient local d'intensité et élargi à l'interface gris-CSF.
Recon-tout-s # Lancer la reconstruction de la surface corticale
A la fin du processus de reconstruction, vous obtiendrez deux modèles de maillage composé d'environ 150 000 points pour chaque hémisphère: un blanc (gris-blanc de l'interface) et d'une pie-mère (gris-CSF interface) de surface. Il est important de noter que toutes les surfaces et les volumes restent dans l'espace natif, permettant de mesure tels que le volume, surface, l'épaisseur ou l'indice gyrification être mesurée sans déformation.
- Vérifier l'exactitude de ces surfaces reconstruites:
tkmedit T1.mgz? H.pial # la surface blanche est superposée en vert et la surface piales en rouge
où h désigne l'hémisphère:? lh.pial pour l'hémisphère gauche et rh.pial pour l'hémisphère droit. La figure 2 (en 2 versions: une image GIF animée à être inclus dans le film et un statique pour le site) montre un exemple de corriger les reconstructions des surfaces blanches et de la pie-mère "Bert" sujet distribué avec le paquet FreeSurfer. Si vous devez corriger manuellement le résultat du processus de reconstruction, vous trouverez un tutoriel sur le Wiki FreeSurfer ( http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/WhiteMatterEdits ,
2. Calculer l'indice gyrification locales
Lorsque vous êtes satisfait de vos surfaces, calculer l'indice gyrification locale (L IG) en utilisant la commande:
Recon-tout-LGI-s
Cette commande exécute habituellement pendant environ 3 heures pour les deux hémisphères d'une participante à l'étude, en fonction de la puissance de votre poste de travail. Les différentes étapes du processus de l IG sont passés en revue dans la Fig. 1. L'informatique commence avec la création d'une surface extérieure en utilisant opération de fermeture morphologique. Cette surface extérieure, notée? H.pial_outer_smoothed, est en outre illustré dans la Fig. 3. Puis, environ 800 se chevauchant les régions d'intérêt circulaires sont créées sur la surface extérieure. Pour chacune de ces régions, une région correspondant d'intérêt est défini surla surface piales. Le calcul entier finit avec la création d'une carte individuelle contenant un L IG valeur pour chaque point de la surface corticale (soit ~ 150 000 valeurs par hémisphère).
3. Vérifiez le résultat du calcul l IG pour chaque hémisphère
tksurfer ? H pial-overlay / Surf /? H.pial_lgi-fthresh 1
Les valeurs L IG sont superposées sur la surface corticale. Comme je corriger les valeurs IG sont généralement comprises entre 1 et 5, fixant le seuil minimum à 1 (avec le fthresh option) permet une rapide vérification: vous ne devriez pas voir toute la zone grise corticale. Un exemple de résultat individuel correct est indiqué dans la Fig. 4.
4. Comparaisons entre les groupes de la statistique
Le but est de quantifier l'effet de groupe à chaque sommet sur la surface corticale tout en contrôlant l'effet du sexe et deâge. Vous aurez besoin de suivre le même processus que si vous voulez comparer l'épaisseur corticale à chaque sommet, mais de donner? H.pial_lgi au lieu de? H.thickness. Deux options sont possibles pour calculer les comparaisons entre les groupes statistiques: les commandes classiques sont listés en premier, et l'interface graphique (Qdec) est brièvement mentionné par la suite.
- La première option pour comparer les résultats L IG entre les groupes, utilisez les commandes ci-dessous; de plus amples détails peuvent être obtenus à https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/GroupAnalysis .
- D'abord vous aurez besoin pour créer un modèle d'étude spécifique donnant à tous vos sujets en entrée:
make_average_subject - sujets ...
La commande ci-dessus va créer un objet nommé "moyenne". Alternativement, vous pouvez utiliser l'objet "FSAVErage »distribuée comme une partie de la distribution FreeSurfer. - Ensuite, créez le fichier texte contenant la description des sujets impliqués dans votre étude (le «Groupe FreeSurfer descripteur de fichier"). Votre FSGD.txt devrait ressembler à ceci:
GroupDescriptorFile 1
Groupe Control_Male
Groupe Control_Female
Groupe Patient_Male
Groupe Patient_Female
Variables de l'âge
Entrée Patient_Male 20
Entrée Control_Female 23
<...> - Rééchantillonner les données L IG dans l'espace de l'objet en moyenne en utilisant la commande suivante pour chaque hémisphère:
mris_preproc - FSGD.txt FSGD - cible moyenne -? hémi h - AME pial_lgi - hors h.lgi.mgh
- Lisser les données sur la surface corticale de réduire le rapport signal à bruit:
mri_surf2surf -? hémi h - s en moyenne -? sval h.lgi.mgh - FWHM 10 - h.10.lg TVAL?i.mgh
- Calculer la comparaison avec un groupe au niveau de chaque sommet. Pour cela vous devrez créer un fichier texte de contraste (par exemple dans le cas de la FSGD.txt décrit ci-dessus, le "contrast.txt" va contenir les valeurs «1 1 -1 -1 0" pour calculer la différence entre les contrôles et patients tout en contrôlant pour l'âge et le sexe). Lancez enfin la comparaison:
moyenne de surfer h - - mri_glmfit - y h.10.lgi.mgh -? FSGD FSGD.txt Doss -? glmdir h.lgi.glmdir? C contrast.txt
- Visualisez les résultats sur votre sujet en utilisant la moyenne tksurfer:
tksurfer moyenne? h gonflés
Puis charger le fichier en tant que superposition sig.mgh trouve dans le dossier? H.lgi.glmdir / contrast.txt / sig.mgh. En utilisant l'option "configurer overlay" vous pouvez encore modifier le seuil de p ainsi que correcte pour les comparaisons multiples en utilisant le taux de fausses découvertes 10.
- L'option alternative pour le groupe de comparaison est d'utiliser Qdec, une interface utilisateur graphique mis en œuvre dans FreeSurfer. L'utilisation de Qdec avec index gyrification locale implique de pré-lisser le L de données GI:
Recon-tout-Qcache mesure pial_lgi moyenne -S
Avec Qdec, le Groupe FreeSurfer descripteur de fichier est remplacé par une version légèrement différente, le tableau de données (qdec.table.dat) qui comprend la description des différents groupes et d'autres variables confondantes telles que l'âge. Une description détaillée de l'utilisation des Qdec est fourni à http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/QdecGroupAnalysis .
Fait à noter, si l'IG L n'est pas disponible dans la liste des variables dépendantes dans Qdec, vous devez ajouter la ligne suivante dans le fichier Qdecrc situé dans votre répertoire home.:
Mesure1 = pial_lgi 5. Analyse
Alternativement, les analyses statistiques pourraient éventuellement être calculé au niveau de la parcellisation corticale intégré dans FreeSurfer 11. À cette fin, la moyenne des valeurs L IG peuvent être extraites pour les 34 régions gyrale d'intérêt pour chaque hémisphère, et ces mesures peuvent encore être comparées entre les différents groupes. Cette analyse de colis-sage (par opposition à l'analyse de vertex-sage décrit ci-dessus) pourrait être intéressante car elle limite la quantité de comparaisons statistiques. Toutefois, l'IG l à chaque point quantifie l'gyrification dans la zone environnante circulaire, de sorte que l'IG L en moyenne dans une région d'intérêt gyrale reflète aussi dans une certaine mesure l'gyrification dans les régions voisines de l'intérêt.
Enfin, bien que les questions les plus importantes ont été décrites dans ce protocole, une solution aux autres problèmes qui peuvent être rencontrés lors de la FreeSurfer ou l IG de traitement peut être trouvée dans les archives de la mailing list FreeSurfer ( http://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferSupport ).
6. Les résultats représentatifs
Comme décrit dans la section 1c du protocole, vous devriez toujours vérifier soigneusement l'exactitude de la reconstruction des surfaces corticales avant l IG calcul. Pendant le défilement entre le frontal et le lobe occipital, une attention particulière que les navires et la membrane ne sont pas inclus dans la surface piales. Vérifiez également que la surface blanche suit précisément l'interface gris-blanc. Un exemple de reconstruction correcte est fournie dans la figure 2 (voir la figure gif animé pour l'ensemble du volume).
A la fin du calcul l IG, vous aurez également pour vérifier le résultat pour les deux hémisphères de chaque sujet.Il ne doit pas être une zone corticale avec un résultat L IG inférieur à 1. L'article 3 du protocole et de la figure 4 montre comment pour vérifier si la sortie correcte du calcul l IG est correcte.

Figure 1. Aperçu du calcul L IG. Tout d'abord, des modèles tridimensionnels mailles corticales sont reconstruits à partir des images brutes en utilisant le pipeline FreeSurfer standard. Ces algorithmes de reconstruction d'utiliser un volume de matière blanche binaires comme point de départ pour surmonter le problème des sillons enterrés. Les modèles de maillage cortical comprend généralement environ 150 000 sommets et sont classiquement utilisés pour calculer l'épaisseur corticale à chaque point. De même, l'indice gyrification locale (L IG) sera calculé à chaque sommet. À cette fin, une surface extérieure est créé. Puis correspondants régions circulaires d'intérêt sont identifiés sur l'extérieur uned de la surface corticale en utilisant l'algorithme d'appariement. Après environ 800 de générer des zones de chevauchement des intérêts, le processus aboutit à la création de cartes individuelles de l IG. Ces cartes peuvent être facilement interprété: un indice de 5 signifie qu'il ya cinq fois plus de surface corticale invaginée dans le sillon de la région environnante que le montant de la surface corticale visible; un indice de 1 signifie que le cortex est à plat dans la zone environnante . Enfin, les comparaisons statistiques sont calculées du groupe au niveau de chaque sommet, à l'instar des comparaisons épaisseur corticale.
Figure 1B. Individuels carte corticale de L IG. Ce petit film montre une rotation de 360 degrés d'une carte individuelle LGI corticale comme indiqué dans la Fig. 1. Il est frappant de constater que les régions corticales supérieures aux valeurs de l'IG L correspondent à la fois première à être créée durant la vie in utero: la scissure de Sylvius, le sillon temporal supérieur et les ensillon traparietal sur la vue latérale du cerveau, et le sillon pariéto-occipitale sur la vue médiale du cerveau. Voir la vidéo

Figure 2. Exemple de reconstruction de la surface corticale suffisante (une coupe coronale). Après la fin du processus de reconstruction, les surfaces corticales doivent être correctement vérifiée à travers l'ensemble du volume cérébral. La surface intérieure corticale (notée surface blanche, en vert sur l'image) devrait suivre précisément l'interface gris-blanc. La surface corticale externe (c'est à dire gris-CSF interface, notée de surface piales, ici en rouge) ne devrait pas inclure n'importe quel morceau de navire ou de la membrane. Fait à noter, l'exemple présenté ici utilise le "Bert" sujet distribué avec le paquet FreeSurfer.
La figure 2B. Exemple de la surface corticale suffisante reconstruction (volume maximal). Cette image gif animé montre la surface corticale de l'hémisphère gauche du "Bert" sujet sur chaque coupe coronale, tel que vu par le défilement de la plus frontale aux sections les plus occipitale coronale avec FreeSurfer. Voir la vidéo

Figure 3. Exemple de surface extérieure calculée comme une partie du processus de l IG (une coupe coronale). La première étape dans le calcul de l IG est la création d'une surface extérieure enveloppant l'hémisphère. Cette surface (notée? H.pial_outer_smoothed dans FreeSurfer) peut être vérifiée à l'aide tkmedit. Ici, le «Bert» distribué avec l'objet FreeSurfer est utilisé comme un exemple.
La figure 3B. Exemple de surface extérieure calculée comme une partie du processus de l IG (volume maximal). Cette image gif animémontre la surface externe de l'hémisphère gauche sur chaque coupe coronale, tel que vu par le défilement de la plus frontale aux sections les plus occipitale coronale avec tkmedit dans FreeSurfer. movie Voir

Figure 4. Exemple de sortie correcte l IG telle qu'elle apparaît avec FreeSurfer. Orientations différentes de la surface corticale de la "Bert" sujet avec l IG valeurs superposées. Le code couleur est le défaut "chaleur" overlay comme on le voit avec les tksurfer dans FreeSurfer. L'utilisation d'un seuil minimal de 1, tous les sommets doivent être colorés et aucune aire corticale devrait APpoire en gris. Fait à noter, la superposition de couleurs peuvent être modifiées en utilisant l'option "Configurer Overlay" dans tksurfer, où les valeurs minimales et maximales, ainsi que l'histogramme de la répartition globale de l'IG L peut également être vérifiée.