Microstructural Superficial White Matter Covariance
This module samples superficial white matter intensities from a quantitative MRI in native qMRI space. This is achieved by constructing a series of equivolumetric surfaces between the gray matter boundary and 3 mm below it following a laplacian field. By parcellating and cross-correlating nodal intensity profiles, this module generates microstructural profile covariance (MPC) matrices describing similiarity in superficial white matter microstructure across the brain.
Superficial white matter surfaces are generated using scripts from the superficial-white-matter repository, available via GitHub by Jordan De Kraker.
-MPC_SWM
Prerequisites 🖐🏼
You need to run -proc_structural
, -proc_surf
and -post_structural
before this stage
Compute registration from qMRI to surface native space.
Generate 16 equivolumetric surfaces between the gray matter boundary and 3 mm below it. Surfaces closest to gray matter boundary is then discarded to account for partial volume effects, resulting in 15 surfaces used for further analyses.
Perform surface-based registration to fsaverage5, fsLR-32k, fsLR-5k templates
Average intensity profiles within parcels defined on the native surface, while excluding outlier vertices
Perform partial correlation, controlling for parcel-wide mean profile, across all pairs of intensity profiles
Terminal:
$ micapipe -sub <subject_id> -out <outputDirectory> -bids <BIDS-directory> -MPC_SWM <options>
Optional arguments:
-MPC_SWM
has several optional arguments:
Optional argument |
Description |
---|---|
|
Specifies an input quantitative MRI on which to sample intensities for the MPC analysis. By default the pipeline will search this regex anat/*mp2rage*T1map.nii* You must specify this flag with the full path to your qMRI (for example, MTR, MTSAT, T2star, T1w/T2w). |
|
Path to scan which will be register to the surface ( e.g INV1, better anatomical contrast ). It MUST be on the same space as the main microstructural image!!.If it is EMPTY will try to find a T1map from here: anat/mp2rage*T1map.nii. Set to ‘FALSE’ to use microstructural_img for registrations. |
|
Default=qMRI. Provide a string with this this flag to process new quantitative map. ( this will create a new directory here: mpc-swm/acq-<qMRI> ) |
|
Specify this option to perform the registration based on synthseg |
|
Specify this option to perform the NON-Linear registration (e.g. MRI with strong geometric distortions). |
Directories created or populated by -MPC_SWM:
- <outputDirectory>/micapipe_v0.2.0/<sub>/mpc-swm/acq_<qMRI>
- <outputDirectory>/micapipe_v0.2.0/<sub>/anat
- <outputDirectory>/micapipe_v0.2.0/<sub>/xfm
Files generated by -MPC_SWM:
- Microstructural image in native Surfer space: <outputDirectory>/micapipe_v0.2.0/<sub>/anat/<sub>_space-fsnative_<qMRI>.nii.gz - Registration outputs from microstructural image to native FreeSurfer space: <outputDirectory>/micapipe_v0.2.0/<sub>/xfm/ <sub>_from-T1map_to-nativepro_0GenericAffine.mat <sub>_from-T1map_to-nativepro_1Warp.nii.gz (optional) <sub>_from-T1map_to-nativepro_1InverseWarp.nii.gz (optional) - Equivolumetric surface sampling ouputs stored in <outputDirectory>/micapipe_v0.2.0/mpc-swm/<qMRI>: - MPC-SWM json card: <sub>_MPC-SWM-<qMRI>.json - Native intensity profiles and MPC: <sub>_surf-fsnative_desc-intensity_profiles.shape.gii - fsaverage5 intensity profiles: <sub>_surf-fsaverage5_desc-intensity_profiles.shape.gii - fsLR-32k intensity profiles: <sub>_surf-fsLR-32k_desc-intensity_profiles.shape.gii - fsLR-5k intensity profiles and MPC: <sub>_surf-fsLR-5k_desc-intensity_profiles.shape.gii <sub>_surf-fsLR-5k_MPC.shape.gii - Parcellated intensity profiles: <sub>_atlas-<atlas>_desc-intensity_profiles.shape.gii - MPC matrices: <sub>_atlas-<atlas>_desc-MPC.shape.gii