@ltetrel and @yohanchatelain have hit a common issue with fMRIPrep’s and sMRIPrep’s default outputs. Sometimes, minimal variations in the original image’s resolution would end up with slightly different NIfTI orientation prescriptions and matrix, making direct aggregation of results in some standardized space impossible.
The default prescription replicates the resolution of the original image (T1w or BOLD) that is going to be resampled into standard space, and then samples the FoV of the target template with that resolution. This is useful when images have large voxels and it feels unnecessary to impose, e.g., 2mm cubic; or very high-res images that should not be rounded to 1mm or 2mm cubic.
To avoid this, one can specify an existing resolution of the template with, e.g.,
--output-spaces MNI152NLin6Asym:res-2. One side issue is that we don’t know if this works for sMRIPrep outputs ATM.
What would be the most adequate way of coming up with an efficient new convention for this?
Some other considerations to note include that for BOLD datasets it would make sense to try to reduce the FoV of the outputs to a bounding box that contains only the brain (plus some safety padding). This would reduce the size of the output data, esp. for very long datasets.