Better option: Use Normalize intensities of a collection of T1 MRIs + Final deformation applied on volume + Create average template instead of this process.
Compute MRI template
Compute MRI template