Freesurfer April 2014

freesurfer@nmr.mgh.harvard.edu

147 participants
229 discussions

Re: [Freesurfer] [tracula] error running trac-all -prep
by Borsodi Florian 20 Dec '17

20 Dec '17

Hi Anastasia, Yes, i am running this on ubuntu 12.04. Sincerely, Florian >>> Anastasia Yendiki 20.08.13 17.38 Uhr >>> Hi Florian - Are you running this on an ubuntu system by any chance? a.y On Tue, 20 Aug 2013, Borsodi Florian wrote: > Hi Anastasia, > > Thank you for your help. > I changed the commas in "/dmri/bvecs" to dots and reran trac-all. But the command "xfmrot" changed the dots back to > commas. So, I changed our system language from de_AUT.UTF-8 to en_US.UTF-8 and reran trac-all. After changing the > system language everything was genereted correctly and trac-all -prep ran without any errors. > I'm so glad. > > Then, I started the next step "trac-all -bedp" and got a syntax error "/usr/local/freesurfer/bin/bedpostx_mgh: 131: > /usr/local/freesurfer/bin/bedpostx_mgh: Syntax error: "(" unexpected. " I also attached the the logfile. I was reading > about this error in the mail archive. Would I have to run "bedpostx" outside of tracula ? > > Sincerely, > Florian > > >>> Anastasia Yendiki 19.08.13 17.54 Uhr >>> > > Hi Florian - Can you please change the decimal commas to points and try > again? > > a.y > > On Mon, 19 Aug 2013, Borsodi Florian wrote: > > > Hi Anastasia, > > > > Here are the bvec and bval files you asked me about. > > > > Sincerely, > > Florian > > > > >>> Anastasia Yendiki 18.08.13 22.52 Uhr >>> > > > > Hi Florian - Can you please also attach the bvecs/bvals files from the > > dmri directory? > > > > Thanks, > > a.y > > > > On Sun, 18 Aug 2013, Borsodi Florian wrote: > > > > > Hi Anastasia, > > > > > > Thanks for your help. I renamed the directories and reran trac-all. It > > > worked for the moment. > > > However, then I got another error: "bvecs and bvals don't have the same > > > number of entries." > > > I don't know why, because they actually have the same number of entries. > > > > > > I attached the logfile so that you may see what is going on. The error > > > occours at line 1463. > > > > > > Thank you in advance for your help. > > > > > > Sincerely, > > > Florian > > > > > > > > > >>> Anastasia Yendiki 15.08.13 18.01 Uhr >>> > > > > > > Hi Florian - I see what's going on. Can you please rename the directory > > > /data/physics/florian/nii into something that doesn't include "nii", and > > > rerun? Because "nii" is a file extension, this confuses some of the > > > scripts. > > > > > > Thanks, > > > a.y > > > > > > On Thu, 15 Aug 2013, Borsodi Florian wrote: > > > > > > > Hi Anastasia > > > > > > > > I ran the command you sent me and it runs without any errors. > > > > I also attached the command line for you. > > > > > > > > Furthermore i ran data with 12 diffusion directions and i keep getting the > > > > same kind of problems. > > > > > > > > sincerely > > > > Florian > > > > > > > > >>> Anastasia Yendiki 02.08.13 19.59 Uhr > > > >>> > > > > > > > > Hi Florian - Can you please run this command directly on the command line > > > > and see if there are any errors? > > > > > > > > flip4fsl > > > > /data/physics/florian/nii/1020000485074_20080319/dmri/dwi_orig.nii.gz > > > > /data/physics/florian/nii/1020000485074_20080319/dmri/dwi_orig_flip.nii.gz > > > > > > > > Thanks, > > > > a.y > > > > > > > > On Fri, 2 Aug 2013, Borsodi Florian wrote: > > > > > > > > > Hi Anastasia > > > > > > > > > > Thank you for your answer. > > > > > > > > > > Yes, those files exist. > > > > > I also extracted the bvec an bval files with anotther software then > > > > dcm2nii > > > > > and i keep getting the same kind of error. > > > > > Also, about the 6 diffusion direction, those data i sent you, are just > > > > > testdata for running tracula. > > > > > > > > > > sincerely > > > > > Florian > > > > > > > > > > >>> Anastasia Yendiki 02.08.13 0.41 Uhr > > > >>> > > > > > > > > > > Hi Florian - Do these files exist? > > > >>/data/physics/florian/nii/1020000485074_20080319/dmri/dwi_orig.mghdti.bvec > > > s > > > > > > > > > > > > >>/data/physics/florian/nii/1020000485074_20080319/dmri/dwi_orig.mghdti.bval > > > s > > > > > > > > > > > > > > > > > > > The format of the bvecs and bvals that you sent seems fine, but (and > > > this > > > > > has nothing to do with the error you're getting) I have to warn you that > > > > > you may not be able to get good results running probabilistic > > > tractography > > > > > on this data. It looks like you're acquiring the same 6 diffusion > > > > > directions, repeated 8 times. Six directions is the absolute minimum for > > > > > reconstructing the tensor model, but will probably not be sufficient for > > > > > the crossing fiber model that bedpostx reconstructs (and is used by > > > > > tracula). Your scan time is better spent acquiring 48 different > > > > > directions, rather than the same 6 directions 8 times. > > > > > > > > > > a.y > > > > > > > > > > On Mon, 29 Jul 2013, Borsodi Florian wrote: > > > > > > > > > > > Hello, > > > > > > > > > > > > I have been trying to run tracula on my data and i keep getting the > > > same > > > > > > kind of error repeatedly. > > > > > > For my imput I am using dicoms from a Siemens scanner with a stadard > > > > > > diffusion protocol. However, tracula cannot read the bvec and bval > > > data > > > > > out > > > > > > of the dicom header. Therefore, I implemented bvec and bval files into > > > > the > > > > > > config file. > > > > > > After doing so, I got the error which says that tracula cannot find > > > the > > > > > > "dwi_orig_flip.mghdti.bvecs" and "dwi_orig_flip.mghdti.bvals" files. > > > > > > I have tried many options, but still cannot find a solution. > > > > > > I am hoping to get some troubleshooting help. > > > > > > > > > > > > I have atteched my log- and config- files, as well as my bvec and bval > > > > > > files. > > > > > > To uphold patient confidentiality, I have encode the name patient's > > > data > > > > > as > > > > > > "AAAAAA". > > > > > > > > > > > > I believe the problems are the bvec and bval files or their > > > structures. > > > > > The > > > > > > log-file lete me assume this, especially the lines 202 - 205, 264 - > > > 268 > > > > > and > > > > > > 1480 - 1484. > > > > > > > > > > > > Thank you in advance for your help, and your speedy response. > > > > > > > > > > > > With kind regards, > > > > > > > > > > > > Florian Borsodi > > > > > > Department of Neurology > > > > > > Medical University of Graz > > > > > > > > > > > > > > > > > > > > > > > > > > > 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 Partners Compliance > > > > > HelpLine at > > > > > http://www.partners.org/complianceline . If the e-mail was sent to you > > > in > > > > > error > > > > > but does not contain patient information, please contact the sender and > > > > > properly > > > > > dispose of the e-mail. > > > > > > > > > > > > > > _______________________________________________ > > > > Freesurfer mailing list > > > > Freesurfer(a)nmr.mgh.harvard.edu > > > > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > > > > > > > > > > > > > > > > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >

3 2

Raw Cortical thickness values
by Muhammad Naveed Iqbal Qureshi 06 Jul '17

06 Jul '17

Dear Community, I need to know how can I read the numeric cortical thickness values from freesurfer processed files. (as there are too many in the surf folder) Thank you. Best Regards Naveed

2 2

FSFAST mkcontrast-sess with spmhrf 2 question
by Joseph Dien 25 Mar '16

25 Mar '16

Hi, I have a question about how mkcontrast-sess works. I ran an analysis using the mkanalysis-sess option spmhrf 2 so there are three regressors for each predictor, the primary, the latency, and the dispersion. When specifying the contrast weights for mkcontrast-sess, the documentation indicates that they are specified in terms of the conditions as numbered in the paradigm file, not the individual regressors. Furthermore there only appears to be one contrast value output for each contrast, not three. How are the three regressors being handled? I can think of several scenarios: 1) the contrast weights are not actually in terms of conditions (the documentation is incorrect), they are actually in terms of the regressors (so contrasting conditions 1 and 2 could be specified as -a 1 -a 2 -a 3 -c 4 -c 5 -c 6). 2) the latency and dispersion regressors are being ignored (a common practice). The contrast weights should therefore be specified as -a 1 -c 2. 3) The Calhoun et al (2004) approach is being used to combine the three regressors into a "derivative boost" amplitude measure. The contrast weights should therefore be specified as -a 1 -c 2. Thanks for any help you can give me! Joe -------------------------------------------------------------------------------- Joseph Dien, Senior Research Scientist University of Maryland E-mail: jdien07(a)mac.com Phone: 202-297-8117 http://joedien.com//

2 17

DEPRECATED warning in Freesurfers
by Peter Wendorff 13 Oct '15

13 Oct '15

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi. I'm working with freesurfer as a software developer in a cloud project. My task is to enable control of freesurfer via network (e.g. REST interfaces). That works fine so far, but I stumbled over the following pearl deprecated warning: "Use of ?PATTERN? without explicit operator is deprecated at /usr/local/freesurfer//mni/bin/sharpen_volume line 153." Which has been printed to the log files repeatingly. I'm not a pearl expert, but I think, it's due to a more up to date pearl version used here, that deprecated the usage of ? delimiters for patterns without explicit operator. Using / as a delimiter works fine. I changed that locally to (new line 153): ($output_volume =~ /^([\S]+).mnc/) && ($base_name = $1) || die "sharpen_volume failed: output volume does not appear to be" ." a minc volume.\n"; (the original code here was): ($output_volume =~ ?^([\S]+).mnc?) && ($base_name = $1) || die "sharpen_volume failed: output volume does not appear to be" ." a minc volume.\n"; I'm not sure this is the best solution. It's even not strictly necessary to change this - it's only a warning that's thrown; but I think, it's useful. How can I submit this as a patch? Is it enough to report here? What to do? regards Peter Wendorff P.S.: For referene the header parameters of my sharpen_volume for revision comparison: #---------------------------------------------------------------------------- #$RCSfile: sharpen_volume.in,v $ #$Revision: 1.1 $ #$Author: bert $ #$Date: 2003/04/16 14:29:34 $ #$State: Exp $ #--------------------------------------------------------------------------- # ------------------------------ MNI Header - ---------------------------------- #@NAME : sharpen_volume #@INPUT : #@OUTPUT : #@RETURNS : #@DESCRIPTION: modifies intensities so as to produce a sharper histogram #@METHOD : #@GLOBALS : #@CALLS : #@CREATED : February 28, 1996 #@MODIFIED : #----------------------------------------------------------------------------- -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.11 (GNU/Linux) Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/ iEYEARECAAYFAlGTa5kACgkQi8ffXNvWmYSCUQCghMg1EffqfeJK4PSfo20uF8/J H0MAnjJc0WbHPf5gLDUw70K8fnQPjS+L =g+6e -----END PGP SIGNATURE-----

5 10

PhD position in tractography validation, Inserm U930 "Brain and imaging" lab, University François-Rabelais of Tours
by Frédéric Andersson 06 Oct '14

06 Oct '14

"A funded PhD Student position will be available in the “Brain and Imaging” Inserm unit (national institute for health) settled in Tours, France, from October 1st, 2014 to September 30th, 2017. During this time period, the candidate will have: (1) to develop tools and pertinent criteria for comparison of MR-based tractography results to data obtained from dissection; (2) to compare different algorithms used for MR-based tractography; (3) to develop a tool aiming to help decision making for the choice of optimal algorithms and parameters used for tractography. This work will be co-directed by Pr C Destrieux (Anatomist) and Pr G Venturini (Computer Science) in the university of Tours. It will be funded for the whole time period by a grant we obtained from the Ministry of Research (about 1400 € per month). Potential candidates have to send a letter and a CV to christophe.destrieux(a)univ-tours.fr and gilles.venturini(a)univ-tours.fr before the 15th of April. After an internal selection of the applications, they will have to defend their candidateship during an interview organized by the Doctoral School of the university at the end of May 2014."

1 1

Re: [Freesurfer] possible to use freesurfer segmentations to get mean values for anatomic structures from other sequences
by Douglas Greve 25 Sep '14

25 Sep '14

Hi Sal, yes you can. Use bbregister to register the lowb image to the anatomical. This creates a registration matrix. Use mri_vol2vol and the registration to map the ADC map to the anatomical, then use mri_segstats to extract the values doug ps. Please post questions to the FS list instead of us directly so that others can benefit and it gets archived. thanks! On 5/25/13 1:03 PM, Salil Soman wrote: > Dear Dr. Greve, > > Is it possible to register a nonstructural MR sequences (e.g. an ADC > map) with the aseg.mgz file (or aparc+aseg.mgz file, etc) and then, > using the segmentation from the aseg (or aparc+aseg) file to output a > mean ADC value for each anatomic area segmented? > > Thank you for your time and consideration. > > Best wishes, > > Sal > > Salil Soman, MD, MS

3 11

Cortical Thickness Correction Factor
by David Tate 08 Aug '14

08 Aug '14

We are working on a study that includes data from multiple scanners and initial indication suggest that there may be small differences in cortical thickness measures that are scanner related (human and other phantom data were acquired at both sites). Is it possible to include a correction factor adjusting for the scanner differences when using glmfit or qdec? Thanks, David David F. Tate, Ph.D. Research Neuropsychologist Defense and Veterans Brain Injury Center Contractor, General Dynamics Information Technology Brooke Army Medical Center MCHE MDU (DVBIC) 3551 Roger Brooke Drive JSBA Ft. Sam Houston, TX 78234-4504 Office (210) 916-8090

4 4

lgi for longitudinal scans
by Anna Jonsson 23 Jul '14

23 Jul '14

Dear group, If I want conduct gyrification analyses on the long runs that have been through the longitudinal stream, do I just perform the "normal" baseline command on the .long data? thank you, kind wishes Anna

2 5

Choi2012 convert to FSL MNI152 2mm
by Kathy L. Pearson 10 Jul '14

10 Jul '14

To convert an example Choi2012 file to FSL 2mm space in FreeSurfer 5.2, I have tried both of the mri_vol2vol commands below. The results as shown in freeview are offset by several voxels from the original. Is there a transform or an intermediate step that can make this work -- or a different command? Also, is there a recommended way to extract each individual value (here 1 ... 7) to create separate masks in FSL 2mm space? Thanks for any advice, Kathy Pearson #!/bin/bash input=$FREESURFER_HOME/average/Choi_JNeurophysiol12_MNI152/Choi2012_7Networks_MNI152_FreeSurferConformed1mm_TightMask.nii.gz mri_vol2vol --mov $input --targ $FSLDIR/data/standard/MNI152_T1_2mm_brain.nii.gz --reg $FREESURFER_HOME/subjects/fsaverage/mri.2mm/reg.2mm.mni152.dat --o output1.nii.gz --no-save-reg --interp nearest --precision uchar mri_vol2vol --targ $input --mov $FSLDIR/data/standard/MNI152_T1_2mm_brain.nii.gz --inv --reg $FREESURFER_HOME/average/mni152.register.dat --o output2.nii.gz --no-save-reg --interp nearest --precision uchar # Compare results freeview $input output1.nii.gz output2.nii.gz&

4 4

GLM: what does "mris_preproc" do?
by ye tian 18 Jun '14

18 Jun '14

Hello Freesurfers, I understand that the output of "mris_preproc" can be the input of "mris_glmfit". In mris_glmfit, a GLM is performed. Does "mris_preproc" give you a smoothed average of all the subjects? If that's the case, how can you perform a GLM with just the average? Don't you need the data (smoothed or not) of all subjects? Thank you very much! Sincerely, Ye

2 8

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer April 2014