Freesurfer

freesurfer@nmr.mgh.harvard.edu

16 participants
30085 discussions

Right-side bias in thalamic segmentation pipeline?
by Gallagher, Evan (NGG) 30 Sep '20

30 Sep '20

External Email - Use Caution Hi all, I've been working on a project that involves the thalamic segmentation pipeline published by Iglesias et al. in 2018 in NeuroImage ("segmentThalamicNueclei.sh"). While the pipeline itself works great, we have noticed that the right hemisphere VOIs tend to have slightly larger volumes than the left hemisphere VOIs. This has been true in two separate datasets (n=24 and n=25), neither of which has an obvious physiological trait that would lead us to expect these sorts of hemisphere differences. This is my first time posting a question on this list, so I'm not sure if images are OK. If yes, there are a couple posted at the bottom of this message that should visualize what I'm talking about. If not, I can say that the volume differences are most pronounced in the largest VOIs (PuM, VPL, VLp, etc); in these regions, the right side is ~200-300 cmm larger than the left side. Additionally, the average volumes on both sides are reasonable based on Figure 8 of Iglesias et al (in the sense that our average volumes are well within the violin plots in that figure). I'm wondering if anyone has encountered this trend before, and/or if there's a correction I might implement to equalize the volumes. We want to be understand-to the extent possible-the extent to which these differences are driven by biology vs image processing. Thanks so much, Evan Images follow: [cid:image001.png@01D69732.FB8CB570] [cid:image002.png@01D69732.FB8CB570] [cid:image004.png@01D69732.FB8CB570]

1 0

Question about ICV analysis !
by Martin Juneja 30 Sep '20

30 Sep '20

External Email - Use Caution Hello FreeSurfer experts, I have this interesting situation. I understand its important to include ICV as a covariate while conducting cortical volume (CV) analysis. So for my clinical sample (CL) vs healthy-controls (HCs): 1. If I do *network-wise* analysis (i.e., after parcellating the whole brain into 17 networks), I get highly significant differences in CV (MANCOVA: CL < HCs) for a specific network (say N1) - *if I include ICV as a covariate*, but not otherwise. 2. If I do *whole-brain vertex-wise* analysis, I get highly significant differences in CV (CL < HCs) for a specific cluster (which is overlapping with an area of N1) - *if I do not include ICV as a covariate*, but not otherwise. I am not sure how to interpret this i.e., ICV as a covariate for network-wise analysis plays an opposite role i.e., makes my findings stronger compared to ICV as a covariate for whole-brain vertex-wise analysis i.e., makes my findings weaker. Does ICV as a covariate add some kind of noise during whole-brain vertex-wise analysis? I would really appreciate any help in understanding this.

2 1

brainGraph (graph theory analysis in R) major release v3.0.0
by Chris Watson 30 Sep '20

30 Sep '20

External Email - Use Caution Dear Freesurfer users, (apologies for cross-posting) I am pleased to announce the 3nd major release of my R package for graph theory analysis of brain MRI data, "brainGraph". The main page on CRAN is at: https://cran.r-project.org/web/packages/brainGraph/index.html The GitHub repository is at: https://github.com/cwatson/brainGraph *Installation* The latest version is already on CRAN for Linux distributions, and should be for Windows and Mac OS X soon. To install directly from CRAN: > install.packages('brainGraph') In the meantime (and for development versions), you may install from GitHub using the R package "devtools": > devtools::install_github('cwatson/brainGraph') *Getting help* The User Guide can be found at: my GitHub Pages site (PDF warning) <https://cwatson.github.io/files/brainGraph_UserGuide.pdf>. Chapter 4 gives a general Overview, and Chapter 5 has some basic details for getting started. Chapters 6 and 7 provide code and details for structural covariance and single-subject networks, respectively. You may also get help by opening an Issue <https://github.com/cwatson/brainGraph/issues> on the GitHub repository, or you may join the Google Group <https://groups.google.com/forum/?hl=en#%21forum/brainGraph-help>. *New features* You can find the release notes/changelog in NEWS.md, at: https://github.com/cwatson/brainGraph/blob/master/NEWS.md There are too many differences to list here, but the most significant new feature is related to GLM-based analyses (including for Multi-threshold Permutation Correction, Network Based Statistic, and mediation analysis). They are now significantly faster, and (in certain situations) is faster than other solutions in R. Please let me know if you have any issues, bug reports, feature requests, etc. Chris Watson

1 0

Regressing out gender
by WON JONG CHWA 30 Sep '20

30 Sep '20

External Email - Use Caution Hello, I am currently attempting a paired t-test of subjects between two time points, regressing out the effects of age and gender. I am currently running the paired analysis tutorial, and have done the step of computing the differences between pairs using mris_preproc --paired-diff. However, I am having trouble with the FSGD and contrast files. Currently, this is my FSGD file: GroupDescriptorFile 1 Class M Class F Variables Age Input subject1pair M 80 Input subject2pair F 82 ... I am not sure how to make the contrast file for this, I would like to answer the question of " does the mean paired difference differ from 0, regressing out the effect of age and gender?" Any assistance would be greatly appreciated. William

2 1

extract morphology data
by 姜黎黎 29 Sep '20

29 Sep '20

External Email - Use Caution Hello, I am Lili Jiang from Institute of Psychology Chinese Academy of Sciences. I used freesurfer to do MRI data analysis. When I used mri_preproc to extract morphology data, I always got log file writing 'Reading source surface reg ~/sub001/surf/rh.fsaverage5.sphere.reg No such file or directory mri_surf2surf: could not read surface ~/sub001/surf/rh.fsaverage6.sphere.reg No such file or directory'. It seemed that mri_surf2surf loaded by mris_preproc always wanted to load reg file with 'fsaverage5'. But when I extrated morphology data in the resolution of fsaverage, then everything went well. Could you tell me what's wrong with this? ps. I just upgraded freesurfer to v6.0 last summer. Best, Lili

2 4

LGI value extraction
by 祝菲 29 Sep '20

29 Sep '20

External Email - Use Caution Dear expert, I have used freesurfer v6.0 to calculate local Gyrification index ( LGI ). And I have output the file named lh.pial_lgi and extracted the value of LGI by using the following command: > mris_anatomical_stats -a ${SUBJECTS_DIR}/S005/label/lh.aparc.annot -t pial_lgi -f ${SUBJECTS_DIR}/S005/stats/lh.aparc_lgi.stats S005 lh > aparcstats2table --subjects S005 --hemi lh --meas thickness --skip —parc aparc_lgi --tablefile result_lgi_lh.txt However, I found that the header of the output result showed the “thickness" instead of “LGI". Does this “thickness" represent the index of “LGI"? I have attached the result of one sample in the attachment. Thank you again for taking your precious time to solve my problem！ kind regards zhu fei

2 1

preproc-sess - motion correction problem
by Nasr, Shahin,Ph.D. 29 Sep '20

29 Sep '20

Hi Surfers, I have faced an issue with preproc-sess and I appreciate your help. I have a single slice EPI scan (number of TRs > 500; TR=0.2). When I tried to do preproc-sess, mc-afni2 is called and it crashed. When I checked the log, I noticed this: 3dvolreg.afni -verbose -dfile 005/fmcpr.mcdat -1Dmatrix_save 005/fmcpr.mat -base 005/tmp.mc-afni2.4869/tempvol.nii.gz -prefix 005/tmp.mc-afni2.4869/outvol.nii.gz 005/tmp.mc-afni2.4869/invol.nii.gz ++ 3dvolreg: AFNI version=AFNI_2011_12_21_1014 (Aug 29 2012) [64-bit] ++ Authored by: RW Cox *+ WARNING: If you are performing spatial transformations on an oblique dset, such as 005/tmp.mc-afni2.4869/tempvol.nii.gz, or viewing/combining it with volumes of differing obliquity, you should consider running: 3dWarp -deoblique on this and other oblique datasets in the same session. See 3dWarp -help for details. ++ Oblique dataset:005/tmp.mc-afni2.4869/tempvol.nii.gz is 30.100002 degrees from plumb. ++ Reading in base dataset 005/tmp.mc-afni2.4869/tempvol.nii.gz ++ Oblique dataset:005/tmp.mc-afni2.4869/invol.nii.gz is 30.100002 degrees from plumb. ++ Reading input dataset 005/tmp.mc-afni2.4869/invol.nii.gz ++ Edging: x=3 y=3 z=0 ++ Creating mask for -maxdisp + Automask has 847 voxels + 847 voxels left in -maxdisp mask after erosion ++ Initializing alignment base *** mri_3dalign_setup: cannot use nD images! ** Can't initialize base image for alignment ERROR: 3dvolreg.afni Invalid null command. Any suggestions? Please note that, this data is properly registered (manually and based on an auxiliary volume) and there is no problem there. Regards

2 1

voxel wise cortical thickness and subcortical volume
by Zeng, Qi 29 Sep '20

29 Sep '20

External Email - Use Caution Hi, Is it possible to use Freesurfer to extract global voxel-wise cortical thickness and subcortical volume? Thank you so much! Best, Qi

2 1

[WARNING: UNSCANNABLE EXTRACTION FAILED]RE: [EXTERNAL]Re: a question on Talairach transformation
by Yi Li 29 Sep '20

29 Sep '20

External Email - Use Caution Hi Douglas, I ran mri_convert --apply_transform talairach.xfm brainmask.mgz brainmask.talairach.mgz But brainmask.talairach.mgz is all black. No error was thrown during talairach transformation computation. Could you help to check what's going on here? I attached 4 relevant files in the mail. Thanks and best regards, Yi From: freesurfer-bounces(a)nmr.mgh.harvard.edu<mailto:freesurfer-bounces@nmr.mgh.harvard.edu> [mailto:freesurfer-bounces@nmr.mgh.harvard.edu] On Behalf Of Douglas N. Greve Sent: Tuesday, September 15, 2020 1:33 PM To: freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] [EXTERNAL]Re: a question on Talairach transformation Try mri_convert --apply_transform talairach.xfm brainmask.mgz brainmask.mni305.mgz On 9/15/2020 1:24 PM, Yi Li wrote: External Email - Use Caution Is there command in freesurfer that applies the linear transform (I suppose talairach.xfm) to brainmask.mgz to get the registered image? Best, Yi From: freesurfer-bounces(a)nmr.mgh.harvard.edu<mailto:freesurfer-bounces@nmr.mgh.harvard.edu> [mailto:freesurfer-bounces@nmr.mgh.harvard.edu] On Behalf Of Douglas N. Greve Sent: Tuesday, September 15, 2020 1:20 PM To: freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu> Subject: [EXTERNAL]Re: [Freesurfer] a question on Talairach transformation No, it is in native individual space. We compute a linear transform but do not apply it On 9/15/2020 1:07 PM, Yi Li wrote: External Email - Use Caution Hi, I am a beginner of Freesurfer. I would be grateful if you could help me understand the excellent software. I got brainmask.mgz at the end of autorecon1 and brain.mgz in the middle of autorecon2. Are these the images registered to MNI305 template? In other words, is my understanding correct that brainmask.mgz is the image after orig.mgz is linearly registered to MNI305 (plus normalization and skullstripping), while brain.mgz is the image after normalized image is nonlinearly registered to MNI305 (plus cropping based on brainmask.mgz and further normalization)? Thank you so much for your help in advance. Best regards Yi Li --- The information in this email, including attachments, may be confidential and is intended solely for the addressee(s). If you believe you received this email by mistake, please notify the sender by return email as soon as possible. _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer --- The information in this email, including attachments, may be confidential and is intended solely for the addressee(s). If you believe you received this email by mistake, please notify the sender by return email as soon as possible. _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer --- The information in this email, including attachments, may be confidential and is intended solely for the addressee(s). If you believe you received this email by mistake, please notify the sender by return email as soon as possible. --- The information in this email, including attachments, may be confidential and is intended solely for the addressee(s). If you believe you received this email by mistake, please notify the sender by return email as soon as possible.

2 1

mri_surf2surf --projfrac error: cannot spec both --sval-xyz and --sval
by 谢娆 29 Sep '20

29 Sep '20

External Email - Use Caution Hi, To begin with, I want clarify my goal here: I want to resample my resting-state functional data in fsaverage5 and then do some preprocessing before statistical analysis, and I get stuck in smooth now. I tried the following command: mri_surf2surf --srcsubject sub-40 --srcsurfval /home/disk2/Data/HC/HCprep/derivatives/fmriprep/sub-40/ses-1/func/sub-40_ses-1_task-rest_run-1_space-fsnative_hemi-L_bold.func.gii --trgsubject fsaverage5 --trgicoorder 5 --trgsurfval /home/disk2/Data/HC/HCprep/sub-40_smooth6_proj.mgh --hemi lh --fwhm-trg 6 --projfrac white 0.5 --srcsubject: this is the freesurfer output folder from fmriprep (in ...derivatives/freesurfer/...) --srcsurfval: I used fmriprep --output-space fsnative to resample the resting data to subject native space. --trgsubject/--trgicoorder: I want to resample the data to a fsaverage5 as I mentioned before --projfrac white 0.5: I followed your advice to smooth use midthickness coordinates, but there is an error: ERROR: cannot spec both --sval-xyz and --sval then I tried to add these two parameters: mri_surf2surf --srcsubject sub-40 --srcsurfval /home/disk2/Data/HC/HCprep/derivatives/fmriprep/sub-40/ses-1/func/sub-40_ses-1_task-rest_run-1_space-fsnative_hemi-L_bold.func.gii --trgsubject fsaverage5 --trgicoorder 5 --trgsurfval /home/disk2/Data/HC/HCprep/sub-40_smooth6_proj.mgh --hemi lh --fwhm-trg 6 --projfrac white 0.5 --sval-xyz white --sval /home/disk2/Data/HC/HCprep/derivatives/fmriprep/sub-40/ses-1/func/ but the error still exist: ERROR: cannot spec both --sval-xyz and --sval I searched a lot but cannot found any useful help with it. I also have another question about smooth: I can also smooth my data with following command: mri_surf2surf --srcsubject fsaverage5 --srcsurfval /home/disk2/Data/HC/HCprep/derivatives/fmriprep/sub-15/ses-1/func/sub-15_ses-1_task-rest_run-1_space-fsaverage5_hemi-L_bold.func.gii --trgsubject fsaverage5 --trgsurfval /home/disk2/Data/HC/HCprep/sub-40_smooth6_proj.mgh --hemi lh --fwhm-trg 6 --projfrac white 0.5 (--projfrac didn't work here too) --srcsurfval: I used fmriprep --output-space fsaverage5 to resample the resting data to fsaverage5 in here With --srcsubject fsaverage5 specified here, will every subject use the same coordinate/distance info when performing smooth? Will the command I used above performing smooth using every subjects' own coordinate/distance info? And which command is better? Any help will be appreciated, Thanks a lot! rao

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