Hi all,

I am processing resting state data using AFNI version 21.2.04 ('Nerva'). Prior to preprocessing, I aligned the centers of the anatomical and functional images to the MNI*2009_SSW template. I ran the centered anatomical images through SSwarper and recon-all and used the following proc.py code to process the functional data:

afni_proc.py                                                                                                            \
                 -subj_id ${sid} -script proc.${sid} -scr_overwrite -out_dir ${sid}.results                             \
                 -blocks despike tshift align tlrc volreg blur mask scale                                               \
                     regress                                                                                            \
                 -radial_correlate_blocks tcat volreg                                                                   \
                 -copy_anat ${ss_dir}/anatSS.${sid}.nii                                                                 \
                 -anat_has_skull no                                                                                     \
                 -anat_follower anat_w_skull anat  ${ss_dir}/anatU.${sid}.nii                                           \
                 -anat_follower_ROI aaseg anat ${fs_dir}/aparc.a2009s+aseg.nii.gz                                       \
                 -anat_follower_ROI aeseg epi ${fs_dir}/aparc.a2009s+aseg.nii.gz                                        \
                 -anat_follower_ROI FSvent epi ${fs_dir}/fs_ap_latvent.nii.gz                                           \
                 -anat_follower_ROI FSWe epi ${fs_dir}/fs_ap_wm.nii.gz                                                  \
                 -anat_follower_erode FSvent FSWe                                                                       \
                 -dsets ${raw_epi}                                                                                      \
                 -tcat_remove_first_trs 4                                                                               \
                 -tshift_opts_ts -tpattern ${tpattern}                                                                      \
                 -align_opts_aea -cost lpc+ZZ -giant_move -check_flip                                                   \
                 -tlrc_base MNI152_2009_template_SSW.nii.gz                                                             \
                 -tlrc_NL_warp                                                                                          \
                 -tlrc_NL_warped_dsets ${ss_dir}/anatQQ.${sid}.nii ${ss_dir}/anatQQ.${sid}.aff12.1D                     \
                     ${ss_dir}/anatQQ.${sid}_WARP.nii                                                                   \
                 -volreg_align_to MIN_OUTLIER                                                                           \
                 -volreg_align_e2a                                                                                      \
                 -volreg_tlrc_warp                                                                                      \
                 -blur_size 4                                                                                           \
                 -mask_epi_anat yes                                                                                     \
                 -regress_motion_per_run                                                                                \
                 -regress_ROI_PC FSvent 3                                                                               \
                 -regress_ROI_PC_per_run FSvent                                                                         \
                 -regress_make_corr_vols aeseg FSvent                                                                   \
                 -regress_anaticor_fast                                                                                 \
                 -regress_anaticor_label FSWe                                                                           \
                 -regress_censor_motion 0.3										\
 	      	 -regress_censor_outliers 0.1										\
                 -regress_apply_mot_types demean deriv                                                                  \
                 -regress_est_blur_epits                                                                                \
                 -regress_est_blur_errts                                                                                \
                 -html_review_style pythonic

I am having two issues with alignment:
1. These data are from older participants (60+ years). As you can see in the attached screenshot (csf.png) showing the anat_final edges overlayed on the final_epi_vr_base_min_outlier, there is quite a bit of bright data (CSF) outside of the edges of the anat (e.g., 71S, 21L, 19R). The edges of the gray matter of the anat seem to align well with the edges of the gray matter in the EPI. I wanted to check in to see if you had comments on these bright distortions and if this is problematic for analysis (e.g. extracting mean timeseries from 90 ROIs in MNI space to form functional connectivity matrix).
2. Some images in the dataset are having issues with alignment of the cerebellum/brainstem between the anat and EPI. For example, in slice 1L, you can see that the brain stem is slightly right shifted compared to the anat_final edge overlay. I have attempted to run different cost functions (e.g., lpa, nmi) and have changed or removed the move option (e.g., big_move) but this has not alleviated this issue. Do you have any suggestions on how to adjust align_epi_anat to fix this issue?

Thank you!
Jenna