Dilating of brain mask and some other general queries

Dear @Michiko,

I am working with HCP pipelines, but I cannot give answer specific to HCP data since I am not working with them (probably the dwi and structural HCP data are already in register?).

But generally, you have to assure that the structural and diffusion data are in register. This, however, does not mean that their voxel size has to be identical. Registering and regridding (reslicing) are two distinct things. You can bring the data to common coordinate systems without reslicing by adjusting their vox2RAS matrix in image header.

Registering structural to dwi (without reslicing structural, only modifying vox2ras matrix) I would do as follows:

Extract b=0 images and averaging them:
dwiextract -bzero dwi.mif - | mrmath - mean -axis 3 dwi_b0_mean.nii.gz

Register dwi to structural:
bbregister --s $subject_id --mov dwi_b0_mean.nii.gz --init-coreg --lta dwi_b0_mean2fs.lta --dti

For 5ttgen I personally use norm.mgz, optionally with refined skullstrip by masking by PreFreeSurfer brain mask (maybe such masking is not necessary with HCP data)

mri_mask $SUBJECTS_DIR/$subject_id/mri/norm.mgz ${hcp_rootdir}/${subject_ID}/T1w/T1w_acpc_brain_mask.nii.gz norm_preFS_brainmasked.nii.gz

Getting subject’s norm.mgz in register to dwi (without reslicing)
mri_vol2vol --mov norm_preFS_brainmasked.nii.gz --targ dwi_b0_mean.nii.gz --lta-inv dwi_b0_mean2fs.lta --o norm2dwi.mgz --no-resample

A note: The structural PreFreeSurfer data are AC-PC aligned, if the same is not applied for dwi, then there is optional step to bring dwi to AC-PC alignment before running bbregister to avoid possible failure in registration. There should be acpc.mat file in the HCP data directory which can be used for that. But then the steps are more complex.

A note 2: I personally currently prefer (maybe in contrast to common recommendations, I would appreciate a comment from other experts on this!) to do opposite, i.e. register dwi to structural mainly for following reason:

The structural data processed by HCP pipeline are AC-PC aligned. This is advantageous for all steps which require registration to MNI since the registration have good initial conditions. This truthfully matters for FSL-first (as part of 5ttgen and labelsgmfix), where I found that changing initial head position can have profound effect on subcortical segmentation result.
I think that the only pitfall in (rigid body!) transforming dwi data after they have been preprocessed by dwipreproc, is to pay attention to rotate your .bvec as well. However, when using mif format with dwi data, you no longer have to worry about it since mrtransform does it seamlessly.
You could also do 5ttgen in AC-PC aligned space and transform to dwi space resulting 5tt image (and also FreeSurfer parcellation image after labelsgmfix) but this seems to me less convenient.

Antonin