Hello,
I am trying to trouble shoot my eddy correction script. The issue is that some of my scans are processed without issue however about 1/4 I receive an error for. I underwent similar preprocessing points up until this point for all subjects. Also, when running eddy without the field map, all subjects run fine (which leads me to believe it is an error with the field map input itself).
I am wondering if someone might be able to tell me what the error given for the following command is referring to. I think understanding where eddy is getting tripped up might help me better figure out the issue.
INPUT
eddy --imain=/Volumes/DANIEL/AES133/AES133.Ex/DTI_analysis.nii.gz --mask=/Volumes/DANIEL/AES133/AES133.Ex/COPY_b0_brain_0.4_mask_eroded.nii.gz --acqp=acqp.txt --index=index.txt --bvecs=bvec.txt --bvals=bvalue.txt --repol --field=/Volumes/DANIEL/AES133/AES133.Ex/field_map/FieldMap2dti --out=eddy_corrected_data
OUTPUT (occurring during the end of the command running)
Segmentation violation, Address not mapped, Offending address = 0x509758000
eddy 0x00000001001dc733 SPLINTERPOLATOR::Splinterpolator<float>::value_at(double const*) const
eddy 0x000000010014c80c SPLINTERPOLATOR::Splinterpolator<float>::operator()(double, double, double, double) const
eddy 0x000000010014bfe5 NEWIMAGE::volume<float>::splineinterpolate(float, float, float) const
eddy 0x0000000100155e49 NEWIMAGE::volume<float>::interpolate(float, float, float) const
eddy 0x00000001000d274a void NEWIMAGE::raw_general_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume4D<float> const&, std::__1::vector<int, std::__1::allocator<int> > const&, std::__1::vector<int, std::__1::allocator<int> > const&, std::__1::vector<unsigned int, std::__1::allocator<unsigned int> >, NEWMAT::Matrix const*, NEWMAT::Matrix const*, NEWIMAGE::volume<float>&, NEWIMAGE::volume4D<float>&, NEWIMAGE::volume<char>*)
eddy 0x00000001000d3511 void NEWIMAGE::raw_general_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume4D<float> const&, std::__1::vector<int, std::__1::allocator<int> > const&, std::__1::vector<int, std::__1::allocator<int> > const&, NEWIMAGE::volume<float>&, NEWIMAGE::volume4D<float>&, NEWIMAGE::volume<char>&)
eddy 0x00000001000ca931 void NEWIMAGE::affine_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume<float>&, NEWIMAGE::volume<char>&)
eddy 0x00000001000bd811 EDDY::ECScan::field_for_model_to_scan_transform(std::__1::shared_ptr<NEWIMAGE::volume<float> const>, NEWIMAGE::volume<float>*, NEWIMAGE::volume<float>*) const
eddy 0x00000001000d7e9c EDDY::EddyUtils::transform_model_to_scan_space(NEWIMAGE::volume<float> const&, EDDY::ECScan const&, std::__1::shared_ptr<NEWIMAGE::volume<float> const>, bool, NEWIMAGE::volume<float>&, NEWIMAGE::volume<float>*, NEWIMAGE::volume4D<float>*)
eddy 0x00000001000ca52b EDDY::EddyUtils::TransformModelToScanSpace(NEWIMAGE::volume<float> const&, EDDY::ECScan const&, std::__1::shared_ptr<NEWIMAGE::volume<float> const>)
eddy 0x00000001000d9ac2 EDDY::EddyUtils::param_update(NEWIMAGE::volume<float> const&, std::__1::shared_ptr<NEWIMAGE::volume<float> const>, NEWIMAGE::volume<float> const&, EDDY::Parameters, bool, float, unsigned int, unsigned int, unsigned int, EDDY::ECScan&, NEWMAT::ColumnVector*)
eddy 0x000000010005af58 EDDY::Register(EDDY::EddyCommandLineOptions const&, EDDY::ScanType, unsigned int, std::__1::vector<float, std::__1::allocator<float> > const&, EDDY::SecondLevelECModel, bool, EDDY::ECScanManager&, EDDY::ReplacementManager*, NEWMAT::Matrix&, NEWMAT::Matrix&)
eddy 0x0000000100055380 EDDY::DoVolumeToVolumeRegistration(EDDY::EddyCommandLineOptions const&, EDDY::ECScanManager&)
eddy 0x000000010004f248 main
libdyld.dylib 0x00007fff97589235 start
end
########################################################################
To unsubscribe from the FSL list, click the following link:
https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=FSL&A=1
|