Hello everyone
I'm trying to use eddy after topup to correct susceptibility and eddy current distorsions on DTI images (with 32 directions)
topup is working very well but when i want to use eddy i get a weird error.
At first, i thought it was a memory problem, but i tested it with another computer with the same results.
my command line is :
eddy --imain=MAP-VS-287_DTI_32_dir_SENSEmcverter.nii --mask=b0_Pbrain_mask.nii.gz --index=index.txt --acqp=acqparams.txt --bvecs=MAP-VS-287_DTI_32_dir.bvec --bvals=MAP-VS-287_DTI_32_dir.bval --fwhm=0 --topup=my_topup_results --flm=quadratic --out=DTI_ec.nii
and i get this :
Running Register
Segmentation violation, Address not mapped, Offending address = 0x512ae0000
eddy 0x000000010010ad86 SPLINTERPOLATOR::Splinterpolator<float>::value_at(double const*) const
eddy 0x00000001000c77d0 NEWIMAGE::volume<float>::splineinterpolate(float, float, float) const
eddy 0x00000001000e5f20 NEWIMAGE::volume<float>::interpolate(float, float, float) const
eddy 0x000000010002a084 void NEWIMAGE::raw_general_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume4D<float> const&, std::vector<int, std::allocator<int> > const&, std::vector<int, std::allocator<int> > const&, NEWMAT::Matrix const*, NEWMAT::Matrix const*, NEWIMAGE::volume<float>&, NEWIMAGE::volume4D<float>&, NEWIMAGE::volume<char>*)
eddy 0x000000010002b3d4 void NEWIMAGE::affine_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume<float>&, NEWIMAGE::volume<char>&)
eddy 0x0000000100019269 EDDY::ECScan::field_for_model_to_scan_transform(boost::shared_ptr<NEWIMAGE::volume<float> const>, NEWIMAGE::volume<float>*, NEWIMAGE::volume<float>*) const
eddy 0x0000000100030f2c EDDY::EddyUtils::transform_model_to_scan_space(NEWIMAGE::volume<float> const&, EDDY::ECScan const&, boost::shared_ptr<NEWIMAGE::volume<float> const>, bool, NEWIMAGE::volume<float>&, NEWIMAGE::volume<float>*, NEWIMAGE::volume4D<float>*)
eddy 0x0000000100036044 EDDY::EddyUtils::param_update(NEWIMAGE::volume<float> const&, boost::shared_ptr<NEWIMAGE::volume<float> const>, NEWIMAGE::volume<float> const&, EDDY::Parameters, bool, EDDY::ECScan&, NEWMAT::ColumnVector*)
eddy 0x0000000100004cfc Register(EDDY::EddyCommandLineOptions const&, EDDY::ScanType, EDDY::ECScanManager&, unsigned int, NEWMAT::Matrix&, NEWMAT::Matrix&)
eddy 0x0000000100005e7d main
eddy 0x0000000100001ca8 start
??? 0x000000000000000b 0x0
The "offending Adress" is almost the same each time 0x512******
Please help me, i dont' know what to do to make it work
|