Freesurfer

freesurfer@nmr.mgh.harvard.edu

30132 discussions

Raw Data for Clusters
by WON JONG CHWA 17 Oct '18

17 Oct '18

External Email - Use Caution Dear FreeSurfer Experts, I am running Qdec 1.4 and I want to extract raw data (in this case, volume) from the significant clusters after Monte-Carlo simulation. In this particular case, the summary file shows 5 significant clusters after correction: # Cluster Growing Summary (mri_surfcluster) # $Id: mri_surfcluster.c,v 1.51.2.3 2012/05/31 22:10:05 greve Exp $ # $Id: mrisurf.c,v 1.693.2.7 2013/05/12 22:28:01 nicks Exp $ # CreationTime 2018/10/16-23:28:21-GMT # cmdline mri_surfcluster --in /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/sig.mgh --csd /Users/user/MyInstall/freesurfer/average/mult-comp-cor/fsaverage/rh/cortex/fwhm18/abs/th13/mc-z.csd --mask /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/mask.mgh --cwsig /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.sig.cluster.mgh --vwsig /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.sig.vertex.mgh --sum /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.sig.cluster.summary --ocn /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.sig.ocn.mgh --oannot /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.sig.ocn.annot --csdpdf /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.pdf.dat --annot aparc --cwpvalthresh 0.05 --o /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/mc-z.abs.th13.masked.mgh --surf white # cwd /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only # sysname Darwin # hostname N-Volunteer-iMac.local # machine x86_64 # FixVertexAreaFlag 1 # FixSurfClusterArea 1 # # Input /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only/qdec/Verification1/rh-Diff-SZ-HC-Intercept-volume/sig.mgh # Frame Number 0 # srcsubj fsaverage # hemi rh # surface white # annot aparc # SUBJECTS_DIR /Volumes/TRIALS/SCZ/TRIAL3.5_Males_only # SearchSpace_mm2 64328 # SearchSpace_vtx 149926 # Bonferroni 0 # Minimum Threshold 1.3 # Maximum Threshold infinity # Threshold Sign abs # AdjustThreshWhenOneTail 1 # CW PValue Threshold: 0.05 # Area Threshold 0 mm^2 # CSD thresh 1.300000 # CSD nreps 10000 # CSD simtype null-z # CSD contrast NA # CSD confint 90.000000 # Overall max 2.06685 at vertex 34785 # Overall min -4.82491 at vertex 356 # NClusters 5 # Total Cortical Surface Area 64328 (mm^2) # FixMNI = 1 # # ClusterNo Max VtxMax Size(mm^2) TalX TalY TalZ CWP CWPLow CWPHi NVtxs Annot 1 -4.825 356 2272.51 60.8 -9.4 -1.7 0.00010 0.00000 0.00020 3745 superiortemporal 2 -4.419 146292 2658.98 19.8 52.0 -5.7 0.00010 0.00000 0.00020 3104 rostralmiddlefrontal 3 -4.267 138330 1525.99 35.6 35.2 7.0 0.00460 0.00370 0.00550 1936 rostralmiddlefrontal 4 -3.304 140353 1712.65 34.9 -30.2 56.6 0.00220 0.00160 0.00280 3540 postcentral 5 -3.050 28351 1096.87 19.6 -34.8 -8.7 0.04530 0.04260 0.04800 1687 parahippocampal ------ However, when I run the command to extract the raw data from these clusters using this line: mri_segstats --i y.mgh --seg mc-z.abs.th13.sig.cluster.mgh --avgwf mc-z.abs.th13.cluster.ocn.dat --excludeid 0 The terminal window outputs this: rh-Diff-SZ-HC-Intercept-volume user$ mri_segstats --i y.mgh --seg mc-z.abs.th13.sig.cluster.mgh --avgwf mc-z.abs.th13.cluster.ocn.dat --excludeid 0 $Id: mri_segstats.c,v 1.75.2.9 2013/02/16 00:09:33 greve Exp $ cwd cmdline mri_segstats --i y.mgh --seg mc-z.abs.th13.sig.cluster.mgh --avgwf mc-z.abs.th13.cluster.ocn.dat --excludeid 0 sysname Darwin hostname N-Volunteer-iMac.local machine x86_64 user user UseRobust 0 Loading mc-z.abs.th13.sig.cluster.mgh Loading y.mgh Voxel Volume is 1 mm^3 Generating list of segmentation ids Found 4 segmentations Computing statistics for each segmentation 0 -4 6849 6849.000 1 -2 5476 5476.000 2 -1 1687 1687.000 Reporting on 3 segmentations Computing spatial average of each frame 0 1 2 Writing to mc-z.abs.th13.cluster.ocn.dat mri_segstats done ----- So I get a .dat file that only has 3 columns, which does not match up to the 5 that I originally observed after the correction. If you could let me know what the problem is that would be greatly appreciated. Thank you. William

2 1

Using a .asc file as overlay in freeview/tksurfer
by Worker, Amanda 17 Oct '18

17 Oct '18

External Email - Use Caution Hello, I am trying to overlay a .asc file in either freeview or tksurfer, but so far having no luck. As practice and to determine the right format, I am simply using lh.thickness, which I've converted to .asc using the following command: % mris_convert -c lh.thickness lh.white lh.thickness.asc I am trying to overlay this onto lh.pial for one subject, however I cannot load the file and I get the following error message: % mri_read(): couldn't determine type of file /data/project/biobank/normative_modelling/DATA/03129b39-45a0-4491-9fdf-c7062a631d19/surf/lh.thickness.asc surfer: couldn't load /data/project/biobank/normative_modelling/DATA/03129b39-45a0-4491-9fdf-c7062a631d19/surf/lh.thickness.asc. If you were trying to load a functional volume, make sure you selected the right registration method, and if necessary, that the registration file exists. If you were trying to load a volume-encoded value file, make sure it has the same number of values as this surface does vertices (128970). sclv_read_from_volume: error in FunD_New Interestingly, I can load the file as a curvature file, but then can't play around with the threshold. Do you have any idea what I'm doing wrong? Any help would be appreciated! Thanks, Amanda

3 7

license format error
by Georg Schramm 17 Oct '18

17 Oct '18

External Email - Use Caution Dear FreeSurfer Developers, I have a installed freesurfer (6.0.0-centos6_x86_64) on 2 of our linux servers. One the first server: platform: ubuntu 16.04 uname -a: Linux fermi 4.4.0-133-generic #159-Ubuntu SMP Fri Aug 10 07:31:43 UTC 2018 x86_64 x86_64 x86_64 GNU/Linux It runs fine. On the second server: platform: CentOS 7.4.1708 uname -a: Linux r23g36 3.10.0-862.11.6.el7.x86_64 #1 SMP Tue Aug 14 21:49:04 UTC 2018 x86_64 x86_64 x86_64 GNU/Linux I get the following license format error when running recon-all: -------------------------------------------------------------------------- GNU libc version: 2.17 ERROR: Systems running GNU glibc version greater than 2.15 ... -------------------------------------------------------------------------- The license files are new (downloaded this week) and exactly the same on both servers (I compared the hexdumps). I also checked that I have read permissions to the license file in $FREESURFER_HOME on the 2nd server. Do you have any idea what could go wrong on the 2nd server? Bests, Georg

1 0

create label from cluster obtained with mri_glm-fit-sim
by Ting Li 16 Oct '18

16 Oct '18

External Email - Use Caution Dear Freesurfer Experts, I really need your attention and help. Thanks a lot. We want to create label through the mri_glm-fit-sim cluster results. We have found the same topic in 2013 and tried the method you have provided. However, it doesn’t work. If you think it is not the right way, could you give me an example to do it right. I really appreciate your great help. Thanks a lot. Here is my code: mri_annotation2label --subject F999 --hemi lh --annotation cache.th20.neg.sig.ocn.annot --labelbase lh.test_label Here is my error information: subject = F999 annotation = cache.th20.neg.sig.ocn.annot hemi = lh labelbase = lh.test_label surface = white Reading surface /media/set2iscsi/ting/origin_gz_data/F999/surf/lh.white Loading annotations from /media/set2iscsi/ting/origin_gz_data/F999/label/lh.cache.th20.neg.sig.ocn.annot.annot could not read annot file /media/set2iscsi/ting/origin_gz_data/F999/label/lh.cache.th20.neg.sig.ocn.annot.annot No such file or directory INFO: could not load from /media/set2iscsi/ting/origin_gz_data/F999/label/lh.cache.th20.neg.sig.ocn.annot.annot, trying /media/set2iscsi/ting/origin_gz_data/F999/label/lh_cache.th20.neg.sig.ocn.annot.annot could not read annot file /media/set2iscsi/ting/origin_gz_data/F999/label/lh_cache.th20.neg.sig.ocn.annot.annot No such file or directory ERROR: MRISreadAnnotation() failed This is the annot files in my directory F999: ls F999/label/*.annot F999/label/cache.th20.neg.sig.ocn.annot F999/label/lh.BA_exvivo.annot F999/label/rh.aparc.DKTatlas.annot F999/label/lh.aparc.a2009s.annot F999/label/lh.BA_exvivo.thresh.annot F999/label/rh.BA_exvivo.annot F999/label/lh.aparc.annot F999/label/rh.aparc.a2009s.annot F999/label/rh.BA_exvivo.thresh.annot F999/label/lh.aparc.DKTatlas.annot F999/label/rh.aparc.annot

2 12

FreeSurfer Errors (2)
by Sparsh Jain 16 Oct '18

16 Oct '18

External Email - Use Caution Hey! I sent two responses on Friday. The second response was bigger than 1 MB so it awaited moderator approval. I don't know if you received it eventually. Following up to that mail, my reconstruction on the subject with 176 images crashed after about 32 hours of execution with the error 'could not read label' and 'too many symbolic links' in fsaverage. I have attached the entire command line output. Please help me understand what I'm doing wrong ! Thanks Sparsh

2 2

fsfast preproc-sess: viewing preproc surface outputs
by Jacob Matthews 16 Oct '18

16 Oct '18

External Email - Use Caution Hi FS Team, We are running our first fsfast pipeline, and starting with a single subject to test everything. We successfully ran preproc-sess with the following command: preproc-sess -s 11987_20180731 -surface fsaverage lhrh -mni305 -fwhm 5 -per-run -sliceorder odd -fsd bold We wanted to check the outputs. It was clear how to view most intermediate outputs, but we didn't know how to view the surface space outputs: fmcpr.odd.sm5.fsaverage.rh.nii.gz fmcpr.odd.sm5.fsaverage.lh.nii.gz Loading them into freeview as surfaces fails. We thought perhaps they needed to be loaded on top of an fsaverage surface, but didn't know where to start this process. Thanks, Jacob Matthews MRI Specialist Rotman Research Institute Baycrest Hospital 416-785-2500 ext. 3322

2 1

Midline errors
by Maxime Perron 16 Oct '18

16 Oct '18

3 15

qdec analysis with more than 2 covariates
by vittal korann 16 Oct '18

16 Oct '18

External Email - Use Caution Hi Doug, Recently I started working on freesurfer for volume based analysis. For QDEC, I need to include more than 2 variables. For example if I include patients who hail from different places like statutory town, census town and rural. In present scenario, can I include all three as covariates. Here is the error that qdec showing while running GLM: "ERROR: Subject 2 seems to have a mismatched number of discrete factors (expected 4, found 1)" "Reading source surface /usr/local/freesurfer/subjects/fsaverage/surf/lh.white ERROR: no contrasts specified." "Error in Analyze: command fai;ed mri_glmfit --y /usr/local/freesurfer/subjects/qdec/qdecAnalys/y.mgh --fsgd /usr/local/freesurfer/subjects/qdec/qdecAnalys/qdec.fsgd dods -glmdir /usr/local/freesurfer/subjects/qdec/qdecAnalys/ -surf fsaverage lh -label /usr/local/freesurfer/subjects/fsaverage/label/lh.aparc.label" Could you help us out in this regard. Would appreciate any help!! Korann

2 3

Postdoctoral positions available
by Chinappen, Dhinakaran Maheswaran 16 Oct '18

16 Oct '18

Three, NIH-supported postdoctoral positions in computational neuroscience are available immediately at Boston University in the labs of Uri Eden and Mark Kramer. The positions center on two projects: Project: Measure, model, and modulate cross-frequency coupling. Team: An interdisciplinary team of researchers at BU (Mark Kramer and Uri Eden) and the University of Minnesota (Alik Widge). Goals: (i) Develop new statistical modeling methods to assess cross-frequency coupling (CFC) in real-time, and deploy these tools in in vivo animal experiments. (ii) Develop simple, biophysical models to identify candidate mechanisms of CFC. Project: Assess human seizure dynamics across spatial and temporal scales. Team: An interdisciplinary team of researchers at BU (Mark Kramer and Uri Eden) Harvard/MGH (Sydney Cash), and Columbia University (Catherine Schevon). Goals: (i) Analyze human seizure activity in a unique data set, consisting of invasive brain voltage recordings across spatial scales - from single units to large brain areas, and temporal scales - from high-frequency oscillations to slow brain rhythms. (ii) Build biophysical models to probe how human seizures start, propagate, and stop. For either project, the successful candidate will have a strong background in statistics and computer programming; experience in interdisciplinary neuroscience research and biophysical modeling is preferred, although not essential. The positions are based at Boston University in the labs of Uri Eden and Mark Kramer. Both projects are highly interdisciplinary and will involve frequent and extensive interactions with neuroscientists at UMN (Alik Widge), or clinicians at Harvard/MGH (Sydney Cash) and Columbia University (Catherine Schevon). To apply, or for more information about the position, please send an email to Mark Kramer (mak(a)math.bu.edu) Please include an updated CV, names and contact information of three references, and a brief statement of your interest in the position and relevant qualifications. Boston University offers a vibrant community of statistical, mathematical, and neuroscience faculty and students, and many opportunities for collaboration within and outside the University. Dhinakaran M. Chinappen Clinical Research Coordinator Neurology Department Massachusetts General Hospital m: 1-267-257-1224

1 0

Postdoctoral positions available
by Manish Saggar 16 Oct '18

16 Oct '18

External Email - Use Caution [https://tr.cloudmagic.com/h/v6/emailtag/tag/2.0/1539198474/b66430ca7e69eb44…] [Apologies for cross-posting] The Brain Dynamics Lab (http://bdl.stanford.edu) at Stanford is looking to hire 1-2 postdocs with expertise in neuroimaging data analysis to develop/apply computational methods that can quantify transitions in brain activity (or connectivity) in health and disease (specifically, major depression and ADHD). We use a number of neuroimaging modalities including simultaneous EEG-fMRI and EEG-fNIRS, structural/functional connectivity, and independent fMRI/fNIRS/EEG. We also use computational modeling (e.g., The Virtual Brain) for realistic bio-physical network model simulations and neuromodulation (e.g., TMS) to gain mechanistic understandings for the observed aberrant brain dynamics in patients with psychiatric illness. Other computational tools used in the lab include Topological Data Analysis (TDA) - including Mapper and Persistent Homology - and Deep Learning (DL). Stanford University is located in sunny northern California, deep in the heart of Silicon Valley and a short train ride to San Francisco. The postdoc will be able to contribute to a number of ongoing research projects in our lab and collaborating labs by analyzing neuroimaging data from those studies. In addition, the postdoc will have the opportunity to develop their own research program with the domain of scientific or methodological problems currently being investigated within the lab (http://bdl.stanford.edu) The successful candidate should: • Have a PhD in a relevant field of study and demonstrated experience with the processing and statistical analysis of neuroimaging data. • Have strong and demonstrated computer programming skills, including Python. Experience with Nipype is preferred. • Have strong skills in statistics and machine learning (esp. TDA, DL, etc.) • Have solid experience with computational data analysis, including strong UNIX/Linux skills and experience using high performance computing clusters. Our lab is strongly committed to diversity in science, and we particularly welcome applications from members of traditionally underrepresented groups. Interested candidates should email the following to saggar(a)stanford.edu: • A cover letter describing your background and research interests • CV • Names and email addresses of two references The initial position is for one year, with the option to renew for 2-3 years, given satisfactory performance and available funding. The positions are available immediately and applicants will be considered on a rolling basis. m- Manish Saggar, Ph.D. Assistant Professor, Psychiatry Stanford University http://bdl.stanford.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