Hi Saad, Thanks for talking a look! Line 59 is modified from the original script. It basically allowed you to run the script from outside the directory by splitting the string at the "/". Thats how I had my scripts set up.. For example: rotate_bvecs2 testdir/test.ecclog testdir/test.bvec; To do it within a directory, line 59 does need to be altered to: matfile=$(remove_ext $(echo ${line} | grep processing | awk '{print $2}')); For example: rotate_bvecs2 test.ecclog test.bvec; I have a attached a zip file with: 1. the corrected rotate_bvecs2 script 2. a text file with sample output (what the terminal prints when all goes well), as well as the command syntax 3. a working bvec file 4. a working ecclog file 5. my corrected "test_rot.bvec" file Please let me know if this works on your end. Thanks for catching this error! Best, Dave Deriso On Thu, Mar 11, 2010 at 1:48 AM, Saad Jbabdi <[log in to unmask]> wrote: > Hi Dave, > Thanks for your suggestions. > I agree with you that it is best to name the output <blah>_rot > However, your script does not work - in particular I am wondering about > line 59 - what is that for? That line is causing the script to fail from my > side. > > So I will make the relevant changes re the ouput name, but will keep the > rest as it is... > > Cheers, and thanks again. > Saad > > > On 9 Mar 2010, at 08:54, Dave Deriso wrote: > > Hello, > > I'm new here. My name is Dave Deriso and I am an undergrad researcher at > UCSD Psychiatry/Psychology under Dr. Martin Paulus and Dr. Ramachandran. I > was following the DTI pipeline on the FSL site and everything was going well > until I tried to run rotate_bvecs. For some reason, it just wouldn't work. I > doubt my OS has anything to do with it, but I am using a mac air laptop > running OSX Snow Leopard (10.6.2). > > Anyway, I decided to go ahead and rewrite the problematic portions. The new > script (attached) runs without problems. In addition, I added: > > - safety measure: make sure you are running an ecclog file instead of the > ecc nifti. > - changed the output to "myBvec*_rot.*bvec", this basically adds a "_rot" > to your existing file name. The original version simply made a backup copy > of your original bvec file, and named it bvec.bvec_old. bvec_rot is better > because it indicates that this file has been transformed.. > > Please have a look and make sure I didn't make an error somewhere! > > Best, > > Dave Deriso > UCSD Psychiatry > <rotate_bvecs2> > > > -- > Saad Jbabdi > University of Oxford, FMRIB Centre > > JR Hospital, Headington, OX3 9DU, UK > (+44)1865-222466 (fax 717) > www.fmrib.ox.ac.uk/~saad <http://www.fmrib.ox.ac.uk/%7Esaad> > > > > > > > > > > >