Freesurfer July 2021

freesurfer@nmr.mgh.harvard.edu

91 participants
121 discussions

Freesurfer longitudinal stream error
by Susanne Gellißen 09 Jul '21

09 Jul '21

External Email - Use Caution Dear FreeSurfer Developers, I am processing MRI data using the freesurfer processing stream with two timepoints. I already executed the individual freesurfer-processing for each subject and timepoint and also the recon-all -base step. All finished without errors. When running the recon-all -long command, the processing stream exits at the same step in the script for all subjects and timepoints. I cannot identify a particular error message (ony at the end the script says exited with errors and no stats files are produced). I am using freesurfer 7.1.1 on an iMAc with macOS Catalina. Please find attached the recon-all log files for the longitudinal streams (2 folders with “long” in the foldername) as well as for the preparation steps. Thank you. --------------------------------------------------- PD Dr. med. Susanne Gellißen Department of Diagnostic and Interventional Neuroradiology University Medical Center Hamburg-Eppendorf Martinistraße 52, 20246 Hamburg E-mail: s.siemonsen(a)uke.de -- _____________________________________________________________________ Universitätsklinikum Hamburg-Eppendorf; Körperschaft des öffentlichen Rechts; Gerichtsstand: Hamburg | http://secure-web.cisco.com/1HvVVbdMRr6qAu-L5y0c3HUFOZFKnzwzd2mFmsc4xoj1PMv… Vorstandsmitglieder: Prof. Dr. Burkhard Göke (Vorsitzender), Joachim Prölß, Prof. Dr. Blanche Schwappach-Pignataro, Marya Verdel _____________________________________________________________________ SAVE PAPER - THINK BEFORE PRINTING

2 1

Visualizing specific cortical parcels in Desikan Killiany atlas
by Ellen Ji 09 Jul '21

09 Jul '21

External Email - Use Caution Dear experts, I have a basic visualization question. I would like to create a surface brain figure with specific cortical parcels colored. Similar to https://secure-web.cisco.com/1DryTPj2WrD-LULa8TF6HVh2xwKTLAaWDFvIpAUrs4z6O2… but not with ALL parcels in the atlas selected. I know that I can go into a subject folder and load this freeview -f surf/lh.pial:annot=aparc.annot in order to see the entire Desikan Killiany atlas. However, I would like to select 4-5 rois. What's the best way to do this? It won't be related to any results and doesn't need to be from a particular subject's brain (can also be fsaverage). It's purely to highlight some cortical parcels for visualisation purposes. Thanks, Ellen

2 1

Transformation Applied to T1.mgz during recon-all
by Alexander Rockhill 09 Jul '21

09 Jul '21

External Email - Use Caution -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA256 Hello, I am trying to replicate the processing steps that occur between when a T1 is given as input to recon-all and the T1.mgz that is stored in the $SUBJECTS_DIR/$SUBJECT/mri directory for mne-bids and dealing with mne's internal dependency on the surface RAS coordinate frame. It is my impression that 1) the image is resampled to 256 x 256 x 256 and then 2) the image is reoriented to head-first supine aka LIA aka scanner space. But when I replicate these steps using nibabel, the affine is off by a rotation and translation: T1.mgz ------ [[ -1. , -0.00000015, -0.00000003, 132.77600098], [ 0.00000003, 0.00000001, 1. , -94.00223541], [ -0.00000002, -1.00000012, -0.00000001, 149.55487061], [ 0. , 0. , 0. , 1. ]] orig resampled and translated ----------------------------- [[ -0.99995208, -0.00905184, 0.00373623, 164.44344257], [ 0.00337327, 0.03978604, 0.99920255, -98.49012464], [ 0.00919327, -0.9991672 , 0.03975359, 141.93864847], [ 0. , 0. , 0. , 1. ]] orig to T1.mgz transform affine ------------------------------- [[ 0.99995203, -0.00337331, -0.00919312, -30.68693093], [ 0.0037362 , 0.99920251, 0.03975359, -1.8476189 ], [ 0.00905186, -0.03978605, 0.99916736, 2.32735433], [ 0. , 0. , 0. , 1. ]] Since there are off-diagonal entries in the last affine, I assume there is some kind of rotation, and since there are non-zero entries in the last column, I assume there is some kind of translation as well. I am not sure what either of those would be from so any advice or help would be appreciated. Alex Rockhill PhD Candidate University of Oregon -----BEGIN PGP SIGNATURE----- Version: BCPG v1.65 iQEcBAABCgAGBQJg4hyMAAoJEOCF7T6SI2rb0VcH/2S6Bs16Bg0hreKYLugx4H2e /gr6eZ06VQa/9+JJZSmX19/eRyQLMfnK823LfUoEl9XqU7HN2PVDi7/5h+cbrw77 2qofJU2r70lADXlP9H82kf90N4enym1Yqem+q1K02QQrNqgs9kjkolQ9npiCkOv7 gqnS6DKf0zWOvptxGErzSqZxQAsyIYV2VBIctWzcUa+Sr5l49hYWg67HTz4W0C1k iLUSXnL40jk6xCGbG56y3nv44O1jDyxWQZ+5mOaDK9wANu48Ym60akTnNXY8BktQ p8KVlzw6OVC2mT89ReXwYoMtUPUIqrpMpVXxdpifqhBPwOFj/ICiaMeE/J/bMKc= =ijvn -----END PGP SIGNATURE-----

2 1

setwdelay error 1 weight
by Marilena Wilding 09 Jul '21

09 Jul '21

External Email - Use Caution Dear Freesurfers, for my FIR model I’d like to use the -setwdelay option in the mkcontrast-sess command. My command looks like this: mkcontrast-sess \ -analysis FIR_group_20-lh -contrast global-v-local \ -a 2 -c 1 \ -setwdelay When I run this the promt appears for the delay weights. But when I enter them, I always get an error: mkcontrast-sess -analysis FIR_group_20.lh -contrast global-v-local -a 2 -c 1 -setwdelay INFO: Found 2 Non-Null Conditions INFO: Found 29 Delays Condition Weights: -1.00000000000000000000 1.00000000000000000000 ################################### ### Enter 29 Delay Weights (Return) 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ERROR: entered 1 weights, needed 29 I'm entering 29 values (I want to include only the first seven TPs) but apparently only 1 value is detected. Could you please help me? I'm really out of ideas here. Thanks a lot in advance! Marilena

2 1

What is the relationship between qdec results and aparc.stats?
by Satoko KINA 09 Jul '21

09 Jul '21

External Email - Use Caution Hello Freesurfer Developers, As a result of whole-brain analysis using qdec, specific regions with significant differences between the two groups were depicted (e.g., paracentral cortical thickness). Based on this result, we compared the paracentral region on the Desikan-killiany atlas from aparc.stats between the two groups, and found significant differences not only in CT but also in items rejected by qdec (e.g., curv). How should we interpret this result? Which result is more reliable? Can we consider the result of the analysis focused on a narrow region of each ROI as a more detailed result because qdec is a whole brain analysis? I am not familiar with statistics, so please forgive me if this is a misguided question. Thank you for help! Satoko

2 1

Edit pial surface
by Rosalia Dacosta Aguayo 09 Jul '21

09 Jul '21

External Email - Use Caution Dear Freesurfer's experts, I wonder if there is any chance to correct the pial surface when it does not cover entirely certain regions of GM, such for example, orbital frontal cortex Kind regards, Rosie

2 1

Question about Extracting QDEC Results
by Sonya Gupta 09 Jul '21

09 Jul '21

External Email - Use Caution Hi, Hope you are having a wonderful evening! I was wondering, is there a way to extract results from QDEC as a NIFTI file? If not, in what format can you the extract QDEC results? Thank you! Sincerely, Sonya

2 1

Help segmenting all subjects in a folder
by g.toniol＠libero.it 09 Jul '21

09 Jul '21

External Email - Use Caution Hello FreeSurfer Developers, Thank you for your work ! I successfully set up Freesurfer and ran a test on my computer, and I am now attempting to run recon-all on all subjects in a folder ("/home/users/s/saccaro/Groupe_0" in the example script below). I run the analysis on a HPC (SLURM system). I found that, concerning the recon-all command, "If no name is given after --s, the code will be run on ALL subjects in $SUBJECTS_DIR https://secure-web.cisco.com/1cAYZmFGVfas4DCf91-Y9WZa9anAVybKAgr13rFXfcD8Bm… ." but I do not find information on how this could work also for the previous step, mri_convert: is there a way to submit a job to segment all subjects in the folder "Groupe_0", and not just S001 (in the example below)? All folders are structured as follows: /S@@@/scans/struct/@@@@@@.nii . I am new to Freesurfer and I apologize if my question is naive but I cannot find a solution in any of the online resources. Thank you very much for your time! Best, Giorgio - Script, without SBATCH settings (I join the full version): module load FreeSurfer/6.0.0-centos6_x86_64 source /home/users/s/saccaro/SetUpFreeSurfer.sh FSDIR="/home/users/s/saccaro/Groupe_0" mksubjdris ${FSDIR}/S001; mri_convert -i ${FSDIR}/S001/scans/struct/s44891-153636-00001-00192-1.nii -o ${FSDIR}/S001/scans/struct/s44891-153636-00001-00192-1.mgz recon-all -all -s S001

2 1

fcseed-config
by Nasiriavanaki, Zahra 09 Jul '21

09 Jul '21

Hello Freesurferers, I’m trying to use a functional label (extracted from task functional BOLD data), as a seed for functional connectivity. I tried to use the fcseed-config with -roi , but when I ran the fcseed-sess, I got this error: “If: expression syntax error”. Could you please let me know what the problem is? fcseed-config -roi ${label}.label -fcname ${label}.dat -fsd rest -rlf restruns -mean -cfg ${label}.config -force fcseed-sess -s ${subj} -d ../ -cfg /${label}.config Thanks, Mona Zahra (Mona) Nasiriavanaki, M.D. Postdoctoral Research Fellow Department of Psychiatry Massachusetts General Hospital 149 13th Street, 149-2615 Boston, MA, USA, 02129

2 1

Can I extract GLM values?
by Thalhammer, Melissa 09 Jul '21

09 Jul '21

External Email - Use Caution Dear FS experts, I have run mris_preproc and mri_glmfit for group analysis as described in the docs. My current project deals with the analysis of cortical thickness of preterm- vs. term-born adults and I want to correlate the difference in cortical thickness between those two groups with gene expression data. Thus, I would need the group difference ROI-wise as values in a spreadsheet/text file to perform correlation lateron. I have used the following code with the --table flag and have used cortical thickness data extracted with aparcstats2table as input. Is there a way to extract the ROI-wise group difference, corrected for the covariates included into my GLM model? Best, Melissa Thalhammer

2 1

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

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer July 2021