Overview of the intermediate steps of the image data processing pipeline: each whole brain T1 3D datasets was first bias-corrected (1) and, using CAT12, segmented into GM (2) and WM (not shown). Then, an iterative non-linear registration (3–8) to the respective group mean tissue map was applied, resulting in ever crisper tissue maps (upper row) and corresponding deformation fields (lower row, illustrated here by their Jacobian determinants). This results in conventional DARTEL templates (see Fig. 2, Fig. 4, top rows). The tissue maps were also submitted to the CerebroMatic toolbox, resulting in synthetic DARTEL templates (see Fig. 2, Fig. 4, bottom rows).