Ok, I think that I found the problem.
What a mess, but let me explain, I now remember another detail that I did here.
I extracted the voxel-based time series with 3dmaskdump using the option "-no ijk" so that the resulting .1D files do NOT include the voxel coordinates.
Why? Because I wanted to run interpolation in Python, and it was easier for me to script it this way.
Subsequently, I also extracted the voxel coordinates from the atlas file that I created.
Now, comparing the voxel coordinates between the atlas and the 3dmaskdump extracted errts files shows me that the voxel coordinates are NOT the same.
Sorry! I feel like an idiot.
I thought that both voxel coordinates should match, since the preprocessed errts files used the MNI152 template.
I guess that is the core of the problem.
Edit: Paul, forget what I wrote above. If I extract all voxel coordinates from the preprocessed errts. file, of course I will get many voxels that do not match the atlas, since the errts file includes all voxels, whereas the atlas/mask/roi does not. Furthermore, 3dinfo showed that "Resample_Yeo7N_1000ribbon.nii" (my atlas file) an the errts file do match.I have no idea anymore where the problem comes from, I am fully lost now.
Edited 2 time(s). Last edit at 04/10/2023 11:43PM by Philipp.