Freesurfer

freesurfer@nmr.mgh.harvard.edu

1 participants
30133 discussions

Is the average surface in MNI space?
by Ling Zhao 10 Jan '17

10 Jan '17

Dear experts, 1. From the FreeSufer wiki we all knows all of the processing and the generated files are in the Talairach (MNI305) coordinate. The FreeSurfer operates in a conformed space, which is different from the original structural image space that it has received as an input. If I want to use the generated surfaces or labels as masks (which is in Talairach space) for DTI probtrackx (which is need in diffusion space). How should I do? 2. In your wiki I see the information as following: “Average Surface Construction The averagesurface is constructed by computing the Talairach (MNI305) coordinate at each vertex for each subject (based on the talairach.xfm matrix for that subject). These coordinates are then mapped into the common surface space and averaged together so that the coordinate of a vertex in the common surface space is the average of the Talairach coordinates from the corresponding surface locations of the individuals. This results in the averagesurface overlaying nicely on the average (MNI305) volume. ” Is the MNI305 in the MNI space? Can you tell me the meaning of the MNI305 in “Talairach (MNI305) coordinate”? After the averagesurface overlaying nicely on the average (MNI305) volume. Is the averagesurface in the MNI space or in the Talairach space / common surface space? Best wishes Ling Zhao -- Center for the Study of Applied Psychology, Key Laboratory of Mental Health and Cognitive Science of Guangdong Province, School of Psychology, South China Normal University, Guangzhou 510631, P. R. China

2 1

reposition surface for pial surface
by Antonin Skoch 10 Jan '17

10 Jan '17

Dear Paul, I could chime in with my current limited experiences with this tool (using V6beta version of freeview). Maybe they could be helpful for you. Most straigthforward is using "edit coordinate": Select the surface you want to edit. Select tools-reposition surface - edit coordinate. With shift + left mouse button click on the surface to pick particular vertex. Manually enter new coordinates of the vertex. Save surface to store edits. This is prohibitively time consuming in case of large error, so it is usable only in case of small errors. However, you can try to use combination of "edit coordinate" and "smooth". By smoothing you can nudge some area of adjacent vertices to desired position. I was partially successful also with using reposition surface - volume based - coordinate: With shift + left mouse button select particular vertex In "volumes", load volume according to what you want to reposition the surface. With shift + ctrl + left mouse button click at the coordinate you want to reposition the surface to. Save surface to store edits. Regarding the parameters: The "size" textfield determine the extent of the vertex neighbourhood which would be repositioned. Do not use too large value since then the behaviour is unpredictable, in my case the maximal resonable value is around 5. I am not sure about "sigma" parameter, maybe it has to do something with the smoothing (of the underlying volume?), however in my case altering this parameter did not change behaviour much. Also, the "gradient direction" in my case did not change much, in many cases it produced identical result regardless of the gradient direction selected. In most cases, I was not satisfied with the result of the reposition, the particular vertex moved to desired position, however position of adjacent vertices did not make sense, mostly the surface got very jagged and irregular, with worsening of the position in areas where the position was good previously, etc. However, subsequent using of "smooth" usually improved the shape and after many iterations I finally arrived to desired shape of surface. I also tried to study the relevant source code to learn the behavior of this tool and the significance of the parameters, but it is very complicated to me to comprehend (so far..). Note that after using this tool you should rerun remaining code of recon-all ReconAllDevTable https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllDevTable (for development version) would tell you which particular steps should be rerun. Please note that in case of rerun of recon-all from start you will lose all modifications of surfaces you achieved using this tool. And, particularly, many surface errors can be corrected by editing the volume (particularly wm.mgz, brainmask.mgz, aseg.mgz) or adding control points to the white matter. Large errors of ?h.white (and subsequently ?h.pial) can be caused of wrongly corrected topological defect (which you can diagnose by inspecting ?h.orig.nofix and ?h.orig). In that case editing of wm.mgz or aseg.mgz would help to correct the issue. I also came across the error when ?h.white was cut off the particular gyrus near midline, which was caused by erroneous labeling of cortex in aseg.mgz (gyrus was attributed to contralateral hemisphere). Editing of aseg.mgz by assignment of the gyrus to correct hemisphere corrected the error. Also, altering several of global parameters of mris_make_surfaces, as listed for example here http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg37990.html could help you to get satisfactory results (however, the altering of these parameters should be done in my opinion for all subjects in the study to avoid bias). Regards, Antonin Thanks bruce for your response. Since there is no documentation, how can i use it because I seen trends on the forum where it was suggested to use reposition surface ( e.g. https://mail.nmr.mgh.harvard.edu/pipermail/freesurfer/2012-October/026117.h…) Use = directions and parameters to input. Thanks Best, PaulOn Tue, Jan 10, 2017 at 12:42 PM, Bruce Fischl <fis...(a)nmr.mgh.harvard.edu> wrote: > Hi Paul > > sorry, this is something that Ruopeng put together for me to mess around > with - there isn't any documentation > > cheers > Bruce > > On Tue, 10 Jan 2017, miracle ozzoude wrote: > > Hello Freesurfer, I would like to learn how to use the reposition surface >> tool in freeview. Specifically, what parameters to enter in the gui (see >> screenshots). Is there any documentation on how to use it if not, can >> someone with experience work me through it. Thank you. >> best, >> Paul

2 1

Help - Error: cannot unpack mosaics without ASCII header
by Caleb Comoglio 10 Jan '17

10 Jan '17

Hello, I am having trouble with the MRI_convert portion of recon-all. I get the error below when I run recon-all and when I run MRI_convert independently of recon-all. WARNING: file Desktop/3T_20161201_131733544/MR000004 does not contain a Siemens ASCII header has this file been anonymized? ERROR: cannot unpack mosiacs without ASCII header Prior to getting the error I get the following message: Scanning Directory INFO: Found 2316 files in Desktop/3T_20161201_131733544 INFO: Scanning for Series Number 2 INFO: found 192 files in series INFO: loading series header info. INFO: sorting. RunNo = 1 WARNING: Run 1 appears to be truncated Files Found: 192, Files Expected (lRep+1): 1 Do you have any suggestions? Additional system info is below: 1) OS: MacOS Sierra (10.12.2) 2) FS software version: freesurfer-Darwin-lion-stable-pub-v5.3.0 3) Data is MPRAGE directly from a siemens machine *Caleb C. Comoglio, BSBME* Graduate Research Fellow Bioelectronics Laboratory, IUPUI ccomogli(a)iupui.edu (574) 727-1620

2 1

Multiple comparison correction on subcortical volume
by neuroimage analyst 10 Jan '17

10 Jan '17

Hi, I extracted thalamic, caudate, and putamen volume (left hemisphere), and cingulate thickness (right hemisphere) from freesurfer's parcellation. Then I compare these volumes and thickness across two groups using R and found significant differences in these measurements. Do I have to correct for multiple comparisons? If yes, how do I do that? Does this make sense to divide the p-value/4 and if the p-value < p-value/4 then the test is Bonferroni corrected? If no, what do I say in my results that the test statistic is not inflated? Thanks Regards --VM

2 1

autorecon1 crash (with nu_correct).
by Jamie Hanson 10 Jan '17

10 Jan '17

Hello Freesurfer List, I was running into issues processing subjects, with errors just getting through the autorecon1 step of the pipeline. My terminal outputs the following: ==== /home/jamielh/Volumes/Hanson/Pardini_Projects/K_award/10007/mri /usr/local/pkg/freesurfer-dev//bin/mri_nu_correct.mni --no-rescale --i orig.mgz --o orig_nu.mgz --n 1 --proto-iters 1000 --distance 50 nIters 1 $Id: mri_nu_correct.mni,v 1.27 2016/02/26 16:19:49 mreuter Exp $ Linux pfc 4.4.0-53-generic #74-Ubuntu SMP Fri Dec 2 15:59:10 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux Tue Jan 10 12:53:44 EST 2017 Program nu_correct, built from: Package MNI N3, version 1.10, compiled by nicks@minerva (x86_64-unknown-linux-gnu) on 2010-02-20 at 17:32:37 /usr/bin/bc tmpdir is ./tmp.mri_nu_correct.mni.129237 /home/jamielh/Volumes/Hanson/Pardini_Projects/K_award/10007/mri mri_convert orig.mgz ./tmp.mri_nu_correct.mni.129237/nu0.mnc -odt float mri_convert.bin orig.mgz ./tmp.mri_nu_correct.mni.129237/nu0.mnc -odt float $Id: mri_convert.c,v 1.226 2016/02/26 16:15:24 mreuter Exp $ reading from orig.mgz... TR=0.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, -1.32349e-23, -1.32349e-23) j_ras = (0, 1.49012e-08, -1) k_ras = (0, 1, 1.49012e-08) changing data type from uchar to float (noscale = 0)... writing to ./tmp.mri_nu_correct.mni.129237/nu0.mnc... -------------------------------------------------------- Iteration 1 Tue Jan 10 12:53:47 EST 2017 nu_correct -clobber ./tmp.mri_nu_correct.mni.129237/nu0.mnc ./tmp.mri_nu_correct.mni.129237/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.129237/0/ -iterations 1000 -distance 50 Can't use 'defined(@array)' (Maybe you should just omit the defined()?) at /usr/local/pkg/freesurfer//mni/bin/nu_estimate_np_and_em line 165. nu_correct: crashed while running nu_estimate_np_and_em (termination status=65280) ERROR: nu_correct Linux pfc 4.4.0-53-generic #74-Ubuntu SMP Fri Dec 2 15:59:10 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux recon-all -s 10007 exited with ERRORS at Tue Jan 10 12:53:47 EST 2017 For more details, see the log file /home/jamielh/Volumes/Hanson/ Pardini_Projects/K_award/10007/scripts/recon-all.log To report a problem, see http://surfer.nmr.mgh.harvard. edu/fswiki/BugReporting ==== Searching the FS list, I saw that I originally thought this was an incompatibility with the current stable Freesurfer version and the most recent MNI tools. To troubleshoot that issue, I installed the dev version of Freesurfer (as noted by this thread: https://mail.nmr.mgh. harvard.edu/pipermail///freesurfer/2016-September/047734.html) The error however still persists (even using the dev version). Anyone else run into this issue or have thoughts about troubleshooting it? Thanks much! Best, Jamie. -- Jamie Hanson Assistant Professor, Psychology Research Scientist, Learning Research and Development Center University of Pittsburgh Personal website: jamiehanson.org

2 1

Re: [Freesurfer] Brighter image at baseline vs. follow up using longitudinal processing stream
by Tamara Tavares 10 Jan '17

10 Jan '17

Hi Dr. Reuter, Thank you for the clarification. The differences I am seeing are also in the T1.mgz file in the longitudinal run folder. I assume it is okay to use the T1.mgz from the longitudinal outputs as the original image as I can easily compare the same slice across different time points? Also, on a related note, I have two scans (time 1 and time 2) that doesn't seem to align perfectly with the template or with each other. When I use the opacity slider to transition from the template to time 1 it looks like the brain is 'moving in'. When I transition from template to time 2 the brain seems to 'pop out'. Similarly, when I transition from time 1 to time 2, it looks like the brain is popping out. With my other images, there is a smooth transition from time 1 to time 2 and a relatively smooth transition from the template to the time points. Do you think this is problematic? Thank you again for your help. Best, Tamara Message: 4 Date: Wed, 28 Dec 2016 11:20:36 +0100 From: Martin Reuter <mreuter(a)nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] Brighter image at baseline vs. follow up using longitudinal processing stream To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu> Message-ID: <F2352B06-D00A-409A-9FB3-7AE1FD6EA20B(a)nmr.mgh.harvard.edu> Content-Type: text/plain; charset=us-ascii Hi Tamara, three out of how many? I would probably remove those three as whatever you measure in them could be driven by the acquisition differences, but it depends on what acquisition differences you have. If it is only small differences, like tiny deviations in the voxel size etc, it should be OK. I have no idea where the other intensity differences could come from. If you see them in the original inputs, maybe an MR physicist could take a look and advise. If you have a concrete case where the inputs are OK, but you think that in FS processed images show different contrast or intensities, I would need to have those images to see if I can replicate that. Best, Martin

2 1

reposition surface for pial surface
by miracle ozzoude 10 Jan '17

10 Jan '17

Hello Freesurfer, I would like to learn how to use the reposition surface tool in freeview. Specifically, what parameters to enter in the gui (see screenshots). Is there any documentation on how to use it if not, can someone with experience work me through it. Thank you. best, Paul

2 2

ERROR: fast_selxavg3() failed\n
by Michael Schirner 10 Jan '17

10 Jan '17

Dear experts, I'm attempting to perform step 6 in the FsFastTutorialV5.1 (https://surfer.nmr.mgh.harvard.edu/fswiki/FsFastTutorialV5.1/FsFastFirstLev…; "First-level Analysis using selxavg3-sess"), but get the following error when I run the selxavg3-sess command: ------- octave output -------------------- LD_LIBRARY_PATH /homec/hbu23/hbu231/Software/Octave/3.8.2/lib:/homec/hbu23/hbu231/Software/Libraries/octave /homec/hbu23/hbu231/Software/Octave/3.8.2/bin: Permission denied. ------------------------------------------ ERROR: fast_selxavg3() failed\n I've searched the list and a similar error has been reported but the proposed solution did not work. (https://mail.nmr.mgh.harvard.edu/pipermail//freesurfer/2016-April/044814.ht…) Specifically, it does not work despite the named directory and several enclosing directories received full read, execute and for tests even write permissions (also, as a test, the freesurfer folders received full permissions): drwxrwxrwx 2 hbu231 hbu23 512 Jan 9 16:52 bin As I’m attempting to run it on a supercomputer I cannot execute it as a superuser. Does anyone have any thoughts on how to solve this? Thanks! Michael ----------- 1) FreeSurfer version: freesurfer-Linux-centos6_x86_64-stable-pub-v5.3.0 2) Platform: CentOS Linux release 7.2.1511 3) uname -a: Linux jrl11 3.10.0-327.36.3.el7.x86_64 #1 SMP Mon Oct 24 16:09:20 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux

2 1

automatic segmentation of hippocampal subregions in young children
by Tracy Riggins 10 Jan '17

10 Jan '17

Dr. Inglesias I am a developmental psychologist studying memory and brain development in early childhood at the University of Maryland. My lab has acquired whole-brain standard resolution T1-weighted scans and ultra-high resolution T2-weighted scans from 125 child participants aged 4-8 years. We are using your automatic segmentation for hippocampal subregions to obtain volumes for assessment in our study. As we write up the results for publication, I began to wonder if reviewers will ask for verification that the automatic segmentation is valid for our young age group. I know you tested Alzheimer's and an elderly sample in your paper, so I suspect it is valid for a variety of populations. But I am writing to inquire your thoughts on this issue. Do you think verification is needed? If so, would you suggest we hand trace a subset of our data to verify the automatic methods? I would appreciate any advise as to how to best approach this. Thank you, in advance, for your consideration of my inquiry. Sincerely, Tracy Tracy Riggins, Ph.D. Associate Professor University of Maryland

2 1

dice coefficient computing
by N Saf 10 Jan '17

10 Jan '17

Dear freesurfers, I analysis my T1 MRI with FreeSurfer (aseg.mgz 256x256x256) and I have manual segmentation for hippocampus in my subject space(208x176x176)as well. I want to compute the dice coefficient with one simple Matlab code. how can I do it correctly ?if I want to bring the aseg to subject space , which command with which options is the best way to do it correctly ?any help can be appreciated.(I did it for the FSL segmentation result but it was simple because the result is in my subject space). Best Regards, Nazanin

2 4

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