Hey Thijs,
Thank you so much for such a quick response and great fix, and actually checking the data all the way through the CSD! Such help is unparalleled elsewhere in my professional world, Thijs!
Yes, I had seen that warning about two volumes not being assigned to a shell and after googling around, read Rob’s post here and assumed that it would have be taken care of automatically, thinking they would be internally rejected by the processing. Thanks for also pinpointing where the crash actually originated from! But yeah, I think such care should be taken at the input level.
We actually have more volumes with b=350 but the reason there are only two here is because of the outlier based volume filtering from eddy_quad in FSL. I will be careful to check that there are enough volumes left at each shell with this filtering process and if there are too few (say less than six) to form a shell, we will completely remove that b-value data.
Thanks again for such an invaluable help, my friend, and hope all is well with you and your loved ones too!
Sincerely,
Nagesh