Hi Rick,

Thanks for your suggestion.

I have tried the method you suggested. But it seemed didn't work. Since we cannot have the final grid of the processed data until we finish processing, I just used the processed data I did before for other subjects as the master dataset. The used data were following example 11 which used FreeSurfer to segment tissues. And the voxel resolution of the processed data is 2.5 mm^3.

I then used this dataset as the mater to resample the template ventricle using the following code:
3dresample -master errts+tlrc -prefix template_ventricle -input aparc+aseg_REN_vent.nii.gz
(But actually the final processed dataset is exactly 2.5 mm^3, they should be matched with the former one using 3dresample -dxyz. Am I right?)

However, when running the script generated from proc.py, error occured as follows:
3dToutcount -automask -fraction -polort 3 -legendre pb00.koon_pre.r01.tcat
++ 3dToutcount: AFNI version=AFNI_17.1.09 (Jun 6 2017) [64-bit]
** Can't open dataset pb00.koon_pre.r01.tcat

I cannot come up with the reason of this error. Could you please help clarify? Thanks!

Tracy