Freesurfer May 2023

freesurfer@nmr.mgh.harvard.edu

75 participants
109 discussions

Unit of the Numbers in cache.th40.neg.y.ocn.dat
by Na, Xiaoxu 17 Feb '24

17 Feb '24

External Email - Use Caution Dear Freesurfer Experts, As I'd like to extract the average value of each subject in the cluster after running group analysis, may I know the unit of these numbers? Since the one I'm running is surface area, and the cluster is 444mm^2, all the numbers in the file cache.th40.neg.y.ocn.dat is in the range of [0.55 1.05], may I know if this is correct and if so, how the unit is for those? Thanks a lot! ClusterNo Max VtxMax Size(mm^2) MNIX MNIY MNIZ CWP CWPLow CWPHi NVtxs WghtVtx 1 -4.8321 40067 444.11 -58.1 -16.6 -23.9 0.0002 0 0.0004 670 -2784.11 Regards, Xiaoxu ---------------------------------------------------------------------- Confidentiality Notice: This e-mail message, including any attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by reply e-mail and destroy all copies of the original message.

4 11

SAMSEG choice of the input
by Giulio Siracusano 10 Sep '23

10 Sep '23

External Email - Use Caution Dear Experts, I'm trying to understand what is the best input choice for running SAMSEG --recon in order to have the closest results as possible when compared against ASEG after executing (separately) an instance of recon-all all (and the raw input files in NIFTI format). I'm very well aware that ASEG and SAMSEG are different methods. I'm using the aseg.stats file as a results of the recon-all all command (here renamed as aseg.recon-all.stats for an easier identification). These stats should be my reference. Currently I've run the samseg command using the following input variants: 1) raw input file (brain MP-RAGE sequence) converted into a mgz archive (Raw.mgz). samseg --t1w ./Raw.mgz --s 002_S_0685_DIRECT --recon --refmode t1w The corresponding results are attached into (samseg.direct.stats) 2) This time I've used the file T1.mgz which is one of the intermediate results generated after running recon-all -autorecon 1. Specifically, this is the output of (mri_normalize <https://secure-web.cisco.com/1JV3nYNrViePm1qpOtzuMUtCndnVtCut-gp7OyXxNuXIOk…> -g 1 -mprage nu.mgz T1.mgz). samseg --t1w ./T1W/T1.mgz --s 002_S_0685_DIRECT_T1W --recon --refmode t1w The corresponding samseg.stats file is attached and called (samseg.direct.t1.stats). 3) this time I've used the file brainmask.mgz which is the last intermediate result generated after running recon-all -autorecon 1. Specifically, this is the output of mri_watershed <https://secure-web.cisco.com/1DjG5BDGTEbHYOWKmd4lc63SOmX8IrmziV9_NRWhpCskwO…>. samseg --t1w ./brainmask.mgz --s 002_S_0685_DIRECT_BRAIN --recon --refmode t1w The corresponding samseg.stats file is attached and called (samseg.brain.stats). I've noticed that the SAMSEG stats file (samseg.brain.stats) as obtained when using the brainmask.mgz as input file is the closest with the ASEG statistics (aseg.recon-all.stats). I also observe that the unknown volume is the lowest (# Measure Unknown, 5329.102436, mm^3) when using brainmask.mgz as input for SAMSEG, whereas is higher (Measure Unknown, 325302.894793, mm^3) when using T1.mgz as input file and is the highest (Measure Unknown, 946924.354060, mm^3) for the raw input. I know that for running recon-all we use the raw input files (for example, in NIFTI format), but I wanted to know what is the best possible configuration for SAMSEG in terms of input file in order to have the closest match between SAMSEG and ASEG statistics. Could it be possible to use the intermediate result T1.mgz or brainmask.mgz as input file for running samseg command? Is it a technically correct procedure? Do we need to use only the raw data as input for running samseg --recon command? Hope this is clear Regards GS

3 5

How is cache.th40.pos.sig.voxel.mgh computed, exactly?
by Mark Ebden 04 Aug '23

04 Aug '23

External Email - Use Caution Hi, I was wondering if anyone knows how the numbers in cache.th40.pos.sig.voxel.mgh are computed, or where I'd look to learn this. There is a one-liner in the documentation<https://secure-web.cisco.com/1R2vQZUhHj-zAreoIH8kMlOuUPh5jYx8TyPmfmDzbM1Dl-…> but I'd like the mathematical details of the correction etc. so I can determine what the final numbers (from 0 to 1) represent. I've noticed from browsing the plot in FreeView that significant areas have values close to 1, which is puzzling. Thanks if you could point me in the right direction. Kind regards, Mark Mar [cid:b7c6e65e-f85e-4206-aaa0-1ec00aa2fde2] Mark Ebden, D.Phil (Pronouns: He/Him) Research Analyst The Hospital for Sick Children 555 University Ave, Toronto ON M5G 1X8 P 416.888.0841 [cid:3b077a75-32bb-46ed-b1d2-b74379e6a8d1]<https://secure-web.cisco.com/17-WtJwO2uIWyZ7umzDRCiH0I5Q1NZ-tfYcx_zx9tznEDb…> ________________________________ This e-mail may contain confidential, personal and/or health information(information which may be subject to legal restrictions on use, retention and/or disclosure) for the sole use of the intended recipient. Any review or distribution by anyone other than the person for whom it was originally intended is strictly prohibited. If you have received this e-mail in error, please contact the sender and delete all copies.

2 6

Asymmetry analysis
by Carlos Alcocer 20 Jul '23

20 Jul '23

External Email - Use Caution Hello, I am trying to do an asymmetry analysis using steps detailed in https://secure-web.cisco.com/1E9BKFPFixB_mpRwvgfFNlEN0qYfJdo8MVLo6aeqDf-afy… However, step 1.3 starts with running: mris_apply_reg --src lh.map.mgh --trg lh.map.lh.fsaverage_sym.mgh --streg $SUBJECTS_DIR/subject/surf/lh.fsaverage_sym.sphere.reg $FREESURFER_HOME/subjects/fsaverage_sym/surf/lh.sphere.reg What is this lh.map.mgh file? and where can I find it/how can i make it? Thank you, Carlos

2 3

Use synthetic in recon-all pipeline
by Ruifeng Dong 20 Jul '23

20 Jul '23

External Email - Use Caution Hi all, I am running freesurfer 7.3.2 for segmentation of T1 data. It seems like still using the old skull stripping tools by default. How should we change it to using synthstrip? (Question 1) Also, is it a good idea to also turn on the --no-csf flag for cortex and subcortical segmentation? (Question 2) If using synthstrip, do we still need the -T2 option for better treatment of dura regions? (Question 3) Thank you. Best, Ruifeng

3 4

Generate Stats File for fornix.mgz in T1 weighted MRI ?
by Ahsan Ali 20 Jul '23

20 Jul '23

External Email - Use Caution Dear Sir/Madam, I hope you are doing well. I'm having a ADNI T1 weighted MRI database. I have run the recon-all command successfully and got a couple of default stats files in the stats folder. Then fornix.mgz is estimated using command "mri_cc -aseg aseg.auto_noCCseg.mgz -o aseg.auto_CCseg.mgz -f folder_name" successfully. Now what can be done to generate the stats for only fornix.mgz? I've been struggling with this for a long time, please help me with this. I'm using linux-ubuntu. [image: IMG-20230526-WA0000.jpg] With Regards Ahsan --------------------------------------------- Ahsan Ali Ph.D Scholar Department of Applied Mechanics and Biomedical Engineering (AMBE) Indian Institute of Technology Madras Chennai-600036, Tamil Nadu, INDIA Mob.:- +91-8899605552

2 3

Question about mri_convert --conform
by ruifeng14＠gmail.com 20 Jul '23

20 Jul '23

External Email - Use Caution Hi all, My input 001.mgz has data type short, but converted to unsigned byte. The plot shows the intensity of thee two images. How is the data conversion done? Thanks. Best, Ruifeng

4 5

Tracula statistics outputs for version 7.2
by King-Wai Chu 18 Jul '23

18 Jul '23

External Email - Use Caution Hi, Upon further examination of pathstats.byvoxel.txt and pathstats.overall.txt, I have 2 more questions. Since my previous 3 questions still have not been answered, I'm including them here so that they are all in one place. The 2 new questions are 4) and 5), my previous 3 questions are appended after. 4) I have noticed that the averages of AD, RD, MD and FA in pathstats.byvoxel.txt are exactly the same as the *_Avg_Center in pathstats.overall.txt. According to the tutorial, *_Avg_Center in pathstats.overall.txt is the "average over highest-probability path only" and the AD, RD, MD and FA in pathstats.byvoxel.txt is the value at the "position on the maximum a posteriori path". Does it mean that "highest-probability path only" and "position on the maximum a posteriori path" are the same or one of them should be changed to the other one? 5) Are the averages of AD_Avg, RD_Avg, MD_Avg and FA_Avg in pathstats.byvoxel.txt suppose to match *_Avg(average over the entire support of the path distribution) in pathstats.overall.txt? When I removed the duplicate points in pathstats.byvoxel.txt, the averages of AD_Avg, RD_Avg, MD_Avg and FA_Avg come very close to the value of *_Avg in pathstats.overall.txt I don't think I can calculate the averages by including the duplicate points in 4) but removing them in 5). Either including them in both or not at all. Questions from my previous email. 1) There are duplicate points in pathstats.byvoxel.txt, is it normal? I didn't try it on all 42 tractcs. For the few I've looked at, they all have duplicate points. I used the default setting in the config file to run trac-all. 2) When I do along-tract statistical analysis, sig.mgh and z.mgh are 1D files. When embed the sig.mgh and z.mgh statistical maps from the analysis into the template .trk (*.display.trk, which is 3D), it seems to me that sig.mgh and z.mgh don't have enough information to give *.display.trk the voxel values for the *.overlay.trk. There is no one-to-one voxel corresponding between sig.mgh and *.display.trk. What am I missing? Or what's the embedding process do? 3) When I convert .trk(*.overlay.trk, for example) to nii.gz or ascii using dmri_trk2trk or mri_convert the voxel values of the .trk file don't get converted into the .nii.gz or ascii file, only the x,y,z coordinates are converted. Is there a way to get the voxels values of the .trk file when converting to .nii.gz or ascii format? Thanks, King-Wai -- King-Wai Chu, Ph.D. Program Analyst Mental Illness Research, Education and Clinical Center (MIRECC, VISN3) James J. Peters Dept. of Veterans Administration Medical Center (OOMH) 130 W. Kingsbridge Road Bronx, NY 10468

3 3

recon-all-clinical Longitudinal workflow
by Andrew Dedinas Geoly 10 Jul '23

10 Jul '23

External Email - Use Caution Hello Freesurfer experts! I am conducting a longitudinal GMV and CT analysis of mTBI participants who received an experimental single-dose pharmacological treatment. We collected GE Bravo (0.9mm isotropic) T1w data pre, immediately after (~ 4 days after intervention) and 1-month following the intervention. I have been working with the recon-all longitudinal pipeline stream (7.1.1) but have noticed pretty significant errors in terms of skull stripping and segmentation (participants are veterans with somewhat random lesion / atrophy patterns). I have tested your new recon-all-clinical program (7.4.0) on some of my worst cases and it is performing remarkably well. I am doing some more digging into the various synth segmentation packages available in FS newest versions and was wondering what the best course of action would be for integrating recon-all-clinical or synthseg into a longitudinal workflow. 1.Could I use the recon-all-clinical outputs as the "cross" for later base and long template construction? 2. SynthSeg and SCLimbic seem to have very good test-retest reliability-- would it be feasible to use the independent segmentation extractions from those and bypass the longitudinal stream entirely (not saying I would want to but I am curious)? Any and all help here would be greatly appreciated, and I thank you very much for your time. Warmest regards, Andrew --- Andrew D. Geoly, MS Research Data Analyst 1 Brain Stimulation Lab Department of Psychiatry and Behavioral Sciences Stanford University School of Medicine http://secure-web.cisco.com/1QSBlrn3-LhYcwEz_Sjst4ZqWKpoVeamcijnxImhD8Wxefq… [cid:40f0dd9c-7a7e-4ec0-83d1-c42f63d1a100]

2 4

recon-all-clinical
by Andrew Dedinas Geoly 20 Jun '23

20 Jun '23

3 2

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer May 2023