Freesurfer

freesurfer@nmr.mgh.harvard.edu

16 participants
30084 discussions

Statistical analyses with QDEC
by Francesco Cardinale 27 Jan '16

27 Jan '16

Dear Freesurfer experts, I have some questions about statistical analyses. How to adjust for sex and age? I use the precompiled Freesurfer 5.3.0 MacOSX version, and its tool QDEC to compare cortical thickness between two groups. I would like to use age and sex as nuisance factors. No problem with the age, but I don’t see the sex in the list of usable variables if I code it as a categorical variable (male / female). On the opposite, If I code it as numbers (male = 1, female = 2) I can use it, but how can I declare that it is a factor and not a continuous numerical variable? I don’t want that “female” affects the fitting of the GLM with a value that means “male” * 2. What about the use of other morphometric properties as nuisance factors? For example, how could I use sulcality or LGI as nuisance factors? Is it possible to use more than one scalar simultaneously? Similarly to question 2, is it possible to use home-made *.mgh scalars as a nuisance factors? How? Do the command-line tools allow such analyses better than QDEC? Related to my questions, which features do they add? In this field of interest, are there any major advantages coming with Freesurfer 6.0? Thanks a lot and best regards Cico Francesco Cardinale, MD Neurosurgeon "Claudio Munari" Centre for Epilepsy and Parkinson Surgery - Ospedale Niguarda "Ca' Granda" Piazza dell'Ospedale Maggiore, 3 - 20162 - Milano - Italia <x-apple-data-detectors://0/0> phone 0039 02 64442917 <tel:0039%2002%2064442917> fax 0039 02 64442868 <tel:0039%2002%2064442868> e-mail francesco.cardinale(a)ospedaleniguarda.it <mailto:francesco.cardinale@ospedaleniguarda.it>

1 0

Re: [Freesurfer] Reg: CORRECTING DEFECT 126 (vertices=61161, convex hull=14255)
by Bruce Fischl 27 Jan '16

27 Jan '16

Hi Jesu can you please post this type of question to the FS list so others can reply and/or learn from the answer? This defect is too big - the defect correction grows quadratically in time with the number of vertices in the defect. vertices=61161 means that you have a defect that is about half the size of the total surface. You need to look at the ?h.orig.nofix and/or the ?h.inflated.nofix and figure out what the defect is and correct it by either editing the wm.mgz or adding control points. If you can't figure this out you can upload the entire subject dir to our ftp site and one of us will take a look. cheers Bruce On Wed, 27 Jan 2016, Jesu Kiran Spurgen wrote: > Dear Prof. Fischl, > Good Afternoon! I have been running freesurfer for a patient data, and have > to find this retesselalation of defect. I have successful with some data but > for some others it looks frozen even after 4 days. Could you please kindly > suggest what I could do here? How could I get some results overriding this > particular step? > > I would extremely grateful for your kind reply and help. > > Warm regards, > > > Jesu Kiran Spurgen > Postgraduate Student, Advanced Medical Imaging > KU Leuven, Belgium > >

1 0

Re mri_binarize --i sig.cluster.mgh --min 1.3 --o thresholded.mgh 1.3 = -log10(.05)
by Erik Lindberg 27 Jan '16

27 Jan '16

Dear Bruce and Freesufers, I tried to run the command and then to load it as an overlay - however I get nothing on the surface of the fsaverage brain. The program is telling me that it is including x-number of voxels into the mask but when I check I see nothing on the surface. What is going wrong here? Further the approach you proposed doing fMRI on the surface sounds interesting – but I have never used this approach. Could you give me a hint on the steps I need to take in order to do this? Thanks in advance! Eric yes, but do you want to do the fMRI analysis on the surface or in the volume? There are lots of advantages to doing it on the surface. If your sig.cluster.mgh was a surface overlay then the threholded.mgh might be exactly what you want. Load it on an inflated surface and take a look cheers Bruce

2 1

3 years Post Doc position at the Neurofunctional Imaging Group in Bordeaux - France
by Fabrice Crivello 27 Jan '16

27 Jan '16

Dear all, Please find enclose a 3 years Post Doctoral position that may be of interest. Post Doc in data processing and databasing of large neuroimaging datasets. The Neurofunctional Imaging Group (GIN) is a CNRS-CEA joint research team of the Neurodegenerative Disease Institute, UMR 5293 at Bordeaux University. The GIN is a multidisciplinary research team gathering scientists with training in instrumentation, cognitive neurosciences, signal processing and databasing. The GIN is a partnering member of the MULTI-LATERAL (Multi-level Integrative Analysis of Brain Lateralization for Language, http://flagera.eu/?q=FLAG-ERA-call-2015-projects <http://flagera.eu/?q=FLAG-ERA-call-2015-projects>) project funded by the FLAG-ERA European consortium set up to contribute to the construction of the Human Brain Project Flagship project. This partnership aims to identify the anatomical, functional and genetic determinants of brain lateralization for language functions. Within this context, the GIN offers a 3 years Post Doc / Research Engineer position immediately available. The core work of the Post Doc will to extract accurate structural and intrinsic connectivity asymmetry phenotypes across a range of large scale imaging datasets: First, the post-doc will contribute to develop improved methods and dedicated software to reliably and automatically measure individual differences in lateralization for anatomy and resting state intrinsic connectivity. In a second phase, the Post Doc will apply the methods in brain imaging datasets (about 10 000 subjects, including the BIL&GIN, BIG, I-Share and UK-biobank neuroimaging databases) having genetic data available, for the purposes of association and rare variant analysis followed by integrated genome-level analysis with transcriptomic (lateralized gene expression) dataand genomic gene-set analysis. The applicant will benefit from the support of existing research in this topic area and preliminary works in the accurate definition of both structural and resting state asymmetry characterization together with data management of large neuroimaging cohorts. The qualified applicant should have a PhD (or equivalent) in neuro-image analysis, computer science (signal and/or image processing) or related field and a strong background / experience in neuroimaging data processing. Applicants should have experience with tools for analysis of neuroimaging data (Freesurfer, SPM,…) and a relevant programming experience (Unix, C/C++ / python, MATLAB). For further information, please contact Dr. Fabrice Crivello (fabrice.crivello(a)u-bordeaux.fr <mailto:fabrice.crivello@u-bordeaux.fr>). Salary is according to Bordeaux University salary grid depending on applicant experience. The position is limited for 3 years. Job location is in Bordeaux - France. Interested candidates should send their applications including CV with publication list and recommendation letters (names and contact information) and a written summary of research interests to fabrice.crivello(a)u-bordeaux.fr <mailto:fabrice.crivello@u-bordeaux.fr>. Dr Fabrice CRIVELLO Mob : +33 (0)6 81 06 47 53 Tel : +33 (0)5 47 30 44 03 Fax : +33 (0)5 47 30 43 94 https://www.researchgate.net/profile/Fabrice_Crivello IMN, Institut des Maladies Neurodégénératives, UMR 5293 Equipe 5 : GIN, Groupe d’Imagerie Neurofonctionnelle, CEA - CNRS - Université de Bordeaux Université de Bordeaux 146 rue Léo Saignat - CS 61292 - Case 28 33 076 Bordeaux cedex http://www.imn-bordeaux.org/

1 0

recon-all: core dumped error
by Fengji Geng 26 Jan '16

26 Jan '16

Hello fs experts, We ran a command like ´recon-all -all -s HMN002 -hippocampal-subfields-T1T2 HMN002/T2/*.mgz T1T2´. Then, there was an error as ´Abort (core dumped)´ (which is also attached below). The error log file is also attached here. Error msg: ----------------------------------------------------------------------------------------------------- setting seed for random number generator to 1234 using ../mri/aseg.presurf.mgz aseg volume to correct midline $Id: mris_ca_label.c,v 1.37 2014/02/04 17:46:42 fischl Exp $ $Id: mrisurf.c,v 1.775 2015/12/10 21:16:09 fischl Exp $ reading atlas from /cell_root/software/freesurfer/dev/sys/average/lh.curvature.buckner40.filled.desikan_killiany.2010-03-25.gcs... reading color table from GCSA file.... average std = 1.0 using min determinant for regularization = 0.011 0 singular and 384 ill-conditioned covariance matrices regularized input 0: MEAN CURVATURE, flags 0, avgs 5, name mean_curvature labeling surface... 1592 labels changed using aseg relabeling using gibbs priors... 000: 2975 changed, 134786 examined... 001: 698 changed, 12707 examined... 002: 152 changed, 3873 examined... 003: 50 changed, 880 examined... 004: 29 changed, 284 examined... 005: 16 changed, 164 examined... 006: 7 changed, 87 examined... 007: 1 changed, 42 examined... 008: 0 changed, 4 examined... 206 labels changed using aseg 000: 105 total segments, 64 labels (264 vertices) changed 001: 43 total segments, 4 labels (14 vertices) changed 002: 39 total segments, 0 labels (0 vertices) changed 1 filter iterations complete (10 requested, 0 changed) Abort (core dumped) Linux hippocampus.umd.edu 2.6.32-573.12.1.el6.x86_64 #1 SMP Mon Nov 23 12:55:32 EST 2015 x86_64 x86_64 x86_64 GNU/Linux recon-all -s HMN002 exited with ERRORS at Tue Jan 26 00:23:40 EST 2016 Our freesurfer information: -------- freesurfer-Linux-centos6_x86_64-dev-20151215 -------- Setting up environment for FreeSurfer/FS-FAST (and FSL) FREESURFER_HOME /cell_root/software/freesurfer/dev/sys FSFAST_HOME /cell_root/software/freesurfer/dev/sys/fsfast FSF_OUTPUT_FORMAT nii.gz SUBJECTS_DIR /export/data/hippocampus/HMN/data MNI_DIR /cell_root/software/freesurfer/dev/sys/mni FSL_DIR /cell_root/software/fsl/new/sys Looking forward to hearing how we can resolve this error. Also, thanks in advance. Best, ----- Fengji

2 1

How to register the workshop at Boston
by Yun Wang 26 Jan '16

26 Jan '16

Dear All, I am new to Freesurfer, so I want to take the Freesurfer Workshops at Boston this April. But It appears I have passed the deadline to register. Does anyone know how to make it happen? Thanks! Yun Wang

1 0

DESPOT T1 map
by Worker, Amanda 26 Jan '16

26 Jan '16

Hi All, I have a T1 map which has been synthesised from a DESPOT imaging sequence and I am hoping to use this as input to recon-all. Do you know if this will be possible at all, or if there is anything I can do to prepare the T1 map to run through recon-all? Thanks in advance, Amanda

2 1

mri_aparc2aseg BA.annot
by Qiyuan Tian 26 Jan '16

26 Jan '16

Hi, I was trying to use mri_aparc2aseg to convert the BA.annot (or BA.thresh.annot) to segmentation volume, that will be then converted to EPI space for some diffusion data. I tried both the BA.annot for the individual subject as well as the PALS_B12_Brodmann.annot in fsaverage. The command line I used was: mri_aparc2aseg --s s1 --annot BA.annot --o BAseg.mgz. It seems that this command is working, but when I view the BAseg.mgz file in freeview, the index for the BA areas seems wrong. For example, for BA1, the index says 2004 for ctx-rh-corpuscallosum. It seems that the index for BA areas in FreeSurferColorLUT are from 400 to 414. Is there a way to make the index correct? Or mri_aparc2aseg can only be used for aparc.a2009s.annot or aparc.annot? Thanks a lot. Best QT

1 0

Qdec ROI Analysis
by Larissa McKetton 26 Jan '16

26 Jan '16

Dear Freesurfer community, I am working on an ROI analysis using qdec. I have created all aseg and aparc stats files and would like to determine if there are any differences in volume, surface area, and thickness for certain structures between two groups. I was wondering with the analysis, would it be okay to use the resulting aseg/aparc output stats from either the Desikan/Destrieux atlases and run these values through SPSS or is it better to draw ROI labels on the fsaverage brain, apply the label to each subject, and run the analyses this way. I am asking because the following tutorial states: *"Note*: If you want to use a pre-existing label and register it to fsaverage, be aware that this might involve two instances of resampling and the results might not be as accurate as they would be if you drew the label on fsaverage. Please contact the FreeSurfer team to get more details on this process if you have any concerns." http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/AnatomicalROI#Tableofth… If drawing labels on the fsaverage brain and applying them to each subject is the way to go- would I need to correct for multiple comparisons? Many thanks, Larissa

2 1

-hippocampal-subfields-T1 fails in freesurfer-Darwin-OSX-dev-20160109
by Carl Anderson 26 Jan '16

26 Jan '16

Hi Everyone, Subfields seems to fail to run when part of the recon-all -all command line and also fails when run on its own (but gives no error). Here is the output and recon-all log attached: Thanks! Carl Fitting mesh to synthetic image from ASEG took 468.3773 seconds Transforming points Transforming points 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/local/lib/libgcc_s.1.dylib dyld: Symbol not found: ___emutls_get_address Referenced from: /Applications/freesurfer/bin//../lib/gcc/lib/libgomp.1.dylib Expected in: /usr/local/lib/libgcc_s.1.dylib /bin/bash: line 1: 22256 Trace/BPT trap: 5 /Applications/freesurfer/bin//mri_convert tempT1.mgz T1resampled.mgz -odt float -rt cubic -vs 0.33333 0.33333 0.33333 > /dev/null /Applications/freesurfer/bin//mri_convert tempT1.mgz T1resampled.mgz -odt float -rt cubic -vs 0.33333 0.33333 0.33333 >/dev/null: Trace/breakpoint trap mghRead(/Applications/freesurfer/subjects/ParallelEveryone/CAT104/tmp/hippoSF_T1_v10_right/T1resampled.mgz, -1): could not open file error reading from volume T1resampled.mgz threshold mask volume at 0.5 mghRead(/Applications/freesurfer/subjects/ParallelEveryone/CAT104/tmp/hippoSF_T1_v10_right/T1resampled.mgz, -1): could not open file /Applications/freesurfer/bin//mri_mask: could not read source volume T1resampled.mgz /Applications/freesurfer/bin//mri_mask -T 0.5 T1resampled.mgz asegModBinDilatedResampled.mgz T1resampled.mgz >/dev/null: Signal 124 mghRead(/Applications/freesurfer/subjects/ParallelEveryone/CAT104/tmp/hippoSF_T1_v10_right/T1resampled.mgz, -1): could not open file error reading from volume T1resampled.mgz mghRead(/Applications/freesurfer/subjects/ParallelEveryone/CAT104/tmp/hippoSF_T1_v10_right/hippoMaskResampled.mgz, -1): could not open file threshold mask volume at 0.5 mghRead(/Applications/freesurfer/subjects/ParallelEveryone/CAT104/tmp/hippoSF_T1_v10_right/T1resampled.mgz, -1): could not open file /Applications/freesurfer/bin//mri_mask: could not read source volume T1resampled.mgz /Applications/freesurfer/bin//mri_mask -T 0.5 T1resampled.mgz hippoMaskResampledDilated.mgz T1resampled.mgz >/dev/null: Signal 124 Error using kvlAtlas3D_Matlab /autofs/space/panamint_005/users/iglesias/software/MAC_ATLAS_3D/Atlas3D_131025/itkMGHImageIO.cxx:216: itk::ERROR: MGHImageIO(0x103507280): Can't find/open file: T1resampled.mgz Error in kvlReadCroppedImage (line 11) Error in segmentSubjectT1_autoEstimateAlveusML (line 605) #------------------------------------------ Started at Tue Jan 26 09:17:04 EST 2016 Ended at Tue Jan 26 09:34:33 EST 2016 #@#%# recon-all-run-time-hours 0.291 recon-all -s CAT104 finished without error at Tue Jan 26 09:34:33 EST 2016 -- <<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>> Carl M. Anderson, Ph.D. (mail stop 204) Assistant Professor of Psychiatry, Harvard Medical School & The Neuroimaging Center, McLean Hospital 115 Mill St. Belmont, MA 02478 Ph: 617-855-2972; Fax: 617-855-2770 e-mail: remfract(a)gmail.com; carl_anderson(a)hms.harvard.edu <<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>~<>~<<•>>

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