Freesurfer March 2019

freesurfer@nmr.mgh.harvard.edu

114 participants
172 discussions

skip subjects already processed
by Akshay Patel 18 Mar '19

18 Mar '19

External Email - Use Caution Hello, I have a subset of subjects that were processed with an older version of FreeSurfer and would like to run FreeSurfer in the subjects directory such that it skips all the subjects already processed by my version of freesurfer and only reruns recon-all ... for those aforementioned subjects. When I do this, will the file structure and files for the subject be overwritten? How does FreeSurfer behave in this case? How does it behave for subjects directories that have already been successfully processed with the same command (including options) attempting to process it at the moment? Would appreciate any help thanks. Best, Akshay

1 0

Cortical thickness values extraction issue
by Giuliana Klencklen 18 Mar '19

18 Mar '19

External Email - Use Caution Hi FS experts, I did group-level, surface-based, vertex-wise analysis for baseline and longitudinal data. I used Qdec and do the same work with the fsgd version (mri_glmfit-sim command) to double-check the data. I created label files with tksurfer for each of the clusters showing a significant between-group difference. Then, I used the following command stream to extract the cortical thickness values for each subject and cluster: mri_label2label, mris_anatomical_stats, and aparcstats2table. E.g., mri_label2label --srcsubject fsaverage --srclabel /home/jagust/gklenck/Long_MRI/lh.ac-baseline-rostralmiddlefrontal.label --trgsubject ${s}_tp1 --trglabel ${s}_tp1/label/lh.ac-baseline-rostralmiddlefrontal.label --hemi lh --regmethod surface mris_anatomical_stats -l lh.ac-baseline-rostralmiddlefrontal.label -t lh.thickness -b -f ${s}_tp1/stats/lh.ac-baseline-rostralmiddlefrontal.stats ${s}_tp1 lh aparcstats2table --subjects ${s}_tp1 --hemi lh --parc ac-baseline-rostralmiddlefrontal --meas thickness --tablefile ${out_dir}/lh.ac-baseline-rostralmiddlefrontal.aparc_stats.txt Subsequently, when I conducted between-group comparisons for each cluster, a couple of clusters (7 out of 16) do not show a significant difference - which does not make sense. This problem seems to appear randomly. Help to the FS archives, I tried to use mri_segstats command but was not able to find a correct combination of arguments. If this is the right track to solve this problem, what combination of arguments should I use? And if this is not, do you have any idea how I can solve this problem? Many thanks in advance. Regards, Giuliana -- Giuliana Klencklen, Ph.D. Helen Wills Neuroscience Institute University of California, Berkeley 118 Barker Hall Berkeley, CA 94720-3190 510-395-0040 giuliana.klencklen(a)berkeley.edu

2 8

回复：回复： Error with mri_label2label(Segmentation fault(coredumped))
by 310913949 18 Mar '19

18 Mar '19

External Email - Use Caution Dear experts, I think I need to send you the information of my ROI,i am not sure whether it is something wrong. This is the part of the 4gt.label I drew.And the attachment is the picture of the 4gt.label. Thanks again. #!ascii label , from subject vox2ras=scanner 59455 -1 43.255 58.553 -24.301 0.0000000000 -1 43.255 58.553 -24.301 0.0000000000 -1 42.058 58.629 -24.306 0.0000000000 -1 40.860 58.704 -24.310 0.0000000000 -1 40.860 58.704 -24.310 0.0000000000 -1 40.799 57.729 -24.310 0.0000000000 -1 40.799 57.729 -24.310 0.0000000000 -1 39.601 57.805 -24.314 0.0000000000 -1 39.601 57.805 -24.314 0.0000000000 -1 38.404 57.880 -24.318 0.0000000000 -1 38.404 57.880 -24.318 0.0000000000 -1 37.206 57.955 -24.322 0.0000000000 -1 37.206 57.955 -24.322 0.0000000000 -1 35.947 57.056 -24.326 0.0000000000 -1 35.947 57.056 -24.326 0.0000000000 -1 34.749 57.131 -24.330 0.0000000000 -1 34.749 57.131 -24.330 0.0000000000 -1 33.552 57.207 -24.334 0.0000000000 -1 33.552 57.207 -24.334 0.0000000000 -1 32.293 56.307 -24.338 0.0000000000 -1 32.293 56.307 -24.338 0.0000000000 -1 32.293 56.307 -24.338 0.0000000000 -1 32.232 55.333 -24.338 0.0000000000 -1 32.232 55.333 -24.338 0.0000000000 -1 31.034 55.408 -24.342 0.0000000000 -1 29.836 55.483 -24.346 0.0000000000 -1 29.836 55.483 -24.346 0.0000000000 -1 29.836 55.483 -24.346 0.0000000000 -1 28.639 55.558 -24.350 0.0000000000 -1 28.639 55.558 -24.350 0.0000000000 -1 28.577 54.584 -24.350 0.0000000000 -1 27.318 53.685 -24.354 0.0000000000 -1 27.318 53.685 -24.354 0.0000000000 -1 27.257 52.710 -24.354 0.0000000000 -1 27.257 52.710 -24.354 0.0000000000 -1 27.257 52.710 -24.354 0.0000000000 -1 27.196 51.735 -24.354 0.0000000000 -1 27.196 51.735 -24.354 0.0000000000 -1 27.135 50.761 -24.354 0.0000000000 -1 27.135 50.761 -24.354 0.0000000000 ------------------ 原始邮件 ------------------ 发件人: "邓刊"<310913949(a)qq.com>; 发送时间: 2019年3月14日(星期四) 上午10:03 收件人: "Freesurfer support list"<freesurfer(a)nmr.mgh.harvard.edu>; 主题: [Freesurfer] 回复： Error with mri_label2label(Segmentation fault(coredumped)) External Email - Use Caution Sure.This is the full terminal output: root@zourushi-All-Series:/usr/local/freesurfer/subjects/PT1_4# mri_label2label --srclabel label/4gt.label --s PT1_4 --trglabel 4gtreg.label --hemi rh --paint 30 pial --trgsurf pial Could not set locale No such file or directory Could not set locale srclabel = label/4gt.label srcsubject = PT1_4 trgsubject = PT1_4 trglabel = 4gtreg.label regmethod = surface srchemi = rh trghemi = rh trgsurface = pial srcsurfreg = sphere.reg trgsurfreg = sphere.reg usehash = 1 Use ProjAbs = 0, 0 Use ProjFrac = 0, 0 DoPaint 1 PaintMax 30.000000 SUBJECTS_DIR /usr/local/freesurfer/subjects FREESURFER_HOME /usr/local/freesurfer Loading source label. Found 59455 points in source label. Starting surface-based mapping Reading source registration /usr/local/freesurfer/subjects/PT1_4/surf/rh.sphere.reg 59455 unassigned vertices in label - building spatial LUT... assigning vertex numbers to label... /usr/local/freesurfer/bin/mri_label2label: line 3: 4446 Segmentation fault (core dumped) mri_label2label.bin "$@" ------------------ 原始邮件 ------------------ 发件人: "Greve, Douglas N.,Ph.D."<DGREVE(a)mgh.harvard.edu>; 发送时间: 2019年3月13日(星期三) 凌晨0:06 收件人: "freesurfer(a)nmr.mgh.harvard.edu"<freesurfer(a)nmr.mgh.harvard.edu>; 主题: Re: [Freesurfer] Error with mri_label2label(Segmentation fault(coredumped)) Can you send the full terminal output? On 3/12/19 10:36 AM, 310913949 wrote: External Email - Use Caution Hi Bruce, I drawed a 3D ROI on the original image using freeview,and saved it as an label file. Then I wanted to sample it to my subject and to found the closest surface vertices so that I can extract the cortical thickness of ROI. I used the command line: mri_label2label --srclabel 4gt.label --s PT1_4 --trglabel 4reg.on-pial.label --hemi rh --paint 30 pial --trgsurf pial. But it reported error: Line 3:4445 Segmentation fault(core dumped), mri_label2label.bin "$@". Can you tell me how to solve this problem?Is it wrong to use the mri_label2label? Thanks. Best wishes， Kan Deng _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

1 0

Cortical Thickness Extraction
by Bunting, Emma 18 Mar '19

18 Mar '19

External Email - Use Caution Hi Freesurfer experts, I used a script created by ENIGMA to extract the cortical thickness values for 129 subjects that had undergone reconall processing. We are using these values in VBM analysis. From reading other papers, we notice that all regional volumes are usually normalised by dividing by total intracranial volumes for each subject. I’m just wondering whether the values we are extracting from our freesurfer output files have already had this normalisation applied? The script we’ve used as follows: #!/bin/bash echo 'SubjID,L_bankssts_thickavg,L_caudalanteriorcingulate_thickavg,L_caudalmiddlefrontal_thickavg,L_cuneus_thickavg,L_entorhinal_thickavg,L_fusiform_thickavg,L_inferiorparietal_thickavg,L_inferiortemporal_thickavg,L_isthmuscingulate_thickavg,L_lateraloccipital_thickavg,L_lateralorbitofrontal_thickavg,L_lingual_thickavg,L_medialorbitofrontal_thickavg,L_middletemporal_thickavg,L_parahippocampal_thickavg,L_paracentral_thickavg,L_parsopercularis_thickavg,L_parsorbitalis_thickavg,L_parstriangularis_thickavg,L_pericalcarine_thickavg,L_postcentral_thickavg,L_posteriorcingulate_thickavg,L_precentral_thickavg,L_precuneus_thickavg,L_rostralanteriorcingulate_thickavg,L_rostralmiddlefrontal_thickavg,L_superiorfrontal_thickavg,L_superiorparietal_thickavg,L_superiortemporal_thickavg,L_supramarginal_thickavg,L_frontalpole_thickavg,L_temporalpole_thickavg,L_transversetemporal_thickavg,L_insula_thickavg,R_bankssts_thickavg,R_caudalanteriorcingulate_thickavg,R_caudalmiddlefrontal_thickavg,R_cuneus_thickavg,R_entorhinal_thickavg,R_fusiform_thickavg,R_inferiorparietal_thickavg,R_inferiortemporal_thickavg,R_isthmuscingulate_thickavg,R_lateraloccipital_thickavg,R_lateralorbitofrontal_thickavg,R_lingual_thickavg,R_medialorbitofrontal_thickavg,R_middletemporal_thickavg,R_parahippocampal_thickavg,R_paracentral_thickavg,R_parsopercularis_thickavg,R_parsorbitalis_thickavg,R_parstriangularis_thickavg,R_pericalcarine_thickavg,R_postcentral_thickavg,R_posteriorcingulate_thickavg,R_precentral_thickavg,R_precuneus_thickavg,R_rostralanteriorcingulate_thickavg,R_rostralmiddlefrontal_thickavg,R_superiorfrontal_thickavg,R_superiorparietal_thickavg,R_superiortemporal_thickavg,R_supramarginal_thickavg,R_frontalpole_thickavg,R_temporalpole_thickavg,R_transversetemporal_thickavg,R_insula_thickavg,LThickness,RThickness,LSurfArea,RSurfArea,ICV' > CorticalMeasuresENIGMA_ThickAvg.csv echo 'SubjID,L_bankssts_surfavg,L_caudalanteriorcingulate_surfavg,L_caudalmiddlefrontal_surfavg,L_cuneus_surfavg,L_entorhinal_surfavg,L_fusiform_surfavg,L_inferiorparietal_surfavg,L_inferiortemporal_surfavg,L_isthmuscingulate_surfavg,L_lateraloccipital_surfavg,L_lateralorbitofrontal_surfavg,L_lingual_surfavg,L_medialorbitofrontal_surfavg,L_middletemporal_surfavg,L_parahippocampal_surfavg,L_paracentral_surfavg,L_parsopercularis_surfavg,L_parsorbitalis_surfavg,L_parstriangularis_surfavg,L_pericalcarine_surfavg,L_postcentral_surfavg,L_posteriorcingulate_surfavg,L_precentral_surfavg,L_precuneus_surfavg,L_rostralanteriorcingulate_surfavg,L_rostralmiddlefrontal_surfavg,L_superiorfrontal_surfavg,L_superiorparietal_surfavg,L_superiortemporal_surfavg,L_supramarginal_surfavg,L_frontalpole_surfavg,L_temporalpole_surfavg,L_transversetemporal_surfavg,L_insula_surfavg,R_bankssts_surfavg,R_caudalanteriorcingulate_surfavg,R_caudalmiddlefrontal_surfavg,R_cuneus_surfavg,R_entorhinal_surfavg,R_fusiform_surfavg,R_inferiorparietal_surfavg,R_inferiortemporal_surfavg,R_isthmuscingulate_surfavg,R_lateraloccipital_surfavg,R_lateralorbitofrontal_surfavg,R_lingual_surfavg,R_medialorbitofrontal_surfavg,R_middletemporal_surfavg,R_parahippocampal_surfavg,R_paracentral_surfavg,R_parsopercularis_surfavg,R_parsorbitalis_surfavg,R_parstriangularis_surfavg,R_pericalcarine_surfavg,R_postcentral_surfavg,R_posteriorcingulate_surfavg,R_precentral_surfavg,R_precuneus_surfavg,R_rostralanteriorcingulate_surfavg,R_rostralmiddlefrontal_surfavg,R_superiorfrontal_surfavg,R_superiorparietal_surfavg,R_superiortemporal_surfavg,R_supramarginal_surfavg,R_frontalpole_surfavg,R_temporalpole_surfavg,R_transversetemporal_surfavg,R_insula_surfavg,LThickness,RThickness,LSurfArea,RSurfArea,ICV' > CorticalMeasuresENIGMA_SurfAvg.csv for subj_id in $(ls -d Subj*); do #may need to change this so that it selects subjects with FS output printf "%s," "${subj_id}" >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%s," "${subj_id}" >> CorticalMeasuresENIGMA_SurfAvg.csv for side in lh.aparc.stats rh.aparc.stats; do for x in bankssts caudalanteriorcingulate caudalmiddlefrontal cuneus entorhinal fusiform inferiorparietal inferiortemporal isthmuscingulate lateraloccipital lateralorbitofrontal lingual medialorbitofrontal middletemporal parahippocampal paracentral parsopercularis parsorbitalis parstriangularis pericalcarine postcentral posteriorcingulate precentral precuneus rostralanteriorcingulate rostralmiddlefrontal superiorfrontal superiorparietal superiortemporal supramarginal frontalpole temporalpole transversetemporal insula; do printf "%g," `grep -w ${x} ${subj_id}/stats/${side} | awk '{print $5}'` >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%g," `grep -w ${x} ${subj_id}/stats/${side} | awk '{print $3}'` >> CorticalMeasuresENIGMA_SurfAvg.csv done done printf "%g," `cat ${subj_id}/stats/lh.aparc.stats | grep MeanThickness | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%g," `cat ${subj_id}/stats/rh.aparc.stats | grep MeanThickness | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%g," `cat ${subj_id}/stats/lh.aparc.stats | grep MeanThickness | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_SurfAvg.csv printf "%g," `cat ${subj_id}/stats/rh.aparc.stats | grep MeanThickness | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_SurfAvg.csv printf "%g," `cat ${subj_id}/stats/lh.aparc.stats | grep WhiteSurfArea | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%g," `cat ${subj_id}/stats/rh.aparc.stats | grep WhiteSurfArea | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%g," `cat ${subj_id}/stats/lh.aparc.stats | grep WhiteSurfArea | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_SurfAvg.csv printf "%g," `cat ${subj_id}/stats/rh.aparc.stats | grep WhiteSurfArea | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_SurfAvg.csv printf "%g" `cat ${subj_id}/stats/aseg.stats | grep IntraCranialVol | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_ThickAvg.csv printf "%g" `cat ${subj_id}/stats/aseg.stats | grep IntraCranialVol | awk -F, '{print $4}'` >> CorticalMeasuresENIGMA_SurfAvg.csv echo "" >> CorticalMeasuresENIGMA_ThickAvg.csv echo "" >> CorticalMeasuresENIGMA_SurfAvg.csv done Please let me know if you need any other info, Best wishes, Emma Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

2 1

Video of brain inflating from the anatomy
by Mageshwar Selvakumar 18 Mar '19

18 Mar '19

External Email - Use Caution Dear Freesurfer Community, We are working on a publication on chronic pain and for explanation purposes, we would like to show how the brain inflates from the anatomy. For this purpose, we thought a video on how inflation takes place could be interesting. I would like to know if there are some videos already made in this perspective and where it can be obtained. Looking forward to hearing from you. Thank you, Kind regards Mageshwar Selvakumar, M.Sc Department of Pharmakologische Bildgebung und Bildanalyse Friedrich Alexander Universität Erlangen Nürnberg Office: Fahrstr. 17, 91054 Erlangen Phone: + 49-9131-8522898

2 1

Longitudinal QDEC
by AKUDJEDU, THEOPHILUS 16 Mar '19

16 Mar '19

External Email - Use Caution Dear All, I am running a longitudinal (2-time points) group cortical analysis in Freesurfer 5.3 with Qdec. After shrinking the longitudinal QDEC table into a cross-sectional form, I added a text file containing the lines MEASURE1 = long.thickness-avg MEASURE2 = long.thickness-rate MEASURE3 = long.thickness-pc1 MEASURE4 = long.thickness-spc in the qdec directory to guide qdec to know about the newly created files (e.g. lh.long.thickness-spc...). However, it doesn't seem to be recognizing/picking up the new files or not reading the text file named Qdecrc.txt (which has the lines above. Will be glad if you could offer any help. Thanks Theo

1 0

One Group GLM Command Line Stream: correlating thickness with cognitive scores
by Arsenije Subotic 15 Mar '19

15 Mar '19

External Email - Use Caution Dear experts, I was hoping on examining the correlation between thickness and cognition within one group only and visualizing this using the command line stream. For my analysis I was planning on using cognitive scores that are adjusted according to published norms for age, sex, and education where possible. I do not want to use the scores as a covariate however, but simply find an association if possible. I'm assuming of course that I would use one of the one group FSGD formats https://surfer.nmr.mgh.harvard.edu/fswiki/FsgdExamples? Thanks! Arsenije --- Arsenije Subotic | MSc Student Department of Clinical Neurosciences Hotchkiss Brain Institute, University of Calgary | T2N 4N1, Canada Tel: (403) 918-6970 arsenije.subotic1(a)ucalgary.ca

2 1

pial surface and eTIV
by Ryan Wales 15 Mar '19

15 Mar '19

External Email - Use Caution Hello, I'm using freesurfer ver 6.0.0. If there are defects with the pial surface (i.e. the pial surface does not cover all the gray matter that it should), can you also assume that the eTIV has not been estimated incorrectly, or do pial surface inaccuracies not indicate eTIV inaccuracy? I know eTIV is estimated based on the linear transform to the talairach.xfm, but I'm not sure if there's any way to determine eTIV accuracy directly. I'm quality checking my images and am curious if I should remove subjects based on poor pial surfaces or if that does not matter. Important to note, at the moment, I'm only interested in pulling hippocampal volumes and eTIV (to get relative hippocampal volumes). Similarly, would poor pial segmentation affect the reported hippocampal volumes at all? What do you recommend? Thanks, Ryan -- Ryan Wales Graduate Student Cognition and Motor Control Neuroscience Laboratory Integrative Neuroscience Psychology Department Stony Brook University E-mail: ryan.wales(a)stonybrook.edu <ryanwales(a)stonybrook.edu>

2 1

Re: [Freesurfer] Welch's t-test for HOV violation
by Hannah CK 15 Mar '19

15 Mar '19

External Email - Use Caution Hi Doug, We actually do want the voxel-by-voxel results. Could you please give details on how we would do that? Specifically, would we need to use the command line to run the Welch's T/F or would we extract values for each voxel and put them into another software to conduct the Welch's T/F? Thank you again. Greve, Douglas N.,Ph.D. <https://www.mail-archive.com/search?l=freesurfer@nmr.mgh.harvard.edu&q=from…> Thu, 14 Mar 2019 08:45:56 -0700 <https://www.mail-archive.com/search?l=freesurfer@nmr.mgh.harvard.edu&q=date…> You can use mri_segstats to get whole-cortex means (use --slabel and specify the lh.cortex.label and add --id 1). You can do the same thing with the variances. The problem with doing it after mri_glmfit is that the variances will not reflect the spatial averaging over cortex. If you really just want one number for cortex, you should use mri_segstats to extract the mean values from the stack (ie, the --y input), using --avgwf. This will give you a column of numbers, one mean for each subject. Then you can process this however you want. When I first responded, I thought that you want to do this on a voxel-wise basis. doug On Thu, Mar 14, 2019 at 10:10 AM Hannah CK <hckalman(a)gmail.com> wrote: > Hi Dr. Greve, > > Sorry about that - I had followed the "Reply via email" link on the archive page instead of posting to the forum directly. > > Thank you for this response. Can extracting the means be done at the whole-brain level? If we get the means and variances for each group, I'm assuming we would need to use ROIs. Could you please clarify whether it's possible to extract means across the whole brain? > > Thank you, > Hannah > > > > > Re: [Freesurfer] Welch's t-test for HOV violation > 2019-03-13 Thread Greve, Douglas N.,Ph.D. > Hi Hannah, please include the previous correspondence so that we have context. > Also, please remember to post to the list and not to us personally. > thanks! > doug > > > > > On 3/13/19 10:52 PM, Hannah CK wrote: > > External Email - Use Caution > > Hi Dr. Greve, > > Thank you for this response. Can extracting the means be done at the > whole-brain level? If we get the means and variances for each group, I'm > assuming we would need to use ROIs. Could you please clarify whether it's > possible to extract means across the whole brain? > > Thank you, > Hannah > > > > ___ > Freesurfer mailing listFreesurfer@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > Re: [Freesurfer] Welch's t-test for HOV violation > 2019-03-11 Thread Greve, Douglas N.,Ph.D. > > > I used to have this about 20 years ago, but I stopped supporting it when it did not fit cleanly into the GLM ... The shortest route is probably to do two separate analyses with mri_glmfit, one for each group. This will output means and variances for each group. Then use fscalc to compute maps of the Welch's t. > > > On 3/11/19 1:28 PM, Hannah CK wrote: > > > > External Email - Use Caution > > > > We are conducting t-test comparisons across groups in which there is a > > large discrepancy in group sizes. HOV is violated. I've searched for > > ways to run Welch's t-test in FreeSurfer (or a similar analysis that > > does not assume homogeneity of variance) but am not finding one. > > > > Could anyone please advise on how to do this? > > > > Thank you. > >

2 1

Color scale is off
by Nillo, Ryan Michael R 15 Mar '19

15 Mar '19

External Email - Use Caution Hello FreeSurfer users, I am using FreeSurfer version 6 with the latest dev version of freeview on a Linux machine. When I load a surface into freeview with a color scale, I notice the scale does not correlate with the values in the surface. The min and max values are always some small number (e.g. 3.0 x 10^-6). This happens even if I set the threshold in the command line. To fix it, I would need to configure the surface and click “apply”. In the past, the color scale would format itself correctly without me having to go through that configure step. Example command line: freeview -f fsaverage/surf/ph.pial:overlay=lh.thickness.betas:overlay_threshold=0,0.05:overlay_method=linear -colorscale -ss lh.lateral.thickness.beta.tiff Thanks in advance, Ryan Nillo

1 0

← Newer
1
...
7
8
9
10
11
12
13
...
18
Older →

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 March 2019