DWMRI:DSIstudio
⭐ Se pueden descargar los datos desde esta liga. 📺 El video del tutorial está en youtube ℹ️ Los datos fueron adquiridos en un resonador GE MR Discovery 750 de 3 teslas.
-
dwis.nii.gz: Imágenes de difusión -
dwis.bvec: Vectores del gradiente de difusión -
dwis.bval: b-valores -
revpe.nii.gz: Dos imágenes b=0 con codificación de fase en sentido inverso. Se usará para pre-procesamiento dedwis.nii.gz -
revpe.bvec: vectores de difusión -
revpe.bval: b-valores.
Todos estos archivos fueron convertidos usando dcm2niix.
Las imágenes dwis.nii.gz tienen los ojos alargados, porque fueron adquiridas con EPI positive blips (P>>A); son muchos volúmenes (106). En cambio, las imágenes revpe (4.nii.gz) tienen negative blips y son solamente 2 volúmenes (una b=0 y una DWI).
No entraremos en detalles de preprocesamiento en este tutorial. Se recomienda leer la documentación de eddy en el sitio de fsl. Se sugiere usar una máquina con CUDA y fsl versión 6.0.2 o superior. No correrlo con CUDA implica varias horas de preprocesamiento.
Al terminar veremos que las imágenes están mucho menos distorsionadas:
Las imágenes preprocesadas de este tutorial, pre-cocinadas para fines prácticos, se llaman dwi_preproc_corrected.{bvec,bval,nii}.
Revisamos que todo esté en orden. Usaré mrtrix3 para generar el tensor de difusión. Me gusta el formato .mif, pero bien podría escribir .nii.gz si me diera la gana.
Primero genero una máscara:
dwi2mask -fslgrad dwi_preproc_corrected.{bvec,bval,nii} dwi_preproc_corrected_mask.nii
Ahora estimo el tensor de difusión:
dwi2tensor -fslgrad dwi_preproc_corrected.{bvec,bval,nii} -mask dwi_preproc_corrected_mask.nii dt.mif
extraigo métricas del tensor:
tensor2metric -fa fa.mif -vector v1.mif -adc adc.mif dt.mif
Y para visualizar los tensores sin que me distraigan los tensores gigantes del LCR, le pongo una máscara basada en ADC:
mrcalc adc.mif 0.001 -lt dt.mif -mul dt_masked.mif
Para finalmente visualizar:
Aunque yo soy muy feliz con mrtrix, en este documento utilizaré dsi-studio para hacer tractografía de manera automática. Existen versiones para todos los sistemas operativos. El programa es muy lógico y nos muestra el orden de los pasos que haremos. En este tutorial nos quedaremos en la creación de tractografía, y no lo llevaremos al análisis de conectividad (connectometry).
Empezamos cargando los datos. Primero que nada, nos aseguramos que tenemos archivos .bvec y .bval con el mismo nombre que nuestro archivo .nii.gz. En mi caso tenemos dwi_preproc_corrected.nii, dwi_preproc_corrected.bvec, y dwi_preproc_corrected.bval.
Al terminar de cargar (unos cuantos segundos), nos mostrará la B-table, una tabla que nos muestra los valores
La parte inferior de esta tabla me pregunta dónde guardar el archivo que se generará, y con qué nombre. Este archivo contendrá toda la información de las imágenes y los bvecs y bvals.
En la pantalla principal de dsi-studio, ahora vemos el archivo que creamos (terminación .src.gz). Le damos doble clic para comenzar a procesar.
En la primer pestaña podemos hacer control de calidad de nuestros datos. Dando clic en cada columna de la B-table podemos ver la imagen correspondiente, y cambiar rebanadas. Los botones nos permiten identificar rebanadas o volúmenes como datos a ignorar (en caso de que tengan artefactos).
En el siguiente paso generamos una máscara binaria para únicamente procesar donde hay tejido. También seleccionamos el modelo de difusión a ajustar (DTI, GQI, o QSDR). Seleccionaré GQI, y doy clic a Run Reconstruction. Como podemos ver, podríamos seleccionar únicamente DTI si nos interesara, pero GQI de todas maneras va a hacer DTI, así que tenemos ambos con mínimo esfuerzo.
Esto generará un archivo .fib en la misma carpeta donde habíamos creado el archivo .src.gz, y podremos encontrarlo en la pantalla principal de dsi-studio. Le damos doble clic.
Esto abrirá una nueva pantalla para todo lo relacionado a tractografía.
Agreguemos la imagen T1 para que podamos hacer un mejor registro de imágenes hacia el atlas y que la tractografía automática funcione mejor. Esto lo hacemos en Slices -> Insert other images. Y seleccionamos nuestra imagen T1 en formato .nii.gz. En Slices-->Adjust registration podemos ver qué tan bien o mal fusionadas están.
Ahora habilitamos autotrack.
Y en Target Podemos seleccionar el fascículo que queramos ver, y le damos Fiber tracking
Por ejemplo, el Arcuate_Fasciculus_L :
Esto aún lo podemos editar, usando regiones de interés que dibujamos en el panel inferior izquierdo, donde pueden apreciarse los streamlines que intersectan con el plano de imagen.
Por ejemplo, vemos unos streamlines que se van demasiado inferiores en el lóbulo temporal, y queremos quitarlos:
Dibujamos una región donde intersectan, y la convertimos en ROA (region of avoidance).
Damos clic derecho a nuestro tracto Arcuate_Fasciculus_L y seleccionamos Filter Tracks by ROI/ROA/END.
Veremos que esos streamlines desaparecen.
