I’m trying to incorporate T2 image along with a T1 to generate 5tt. I am using FLIRT / transformconvert / mrtransform to register images. FLIRT -in T1 -ref B0 2 followed by transformconvert and mrtransform (inverse) to get coregistered high resolution T1. Now, I’ve tried to bring my T2 “to condition” by registering it to: 1.initial T1 and 2. T1-coreg (which is a result form step 1). In both cases FLIRT is followed by transformconvert / mrtransform. Results of both visually are very similar. Neither result 1 or 2 gets through giving error: Provided T2 image does not match input T1 image. Please advise!
The way that the
5ttgen fsl script incorporates a T2 image into the segmentation process is by passing both to FSL’s
fast command, where both image modalities are used within the hMRF model. This requires that both images be defined on the same image grid. Therefore our usual advice in the context of registration of altering just the image header transformation matrix and specifically not re-gridding the image data does not apply in this instance: the T2-weighted image must be explicitly re-sampled to the voxel grid of the T1-weighted image.
Thank you, Rob! By now I’ve completed 5tt generation without T2, but will try in a few days with the T2 - curious to see the difference.