Hi Jesper,
Thanks so much for the suggestion.
I re-ran the commands with --ol_ec=2 as you suggested, and they completed successfully.
I found a very brief explanation of this otherwise-hidden eddy option in an earlier email from you:
https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=fsl;12467485.1811
It sounds like this flag shouldn't be used routinely, but kept as an alternative if this problem occurs with datasets that have only a small number of volumes with the lowest b-value?
I'm surprised that you wrote "not being able to find a single volume in the b=500 shell that doesn’t have at least one outlier slice" - so, ideally, its better to have few outlier slices for this procedure? That seems non-intuitive - or are you expecting a certain level of motion? (At least, it might upset our "professional" subjects who pride themselves on not moving in the scanner!)
On a somewhat related note - I've also noticed that the cnr_maps file being created by eddy with s2v are ending up with NaNs in the file; while when the same data is processed without s2v correction but otherwise identical parameters, this is not the case.
I get NaNs if I use fslstats on the files, and eddy_quad is logging:
!!!Warning!!! NaNs detected in the CNR maps!!!
although it goes on to complete the QC report. The SNR and CNR maps in the QC report look almost identical whether s2v was used or not - so its not immediately obvious where the NaN voxels are, and whether they can be safely ignored.
Thanks,
Ross
------------------------------
Date: Mon, 9 Dec 2019 14:22:27 +0000
From: Jesper Andersson <[log in to unmask]>
Subject: Re: EDDY::: Unable to find volume with no outliers in shell 0 (was: Re: eddy: msg=ECScanManager::set_slice_to_vol_reference: ref index out of bounds)
Hi again Ross,
> I am getting a very similar error message to the last poster (Nikitas); ie:
>
> EDDY::: DoVolumeToVolumeRegistration: Unable to find volume with no outliers in shell 0 with b-value=49?
> EDDY::: Eddy failed with message EDDY::: eddy.cpp::: EDDY::ReplacementManager*
> EDDY::DoVolumeToVolumeRegistration(const EDDY::EddyCommandLineOptions&, EDDY::ECScanManager&): Exception
> thrown
>
> This is occurring when I run eddy_cuda with slice-to-volume correction; however when I run eddy_cuda without s2v
> correction, there are no problems.
>
> In this dataset, I also have a small number of volumes - 6 - with low-b-value ~ 500. Does this mess with the s2v process? > (No volumes with b ~ 50 - so not sure where it found that number!)
as I suspected this is an issue with not being able to find a single volume in the b=500 shell that doesn’t have at least one outlier slice. It seems that for some reason the final number it the error message is overwritten with a ?.
I can’t re-run your command here because I don’t have the "slspec_ABCD102.txt” file. Can I ask you to please re-run the s2v command with --ol_ec=2 and let me know what you get?
Jesper
########################################################################
To unsubscribe from the FSL list, click the following link:
https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=FSL&A=1
|