Hello,
After running this code for eddy_correction.....
@ linenum = 1
foreach subj (`ls -d /Volumes/DANIEL/AES13* | sed 's/\// /g' | awk '{print $3}'`)
foreach cond (Ex Rest)
echo "subj = " ${subj} ", cond = " ${cond}
cd /Volumes/DANIEL/${subj}/${subj}.${cond}
#eddy correct function by fsl (on other computers probably needs to be eddy_openmp)
#input origional dti dataset, binary mask that has been automatically and manually edited, acqp (text file/get from mri tech), index (text file), bvecs, and bvals.
#replace outliers (--repol), output as eddy_corrected_data
eddy --imain=DTI_analysis.nii.gz --mask=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/${subj}/${subj}.${cond}/field_map/FieldMap2dti* --out=eddy_corrected_data
end
end
All my subjects were processed except 1. For one of my subjects I received the error.....
Segmentation violation, Address not mapped, Offending address = 0x510477000
eddy 0x0000000106dd7733 SPLINTERPOLATOR::Splinterpolator<float>::value_at(double const*) const
eddy 0x0000000106d4780c SPLINTERPOLATOR::Splinterpolator<float>::operator()(double, double, double, double) const
eddy 0x0000000106d46fe5 NEWIMAGE::volume<float>::splineinterpolate(float, float, float) const
eddy 0x0000000106d50e49 NEWIMAGE::volume<float>::interpolate(float, float, float) const
eddy 0x0000000106ccd74a 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 0x0000000106cce511 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 0x0000000106cc5931 void NEWIMAGE::affine_transform<float>(NEWIMAGE::volume<float> const&, NEWMAT::Matrix const&, NEWIMAGE::volume<float>&, NEWIMAGE::volume<char>&)
eddy 0x0000000106cb8811 EDDY::ECScan::field_for_model_to_scan_transform(std::__1::shared_ptr<NEWIMAGE::volume<float> const>, NEWIMAGE::volume<float>*, NEWIMAGE::volume<float>*) const
eddy 0x0000000106cd2e9c 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 0x0000000106cc552b EDDY::EddyUtils::TransformModelToScanSpace(NEWIMAGE::volume<float> const&, EDDY::ECScan const&, std::__1::shared_ptr<NEWIMAGE::volume<float> const>)
eddy 0x0000000106cd4ac2 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 0x0000000106c55f58 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 0x0000000106c50380 EDDY::DoVolumeToVolumeRegistration(EDDY::EddyCommandLineOptions const&, EDDY::ECScanManager&)
eddy 0x0000000106c4a248 main
libdyld.dylib 0x00007fffdbd02235 start
??? 0x000000000000000a 0x0
Any thoughts on what the issue is and how I might go about resolving it?
Thank you!
########################################################################
To unsubscribe from the FSL list, click the following link:
https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=FSL&A=1
|