Freesurfer

freesurfer@nmr.mgh.harvard.edu

3 participants
30135 discussions

Cannot allocate memory
by Mayer Kristina 15 May '14

15 May '14

Hi there, I am using freesurfer on MRI data of children. In one dataset I have following error appearing, I did not find described before. MRISalloc(41043678, 76): could not allocate faces Cannot allocate memory Linux cs2.kispi.int 3.13.5-103.fc19.x86_64 #1 SMP Mon Mar 3 18:46:36 UTC 2014 x86_64 x86_64 x86_64 GNU/Linux Waht do I have to check or edit manually? Thanks! Best, Kristina

2 4

Frontal WM volume
by Yang, Daniel 14 May '14

14 May '14

Dear FreeSurfer group, If I want to extract frontal WM volume, is it something doable in FreeSurfer? Much obliged, Daniel

2 3

Measuring cortical thickness at a specific location
by pfotiad＠nmr.mgh.harvard.edu 14 May '14

14 May '14

Hi FS experts, I was wondering whether it would be possible to measure the average cortical thickness at a specific slice (either axial, coronal, or sagittal). I read in Bruce's paper: "Measuring the thickness of the human cerebral cortex from magnetic resonance images" that "Measuring the thickness from the coronal slice at the point indicated by the green cross would result in an estimate in excess of 1 cm." Is that always the case, or has there been a way to bypass that issue? Thank you for your time, Panos

2 10

Registration and Surface flatten
by Caka 14 May '14

14 May '14

Dear Doug and all, I have already done two things: (1) 'recon-all’ to the segmentation for 140 subjects; (2) surface registration to fsaverage template for 1 subject using: mri_surf2surf --hemi lh --srcsubject sub001 --srcsurfval thickness --src_type curv --trgsubject fsaverage --trgsurfval ./sub001-thickness-lh.mgz Now I have got the 'sub001-thickness-lh.mgz' file. I have two questions: (1) Is this 'sub001-thickness-lh.mgz’ a surface file or a volume file? How can I read it? I used ‘freeview sub001-thickness-lh.mgz’ to read it, but a dark scene was shown, nothing can be found. I also used ‘mris_convert -c $SUBJECTS_DIR/sub001/sub001-thickness-lh.mgz $SUBJECTS_DIR/sub001/surf/lh.orig lh.thickness.asc’, an error occurred like ‘ERROR: number of vertices in /disk250/data_analysis/sub001/sub001-thickness-lh.mgz does not match surface (163842, 132230)’. (2) After the registration, I want to flatten the surface of curvature or thickness into a full plane. The wiki says we should first specify 5 vertices and cut the whole area without the midline region, then ‘kris_flatten’ can be used. Indeed, this procedure requires quite a lot manually operation. I am wondering if there are some automatically commands that can do this flatten projects for 140 subject? Any reply would be really appreciated! Qiongmin

3 2

Coomand not found
by Francesca Strappini 14 May '14

14 May '14

Hi, I'm just starting out with freesurfer. I downloaded CentOS6_x86_64 and I have Ubuntu 12.04 LTS. I have set up the licence file and configured freesurfer: tcsh setenv FREESURFER_HOME /media/Shared_Data/freesurfer source $FREESURFER_HOME/SetUpFreeSurfer.csh -------- freesurfer-Linux-centos6_x86_64-stable-pub-v5.3.0 -------- Setting up environment for FreeSurfer/FS-FAST (and FSL) FREESURFER_HOME /media/Shared_Data/freesurfer FSFAST_HOME /media/Shared_Data/freesurfer/fsfast FSF_OUTPUT_FORMAT nii.gz SUBJECTS_DIR /media/Shared_Data/freesurfer/subjects MNI_DIR /media/Shared_Data/freesurfer/mni When I try to run any type of command, like tkmedit bert orig.mgz, I get this error: command not found. Thanks for the help! Francesca

5 12

problems about extracting specific vertex thickness value and getting a average thickness for a specific mask
by charujing 13 May '14

13 May '14

Hi FS experts and users I have done the each subjects' surface regestration on fsaverage surface, and merge all of them into a *.mgh. Group analysis is also completed. Now I have two questions: 1,I want to extract the VtxMax vertex number thickness for all subjects. The VtxMax vertex number can be got in the *.sig.cluster.summary file. And it represents the most significant vertex in that cluster. I search my question in google and find these scripts which maybe helpful for me: for i in `seq 1 $numsubject` do mri_convert --frame $i y.mgh $SUBJECTS_DIR/fsaverage/surf/lh.$i.mgh mris_convert -c $i.mgh $SUBJECTS_DIR/fsaverage/surf/lh.inflated lh.$i.asc done I got the lh.*.asc file and open it. There are five columes in this file. I guess that: first colume is vertex number, the 2nd-4th columes represent MNI305 coordinates in surface, and the 5th colume is thickness of this vertex. If I want to deal with my 1st question, all I need to get the lh.*.asc file and search the specific vertex number in this file. All of above, is there anything wrong? 2, I created a surface mask, and want to extract the average thickness value in this mask for every subject. How can I complete it? Can I extract all vertex thickness in that mask, and average them in each subject? If so, however, it will be a bot of inconvenience. I want to some simple codes to deal with it. Any reply will be highly appreaciated. Thanks All the best. Rujing Zha 2014-05-10 charujing

3 4

freeview segmentation fault
by Mark D. Strohm 13 May '14

13 May '14

Hello- We’re having a problem with our Freesurfer 5.3 install: freeview exits with a Segmentation fault. No windows open, and there is no other other output. strohm::funcserv1:=>freeview Segmentation fault All I have been able to find is that it appears to be making an inappropriate ioctl call to the freeview script. This is freesurfer-Linux-centos6_x86_64-stable-pub-v5.3.0 on openSUSE 12.3. Any help would be appreciated Thank you. -Mark- 8011 open("/space/raid/fmri/freesurfer/bin/freeview", O_RDONLY) = 3 8011 dup2(3, 16) = 16 8011 close(3) = 0 8011 read(16, "#!/bin/tcsh -ef\nsource $FREESURFER_HOME/bin/tcl_setup\nsource $FREESURFER_HOME/bin/vtk_setup\nsource $FREESURFER_HOME/bin/qt_setup\nif($#argv > 0) then\n if(\"$argv[1]\" == \"-tkmedit\") then\n tkmeditfv $argv[2-$#argv]\n exit 0\n endif\nendif\nfreeview.bin $argv\n", 8192) = 259 8011 lseek(16, 0, SEEK_CUR) = 259 8011 ioctl(16, SNDCTL_TMR_TIMEBASE or SNDRV_TIMER_IOCTL_NEXT_DEVICE or TCGETS, 0x7fffe2b9a2f0) = -1 ENOTTY (Inappropriate ioctl for device) […] 8011 ioctl(16, SNDCTL_TMR_TIMEBASE or SNDRV_TIMER_IOCTL_NEXT_DEVICE or TCGETS, 0x7fffe2b9a480) = -1 ENOTTY (Inappropriate ioctl for device) 8011 write(17, "Segmentation fault\n", 19) = 19 8011 exit_group(139) = ? 8011 +++ exited with 139 +++ -- Mark D. Strohm, Ph.D. System Administrator, UCLA, Semel Institute & Department of Psychology Center for Cognitive Neuroscience (CCN) & Functional Neuroimaging Core (FUNC) 2567D Franz Hall, 1-310-206-5618, mstrohm(a)ucla.edu

1 0

distance comparisons across subjects in fsaverage group sphere
by Milde, Christopher 13 May '14

13 May '14

Dear Freesurfers, In Fischl et al. 1999 it is highlighted, that the topological structure of the original surface is preserved in the spherical space and due to the nature of a sphere, computations (distance, area) measures are easy to make. My question now, concerns the comparability of distance measures between subjects in the group fsaverage sphere. My aim is to compare distances between peak voxels coordinates of two study populations, to figure out if a systematic difference exist between the peak voxel distances of the two study populations. Can I compare the distances between peak voxels (vertices) between different subjects based on the fsaverage group sphere? So, is the comparability of distances between vertexes of individually spheres preserved in group fsaverage space? Ps: I know that the distances in spherical space are not Euclidian but still geodesic Best wishes and thanks in advance, Christopher Milde, M.Sc. Biol. Institute for Cognitive and Clinical Neuroscience Central Institute of Mental Health Square J 5 68159 Mannheim, Germany Christopher Milde, M.Sc. Biol. Institute for Cognitive and Clinical Neuroscience Central Institute of Mental Health Square J 5 68159 Mannheim, Germany Phone: +49-621-1703-6313 E-mail: christopher.milde(a)zi-mannheim.de Homepage: http://www.zi-mannheim.de/ Office: Forschungs- und Verwaltungsgebäude, Room 230

3 4

Freeview/Freesurfer Version Compatibility
by Sarah Blankenship 13 May '14

13 May '14

During the Spring Freesurfer workshop the newest version of Freeview was used as the primary viewing tool. Are there compatibility issues with using the newest version of Freeview while running analyses in an older version of Freesurfer (v5.1.0)? That is, if you run recon-all through v5.1.0, can you view and make manual edits using the newest release of Freeview (11 May 2014 build) or should you only use the Freeview build that is specific to your Freesurfer version (22 May 2011 build for 5.1.0)? Thank you, Sarah Blankenship Graduate Research Assistant Neuroscience and Cognitive Science Program University of Maryland, College Park

2 1

Re: [Freesurfer] nonlinear models in longitudinal analysis
by Jodie Davies-Thompson 13 May '14

13 May '14

Hi Martin, Thanks for your reply. I'm not entirely sure what you mean (sorry - I'm new to this!). Do you mean something like this? For examining an increase between timepoint 1 and 2, but the same between timepoint 2 and 3? fsid fsid-base years years-to-years Day1 base 0 0.4 Day2 base 0.4 0.4 Day3 base 0.8 0 All the best, Jodie Message: 6 Date: Fri, 09 May 2014 16:52:25 -0400 From: Martin Reuter <mreuter(a)nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] nonlinear models in longitudinal analysis To: freesurfer(a)nmr.mgh.harvard.edu Message-ID: <536D4009.7030907(a)nmr.mgh.harvard.edu> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Hi Jodie, 'years' should contain the real time between time points (usually, unless you have test-retest where time is irrelevant, or different task that can be ordered as you like). For what you are trying to do you should use linear mixed effects models with a piecewise linear setup. This will allow you to estimate a slope before and after timepoint 2 and compare those slopes. Your design matrix would probably have the column of 1, years and a column with years-years(at tp2 of that subject) if that value is positive 0 otherwise Testing if the coef of the 3rd column is > 0 is asking if the slope increases in the second interval. Not exactly what you want, but close and of course you can test if the slope in interval 1 is different from zero (second column) and if the slope in interval 2 is different from zero ( 0 1 1, ie the sum of the last two columns). Best, Martin On 05/08/2014 03:48 PM, Jodie Davies-Thompson wrote: > Dear Freesurfer experts, > > I am attempting to do a longitudinal analysis of cortical thickness with 3 timepoints. So far, I have successfully completed the analysis for long_mris_slopes comparing timepoint 1 to timepoint 2, timepoint 2 to timepoint 3, and a linear increase for timepoint 1 to 2 to 3. > > However, I am primarily interested in finding regions which show a) an increase between timepoint 1 and 2, but no difference between 2 and 3, and b) no difference between timepoint 1 and 2, but an increase between 2 and 3. > > Would the correct way of 'modelling' these two in the qdec file be: > > a) > fsid fsid-base years > day1 base 0 > day2 base 3 > day3 base 3 > > and b) > fsid fsid-base years > day1 base 0 > day2 base 0 > day3 base 3 > > All the best, > Jodie > > - > --- > Jodie Davies-Thompson, Postdoctoral Fellow > Department of Ophthalmology & Visual Sciences > UBC/VGH Eye Care Centre > 2550 Willow Street > Vancouver, BC, V5Z 3N9 > Canada > Tel: 604-875-4111 ext 69003 > - --- Jodie Davies-Thompson, Postdoctoral Fellow Department of Ophthalmology & Visual Sciences UBC/VGH Eye Care Centre 2550 Willow Street Vancouver, BC, V5Z 3N9 Canada Tel: 604-875-4111 ext 69003

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