根据fmriprep的xfm.h5变换lesion不匹配MNI空间

根据以下代码转换lesion从个体空间到MNI空间,发现lesion并不能很好的匹配MNI空间的T1w

root_dir=/data/disk1/zhangyi/Prep_LGStroke/Left_stroke
cd $root_dir
for sub in `ls bids_result_left`
do \
echo $sub
prep_t1=$root_dir/bids_result_left/$sub/fMRIPrep/fmriprep/$sub/anat/${sub}_desc-preproc_T1w.nii.gz
matrix=$root_dir/bids_result_left/$sub/fMRIPrep/fmriprep/$sub/anat/${sub}_from-T1w_to-MNI152NLin2009cAsym_mode-image_xfm.h5
mni_t1=$root_dir/bids_result_left/$sub/fMRIPrep/fmriprep/$sub/anat/${sub}_space-MNI152NLin2009cAsym_desc-preproc_T1w.nii.gz
lesion=$root_dir/lesionmask_L/${sub}_ses-1_lesionmask.nii
mni_lesion=$root_dir/lesionmask_L/${sub}_ses-1_space-MNI152NLin2009cAsym_lesionmask.nii
antsApplyTransforms -d 3 \
-i $lesion \
-o $mni_lesion \
-t $matrix \
-r $mni_t1
done 

用的是fmriprep做的预处理,选取的是他输出的xfm.h5作为变换矩阵

你的lesion的图和preproc_T1w.nii.gz是在同一个空间吗?

解决办法参看https://github.com/PennLINC/qsiprep/issues/532