I hope this email finds you well! My name is Sam Sievertsen, the research coordinator for Dr. Jennifer Forsyth at the University of Washington. We are utilizing the recon-all clinical processing stream from FreeSurfer 7.4.1 on Rocky Linux (within our institutions HPC environment) to process subjects with decreased T1 image quality, and for whom previously processing with the standard recon-all pipeline resulted in missed grey matter during segmentation. Per recommendations in the bug reporting page, I have attached an example recon-all clinical log from this processing in case it is helpful. 

Following the recon-all clinical processing of these subjects, we are using an existing (ENIGMA cortical QC) pipeline to take screenshots of the external cortical reconstruction / aparc+aseg segmentation overlaid on internal images. 

We are able to successfully generate external (reconstruction) cortical QC images using the recon-all clinical output, but encounter an error when attempting to generate the internal (T1 with aparc+aseg.mgz overlay) QC images. Specifically, there is a conditional check within the internal QC image generation script linked above that halts the process if the Nx*Ny*Nz dimensions in the T1 image and the Nx*Ny*Nz dimensions in aparc+aseg.mgz do not align, which is where our process is throwing the error. 

Manual examination of the recon-all clinical image output (e.g., "native.mgz", "synthSR.norm.mgz", and "aparc+aseg.mgz") suggests that the aforementioned dimensions are all different in those images. Conversely, the dimensions in the processed images and segmentations (e.g., "orig_nu.mgz", "aparc+aseg.mgz") from the standard recon-all pipeline appear to be uniform (256) for those same subjects we are processing in this cohort. To hopefully assist with clarity, I have attached the script we are using ("func_make_corticalpngs_ENIGMA_QC.m") to construct these internal QC images, and example output files from the same subject using the standard recon-all pipeline output ("internalQC_standard_recon_all.txt") and recon-all clinical pipeline output ("internalQC_recon_all_clinical.txt"), which show a discrepancy in the aforementioned dimensions. 

Following up on the above information, I am wondering if "SynthSeg: to obtain a volumetric segmentation and linear registration to Talairach space"on the recon-all clinical page about registration to talairach space implies that the recon-all clinical output should have uniform dimensions, or if I am misunderstanding that and there is anything that stands out from this description that we may be doing incorrectly/that other users have encountered as well? 

Would also be very appreciative of any other suggestions you may have for attempting to troubleshoot this issue, thank you very much in advance for your time! Happy to answer any resulting questions or follow up with additional information should it be helpful. 

All the best,