Freesurfer February 2014

freesurfer@nmr.mgh.harvard.edu

128 participants
236 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

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

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

v.5.3 Longitudinal Analysis on v.5.0 processed scans
by Elijah Mak 04 Jun '14

04 Jun '14

Greetings Freesurfer Community, I am running a longitudinal study using Freesurfer v.5.3 to compare rate of cortical thinning between 2 groups of subjects. I hope I can seek some expert advice on the following situation: All data (both the baseline and follow-up) has been processed independently using recon-all. This has been done on Freesurfer 5.0. Since then, I have upgraded to Freesurfer v5.3 to utilize the linear mixed effect capability. For the construction of base and longs, I understand that it is important to add -clean flag to my recon-all strings to overwrite brainmask.mgz etc. When must this be done? I am still new to Freesurer, and I'd greatly appreciate any advice regarding workflow. Many Thanks, Elijah -- Elijah Mak, B.A. Doctoral Candidate *|* Psychiatry University of Cambridge *| *Department of Psychiatry Trinity College, Cambridge, CB2 1TQ

2 4

brain.finalsurfs.mgz editing
by Khadka, Sankalpa 02 May '14

02 May '14

Dear freesurfer community, My question is directed to those who have experience with editing brain.finalsurfs.mgz. According to the tutorials, after making edits on brainmask.finalsurfs.mgz (which is normally done when the pial surface is not accurate), we can run, recon-all -autorecon-pial -subjid subj" However, I found an anomaly when I looked at the steps that -autorecon-pial flag goes through. The very first step of autorecon-pial flag is "mri_mask" command. This takes brain.mgz and brainmask.mgz as inputs and gives brain.finalsurfs.mgz as output. Now, if we have made edits on brain.finalsurfs.mgz, it will be overwritten by the new one and hence our edits won't take an effect. I ran few experiments too. First, I made edits on brain.finalsurfs.mgz and ran recon-all with autorecon-pial flag. I ended up with a new brain.finalsurfs.mgz and my edits did not make any changes. Second, I made edits on brain.finalsurfs.mgz and ran with white, smooth2, inflate2, autorecon3 flag (basically starting from mri_make_surfaces command). On doing so, the edits made desirable changes and brain.finalsurfs.mgz remained intact. Third, I made changes on brainmask.mgz and brain.mgz (inputs for mri_mask). It also works as expected. Finally, I also followed the recommended steps from the tutorial by creating brainmask.finalsurfs.manedit.mgz, however, the edits from "manedit" were not carried on to brainmask.finalsurfs.mgz Could anyone please shed some light on this? What would be the best strategy to follow? Thank you, Sankalpa

2 1

V1 label
by krista kelly 24 Apr '14

24 Apr '14

Hello, I am looking at cortical thickness and surface area using Freesurfer. I have transferred label V1 from freesurfer onto the brains of my participants and when I view them on the inflated brains in tksurfer , I notice that there are holes in the label. Is this an issue when calculating mean cortical thickness and total surface area estimations? Should I fill these in somehow? I've attached a picture to show as an example. Please let me know if the picture cannot be viewed. Thanks! Krista

4 11

Re: [Freesurfer] mri_convert Siemens .IMA error
by Douglas N Greve 22 Apr '14

22 Apr '14

Hi Paul, I was finally able to track down the bug. You can get a new version from here: ftp://surfer.nmr.mgh.harvard.edu/pub/dist/freesurfer/5.3.0-patch/mri_conver… copy it to $FREESURFER_HOME/bin/mri_convert doug On 12/08/2013 06:28 PM, prasser wrote: > Hi Doug, > > Sorry for the delay. Please find 1.IMA attached. > > Thanks, > Paul > > > ---- On Fri, 06 Dec 2013 13:15:29 -0800 *Douglas N Greve > <greve(a)nmr.mgh.harvard.edu>* wrote ---- > > > Can you upload a single dicom file to me? One should be enough > > On 12/06/2013 04:02 PM, Paul Rasser wrote: > >> > >> We are using Freesurfer5.3 on os x 10.8. > >> > > > > > > On 6 Dec 2013, at 2:19 pm, Douglas Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > > <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > > >> > >> This may be a known problem. What version of FS are you using and > >> what is your platform? > >> doug > >> > >> > >> > >> > >> On 12/5/13 3:26 PM, prasser wrote: > >>> Thanks, but the orientation shift appears as a skew. Also, we > don't > >>> have any fiducial markers for this group of scans so it would > be a > >>> little tricky to confirm the correct application of > >>> -out_orientation. I would rather let mri_convert handle this from > >>> the .IMAs. > >>> > >>> Is there any additional information I can provide from from > the .IMA > >>> header etc that could shed some light on why mri_convert is > behaving > >>> in this way? > >>> > >>> Thanks, > >>> Paul > >>> > >>> > >>> > >>> ---- On Wed, 04 Dec 2013 22:11:27 -0800 *Mark Plantz > >>> <markplantz2016(a)u.northwestern.edu > <mailto:markplantz2016@u.northwestern.edu>>* wrote ---- > >>> > >>> Not sure if this helps, but you can alter the output orientation > >>> using the flag "--out_orientation" followed by LIA (for > >>> example). Maybe the input file needs to be in a specific > >>> orientation before it can be converted properly. > >>> > >>> > >>> On Wed, Dec 4, 2013 at 9:56 PM, prasser <prasser(a)zoho.com > <mailto:prasser@zoho.com> > >>> <mailto:prasser@zoho.com>> wrote: > >>> > >>> Hi, > >>> > >>> I'm currently having problems converting .IMA files > >>> collected from a Siemens Skyra to .mgz format using > >>> mri_convert from Freesurfer5.3 on os x. > >>> > >>> Using: > >>> > >>> mri_convert -it siemens -ot mgz 1.IMA test.mgz > >>> > >>> I get the error message: > >>> > >>> does not contain a Siemens ASCII header has this file been > >>> anonymised? > >>> > >>> ERROR: cannot unpack mosiacs without ASCII header > >>> > >>> The .IMA files have been taken straight from the scanner > >>> console. > >>> > >>> If I try: > >>> > >>> mri_convert -it dicom -ot mgz 1.IMA test.mgz > >>> > >>> this creates the .mgz, however there is an orientation shift > >>> and a vertical stripe artefact appears in the sagittal view > >>> (please see attached) when compared to the original .IMA. > >>> > >>> Any suggestions on how to convert the Siemens .IMA files > >>> correctly? > >>> > >>> Thanks, > >>> > >>> Paul > >>> > >>> > >>> _______________________________________________ > >>> Freesurfer mailing list > >>> Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > >>> <mailto: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 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 > <mailto:Freesurfer@nmr.mgh.harvard.edu> > >>> <mailto: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 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 > <mailto:Freesurfer@nmr.mgh.harvard.edu> > >>> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >> > >> _______________________________________________ > >> Freesurfer mailing list > >> Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > <mailto: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 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. > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 > Fax: 617-726-7422 > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > Outgoing: > ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

2 1

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 February 2014