El procesamiento se realiza con la version 0.4 del pipeline

Determinar los valores de anisotropia fraccional (FA) y difusividad media(MD) en los tractos que conectan la corteza Fronto-Parietal para los sujetos del protocolo FACEHBI.

La adquisición DTI se ha realizado con una orientación distinta al T1w. Esto por sí solo no seria un problema pero los cortes se han calculado para adquirir la materia blanca en su totalidad y no el cerebro completo.

El resultado es que las imagenes DTI estan cortadas en la parte superior de la cabeza y el corregistro a cualquier espacio es casi imposible.

Se procede en varios pasos par aintentar hacer un registro decente,

1. Se registra (flirt, fnirt) la imagen FA al atlas FMRIB58_FA
2. Usando esta transformacion (applywarp) se registra el B0 a espacio MNI
3. Se aplica una mascara (fslmaths) del B0 registrado sobre el template MNI y se obtiene un espacio MNI recortado al DTI del sujeto
4. Se registra (antsRegistrationSyN.sh) el template MNI a espacio nativo T1
5. Se usa la transformacion (antsApplyTransforms) para pasar el template MNI recortado a espacio nativo
6. Se aplica una mascara (fslmaths) de este template recortado al T1 y se obtiene un T1 recortado
7. Se registra (epi_reg) el DTI a este T1 recortado
8. Se invierte la matriz de la transformacion (convert_xfm) y se convierte a formato ANTS (c3d_affine_tool).Ya tenemos una transformación valida del T1 al B0!
9. Se pasa (antsApplyTransforms) el T1 a espacio nativo DTI usando esta matriz
10. Se pasa el template MNI a espacio nativo T1 y de ahi a espacio nativo DTI
11. Se pasan los Atlas JHU a espacio nativo T1 y de ahi a espacio nativo DTI

Basicamente,

$ dti_reg.pl -chop -time '12:0:0' facehbi

Una vez efectuado el registro correcto pasamos a planificar la tractografia. Obtenemos la corteza Fronto-Temporal a partir de la segmentacion de FS. Las regiones FP son,

ctx-caudalmiddlefrontal
ctx-inferiorparietal
ctx-middletemporal
ctx-parsopercularis
ctx-parstriangularis
ctx-postcentral
ctx-precentral
ctx-superiorfrontal
ctx-superiorparietal
ctx-superiortemporal
ctx-supramarginal

Los FS-LUT (incluyendo rh y lh) han de ser,

1003
1008
1015
1018
1020
1022
1024
1028
1029
1030
1031
2003
2008
2015
2018
2020
2022
2024
2028
2029
2030
2031

Se toma la segmentacion de FS y se lleva a espacio nativo DTI (Ya tenemos la matriz correcta). De aqui se toman los LUT correspodientes.

Usando las ROIs ya en espacio nativo DTI se lanza la tractografia,

Todo esto se agrupa con el comando,

$ dti_track.pl -t1 -time '12:0:0' facehbi

Ahora hay que extraer los valores de FA y MD en los tractos calculados,

$ for x in `ls -d working/*_probtrack_out`; do mv $x `echo $x | sed 's/out/FPCustom/'`;done
$ dti_metrics_tracks.pl -path FPCustom facehbi

Ademas de esta ROI, voy a tomar los nodos predeterminados de redes funcionales conocidas (DMN, LN, SN).El procedimiento es basicamente el mismo pero ahora la tractografia es basada en atlas.

$ dti_track.pl -t1 -time '12:0:0' -uofm DMN facehbi

$ for x in `ls -d working/*_probtrack_out`; do mv $x `echo $x | sed 's/out/DMN/'`;done
$ dti_metrics_tracks.pl -path DMN facehbi