Freesurfer July 2012

freesurfer@nmr.mgh.harvard.edu

128 participants
194 discussions

mri_vol2surf failed, segmentation error
by Qingyang Li 04 Feb '14

04 Feb '14

Hi all, I am new to the freesurfer world. I was running a pipeline where I met an error. Here is the command I run and the error message: Command: mri_vol2surf --hemi lh --noreshape --o /working_dir/SurfaceRegistration/_subject_id_sub/func_sampler_lh/mapflow/_func_sampler_lh4/lh.rest_mc.mgz --reg /working_dir/SurfaceRegistration/_subject_id_sub/func_bbreg/mapflow/_func_bbreg4/rest_mc_bbreg_sub.dat --projfrac 0.500 --mov /working_dir/SurfaceRegistration/_subject_id_sub/func_volreg/mapflow/_func_volreg4/rest_mc.nii.gz Standard output: srcvol = /working_dir/SurfaceRegistration/_subject_id_sub/func_volreg/mapflow/_func_volreg4/rest_mc.nii.gz srcreg = /working_dir/SurfaceRegistration/_subject_id_sub/func_bbreg/mapflow/_func_bbreg4/rest_mc_bbreg_sub.dat srcregold = 0 srcwarp unspecified surf = white hemi = lh ProjFrac = 0.5 thickness = thickness reshape = 0 interp = nearest float2int = round GetProjMax = 0 INFO: float2int code = 0 INFO: changing type to float Done loading volume Reading surface /NKI_recon_subjects/sub/surf/lh.white Done reading source surface Reading thickness /NKI_recon_subjects/sub/surf/lh.thickness Done Mapping Source Volume onto Source Subject Surface 1 0.5 0.5 0.5 using old Done mapping volume to surface Standard error: Segmentation fault Return code: 139 Interface SampleToSurface failed to run. Thank you all in advance! Yang ------ Qingyang Li Balance and focus. On Sun, Mar 18, 2012 at 8:28 AM, <freesurfer-request(a)nmr.mgh.harvard.edu>wrote: > Mailing list subscription confirmation notice for mailing list > Freesurfer > > We have received a request from 67.244.70.75 for subscription of your > email address, "congel(a)gmail.com", to the > freesurfer(a)nmr.mgh.harvard.edu mailing list. To confirm that you want > to be added to this mailing list, simply reply to this message, > keeping the Subject: header intact. Or visit this web page: > > > https://mail.nmr.mgh.harvard.edu/mailman/confirm/freesurfer/8f8fcc67483bbce… > > > Or include the following line -- and only the following line -- in a > message to freesurfer-request(a)nmr.mgh.harvard.edu: > > confirm 8f8fcc67483bbcee3ff0f76888c79e7704b4a03a > > Note that simply sending a `reply' to this message should work from > most mail readers, since that usually leaves the Subject: line in the > right form (additional "Re:" text in the Subject: is okay). > > If you do not wish to be subscribed to this list, please simply > disregard this message. If you think you are being maliciously > subscribed to the list, or have any other questions, send them to > freesurfer-owner(a)nmr.mgh.harvard.edu. > > > The information in this e-mail is intended only for the person to whom it > is > addressed. If you believe this e-mail was sent to you in error and the > e-mail > contains patient information, please contact the Partners Compliance > HelpLine at > http://www.partners.org/complianceline . If the e-mail was sent to you in > error > but does not contain patient information, please contact the sender and > properly > dispose of the e-mail. > >

4 4

Freesurfer - problem with automatic labeling of V1 using the Hinds method in 5.1.0
by Andrew Bock 08 Mar '13

08 Mar '13

Hello, I am having a similar problem to one described on the freesurfer mail archive, where I receive the following error when trying to identify V1 using the Hinds method (i.e. recon-all -s <subject> -label_v1): processing subject V1_average... MRISreadVertexPosition(white): could not open file /usr/local/freesurfer/ subjects/V1_average/surf/rh.white No such file or directory MRISreadOriginalProperties: could not read surface file white No such file or directory could not read surface positions from white No such file or directory As Prof. Krish Singh indicated in July, it appears the problem occurs when trying to open rh.white (or lh.white) from the surf subdirectory in V1_average, which doesn't exist. I have run recon-all -s <subject> -make all, but it returns "make: Nothing to be done for `all'." Any help you can provide would be greatly appreciated! Thanks in advance, Andrew -- Andrew Bock, Ph.D. Department of Psychology University of Washington Seattle, WA 98195 Tel: 206.543.3817 Fax: 206.685.3157 abock(a)u.washington.edu

4 8

entorhinal cortex volume question
by Cat Chong 29 Jan '13

29 Jan '13

Hello experts, I am very new to freesurfer so please excuse my basic question: We want to get measurements of entorhinal cortex volumes on a group of people. I found these results listed in ?h.aparc.stats. I also noticed another ?h.entorhinal_exvivo.stats, with different results. Which should I use for a group analysis of entorhinal volume and surface area, and why are the results very different? very best regards, catherine

3 4

selxavg3 error : NaN in Xtmp.mat file
by Sid Kundu 07 Dec '12

07 Dec '12

Hi Doug, I created a new first level analysis along with new contrasts and event numbers for a group of subjects I had previously done analyses on. All but one of the subjects finished selxavg3 without error, but one subject gave the following error: Saving X matrix to /home/sidkundu/freesurfer/sessions/Heading-TTP/LMV2010_N19_Heading-TTP/bold/TTP-analysis-3/Xtmp.mat Error using svd Input to SVD must not contain NaN or Inf. Error in cond (line 39) s = svd(A); Error in fast_selxavg3 (line 248) XCond = cond(XtX); >> ------------------------------------------ ERROR: fast_selxavg3() failed\n I looked at the generated Xtmp.mat file and it shows the design matrix X has all zeros in 2 columns, resulting in NaNs in the same columns for Xn, and corresponding NaNs in XtX. I'm wondering what would cause this or why it would only occur for only one subject. _________________________________ Sid Kundu Biomedical Engineering Boston University

3 6

flipping surface data
by kelsi 04 Dec '12

04 Dec '12

Hi Freesurfers, Is there a tool to flip the orientation of surface data to the opposite hemisphere? We are working with stroke patients, and I would like to do analysis of the contralesional hemisphere of all the patients. However, this requires putting them on the same hemisphere for ttests, etc. I saw an earlier thread online about a program for flipping in the works, I was just wondering if this is now possible? Thank you! Kelsi

3 3

AAL labels
by MCLAREN, Donald 23 Oct '12

23 Oct '12

Has anyone created an AAL label file for Freesurfer? I've done some ROI analysis using AAL regions and want to display regional values on the freesurfer surface. Best Regards, Donald McLaren ================= D.G. McLaren, Ph.D. Research Fellow, Department of Neurology, Massachusetts General Hospital and Harvard Medical School Postdoctoral Research Fellow, GRECC, Bedford VA Website: http://www.martinos.org/~mclaren Office: (773) 406-2464 ===================== This e-mail contains CONFIDENTIAL INFORMATION which may contain PROTECTED HEALTHCARE INFORMATION and may also be LEGALLY PRIVILEGED and which is intended only for the use of the individual or entity named above. If the reader of the e-mail is not the intended recipient or the employee or agent responsible for delivering it to the intended recipient, you are hereby notified that you are in possession of confidential and privileged information. Any unauthorized use, disclosure, copying or the taking of any action in reliance on the contents of this information is strictly prohibited and may be unlawful. If you have received this e-mail unintentionally, please immediately notify the sender via telephone at (773) 406-2464 or email.

2 1

Freesurfer and FSL
by Ana Arruda 10 Oct '12

10 Oct '12

Dear experts, I'm an MSc student and part of my final project is related to the measurement of the volume of the hippocampus. The last thing I have to do is to compare two softwares, FSL and Freesurfer. On Freesurfer, I've used the following commands: recon-all -i where_my_subject_is -autorecon1 -s output_name then, related to the same subject recon-all -autorecon2 -s output_name After approximately 12 hours, and after comparing the results of Freesurfer and FSL, I found that those obtained with Freesurfer are bigger than those obtained with FSL. According to the literature the results are the reverse. Do you know any idea what I'm doing wrong? Thank you for your attention, Ana Arruda

6 25

Re: [Freesurfer] integrating FSL and Freesurfer
by noam Schneck 09 Oct '12

09 Oct '12

Hi Doug Here are more detailed versions of my first two questions. 1) I used tkmedit to look at the mean_func produced by a gfeat analysis of one subject's three functional runs overlayed onto the anatomical image produced by recon-all on that same subject(bert/mri/orig.mgz). The registration matrix I used was the anat2std.dat that was produced by running reg-feat2anat on all three of the feat analyses from the individual runs. When I did this I saw that the mean_func image was outside the boundaries of the anatomical image. Correspondingly, when I overlayed the zstat.nii.gz from the gfeat analysis onto the anatomical image I had activation outside of the anatomical image brain. 2) Originally, I ran BBR on a subject that I had run Recon-all for. I used BBR to register the example_func.nii.gz to the bert/mri/orig.mgz. I had BBR output an fsl.mat. I then used ApplyXFM in fsl to apply BBR's fsl.mat to register the example_func.nii.gz onto a T1 image that had been created by converting bert/mri/brainmask.mgz into niftii with mri_convert. I tried it again, this time converting the orig.mgz to niftii. In both cases the example_func in T1 space image produced by ApplyXFM is rotated about 90 degrees and also inverted. This is viewed in fslview, The corresponding T1 images created by mri_convert are not. What can I do to allow me to use BBR's fsl.mat within ApplyXFM?

3 4

More than two time-points in longitudinal base = Memory allocation error?
by Liz Bowman 03 Oct '12

03 Oct '12

Hello, I have several patients scans taken over multiple time points (up to nine time points in some patients). I thought I might be able to use these to create a base template for longitudinal processing, however the process fails and exits with errors if I try to create a base with more than two time points from a patient. This is a sample of the error log of the most recent attempt, with three timepoints: Resolution: 1 S( 128 128 128 ) T( 128 128 128 ) Iteration(f): 1 -- diff. to prev. transform: 0.876663 Iteration(f): 2 -- diff. to prev. transform: 0.0148738 Iteration(f): 3 -- diff. to prev. transform: 0.00156813 < 0.01 :-) Resolution: 0 S( 256 256 256 ) T( 256 256 256 ) Iteration(f): 1 -- diff. to prev. transform: 0.263086 Iteration(f): 2mri_robust_template(1035) malloc: *** mmap(size=262144) failed (error code=12)*** error: can't allocate region*** set a breakpoint in malloc_error_break to debugMRIalloc(256, 256, 256): could not allocate 262144 bytes for 68th slice Cannot allocate memoryDarwin 550D-UOM88670.local 10.8.0 Darwin Kernel Version 10.8.0: Tue Jun 7 16:32:41 PDT 2011; root:xnu-1504.15.3~1/RELEASE_X86_64 x86_64 recon-all -s KSx3 exited with ERRORS at Mon Jul 23 10:15:13 EST 2012 For more details, see the log file /Applications/freesurfer/subjects/KSx3/scripts/recon-all.logTo report a problem, see http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting I would be grateful for any advice. Regards, EA Bowman

3 4

Error with mri_volcluster
by Antonella Kis 01 Oct '12

01 Oct '12

Dear Doug, I was running the mri_glmfit-sim \ --glmdir GroupAnalysis_abs.glmdir \ --sim perm 1000 3 permcsd \ --sim-sign abs \ --cwpvalthresh .05 just to get an idea of how my FA Group Analysis data is looking and I got an error before to get the permcsd.sig.cluster.summary file from the mri_volcluster: mri_volcluster --in GroupAnalysis_abs.glmdir/contrast/sig.mgh --mask GroupAnalysis_abs.glmdir/mask.mgh --reg /usr/local/freesurfer/subjects/fsaverage/mri.2mm/reg.2mm.dat --no-fixmni --cwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.mgh --sum GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.summary --ocn GroupAnalysis_abs.glmdir/contrast/permcsd.sig.ocn.mgh --cwpvalthresh .05 --seg fsaverage aparc+aseg.mgz --csdpdf GroupAnalysis_abs.glmdir/contrast/permcsd.pdf.dat --csd GroupAnalysis_abs.glmdir/csd/permcsd.j001-contrast.csd --vwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.voxel.mgh ERROR: cannot find /media/2Tra/Antonella_2011/DTI_FS/GLM_CVS/fsaverage/mri/aparc+aseg.mgz I replaced the fsaverage/mri/aparc+aseg.mgz with the cvs_avg35/mri/aseg.mgz and to get the cluster summary I was running: mri_volcluster --in GroupAnalysis_abs.glmdir/contrast/sig.mgh --mask GroupAnalysis_abs.glmdir/mask.mgh --reg /usr/local/freesurfer/subjects/fsaverage/mri.2mm/reg.2mm.dat --no-fixmni --cwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.mgh --sum GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.summary --ocn GroupAnalysis_abs.glmdir/contrast/permcsd.sig.ocn.mgh --cwpvalthresh .05 --seg cvs_avg35 aseg.mgz --csdpdf GroupAnalysis_abs.glmdir/contrast/permcsd.pdf.dat --csd GroupAnalysis_abs.glmdir/csd/permcsd.j001-contrast.csd --vwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.voxel.mgh but then I am getting the following error: NOT Adjusting threshold for 1-tailed test Using defalt ctab /usr/local/freesurfer/FreeSurferColorLUT.txt INFO: loading mask volume: GroupAnalysis_abs.glmdir/mask.mgh INFO: MRIbinarize01(): nhits = 961260 Found 961260 voxels in mask ERROR: Input volume voxel dimensions do not match those in the registration file. If the input volume is in bshort/bfloat format, check that there is an accompanying bhdr file. Can you please help me to solve this problem. I am not sure what to do. Thank you very much! Antonella PS. Just to remember you that I registered my FA data over the CVS space but I have for each subject from dt_recom the register.dat file so I suppose I need to keeo the fsaverage in the --reg COMMAND-LINE ARGUMENTS.

3 2

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 July 2012