External Email - Use Caution        

Dear All,

I am encountering an issue while running Tracula from the dmrirc file. An
additional line was required to run smoothly, without errors. I added
" set usethalnuc = 0 " and it proceeded without errors, yet when " set
usethalnuc = 1 " was added I encountered an error message.

I might have missed some relevant information on the page, but could you
please clarify what this setting does? Additionally, which preprocessing
step might I have missed if Tracula only works with set usethalnuc = 0?

Thank you for your help in advance!

Kind regards,
Diana

PS: I am also enclosing the test dmrirc file that I was using.
# FreeSurfer SUBJECTS_DIR
# T1 images and FreeSurfer segmentations are expected to be found here
# 

setenv SUBJECTS_DIR /Volumes/Diana/Preprocessed/DTI/Test_script/ADDI_013


# Output directory where trac-all results will be saved
# Default: Same as SUBJECTS_DIR
#
set dtroot = /Volumes/Diana/Preprocessed/DTI/Test_script/ADDI_013

# Subject IDs
#
set subjlist = (ADDI_013_recon_all)

set dcmlist = (data.nii.gz)

# In case you want to analyze only Huey and Louie
# Default: Run analysis on all subjects
#
set runlist = (1)

# Input diffusion DICOMs (file names relative to dcmroot)
# If original DICOMs don't exist, these can be in other image format
# but then the gradient table and b-value table must be specified (see below)
#
set dcmroot = /Volumes/Diana/Preprocessed/DTI/Test_script/ADDI_013

# Diffusion gradient tables (file names can be relative to dcmroot)
# Must be specified only if they cannot be read from the DICOM headers
# The tables must have either three columns, where each row is a gradient vector
# or three rows, where each column is a gradient vector
# There must be as many gradient vectors as volumes in the diffusion data set
# Default: Read from DICOM header
#
set bvecfile = (bvecs)

# Diffusion b-value tables (file names can be relative to dcmroot)
# Must be specified only if they cannot be read from the DICOM headers
# There must be as many b-values as volumes in the diffusion data set
# Default: Read from DICOM header
#
set bvalfile = (bvals)

# Perform registration-based B0-inhomogeneity compensation?
# Default: 0 (no)
#
set dob0 = 0

# Input B0 field map magnitude DICOMs (file names relative to dcmroot)
# Only used if dob0 = 1
# Default: None
#
# set b0mlist = (huey/fmag/XXX-1.dcm dewey/fmag/XXX-1.dcm louie/fmag/XXX-1.dcm)

# Input B0 field map phase DICOMs (file names relative to dcmroot)
# Only used if dob0 = 1
# Default: None
#
# set b0plist = (huey/fphas/XXX-1.dcm \
#               dewey/fphas/XXX-1.dcm \
#               louie/fphas/XXX-1.dcm)
#

# Field mapping TE difference (this is found in the scanner protocol printout)
# Only used if dob0 = 1
# Default: Read from field mapping phase DICOM
#
# set dTE = 0.0025

# Echo spacing for field mapping sequence (from sequence printout)
# Only used if dob0 = 1
# Default: None
#
# set echospacing = 0.7

# Perform registration-based eddy-current compensation?
# Default: 1 (yes)
#
set doeddy = 1

# Rotate diffusion gradient vectors to match eddy-current compensation?
# Only used if doeddy = 1
# Default: 1 (yes)
#
# set dorotbvecs = 1

# Fractional intensity threshold for BET mask extraction from low-b images
# This mask is used only if usemaskanat = 0
# Default: 0.3
#
# set thrbet = 0.5

# Perform diffusion-to-T1 registration by flirt?
# Default: 0 (no)
#
set doregflt = 1

# Perform diffusion-to-T1 registration by bbregister?
# Default: 1 (yes)
#
set doregbbr = 1

# Perform registration of T1 to MNI template?
# Default: 1 (yes)
#
# set doregmni = 1

# MNI template
# Only used if doregmni = 1
# Default: $FSLDIR/data/standard/MNI152_T1_1mm_brain.nii.gz
#
# set mnitemp = /path/to/mni_template.nii.gz

# Perform registration of T1 to CVS template?
# Default: 0 (no)
#
# set doregcvs = 0

# CVS template subject ID
# Only used if doregcvs = 1
# Default: cvs_avg35
#
# set cvstemp = donald

# Parent directory of the CVS template subject
# Only used if doregcvs = 1
# Default: $FREESURFER_HOME/subjects
#
# set cvstempdir = /path/to/cvs/atlases/of/ducks

# Use brain mask extracted from T1 image instead of low-b diffusion image?
# Has no effect if there is no T1 data
# Default: 1 (yes)
#
# set usemaskanat = 1

# Paths to reconstruct
# Default: All paths in the atlas
#
# set pathlist = ( lh.cst_AS rh.cst_AS \
#                 lh.unc_AS rh.unc_AS \
#                 lh.ilf_AS rh.ilf_AS \
#                 fmajor_PP fminor_PP \
#                 lh.atr_PP rh.atr_PP \
#                 lh.ccg_PP rh.ccg_PP \
#                 lh.cab_PP rh.cab_PP \
#                 lh.slfp_PP rh.slfp_PP \
#                 lh.slft_PP rh.slft_PP )

# Number of path control points
# It can be a single number for all paths or a different number for each of the
# paths specified in pathlist
# Default: 7 for the forceps major, 6 for the corticospinal tract,
#          4 for the angular bundle, and 5 for all other paths
#
# set ncpts = (6 6 5 5 5 5 7 5 5 5 5 5 4 4 5 5 5 5)

# List of training subjects
# This text file lists the locations of training subject directories
# Default: $FREESURFER_HOME/trctrain/trainlist.txt
#
# set trainfile = $FREESURFER_HOME/trctrain/trainlist.txt

# Number of "sticks" (anisotropic diffusion compartments) in the bedpostx
# ball-and-stick model
# Default: 2
#
# set nstick = 2

# Number of MCMC burn-in iterations
# (Path samples drawn initially by MCMC algorithm and discarded)
# Default: 200
#
# set nburnin = 200

# Number of MCMC iterations
# (Path samples drawn by MCMC algorithm and used to estimate path distribution)
# Default: 7500
#
# set nsample = 7500

# Frequency with which MCMC path samples are retained for path distribution
# Default: 5 (keep every 5th sample)
#
# set nkeep = 5

# Reinitialize path reconstruction?
# This is an option of last resort, to be used only if one of the reconstructed
# pathway distributions looks like a single curve. This is a sign that the
# initial guess for the pathway was problematic, perhaps due to poor alignment
# between the individual and the atlas. Setting the reinit parameter to 1 and
# rerunning "trac-all -prior" and "trac-all -path", only for the specific
# subjects and pathways that had this problem, will attempt to reconstruct them
# with a different initial guess.
# Default: 0 (do not reinitialize)
#
# set reinit = 0

# set the direction of the scanning 
set pedir = "y"  
# or "-y", "x", "-x", "z", or "-z" depending on your data

# it was advised to use thalamic segmentation, but I have no clue how to do it 
set usethalnuc = 0

_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is 
addressed.  If you believe this e-mail was sent to you in error and the e-mail 
contains patient information, please contact the Mass General Brigham 
Compliance HelpLine at https://www.massgeneralbrigham.org/complianceline 
<https://www.massgeneralbrigham.org/complianceline> .
Please note that this e-mail is not secure (encrypted).  If you do not wish to 
continue communication over unencrypted e-mail, please notify the sender of 
this message immediately.  Continuing to send or respond to e-mail after 
receiving this message means you understand and accept this risk and wish to 
continue to communicate over unencrypted e-mail. 

Reply via email to