Freesurfer

freesurfer@nmr.mgh.harvard.edu

30132 discussions

Weighted Average of Cortical Thickness
by Michael Kranz 24 Apr '17

24 Apr '17

Hi, I had 2 questions relating to measuring global mean cortical thickness with the weighting factor of surface area bh.thickness = ( (lh.thickness * lh.surfarea) + (rh.thickness * rh.surfarea) ) / (lh.surfarea + rh.surfarea) 1.) If this weighting is recommended, then shouldn't it be performed for each vertex in an exploratory whole brain group analyses (e.g., for one vertex in left hemi: lh.thickness*lh.surfarea)/(lh.surfarea). 2.) Are there any references that come to mind to cite in regards to this normalization technique? I see that this following paper uses normalization (but was wondering if you would recommend any others): Vuoksimaa, E., Panizzon, M. S., Chen, C.-H., Fiecas, M., Eyler, L. T., Fennema-Notestine, C., … Kremen, W. S. (2015). The Genetic Association Between Neocortical Volume and General Cognitive Ability Is Driven by Global Surface Area Rather Than Thickness. *Cerebral Cortex (New York, NY)* , *25*(8), 2127–2137. http://doi.org/10.1093/cercor/bhu018 Best, Mike

1 0

volume in native space of .vtk file in mni152 space
by Corinna Bauer 23 Apr '17

23 Apr '17

Dear all, I have a series of .vtk files in MNI152 1mm space from an external thalamus atlas outlining specific thalamic nuclei. I would like to put these into native subject space and calculate the volume from each nuclei. I would imagine that I would do the following: 1. run mni152reg --s --1 on each subject 2. mri_vol2vol --mov mni152space_atlas_file.nii.gz -- targ subject/mri/orig.mgz --reg subject/mri/mni152.orig.mgz.reg --inv --o newoutputfile_subject_space But then I am not certain how to calculate the volume of the new output file. Suggestions? Thanks! Corinna

2 3

Difference between lh.thickness file and lh.thickness.fsaverage.mgh file
by Zhiliang Long 23 Apr '17

23 Apr '17

Dear FS experts: I conducted recon-all analysis on my own subjects, and found that there are several kinds of thickness files, such as the lh.thickness file and lh.thickness.fsaverage file. By using MATLAB scripts implemented in FS, different numbers of vertex were found between the two files. What I am thinking is that the lh.thickness file is the individual thickness map (It has not been transformed into fsaverage space), while the lh.thickness.fsaverage.mgh is the thickness map in fsaverage space. Is that correct? Best, Zhiliang

2 1

R: R: Re: R: Re: Hipposubfields in FreeSuerferColorLUT.txt
by stdp82＠virgilio.it 23 Apr '17

23 Apr '17

Dear list and Eugenio, I have run for each subject mri_label2vol --label SUBJ/label/rh.perirhinal_exvivo.thresh.label --regheader SUBJ/mri/norm.mgz --fill-ribbon --subject SUBJ --hemi rh --o SUBJ/mri/rh.perirhinal.nii.gz next, to use FS-FAST on perihrinal.nii.gz, I have run fcseed-config -seg SUBJ/mri/lh.perirhinal.nii.gz -fcname L_Perirhinal.dat -fsd rest -mean -cfg mean.L_Perirhinal.config This error occurs: ERROR: no segmentations IDs (-segid) or ROI (-roi) specified What's should I do? Thanks Stefano On Mon, 06 Feb 2017 13:59:38 -0800 Do you mean to use them as seeds in FSFAST? You'll need to map them into the conformed space using mri_label2vol with the --regheader option (see the --help); view the volume to make sure it looks ok. Then run fcseed-config specifying the volume with -seg and the code with segid ----Messaggio originale---- Da: stdp82(a)virgilio.it Data: 6-feb-2017 22.27 A: "Freesurfer support list"<freesurfer(a)nmr.mgh.harvard.edu> Ogg: [Freesurfer] R: Re: R: Re: Hipposubfields in FreeSuerferColorLUT.txt Thank you very much.Could you provide me any suggestion to extract CA1 left or CA1 right?And how can I use them in FAST? By codes?Regards,Stefano ----Messaggio originale---- Da: "Iglesias Gonzalez, Eugenio" <e.iglesias(a)ucl.ac.uk> Data: 6-feb-2017 17.42 A: "stdp82(a)virgilio.it"<stdp82(a)virgilio.it>, "Freesurfer support list"<freesurfer(a)nmr.mgh.harvard.edu> Ogg: Re: [Freesurfer] R: Re: Hipposubfields in FreeSuerferColorLUT.txt You would extract CA1 left from lh.hippoSfLabels-T1.v10.mgz and CA1 right from rh.hippoSfLabels-T1.v10.mgz (assuming you’ve used T1s to segment the subfields). Cheers, /E Juan Eugenio Iglesias ERC Senior Research Fellow Translational Imaging Group University College London http://www.jeiglesias.com http://cmictig.cs.ucl.ac.uk/ On 6 Feb 2017, at 16:17, stdp82(a)virgilio.it wrote: Thanks Eugenio for your response. My apologies but this matter is unclear to me. If I would use results for FS-FAST how I can identify for example CA1 right from CA1 left? Additionally, for enthorinal which is the correct code in FAST? Stefano ----Messaggio originale---- Da: "Iglesias Gonzalez, Eugenio" <e.iglesias(a)ucl.ac.uk> Data: 6-feb-2017 15.54 A: "stdp82(a)virgilio.it"<stdp82(a)virgilio.it>, "Freesurfer support list"<freesurfer(a)nmr.mgh.harvard.edu> Ogg: Re: [Freesurfer] Hipposubfields in FreeSuerferColorLUT.txt Ciao Stefano, The results for the left and right hippocampi are in different files, so there’s no ambiguity despite the fact that they indeed share the same label codes. Cheers, /Eugenio Juan Eugenio Iglesias ERC Senior Research Fellow Translational Imaging Group University College London http://www.jeiglesias.com http://cmictig.cs.ucl.ac.uk/ On 6 Feb 2017, at 13:18, stdp82(a)virgilio.it wrote: Hi list, within FreeSuerferColorLUT.txt, do the numbers between 200-226 descrived the hippo-subfield outputs? Additionally, how can I obtain separate right and left masks? Thanks Stefano _______________________________________________ Freesurfer mailing list Freesurfer(a)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 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.

1 0

mris_label2annot not creating annot file
by Das S. 23 Apr '17

23 Apr '17

Hi, I was using the mris_label2annot to create .annot file before but now couldn't figure out why the .annot file is not being created. The mris_label2annot is giving the below output but not creating the label files. My label file names are as node-lh-xx.label and the content of ctab file as below: 0 node-lh-0 25 5 25 0 1 node-lh-1 70 0 180 0 2 node-lh-2 245 245 245 0 3 node-lh-3 205 62 78 0 4 node-lh-4 120 18 134 0 5 node-lh-5 196 58 250 0 6 node-lh-6 0 148 0 0 7 node-lh-7 220 248 164 0 8 node-lh-8 230 148 34 0 9 node-lh-9 0 118 14 0 10 node-lh-10 70 118 14 0 11 node-lh-11 122 186 220 0 12 node-lh-12 236 13 176 0 13 node-lh-13 12 48 255 0 14 node-lh-14 204 182 142 0 *********************************************************** mris_label2annot --s 100307 --h lh --no-unknown --ctab label/eq_parcellation_label/MyColorLUTleft01.txt --a eq40left --ldir label/eq_parcellation_label/left/ Reading ctab label/eq_parcellation_label/MyColorLUTleft01.txt Number of ctab entries 41 INFO: no labels specified, generating from ctab $Id: mris_label2annot.c,v 1.20 2016/01/07 23:28:11 greve Exp $ cwd /home/koushik/Parcellation/100307 cmdline mris_label2annot --s 100307 --h lh --no-unknown --ctab label/eq_parcellation_label/MyColorLUTleft01.txt --a eq40left --ldir label/eq_parcellation_label/left/ sysname Linux hostname koushik-DQ77MK machine x86_64 user koushik subject 100307 hemi lh SUBJECTS_DIR /home/koushik/Parcellation ColorTable label/eq_parcellation_label/MyColorLUTleft01.txt AnnotName eq40left nlables 0 LabelThresh 0 0.000000 Loading /home/koushik/Parcellation/100307/surf/lh.orig Writing annot to /home/koushik/Parcellation/100307/label/lh.eq40left.annot Thanks Sarbani

1 1

aparc files missing a comma
by Mike Schmidt 22 Apr '17

22 Apr '17

Hello FreeSurfer Developers, I'm attempting to parse stats files with regular expressions. Everything seems to work in all of the aseg and wmparc files, but the {lh,rh}.aparc{,a2009s,DKTatlas,pial}.stats files are missing a comma after "CortexVol" on line 25. This causes my regex to miss that measurement on each pass. All other lines match my regex perfectly, only CortexVol misses the comma. This is a new field in version 6; it does not exist in versions 5.2.0 and 5.3.0. This is consistent across all subjects I've viewed, and there are no errors to dig into, but I've attached the recon-all.log file and a few stats files from one of them in case they're of any use. 1) FreeSurfer version: freesurfer-Linux-centos6_x86_ 64-stable-pub-v6.0.0-2beb96c 2) Platform: external linux cluster, unsure exactly 3) uname -a: Linux sequoia 3.12.67-60.64.24-default #1 SMP Fri Dec 9 14:35:11 UTC 2016 (8154761) x86_64 x86_64 x86_64 GNU/Linux 4) recon-all.log: see attached, along with a few aparc files. I have added the comma using sed to all of the aparc files in our system, but I hope this is helpful in removing bugs for future use. Best wishes, Mike Schmidt

1 0

Superior Posterior Insula and DLPFC measurements
by Guilherme Alves 22 Apr '17

22 Apr '17

Hello everyone, I’ve been having a hard time trying to figure out how to get appropriate measures of gray matter volume for both the superior posterior insula and the DLPFC. The atlases included in freesurfer don’t provide such information after the recon-all process is run. Is there a way to get accurate measures without having to manually label these areas for each subject? I’ve seen some answers on the archive defining the posterior insula as G_Ins_lg_and_S_cent_ins, but how can I isolate the superior posterior insula? Is there a way to create a new label based on a combination of pre-existing labels? Thanks in advance, Guilherme.

1 0

Need help to run freesurfer
by Ramesh Babu 22 Apr '17

22 Apr '17

Dear Experts, I have installed in /usr/local/freesurfer. I am trying to follow the command given in the freesurfer web page. $> export FREESURFER_HOME=/usr/local/freesurfer $> source $FREESURFER_HOME/SetUpFreeSurfer.sh But I could not figure out the .sh file in the freesurfer folder. Please see the list of files in the freesurfer folder given below. List of files in the freesurfer folder; /usr/local/freesurfer/average /usr/local/freesurfer/fsafd /usr/local/freesurfer/fsfast /usr/local/freesurfer/lib /usr/local/freesurfer/sessions /usr/local/freesurfer/tktools /usr/local/freesurfer/trctrain /usr/local/freesurfer/DefectLUT.txt /usr/local/freesurfer/FreeSurferColorLUT.txt /usr/local/freesurfer/FreeSurferEnv.csh /usr/local/freesurfer/LICENSE /usr/local/freesurfer/SetUpFreeSurfer.csh /usr/local/freesurfer/sources.csh /usr/local/freesurfer/sources.sh /usr/local/freesurfer/tkmeditParcColorsCMA /usr/local/freesurfer/VERSION /usr/local/freesurfer/WMParcStatsLUT.txt You can notice that there is no SetUpFreeSurfer.sh file. Details about my pc; Ubuntu 16.04 LTS, 64 bit is installed. Running w8.1 via virtual box. Freesurfer is installed in Ubuntu OS. I am using bash shell. pc@RB:~$ echo $0 bash pc@RB:~$ echo $SHELL /bin/bash pc@RB:~$ ps -p "$$" PID TTY TIME CMD 3706 pts/1 00:00:00 bash Need your help to run freesurfer in my pc. Thanks in advance Ramesh

2 3

Multiple comparisons
by Laura Ferrero Montes 22 Apr '17

22 Apr '17

Dear FS team, I would like to correct for multiple comparisons after I run command-line GLM analysis on stats tables. I would like to make a second level analysis between a control group and an autism group. I would like to know if there is a command-line in FreeSurfer or if there is any program or toolkit in Matlab that allows to do that. Thank you, Laura

3 6

mri_glmfit error with aparc table flag
by miracle ozzoude 21 Apr '17

21 Apr '17

Hello freesurfer, I am trying to run mri_glmfit with the --table flag. My aparc_lh_thick_stats.txt contain entorhinal and parahippocampal thickness values for 32 subjects. I ran the glm command and received an error. Please, how can i solve this issue? I have included my command, fsgd file, contrast matrix and error below. Thanks Best, Paul GLM command: mri_glmfit --table aparc_lh_thick_stats.txt --fsgd age.fsgd --C age1.mtx --C age2.mtx --surf fsaverage lh --cortex --glmdir lh.thickness.glmdir ERROR= "gdfRead(): reading age.fsgd INFO: DeMeanFlag keyword not found, DeMeaning will NOT be done. Continuous Variable Means (all subjects) 0 Age 39.9375 17.3312 Class Means of each Continuous Variable 1 Group1 35.6316 2 Group2 46.2308 INFO: gd2mtx_method is dods Reading source surface /Users/MiracleOz/Documents/improvervsdeclinermri/fsaverage/surf/lh.white Number of vertices 163842 Number of faces 327680 Total area 65416.984375 AvgVtxArea 0.399269 AvgVtxDist 0.721953 StdVtxDist 0.195470 $Id: mri_glmfit.c,v 1.241.2.4 2016/12/08 22:02:40 zkaufman Exp $ cwd /Users/MiracleOz/Documents/improvervsdeclinermri cmdline mri_glmfit.bin --table aparc_lh_thick_stats.txt --fsgd age.fsgd --C age1.mtx --C age2.mtx --surf fsaverage lh --cortex --glmdir lh.thickness.glmdir sysname Darwin hostname MacBook-Pro.local machine x86_64 user MiracleOz FixVertexAreaFlag = 1 UseMaskWithSmoothing 1 OneSampleGroupMean 0 y /Users/MiracleOz/Documents/improvervsdeclinermri/aparc_lh_thick_stats.txt logyflag 0 usedti 0 FSGD age.fsgd labelmask /Users/MiracleOz/Documents/improvervsdeclinermri/fsaverage/label/lh.cortex.label maskinv 0 glmdir lh.thickness.glmdir IllCondOK 0 ReScaleX 1 DoFFx 0 Creating output directory lh.thickness.glmdir Loading y from /Users/MiracleOz/Documents/improvervsdeclinermri/aparc_lh_thick_stats.txt Found 66 data colums Found 0 data rows INFO: gd2mtx_method is dods Saving design matrix to lh.thickness.glmdir/Xg.dat Computing normalized matrix Normalized matrix condition is 32.2165 Matrix condition is 20734.7 Found 149955 points in label. ERROR: mri_reshape: number of elements cannot change nv1 = 163842, nv1 = 66 search space = 82219.962936 ERROR: dimension mismatch between y and X. y has 0 inputs, X has 32 rows." FSGD file: GroupDescriptorFile 1 Title MOT Class Group1 Class Group2 Variables Age Input 01053p Group1 23 Input 01054p Group1 43 Input 01061p Group2 39 Input 01062p Group2 19 Input 01074p Group2 28 .............. (I have more participants) CONTRAST MATRIX: 1 -1 0 0

4 15

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