Freesurfer

freesurfer@nmr.mgh.harvard.edu

16 participants
30084 discussions

Native volume to fs_average surface
by Mason Wells 10 May '21

10 May '21

External Email - Use Caution Hi experts, I have some fMRI data that I have pre-processed in Feat. I want to resample a single subject's data to the fs_average surface but unsure of the most optimal way of doing this. At present I am warping the functional data to the MNI brain in the volume (using FSL's applywarp) and then resampling to fs_average using the following command. * mri_vol2surf --src filtered_func_MNI.nii --out filtered_func_surf.mgh --hemi rh --mni152reg My question is whether you think this is the most appropriate approach, and/or if there is way to move directly from native volume to fs_average surface without needing to go via MNI volume? Best wishes, Mason Mason T Wells, MSc PhD student School of Optometry and Vision Sciences & Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology Cardiff University Cardiff CF24 4HQ UK Email: WellsMT1(a)cardiff.ac.uk<mailto:WellsMT1@cardiff.ac.uk> Tel: 02920 879628 Web: Cardiff University webpage<https://secure-web.cisco.com/1qGEoQnnU3besQKg-GQwdKq-WtCCzi1b5tl-HNLywYRQ48…><http://secure-web.cisco.com/1Mj4YBzh-BsbinaViI1QJ6drLVkixPOlMFL2tzAILihwp7N…> Mason T Wells, MSc Myfyriwr PhD Yr Ysgol Optometreg a Gwyddorau’r Golwg & Canolfan Ymchwil Delweddu’r Ymennydd Prifysgol Caerdydd (CUBRIC), Yr Ysgol Seicoleg Prifysgol Caerdydd Heol Maindy Caerdydd CF24 4HQ DU E-bost: WellsMT1(a)caerdydd.ac.uk<mailto:WellsMT1@caerdydd.ac.uk> Ffôn: 02920 879628

2 7

Cortical surface analysis
by Ornaith O Reilly 10 May '21

10 May '21

External Email - Use Caution Hello Freesurfer developers, I am new to Freesurfer and I'm hoping that someone can help me. My goal is to be able to visualise electrodes contacts that are on lying on the cortical surface on a flat map and a spherical map. I've run recon-all and created the flatmap. I'm using Matlab to achieve this and I've loaded the information with the lines lhsphere = fs_read_surf( '/Applications/freesurfer/7.1.1/subjects/IM-0001/surf/lh.sphere'); lhflatpatch=fs_read_patch( '/Applications/freesurfer/7.1.1/subjects/IM-0001/surf/lh.full.flat.patch.3D' ); 1. I want to be able to access the coordinates and show the electrons on the flat map/ sphere and I was wondering how to do this? I understand that the vertex ID number is what relates the coordinates in the original image to these surfaces. I am wondering how to access the vertex ID numbers ? At the moment I think that I only have access to the vertex coordinates and the vertex faces. 2. Additionally, I was wondering if there is a way to display the surface images through Matlab? What would be the best way to show the surface with the electrons overlaid on it? The version of Freesurfer and the platform that I'm using is outlined below: 1) FreeSurfer version: freesurfer-darwin-macOS-7.1.1-20200811-8b40551 2) Platform: macOS High Sierra 10.13.6 (17G7024) I really appreciate any help with this as I'm really stuck. Thanks in advance for any help. Kind Regards, Ornaith O'Reilly

3 10

Cannot select ROI edit tool in Freeview
by Mercy Mazurek 10 May '21

10 May '21

External Email - Use Caution I loaded my brainmask.mgz file generated from a "recon-all -autorecon1" command. When I open up Freeview, I see the icon but it is grayed out and I cannot select it to create a volume ROI on my image. How do I draw an ROI on my scans? Thanks for your help!

2 1

GM Volumes
by Nils B. 10 May '21

10 May '21

External Email - Use Caution Dear Freesurfer developers, I'm hoping you can help me once again. Right now I'm working on how to extract certain measurements out of my data set. I have two groups and MRI data from one time point in elderly people. The data quality is pretty bad but for pilot analysis I want to analyze the gray matter volume difference in some ROI's. I've already found the total GM volume in the aseg-stats as well as the hippocampal volumes but I couldn't find any of these ROI's: subiculum, parahippocampal gyrus,Central sulcus, Cingulate cortex, left frontal cortex, central premotor cortex, ventral premotor cortex. Are any of these regions found in the Freesurfer output or is there a way to generate them? Has Freesurfer even a way to automatically identify these regions? My Set-up: Freesurfer Version: 7.1 Plattform: Ubuntu 18.04 Thank you in advance for the help! Best regards, Nils -- Nils Böer Research Assistant with Bachelor Degree (WHB) *Psychology and Movement Science* University of Paderborn Department of Sports and Health Warburger Str. 100 33098 Paderborn Room: Sp 1.401 Email: nils.tobias94(a)gmail.com / nboeer(a)mail.uni-paderborn.de <nils.tobias94(a)gmail.com> Web: https://secure-web.cisco.com/1WVj0EcABQ4b3chjPYhYGhDOux-whaYVbSMnsWR9otFaBI… https://secure-web.cisco.com/1KCbyE0u2DHoajFIXnQlic2g_jJy0p9EnLmzpN0VFe-d7c… <https://secure-web.cisco.com/1WVj0EcABQ4b3chjPYhYGhDOux-whaYVbSMnsWR9otFaBI…>

2 1

Extract ROI cortical thickness with mri_segstat
by LIN, Yezhe 10 May '21

10 May '21

External Email - Use Caution Dear Freesurfer experts, My name is Yezhe Lin. I am trying to extract the thickness in an area we found significant in group analysis. However, when we did a t-test on the extracted values, the significance was not obtained. I wonder something could go wrong somewhere and what possibly caused the difference. I attached my scripts for the adjusted pipeline FYI. steps #make a 3x3x3 ROI using the peak coordinate obtained in group analysis in standard space (MNI152) by fslview #check ROI location tkmedit -f /usr/local/fsl/data/standard/MNI152_T1_2mm_brain.nii.gz -overlay MNI152_T1_2mm_rh.thickness.CBPHCagesexeTIVmask.nii -fthresh 0.5   # register to fsaverage of the dataset cd ../fsaverage/surf/ (Freesurfer_Home fsaverage) fslregister --s fsaverage --mov /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_brain.nii --reg MNI2mm_to_fsaverage.dat   #To check results, run: tkregister2 --mov /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_brain.nii --reg MNI2mm_to_fsaverage.dat --surf orig   #view your ROI in the fsaverage subject volume to ensure that it is located where you expect it to be: (two pials clear) cd $SUBJECTS_DIR/fsaverage/surf  tkmedit fsaverage T1.mgz \   -overlay /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_rh.thickness.CBPHCagesexeTIVmask.nii \   -overlay-reg MNI2mm_to_fsaverage.dat \   -fthresh 0.5 \   -surface lh.white -aux-surface rh.white #register ROI to fsaverage space cd $SUBJECTS_DIR/fsaverage/surf mri_vol2surf \   --mov /mnt/backup2/restore-2021/CBPFreesurfer/FreesurferResults/TP1/FS/mask/MNI152_T1_2mm_rh.thickness.CBPHCagesexeTIVmask.nii \   --reg MNI2mm_to_fsaverage.dat \ --projdist-max 0 1 0.1 \   --interp nearest \   --hemi rh \   --out rh.fsaverage.thickness.CBPHCagesexeTIVmask.mgh   #view the ROI registration (on the peak center of significant cluster) tksurfer fsaverage rh inflated -overlay rh.fsaverage.thickness.CBPHCagesexeTIVmask.mgh -fthresh 0.5 # extract values from rh.XXX.CBPvsHC_age_sex_eTIV.10.mgh mri_segstats --i ../rh.thickness.CBPvsHC_age_sex_eTIV.10.mgh --seg rh.fsaverage.thickness.CBPHCagesexeTIVmask.mgh --excludeid 0 --avgwf CBPHC_th_results.dat Hope my questions is clear! Appreciate your kind help and time! ------------------ Best regards, Yezhe Yezhe Lin, MD MMed  (she/her/hers) Postdoctoral Associate, Department of Psychiatry University of Rochester Medical Center yezhe.lin.bee(a)gmail.com Twitter: @LinYezhe | LinkedIn | Researchgate  

2 1

average thickness lobe stats
by Agurne Sampedro Calvete 10 May '21

10 May '21

External Email - Use Caution Dear Freesurfer experts, I am trying to obtain the average thickness of the brain lobes. I am particularly interested in the temporal lobe. I have run "mri_annotation2label" (--lobeStrict) and "mris_anatomical_stats" commands and I have obtained the average thickness of the different lobes in each hemisphere. However, I am interested in obtaining data of both hemispheres together, that is, a single thickness average for example for the temporal lobe. Is there a way to calculate it? Could I just sum both averages from each hemisphere? Thank you very much in advance. Kind regards, Agurne

2 1

Question about hippocampal sub-volume extraction
by Wei Shao 07 May '21

07 May '21

External Email - Use Caution Hello,Dear Freesurfer exports I have a question about the pipeline of hippocampal sub-volume extraction: quantifyHAsubregions.sh hippoSf <T1>-<analysisID> <output_file> <OPTIONAL_SUBJECTS_DIR> Let's say I have two different groups of analysisID such as sub-AVP???? and sub-BAY???? I can pass a wildcard to extract the data of T2 scan from all sub-BAY???? subjects by quantifyHAsubregions.sh hippoSf BAY*-T2 result.txt But Is there any way that I can extract the sub-volume from all the two different group subjects at one time? Thanks. Best, Wei

3 4

How white matter voxel is represented on the surface ?
by John Gelburg 07 May '21

07 May '21

External Email - Use Caution Dear list, I am sorry for a newbie question: Suppose I analyze the fMRI data in a volumetric way using SPM. I have a voxel that appears in a white matter. Does this mean that this voxel will not appear on the surface? In other words, there would be no corresponding vertex for my voxel ? I would also appreciate any reading suggestions. Thank you, John

2 1

mri_label2vol error
by Melissa Marie Rundle 07 May '21

07 May '21

External Email - Use Caution Hi, I get the following error message when I run mri_label2vol (freesurfer v 6.0) in a jupyter notebook(6.1.4) with python (v3.8.5): dyld: lazy symbol binding failed: Symbol not found: ___emutls_get_address Referenced from: /Applications/freesurfer/bin/../lib/gcc/lib/libgomp.1.dylib Expected in: /usr/lib/libSystem.B.dylib However, when I run mri_label2vol directly in a terminal window, the command runs just fine. Nor do I have trouble running similar commands such as mri_vol2vol or mri_convert using the jupyter notebook. I tried running "sudo -E fs_update" as suggested in https://secure-web.cisco.com/1fu8uZTsG4VqEpMIJBtTitHsB5L2G2xXvyt_Drs1gEwm3W… but that didn't solve the problem. Any suggestions would be greatly appreciated. Sincerely, Melissa Rundle

2 1

Freesurfer 7 crashes during mris_fix_topology
by Yakupov, Renat /DZNE 07 May '21

07 May '21

External Email - Use Caution Dear all, I have a subject whose Freesurfer processing crashes during topology fixing step. The last lines in `recon-all.log` are: After retessellation of defect 2 (v0=22890), euler #=-15 (110444,329616,219157) : difference with theory (-18) = -3 CORRECTING DEFECT 3 (vertices=8808, convex hull=2020, v0=31369) Segmentation fault (core dumped) Linux compute-0-2.local 2.6.32-504.16.2.el6.x86_64 #1 SMP Wed Apr 22 06:48:29 UTC 2015 x86_64 x86_64 x86_64 GNU/Linux recon-all -s 81489f2bb-M12_T1_01 exited with ERRORS at Tue Apr 27 21:23:59 CEST 2021 For the full `recon-all.log` and `recon-all.error` files please see the attached files. I tried to follow the manual topology fixing instructions from https://secure-web.cisco.com/1CVb3gwDr9GTwil5kr7FQYCZpGmZhLu9ovxhkC6nEa1QhF…, but *.white and *.pial surface files are missing from `surf` directory. It only contains these files: lh.defect_borders lh.defect_labels lh.orig.nofix lh.smoothwm.nofix rh.orig.nofix rh.smoothwm.nofix lh.defect_chull lh.inflated.nofix lh.qsphere.nofix rh.inflated.nofix rh.qsphere.nofix Which files can I use in place of white and pial surface files in order to look for WM segmentation defects? Some additional information about the data: - I am running Freesurfer 7.1 - overall, there are no gross errors in WM segmentation, although the subjects does have significant temporal lobe atrophy - there are two images acquired at different times, and processing fails for both images - interestingly, processing didnt crash with Freesurfer 6. Thank you and have a great weekend, Renat.

2 5

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