Hi Michael <div><br></div><div>Thank you for the quick reply and the advice on an adequate shell for both FS and FSL.</div><div><br></div><div>The commands I had difficulty with stem from the use of FS5.1.0 and the trac-all - prep command that calls a number of FSL routines. The configuration file and the routines that it calls seem to be written for csh/tcsh. </div>
<div><br></div><div>I have encountered errors, and still do even with fslmachtype.sh and the csh apparently working correctly, that imply that my FSL distro is not sourced properly. Here is an example of the output from trac-all -prep -c when it faults before even processing one subject (thank you in advance for any advice):</div>
<div><br></div>
<div>trac-all -prep -c /media/5424CE1A24CDFECC/UofT/data/Tracula/tracula_config.txt <br>INFO: SUBJECTS_DIR is /usr/local/freesurfer/subjects<br>INFO: Diffusion root is /usr/local/freesurfer/subjects<br>Actual FREESURFER_HOME /usr/local/freesurfer<br>
flirt: error while loading shared libraries: libnewimage.so: cannot open shared object file: No such file or directory<br>trac-preproc -c /usr/local/freesurfer/subjects/C001/scripts/dmrirc.local -log /usr/local/freesurfer/subjects/C001/scripts/trac-all.log -cmd /usr/local/freesurfer/subjects/C001/scripts/trac-all.cmd<br>
#-------------------------------------<br>/usr/local/freesurfer/bin/trac-preproc <br>#-------------------------------------<br>#@# Image corrections Wed Jun 8 16:42:47 EDT 2011<br>mri_convert /media/5424CE1A24CDFECC/UofT/data/Tracula/C001_merged.nii.gz /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.nii.gz<br>
mri_convert /media/5424CE1A24CDFECC/UofT/data/Tracula/C001_merged.nii.gz /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.nii.gz <br>$Id: mri_convert.c,v 1.179.2.2 2011/05/16 20:53:47 greve Exp $<br>reading from /media/5424CE1A24CDFECC/UofT/data/Tracula/C001_merged.nii.gz...<br>
TR=17000.00, TE=0.00, TI=0.00, flip angle=0.00<br>i_ras = (-1, 0, 0)<br>j_ras = (0, 1, 0)<br>k_ras = (0, 0, 1)<br>writing to /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.nii.gz...<br>cp /media/5424CE1A24CDFECC/UofT/data/Tracula/bvecs.txt /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.mghdti.bvecs<br>
cp /media/5424CE1A24CDFECC/UofT/data/Tracula/bvals.txt /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.mghdti.bvals<br>flip4fsl /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.nii.gz /usr/local/freesurfer/subjects/C001/dmri/dwi_orig_flip.nii.gz<br>
INFO: input image orientation is LAS<br>INFO: input image determinant is -10<br>fslswapdim /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.nii.gz x y z /usr/local/freesurfer/subjects/C001/dmri/dwi_orig_flip.nii.gz<br>fslorient -forceradiological /usr/local/freesurfer/subjects/C001/dmri/dwi_orig_flip.nii.gz<br>
fslorient: error while loading shared libraries: libnewimage.so: cannot open shared object file: No such file or directory<br>INFO: found /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.mghdti.bvals, converting to FSL format<br>
INFO: found /usr/local/freesurfer/subjects/C001/dmri/dwi_orig.mghdti.bvecs, converting to FSL format<br>mv -f /usr/local/freesurfer/subjects/C001/dmri/dwi_orig_flip.mghdti.bvecs /usr/local/freesurfer/subjects/C001/dmri/bvecs<br>
mv -f /usr/local/freesurfer/subjects/C001/dmri/dwi_orig_flip.mghdti.bvals /usr/local/freesurfer/subjects/C001/dmri/bvals<br>eddy_correct /usr/local/freesurfer/subjects/C001/dmri/dwi_orig_flip.nii.gz /usr/local/freesurfer/subjects/C001/dmri/dwi.nii.gz 0<br>
Input does not exist or is not in a supported format<br>mv -f /usr/local/freesurfer/subjects/C001/dmri/bvecs /usr/local/freesurfer/subjects/C001/dmri/bvecs.norot<br>xfmrot /usr/local/freesurfer/subjects/C001/dmri/dwi.ecclog /usr/local/freesurfer/subjects/C001/dmri/bvecs.norot /usr/local/freesurfer/subjects/C001/dmri/bvecs<br>
grep: /usr/local/freesurfer/subjects/C001/dmri/dwi.ecclog: No such file or directory<br>wc: /usr/local/freesurfer/subjects/C001/dmri/dwi.ecclog: No such file or directory<br>ERROR: Transform file should be an eddy_correct log file or a .mat file<br>
Linux dsb-P35-DS3L 2.6.35-28-generic #50-Ubuntu SMP Fri Mar 18 18:42:20 UTC 2011 x86_64 GNU/Linux<br><br>trac-preproc exited with ERRORS at Wed Jun 8 16:43:10 EDT 2011<br><br>Here are the contents of the FS subject folder after it faults. It is as thought the dwi.ecclog was never written:</div>
<div><br>bvals bvecs.norot dwi_orig.mghdti.bvals dwi_orig.mghdti.bvecs dwi_orig.nii.gz xfms<br><br>Here is my CONFIG FILE:</div><div><br># In case you want to analyze only Huey and Louie<br># Default: Run analysis on all subjects<br>
#<br>set runlist = (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22)<br><br># Input diffusion DICOMs (file names relative to dcmroot)<br># If original DICOMs don't exist, these can be in other image format<br>
# but then bvecfile, bvalfile, and nb0 must be specified (see below)<br>#<br>set dcmroot = /media/5424CE1A24CDFECC/UofT/data/Tracula<br>set dcmlist = (C001_merged.nii.gz \ C002_merged.nii.gz \ C003_merged.nii.gz C004_merged.nii.gz \ C005_merged.nii.gz \ C006_merged.nii.gz \ C007_merged.nii.gz \ C008_merged.nii.gz \ C009_merged.nii.gz \ C010_merged.nii.gz \ C011_merged.nii.gz \ S001_merged.nii.gz \ S002_merged.nii.gz \ S003_merged.nii.gz \ S004_merged.nii.gz \ S005_merged.nii.gz \ S006_merged.nii.gz \ S007_merged.nii.gz \ S008_merged.nii.gz \ S009_merged.nii.gz \ S010_merged.nii.gz \ S011_merged.nii.gz)<br>
<br># Diffusion gradient table<br># Must be specified if inputs are not MGH DICOMs<br># Three-column format, one row for each volume in the diffusion data set<br># Default: Read from DICOM header<br>#<br>set bvecfile = /media/5424CE1A24CDFECC/UofT/data/Tracula/bvecs.txt<br>
<br># Diffusion b-value table<br># Must be specified if inputs are not MGH DICOMs<br># Single-column format, one value for each volume in the diffusion data set<br># Default: Read from DICOM header<br>#<br>set bvalfile = /media/5424CE1A24CDFECC/UofT/data/Tracula/bvals.txt<br>
<br># Number of low-b images<br># Must be specified if inputs are not DICOM<br># Default: Read from DICOM header<br>#<br>set nb0 = 4<br><br># Perform registration-based B0-inhomogeneity compensation?<br># Default: 0 (no)<br>
#<br>set dob0 = 0<br><br># Input B0 field map magnitude DICOMs (file names relative to dcmroot)<br># Only used if dob0 = 1<br># Default: None<br>#<br># set b0mlist = (huey/fmag/XXX-1.dcm dewey/fmag/XXX-1.dcm louie/fmag/XXX-1.dcm)<br>
<br># Input B0 field map phase DICOMs (file names relative to dcmroot)<br># Only used if dob0 = 1<br># Default: None<br>#<br># set b0plist = (huey/fphas/XXX-1.dcm dewey/fphas/XXX-1.dcm louie/fphas/XXX-1.dcm)<br><br># Echo spacing for field mapping sequence (from sequence printout)<br>
# Only used if dob0 = 1<br># Default: None<br>#<br># set echospacing = 0.7<br><br># Perform registration-based eddy-current compensation?<br># Default: 1 (yes)<br>#<br>set doeddy = 1<br><br># Rotate diffusion gradient vectors to match eddy-current compensation?<br>
# Only used if doeddy = 1<br># Default: 1 (yes)<br>#<br>set dorotbvecs = 1<br><br># Fractional intensity threshold for BET mask extraction from low-b images<br># This mask is used only if usemaskanat = 0<br># Default: 0.3<br>
#<br>set thrbet = 0.3<br><br># Perform diffusion-to-T1 registration by flirt?<br># Default: 1 (yes)<br>#<br>set doregflt = 0<br><br># Perform diffusion-to-T1 registration by bbregister?<br># Default: 0 (no)<br>#<br>set doregbbr = 1<br>
<br># MNI template (the only option for inter-subject registration in this version)<br># Default: $FSLDIR/data/standard/MNI152_T1_1mm_brain.nii.gz<br>#<br>set mnitemp = /usr/share/fsl/4.1/data/standard/MNI152_T1_1mm_brain.nii.gz<br>
<br># Use brain mask extracted from T1 image instead of low-b diffusion image?<br># Has no effect if there is no T1 data<br># Default: 1 (yes)<br>#<br>set usemaskanat = 1<br><br># Paths to reconstruct<br># Default: All paths in the atlas<br>
#<br>set pathlist = ( lh.cst_AS rh.cst_AS \<br> lh.unc_AS rh.unc_AS \<br> lh.ilf_AS rh.ilf_AS \<br> fmajor_PP fminor_PP \<br> lh.atr_PP rh.atr_PP \<br> lh.ccg_PP rh.ccg_PP \<br>
lh.cab_PP rh.cab_PP \<br> lh.slfp_PP rh.slfp_PP \<br> lh.slft_PP rh.slft_PP )<br><br># Number of path control points<br># Default: 5<br>#<br>set ncpts = 5<br><br># List of training subjects<br>
# This text file lists the locations of training subject directories<br># Default: $FREESURFER_HOME/trctrain/trainlist.txt<br>#<br>set trainfile = $FREESURFER_HOME/trctrain/trainlist.txt<br><br># Use long (more descriptive) directory hierarchy for saving path distributions?<br>
# By default, paths distributions are saved directly under $subjectname/dpath<br># Default: 0 (no)<br>#<br>set dopathsubdirs = 0<br><br># Number of MCMC burn-in iterations<br># (Path samples drawn initially by MCMC algorithm and discarded)<br>
# Default: 200<br>#<br>set nburnin = 200<br><br># Number of MCMC iterations<br># (Path samples drawn by MCMC algorithm and used to estimate path distribution)<br># Default: 5000<br>#<br>set nsample = 5000<br><br># Frequency with which MCMC path samples are retained for path distribution<br>
# Default: 5 (keep every 5th sample)<br>#<br>set nkeep = 5<br></div><div><br clear="all">Deryk S. Beal, Ph.D., CCC-SLP, S-LP(C), Reg. CASLPO<br>Speech-Language Pathologist<br>C.I.H.R. Post Doctoral Research Fellow<br>Department of Cognitive and Neural Systems Speech Laboratory<br>
Boston University, 677 Beacon Street, Boston, MA 02215 <br><a href="mailto:dsbeal@bu.edu" target="_blank">dsbeal@bu.edu</a><br><a href="http://blogs.bu.edu/dsbeal/about/" target="_blank">http://blogs.bu.edu/dsbeal/about/</a><br>
*************************************************************************<br>This email may contain confidential and/or privileged information for the<br>sole use of the intended recipient. Any review or distribution by others is<br>
strictly prohibited. If you have received this email in error, please<br>contact the sender and delete all copies. Opinions, conclusions or other<br>information expressed or contained in this email are not given or endorsed<br>
by the sender unless otherwise affirmed independently by the sender.<br>
<br><br><div class="gmail_quote">On Wed, Jun 8, 2011 at 10:22 PM, Michael Hanke <span dir="ltr"><<a href="mailto:mih@debian.org" target="_blank">mih@debian.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<div><br>
On Wed, Jun 08, 2011 at 04:42:03PM -0400, Deryk S. Beal, Ph.D. wrote:<br>
> Why is fslmachtype.sh not distributed with the Neurodebian version? One must<br>
> be able to run fsl in tcsh or csh for certain Freesurfer commands.<br>
<br>
</div>I wasn't aware that this is the case. Would you mind sharing what<br>
commands are affected?<br>
<br>
In any case fslmachtype.sh is not shipped, because it is not necessary.<br>
You can safely remove this line from fsl.csh. The Debian package doesn't<br>
provide out-of-the-box csh support, because (not just) in Debian csh is<br>
considered harmful:<br>
<br>
<a href="http://lintian.debian.org/tags/csh-considered-harmful.html" target="_blank">http://lintian.debian.org/tags/csh-considered-harmful.html</a><br>
<br>
However, if there is a problem with Freesurfer we need to work it out --<br>
it is just that we haven't seen problems so far and happily use FSL with<br>
ksh/bash/dash/zsh/...<br>
<br>
Michael<br>
<font color="#888888"><br>
--<br>
Michael Hanke<br>
<a href="http://mih.voxindeserto.de" target="_blank">http://mih.voxindeserto.de</a><br>
</font></blockquote></div><br></div>