Dear FSL users,
We are trying to use FSL eddy correction to fix few DTI files.
Following these commands:
topup --imain=A2P_P2A_b0 --datain=acqparams.txt --config=b02b0.cnf --out=my_topup_results --iout=my_hifi_b0
eddy/eddy --imain=data --mask=my_hifi_b0_brain_mask --acqp=acqparams.txt --index=index.txt --bvecs=bvecs --bvals=bvals --topup=my_topup_results --niter=1 --out=eddy_corrected_data -v
We got this error:
Segmentation violation, Address not mapped, Offending address = 0x7c00000050
/lib64/libc.so.6 ) [0x3e6f73356f] [
eddy/eddy RBD_COMMON::Tracer::AddTrace() [0x61434a]
eddy/eddy NEWMAT::Matrix::operator()(int, int) [0x623be5]
eddy/eddy EDDY::ECScanManager::offset_design_matrix() const [0x4827cd]
eddy/eddy EDDY::ECScanManager::SeparateFieldOffsetFromMovement() [0x488828]
eddy/eddy main [0x470c68]
/lib64/libc.so.6 __libc_start_main [0x3e6f61ed5d]
Segmentation fault (core dumped)
Based on previous posts in this archive - the b02b0.cnf file is ok.
removing the downsampling from the config file also doesn't help.
Some basic debug showed us that the problem is here: (removing this line solves the problem).
// Separate field offset from subject movement in PE direction
if (clo.RegisterDWI()) sm.SeparateFieldOffsetFromMovement();
obviously, we don't want to skip this step.
can anyone advise us? please
Thanks!
|