Freesurfer

freesurfer@nmr.mgh.harvard.edu

19 participants
30088 discussions

Recon editing + FLAIRpial
by Matthieu VANHOUTTE 18 Mar '19

18 Mar '19

External Email - Use Caution Dear Experts, Having launched recon-all process on one subject with supplemental FLAIR data in order to better define pial surface, I still need to follow some recon editing during the QC of this patient. However, when re-launching recon-all to take into account manual editing do I have to precise on more time "-FLAIR FLAIR.nii.gz -FLAIRpial" or will it be implicit? Best regards, Matthieu

2 1

QDEC Longitudinal Analysis: Missing spc option from dropdown menu
by Arsenije Subotic 18 Mar '19

18 Mar '19

External Email - Use Caution Dear Freesurfer experts, I am having issues choosing the spc option for my QDEC two stage longitudinal analysis as it is not there. The spc files are all present in my participants files and the Qdecrc file is in the qdec subdirectory containing the line long.thickness-spc. Is there a way to bypass this problem somehow? Thank you, 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

3 2

Longitudinal QDEC
by AKUDJEDU, THEOPHILUS 18 Mar '19

18 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

3 3

tkregister-sess
by Nasiriavanaki, Zahra 18 Mar '19

18 Mar '19

Dear Freesurfers Hi I was checking the functional on structural data registration quality by using the below command. tkregister-sess -s $subj -fsd bold -per-run -bbr-sum I observed high values (around 0.9) for one of my subjects, and I figured out, that subject was registered to another subject's structural data by mistake. I re-did the recon (in a new folder) and checked the registration again, the values are still the same! However when I check my log files, it seems that the recon is done correctly and the snapshots from QA_tools don't look bad. I appreciate if you could please let me know what you think. Thanks a lot Mona Zahra (Mona) Nasiriavanaki Postdoctoral Research Fellow Martinos Center for Biomedical Imaging Massachusetts General Hospital 149 13th Street, 149-2615 Charlestown, MA, USA, 02129

2 2

Time-series extraction from each vortex in a ROI
by Zhi Li 18 Mar '19

18 Mar '19

External Email - Use Caution Hello Freesurfer experts, I would like to extract time series of pre-processed data from each vortex of ROIs in a annotation file, may I ask how could I achieve this? I tried 'mri_segstats' which only generated the mean time series of each ROI. I also tried Matlab, but 'read_surf' doesn't seem to work in import brain functional data (e.g., fmcpr.up.sm5.fsaverage.lh.nii.gz). If there is a way to import functional data into Matlab? And how could I use the information in the annotation file to get time series of each vortex in a ROI? Thank you very much. Best regards, Lizhi

1 0

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

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