Freesurfer

freesurfer@nmr.mgh.harvard.edu

3 participants
30135 discussions

longitudinal qdec tables
by Andrews, Derek 11 Jul '16

11 Jul '16

Dear FreeSurfer Developers and Community, I have been trying to run some commands using the -qdec-long flag (I.e. asegstats2table & mris_preproc) but am encountering errors. In the past I have used .rtf files to read in my data using the -fsgd flag but am unsure of how to go about this with longitudinal data, therefore I am using the qdec.dat file structure. I believe something must be wrong with the format of my qdec file. I'm working on MacOSx 10.9.3. I created the qdec in text editor and changed the suffix to .dat under get info (1). I also created one in gedit and saved it using a .dat suffix (2). See these qdec files attached. I receive the following errors: asegstats2table --qdec-long long.qdec.table.dat -t ./outdir/aseg.long.table SUBJECTS_DIR : /Users/dandrews/Data/MIAMRC/MIAMRC_SPGR_FS Traceback (most recent call last): File "/Applications/freesurfer/bin/asegstats2table", line 514, in <module> subj_listoftuples = assemble_inputs(options) File "/Applications/freesurfer/bin/asegstats2table", line 320, in assemble_inputs for row in reader: File "/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/csv.py", line 103, in next self.fieldnames File "/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/csv.py", line 90, in fieldnames self._fieldnames = self.reader.next() _csv.Error: new-line character seen in unquoted field - do you need to open the file in universal-newline mode? And mris_preproc --qdec-long long.qdec.table.dat --target fsaverage --hemi -lh --meas thickness --out lh.thickness.mgh ERROR: no subjects specified I imagine something silly with the qdec file must be going on. Any advice would be very appreciated. Best, Derek ________________________________ Derek Sayre Andrews, MSc PhD Candidate The Sackler Institute for Translational Neurodevelopment Department of Forensic and Neurodevelopmental Sciences Institute of Psychiatry, Psychology & Neuroscience, King's College London Telephone: +44 (0)20 7848 5701 Email: Derek.Andrews(a)kcl.ac.uk

2 3

combining -autorecon2-wm and -autorecon2-cp for segmentation correction
by Barbara Kreilkamp 11 Jul '16

11 Jul '16

Dear freesurfers, Could you please have a look at the segmentation results (attached). I ran this subject several times with control points, eventually I combined the edited wm.mgz and the control point step * and I end up with this result, which has improved segmentation, but it is still not optimal in the medial portion of the temporal lobe. * recon-all -autorecon2-wm -autorecon2-cp -autorecon3 -sd . -subjid REF033c -no-isrunning Could you please advise? Thank you so much, Best wishes, Barbara

1 0

right hemisphere specific problem with converting labels from fsaverage to subject space
by Halbertsma, HN (ohk) 11 Jul '16

11 Jul '16

I’m having issues with the conversion of labels drawn in fsaverage to subject’s space (using mri_label2label). I’ve manually drawn labels (V1-bins) in the fsaverage brain in both the left and right hemisphere that I want to convert to a subjects brain. I’ve successfully managed to do this for the left hemisphere of all the subjects, where the V1-bins converted from the fsaverage fall perfectly within subject specific V1. However, for the right hemisphere, the bins seem slightly misaligned for all the subjects (they’re usually a bit bigger than they should be). Any idea what could have caused this?! Many thanks. Hinke Hinke N. Halbertsma, MSc PhD Student University Medical Center Groningen (UMCG) Department of Ophthalmology NeuroImaging Center (Antonius Deusinglaan 2, 9713 AW Groningen), room 124 T. (+31) 050 363 4431 E. h.n.halbertsma(a)umcg.nl<mailto:h.n.halbertsma@umcg.nl> ________________________________ De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht. The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.

2 4

fornix segmentation using recon-all
by Hari Guragain 10 Jul '16

10 Jul '16

Dear Freesurfer Experts, I am interested to get the volume information of the structure "fornix" in the aseg.stats file. I would appreciate if there is any such argument that can be used along with the recon-all to get the volume information of the fornix like below. Till now, I have been using · *recon-all -autorecon-all -subjid **your_subject_name* Any suggestions will be appreciated. # ColHeaders Index SegId NVoxels Volume_mm3 StructName normMean normStdDev normMin normMax normRange 1 4 6489 6958.2 Left-Lateral-Ventricle 34.8514 11.5767 15.0000 82.0000 67.0000 Thank you, Hari Guragain

2 2

Question about write_surf.m
by Dulce Maria Magdaleno Arroyo 10 Jul '16

10 Jul '16

Hi freesurfers! I'm trying to run the function write_surf in matlab, as far as I understand this function is to write a surface to analyze in freesurfer from matlab data. I'm not sure if the extension of the file name has to be different, so I need your help to solve this simple (I guess) problem. The first line is what I'm trying to run and the next lines show me the error. *write_surf('prueba.nofix/', vertices, faces);* *...subtracting 1 from face indices for FreeSurfer compatibility*. *Error using fwriteInvalid file identifier. Use fopen to generate a valid file identifier.Error in fwrite3 (line 29)fwrite(fid, b1, 'uchar') ;Error in write_surf (line 38)fwrite3(fid, TRIANGLE_FILE_MAGIC_NUMBER);* thank you in advance for your help! Dulce

3 3

Surface segmentation for resting state
by ERIK JAHNER 09 Jul '16

09 Jul '16

Dear free surfer experts, I am conducting a resting state fMRI study using a lagged analysis method. To do this method on a volume, I have used 6 mm cubes as rois covering the whole brain but masked for grey matter. I use these values in a correlation map. I would like to do a similar approach here. Is there a way to get measurements ~6mm apart averaging the vertices 3 mm away. basically, I would like to parcel the cortex surface map into 6mm rois evenly distributed over the surface of the cortex. I would like to take advantage of the fact that I can take measurements along surfaces and not across sulci. Is there an automated way or easy way of building this ROI map. Erik Jahner

2 1

Re: FSFAST error: ERROR: fast_selxavg3() failed\n matlab command not found
by ERIK JAHNER 09 Jul '16

09 Jul '16

On Jul 09, 2016, at 09:00 AM, freesurfer-request(a)nmr.mgh.harvard.edu wrote: Send Freesurfer mailing list submissions to freesurfer(a)nmr.mgh.harvard.edu To subscribe or unsubscribe via the World Wide Web, visit https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer or, via email, send a message with subject or body 'help' to freesurfer-request(a)nmr.mgh.harvard.edu You can reach the person managing the list at freesurfer-owner(a)nmr.mgh.harvard.edu When replying, please edit your Subject line so it is more specific than "Re: Contents of Freesurfer digest..." Today's Topics: 1. Longitudinal LME contrast vectors (K) 2. Re: Longitudinal LME contrast vectors (Martin Reuter) 3. Question about NULL event with Optseq2 (Feng Zhou) 4. Cerebral volume without ventricles (Kaiming Yin) 5. Re: FSFAST error: ERROR: fast_selxavg3() failed\n (Douglas Greve) ---------------------------------------------------------------------- Message: 1 Date: Fri, 8 Jul 2016 15:02:58 -0700 From: K <kam.ntnu(a)gmail.com> Subject: [Freesurfer] Longitudinal LME contrast vectors To: freesurfer(a)nmr.mgh.harvard.edu Message-ID: <CAH54DzgqKWSCwm0gu_ecbB3jrQmLc=Y6Nun5j4DujKwaxn78UA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Hello! We are trying to run a longitudinal morphometry analysis with a control group and clinical group scanned twice. We'd like to look at group x time interactions and would like help figuring out the best way to do so. As the time variable we use centered age -- the deviation from the mean age across all subjects and timepoints. To do so, we use the following lme model: intercept + centeredage.controls + centeredage.clinical + group(0,1) + sex(0,1) Contrast vector: [0 1 -1 0 0] where the centeredage column contains only the centered age for that group and zeros for the opposite group. The thinking behind this is to investigate which cortical regions show opposite effects of time between the groups. However, the FreeSurfer lme software requests a specified "time column". Obviously, we cannot add a centered age for all subjects, nor does it seem to make sense to use one of the two group-specific centeredage columns as the time column. Our alternative model for investigating interatction is: intercept + centeredage + group(0,1) + sex(0,1) + groupxcenteredage where centeredage is the time columns and is a common regressor for both groups and thus, groupxcenteredage becomes 0 for control subjects. Contrast vector: [0 0 0 0 1] This would test whether there is an added effect of centered age in the clinical group. Is there a better way to do an group x time interaction analysis here, or is there a solution for the time column? Thanks!

1 0

Re: [Freesurfer] FSFAST error: ERROR: fast_selxavg3() failed\n
by Douglas Greve 09 Jul '16

09 Jul '16

what do you mean nothing happens? After you hit the return key, it should either give you the path to the command or it should say Command not found. doug ps. Please remember to post to the list and not to us individually On 7/8/16 7:16 PM, Erik Jahner wrote: > When I type "which matlab” in the command window, nothing happens. My > $PATH structure is a s follows: > > echo $PATH > /Applications/freesurfer/bin:/Applications/freesurfer/fsfast/bin:/Applications/freesurfer/tktools:/usr/local/fsl/bin:/Applications/freesurfer/mni/bin:/opt/local/bin:/opt/local/sbin:/opt/local/bin:/opt/local/sbin:/usr/local/bin:/usr/bin:/bin:/usr/sbin:/sbin:/opt/X11/bin > > Sorry if my response seems naïve, I’am new to running command > lines although I’m trying learn. > > Erik > > Message: 6 > Date: Fri, 8 Jul 2016 14:30:43 -0400 > From: Douglas N Greve <greve(a)nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>> > Subject: Re: [Freesurfer] FSFAST error: ERROR: fast_selxavg3() > failed\n > To: freesurfer(a)nmr.mgh.harvard.edu <mailto:freesurfer@nmr.mgh.harvard.edu> > Message-ID: <577FF153.1080602(a)nmr.mgh.harvard.edu > <mailto:577FF153.1080602@nmr.mgh.harvard.edu>> > Content-Type: text/plain; charset=windows-1252; format=flowed > > Is matlab in your path? What happens if you type > which matlab > at the command line? > > On 07/08/2016 12:34 AM, ERIK JAHNER wrote: >> >> Hello FreeSurfer Developers, >> >> >> This is Erik Jahner, I?m attempting to complete the final step in the >> tutorial concerning FSFast Functional connectivity walkthrough >> (http://surfer.nmr.mgh.harvard.edu/fswiki/FsFastFunctionalConnectivityWalkth… ), >> >> but for I get the following error when I run the following command: >> >> selxavg3-sess -s subj439rest -a fc.lpccseed.surf.lh >> >> >> Mon Jul 4 11:01:41 PDT 2016 >> >> rawfunc2surf-sess completed >> >> ---------------------------------------- >> >> Started at Mon Jul 4 11:01:40 PDT 2016 >> >> Ended at Mon Jul 4 11:01:41 PDT 2016 >> >> preproc-sess done >> >> ------------------------------------------- >> >> /Applications/freesurfer/LagStructure/subj439rest >> >> Mon Jul 4 11:01:41 PDT 2016 >> >> anadir = >> /Applications/freesurfer/LagStructure/subj439rest/rest/fc.lpccseed.surf.lh >> >> DoGLMFit = 1 >> >> DoContrasts = 1 >> >> UpdateNeeded = 1 >> >> ------------------------------------------ >> >> ------- matlab output -------------------- >> >> matlab: Command not found. >> >> ------------------------------------------ >> >> *ERROR: fast_selxavg3() failed\n* >> >> >> I've searched the list and no similar errors have been reported. I >> have tried to reinstall MATLAB; I have tried following the >> instructions to install the runtime environment; I have tried redoing >> all analysis from scratch. I have noticed that STEP 5 does not work >> unless I change the -cfg from ?L_Posteriorcingulate.config? >> to ?mean.L_Posteriorcingulate.config?. I don?t think this is the >> problem, but I thought I would point it out. >> >> My session directory is ?subj439rest? >> >> Instead of titling the sub directory ?bold?, I titled it ?rest? >> >> >> 1) FreeSurfer version: freesurfer-Darwin-lion-stable-pub-v5.3.0 >> >> 2) Platform: OSX Yousemite >> >> 3) uname -a: Darwin ERIKs-iMac.local 15.5.0 Darwin Kernel Version >> 15.5.0: Tue Apr 19 18:36:36 PDT 2016; >> root:xnu-3248.50.21~8/RELEASE_X86_64 x86_64 >> >> 4) selxavg3-sess-rest-fc.lpccseed.surf.lh-160704110140.log: see attached >> >> ___________________________________________________ >> >> Here are the command lines I put in for steps 4-8, incase this is of >> some use: >> >> STEP 4 >> >> fcseed-config-segid 1010 -fcname L_Posteriorcingulate.dat -fsd rest >> -mean-cfg mean.L_Posteriorcingulate.config >> >> STEP 5 >> >> fcseed-sess-s subj439rest -cfg mean.L_Posteriorcingulate.config >> >> STEP 6 >> >> * fcseed-config -wm -fcname wm.dat -fsd rest -pca -cfg wm.config >> * fcseed-sess -s subj439rest -cfg wm.config >> >> * fcseed-config -vcsf -fcname vcsf.dat -fsd rest -pca -cfg vcsf.config >> * fcseed-sess -s subj439rest -cfg vcsf.config >> >> STEP7 >> >> mkanalysis-sess-analysis fc.lpccseed.surf.lh -surface fsaverage lh >> -fwhm 5 -notask-taskreg L_Posteriorcingulate.dat 1 -nuisreg vcsf.dat 5 >> -nuisreg wm.dat 5 -mcextreg-polyfit 5 -nskip 4 -fsd rest -TR 2 -per-run >> >> STEP 8 >> >> selxavg3-sess-s subj439rest -a fc.lpccseed.surf.lh >> >> >> >> >> _______________________________________________ >> Freesurfer mailing list >> Freesurfer(a)nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> >> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > -- > 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 > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > > > ------------------------------

1 0

Cerebral volume without ventricles
by Kaiming Yin 09 Jul '16

09 Jul '16

Dear FS experts, Could I ask which value in the aseg.stats file represent the whole cerebral volume without ventricles? Thanks, Kaiming

1 0

Longitudinal LME contrast vectors
by K 08 Jul '16

08 Jul '16

Hello! We are trying to run a longitudinal morphometry analysis with a control group and clinical group scanned twice. We'd like to look at group x time interactions and would like help figuring out the best way to do so. As the time variable we use centered age -- the deviation from the mean age across all subjects and timepoints. To do so, we use the following lme model: intercept + centeredage.controls + centeredage.clinical + group(0,1) + sex(0,1) Contrast vector: [0 1 -1 0 0] where the centeredage column contains only the centered age for that group and zeros for the opposite group. The thinking behind this is to investigate which cortical regions show opposite effects of time between the groups. However, the FreeSurfer lme software requests a specified "time column". Obviously, we cannot add a centered age for all subjects, nor does it seem to make sense to use one of the two group-specific centeredage columns as the time column. Our alternative model for investigating interatction is: intercept + centeredage + group(0,1) + sex(0,1) + groupxcenteredage where centeredage is the time columns and is a common regressor for both groups and thus, groupxcenteredage becomes 0 for control subjects. Contrast vector: [0 0 0 0 1] This would test whether there is an added effect of centered age in the clinical group. Is there a better way to do an group x time interaction analysis here, or is there a solution for the time column? Thanks!

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