Hi! I am sometimes getting bad func/anatomical registration which results in low snr on parts of the brain where I wouldn’t expect to get low snr (like the lateral occipital lobe). Here is the fMRIprep command I used. singularity run --cleanenv -B /projects/:/projects/\ -B ${WORK_DIR}:/quest_…

The gray/white contrast in your BOLD is extremely low, so boundary-based registration has little to work with. If you can find some other way to register it and update the header, you can use --bold2t1w-init header --force-no-bbr to tell fMRIPrep not to attempt any registration.

thank you @effigies I do have single band reference files for each echo which all seem to have at least decent gray/white contrast (much better than what is shown in the registration screen shot i included in the last post). Does fmriprep use them? From the documentation it looks like it should bu…

I believe so but can’t check right now. Multi echo has definitely changed since 21.0. Can you try with 23.1?

Hi @effigies I tried it with 23.0.2 on a subset of sessions since we had that on our system already. Registration looked a little better, but the contrast between gray/white matter in boldref images doesn’t look any better. What do you think? How is it making these boldref images? The contrast is…

Here is an example of the SBref (echo 3) versus the boldref image for the same run [image] [image]

Hi @effigies , Just checking. Is your recommendation to register it elsewhere and then do fmriprep with --bold2t1w-init header --force-no-bbr?

Thanks for your patience. Yes, I think that is probably best. I looked at the code again, and we are using the sbrefs for the initial bold reference, which is used for head motion correction, but not for the final, which is used for coregistration. This needs fixing in an updated 23.2 series. That i…

Interesting! On top of that, SBref images used as reference for HMC tend to bias the HMC parameters with the standard cost function used by mcflirt: [image] [mcflirt related behavour] FMRIPREP: Strange motion estimations in confounds_timeseries.tsv (lots of zeros) …

Wow, either I missed this or forgot about it. Have you found any issues in using mutualinfo for BOLD runs without sbrefs? It does seem that we should probably prefer just to motion correct to the first volume and register the first volume to the reference, regardless of whether we use mcflirt or 3d…

More details on this topic here: [mcflirt related behaviour] Wrong motion parameters estimates in confounds_timeseries.tsv (lot of zeros!) opened 08:12PM - 23 Jan 21 UTC (UTC) [julfou81] julfou81…

Poor func/anatomical registration with fmriprep (21.0.2) on multi-echo data set

Software Support

jsein July 12, 2023, 8:13pm 10

Interesting! On top of that, SBref images used as reference for HMC tend to bias the HMC parameters with the standard cost function used by mcflirt: