Freesurfer December 2012

freesurfer@nmr.mgh.harvard.edu

101 participants
161 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

Re: [Freesurfer] Fwd: white matter and pial surface don't track sulcus
by Darren Gitelman 06 Dec '13

06 Dec '13

Yes. That's it. I tried editing the wm.mgz earlier and deleting all voxels labels as ventricle (250) but that didn't solve the problem. Do I edit the aseg.mgz file and just label it as none in this region? Would I use the command tkmedit <subjid> brainmask.mgz -aux brain.finalsurfs.mgz -surfs -segmentation aseg.mgz $FREESURFER_HOME/FreeSurferColorLUT.txt When done do I run recon-all -autorecon2-noaseg -autorecon3? Thanks, Darren On Sun, Sep 16, 2012 at 4:38 PM, Bruce Fischl <fischl(a)nmr.mgh.harvard.edu>wrote: > Check the aseg and see if it was labeled ventricle (and make sure it > really isn't). You may need to edit the aseg > > > > On Sep 16, 2012, at 5:19 PM, Darren Gitelman <d-gitelman(a)northwestern.edu> > wrote: > > Hi > > I have run into an issue that I'm not sure how to correct. What happens > is that in some subjects the surfaces don't seem to properly track sulci > coming off the medial part of temporo-occipital cortex. Please see green > arrows on the attached image. > > In fact the tutorial on topological defects shows an image with a > similar issue, but fixing it is not discussed. > > http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/TopologicalDefect?actio… > > I have checked and these are sulci and not the occipital extent of the > lateral ventricle. > > Perhaps it is not a problem, or if it is then how to approach fixing it? > > Thank you, > Darren > -- > Darren Gitelman, MD > Northwestern University > 710 N. Lake Shore Dr. > Abbott 11th Floor > Chicago, IL 60611 > Ph: (312) 908-8614 > Fax: (312) 908-5073 > > <3809-error.png> > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > 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. > -- Darren Gitelman, MD Northwestern University 710 N. Lake Shore Dr. Abbott 11th Floor Chicago, IL 60611 Ph: (312) 908-8614 Fax: (312) 908-5073

7 8

hippocampal segmentation converting mgz to nii
by lorenzo pasquini 13 Aug '13

13 Aug '13

Dear freesurfer users, I have used freesurfer to segment hippocampal subfields. Then I used freeview to convert the derived .mgz files in .nii files for further analysis. I estimated the hippocampal subfields volumes with the command: kvlQuantifyPosteriorProbabilityImages $ FREESURFER_HOME/ data/ GEMS/ compressionLookupTable.txt posterior_right_* posterior_Right-Hippocampus.mgz I have calculated the number of voxels in the converted .nii files using SPM and realized that they do not fit the estimated volumes through freesurfer Example for CA2-2 estimated voxel number though freesurfer: 6851 voxels found in .nii file: 15812 Did somebody had similar problems? Do you have any suggestion? Are there any incompatibilities between SPM and freesurfer? I did not find any apparent error in my counting script. Thank you very much and regards, Lorenzo Pasquini Can. Ph. D., Medical Life Science and Technology (Technische Universität Muenchen) Neuroimaging Lab, Neuroradiology Klinikum r. d. Isar (Technische Universität) 81675 München Tel.: 0049 (0) 89 4140-7666

3 4

vertices don't match
by Darren Gitelman 11 Jun '13

11 Jun '13

Dear List I am adding to the post on this error from yesterday regarding problems with running qcache. I know this has been posted to the list before, but the suggested fixes and checks don't fix the problem. See the error below (ERROR: number of vertices in /projects/d20277/dek/done/analyzed/3300/surf/lh.volume does not match surface (145332,147148) ERROR: reading curvature file) 1) Running recon-all -s <subjid> -vno-check does not find any problems. It says everything is verified and checks out ok. 2) Running recon-all -s <subjid> -make all returns the message: make: Nothing to be done for `all'. But when I try recon-all -s <subjid> -qcache I get the error on this and several other subjects, but not all subjects. Checking the modification times on some files it looks like the thickness file came after the curv file. [drg@quser03 surf]$ ll *.thickness -rwxrwxrwx 1 dek322 d20277 588607 Sep 11 21:54 lh.thickness -rwxrwxrwx 1 dek322 d20277 593959 Sep 11 22:10 rh.thickness [drg@quser03 surf]$ ll *.curv -rwxrwxrwx 1 dek322 d20277 588607 Sep 11 21:45 lh.curv -rwxrwxrwx 1 dek322 d20277 593959 Sep 11 22:01 rh.curv Any suggestions for fixing the qcache error would be appreciated. Darren -- --------------------------------------------------- #@# 1/1 3300 Wed Sep 19 13:52:44 CDT 2012 -------------- ----------------------- mri_surf2surf --srcsubject 3300 --srchemi lh --srcsurfreg sphere.reg --trgsubject fsaverage --trghemi lh --trgsurfreg sphere.reg --tval ./tmp.mris_preproc.6938/3300.1.mgh --sval /projects/d20277/dek/done/analyzed/3300/surf/lh.volume --sfmt curv --noreshape --no-cortex Source registration surface changed to sphere.reg Target registration surface changed to sphere.reg srcsubject = 3300 srcval = /projects/d20277/dek/done/analyzed/3300/surf/lh.volume srctype = curv trgsubject = fsaverage trgval = ./tmp.mris_preproc.6938/3300.1.mgh trgtype = srcsurfreg = sphere.reg trgsurfreg = sphere.reg srchemi = lh trghemi = lh frame = 0 fwhm-in = 0 fwhm-out = 0 label-src = (null) label-trg = (null) OKToRevFaceOrder = 1 Reading source surface reg /projects/d20277/dek/done/analyzed/3300/surf/lh.sphere.reg Loading source data Reading curvature file /projects/d20277/dek/done/analyzed/3300/surf/lh.volume ERROR: number of vertices in /projects/d20277/dek/done/analyzed/3300/surf/lh.volume does not match surface (145332,147148) ERROR: reading curvature file Linux quser01 2.6.18-128.el5 #1 SMP Wed Dec 17 11:41:38 EST 2008 x86_64 x86_64 x86_64 GNU/Linux recon-all -s 3300 exited with ERRORS at Wed Sep 19 13:52:46 CDT 2012 -- Darren Gitelman, MD Northwestern University 710 N. Lake Shore Dr. Abbott 11th Floor Chicago, IL 60611 Ph: (312) 908-8614 Fax: (312) 908-5073

3 2

Probably activation due to vessels
by Jörg Pfannmöller 07 May '13

07 May '13

Dear Freesurfers, I computed the functional activation in a sensory task using freesurfer's fsfast stream. In the cortex surface line like activated regions are present. Their extension is much larger than the extension of the expected regions. Therefore, I suppose that those regions are due to vessels. Is there a possibility in freesurfer to remove them. Ideally there would be something like a volume based mode which does not sxclude the cortex. In this mode the vessels couls be identified and removed. Any help is highly apprexiated. Respectfully yours pfannmoe

2 13

Regarding cortical thickness
by Bharti Rana 14 Apr '13

14 Apr '13

I am Bharti, a M.Tech./Ph.D student at School of Computer & Systems Sciences, Jawaharlal Nehru University, New Delhi, India. My area of interest is pattern recognition and image processing. In one of my research work, i used freely available freesurfer software to find cortical thickness using recon-all -all command. But, while retrieving cortcal thickness maps using read_curv() function, i am getting different number of vertices for different brains. Are different number of cortical vertices possible for two different brains? I have gone through FAQ provided by Freesurfer but i did not find any related answer. Kindly clarify. Thanks in advance. -- Best Regards Bharti M.Tech./Ph.D. School of Computer & Systems Sciences Jawaharlal Nehru University NEW DELHI, INDIA

2 3

contrast matrix
by xiangbo_2010 28 Mar '13

28 Mar '13

Hi Freesurfer experts, I'm very sorry to bother you, but I am very confused with the following questions: My experimental design includes three discrete factors: factor 1 with three levels (A,B,C ); factor 2 with three levels (O,P,Q); gender (F, M), and one covariate. So I can get 18 classes: FAO, FAP,FAQ,FBO,FBP,FBQ,FCO,FCP,FCQ,MAO, MAP,MAQ,MBO,MBP,MBQ,MCO,MCP,MCQ. I want to perform the interaction between factor 1 and factor 2 regressing out the effect of gender and one covariate, but I don't know the rules for setting the contrasts for the F-test. The contrast matrix I used is: 1 1 1 -1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0 1 -1 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 -1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 1 0 -1 1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 -1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0 1 -1 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 -1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 -1 1 0 -1 1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 is it correct? Any help will be very appreciated. Best wishes, Bo Xiang

6 18

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

Re: [Freesurfer] "nu_estimate_np_and_em" error during recon-all
by Douglas Greve 15 Feb '13

15 Feb '13

Hi Tudor, I don't know how to solve this problem. Maybe Nick knows. doug On 12/19/12 6:59 AM, Tudor Popescu wrote: > Hi Doug, thanks for replying, here is the file attached. In case the > attachment doesn't go through please let me know if I should send it > to a different address. > Tudor > > On 18 December 2012 17:34, Douglas N Greve <greve(a)nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > > Can you send the recon-all.log file (located in the > subject/scripts folder). > doug > > > On 12/18/2012 11:41 AM, Tudor Popescu wrote: > > Hello, > > > > First of all, apologies for asking for help three times today! > THings > > will hopefully become clearer once I get used to FS a bit; I > really do > > appreciate the amazing support provided by the team behind > FreeSurfer! > > > > I ran the following command on one of my structurals: > > > > /recon-all -subject M01 -i M01.nii.gz -all/ > > > > After a few minutes, recon-all ended with the following message: > > > > /Assertion failed at line 827 in file templates/CachedArray.cc/ > > /nu_estimate_np_and_em: crashed while running volume_stats > > (termination status=256)/ > > /nu_correct: crashed while running nu_estimate_np_and_em > (termination > > status=256)/ > > /ERROR: nu_correct/ > > /Linux FreeSurfer 2.6.28-11-generic #42-Ubuntu SMP Fri Apr 17 > 01:57:59 > > UTC 2009 i686 GNU/Linux/ > > > > I apologise for reporting an error that's already been reported > in the > > past, but unfortunately the solutions suggested in those past > threads > > did not work in my case. To address a few of them: > > > > * my PATH (listed below) does not seem to include two MNI > installations > > > //home/virtualuser/freesurfer/bin:/home/virtualuser/freesurfer/fsfast/bin:/home/virtualuser/freesurfer/tktools:/home/virtualuser/freesurfer/mni/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games/ > > * the version of FreeSurfer that I installed is the latest disk > > image available on the website > > (freesurfer-Virtualbox-linux-x86-stable-pub-v5.1.0-full) - I did > > not see any dev version that would presumably include any > missing > > libraries presumably implied by the above error > > * I tried updating the Xubuntu inside my VirtualBox, however it > > reports an insufficient space problem that from I can see > does not > > have a trivial fix > > > > Advance thanks for any help! > > Tudor > > > > > > _______________________________________________ > > Freesurfer mailing list > > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 > Fax: 617-726-7422 > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > Outgoing: > ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > 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. > >

3 19

issue about LME Matlab tools
by ting xu 13 Feb '13

13 Feb '13

Dear all, A few questions come to me when I used LME Matlab toolbox. Then Jorge kindly answered my questions and he also asked me to post our discussion to Freesurfer list so everyone knows what's going on here :) Cheers, Ting ------------------------------------------------------------------- *Subject: *Re: Re: issue about LME Matlab tools *Sent: *Tue, Dec 11, 2012 21:41:26 Hi Ting Please could you post this question to the Freesurfer list, so other people can benefit from this discussion? I'll answer you ASAP. Best -Jorge ------------------------------------------------------ *Subject: *Re: Re: issue about LME Matlab tools *Sent: *Wed, Dec 12, 2012 6:44:36 AM Hi Jorge, Thank you for your explanation in more details. I followed your suggestion, set the convergence epsilon to 10^-5. It truly improved. I am now running the whole brain data, estimations in several voxels were not converge after 50 iterations. Given that iterations takes more computational time, how many iterations do you recommend and what do you usually do if it still not converge? Thank you again:) Warmly regards, Ting ------------------------------------------------------ *Subject:* Re: issue about LME Matlab tools *Sent: *Wed, Dec 12, 2012 00:50:35 Hi Ting You needed more iterations to make the algorithm stop closer to the optimal values. Try: statsFS= lme_fit_FS(X,[1 2],Y,ni,10^-5); or the EM algorithm statsEM = lme_fit_EM(X,[1 2],Y,ni,10^-10); Neither of these algorithms can be guaranty to converge but I have found the FS algorithm to be the most robust and fast. For the mass-univariate setting we limited the number of iterations for the FS algorithm to only 20 due to computational time. Note that for only two repeated measures (as in your data) compound symmetry (a model with a single random effect for the intercept term) likely holds for the covariance matrix among the repeated measures. Although a likelihood ratio test here comparing the model with one random effects against the model with two random effects is barely significant it will not likely survive a multiple comparisons correction. You can not impose structure on D in our toolbox nor it is recommended for general longitudinal data. Let me know any doubt you might have. Best -Jorge ------------------------------------------------------ *Subject: *issue about LME Matlab tools *Sent: *Tue, Dec 11, 2012 14:47:14 Dear Dr. Sabuncu Your recent work "Statistical Analysis of Longitudinal Neuroimage Date with Linear Mixed Effects Models" provided us a great matlab toolbox to apply Mixed Model, especially for imaging data. Now, I am investigating the intraclass correlation based on this toolbox. Sometime I found the function "lme_fit_FS" could not get the accurate estimations. I think the problem may be relevant to the initial value from OLS method, in the case that the within-subject variability estimated from "lme_fit_init" is very close to zero. I attached an example. Y is test-retest imaging data for one voxel, X is design matrix including intercept, head_motion, age, gender. I set random effect for intercept and head_motion. I also attached result from SAS and the model is ok. stats = lme_fit_FS(X,[1,2],Y,ni); Here I'm wondering is there any way to fixed this in matlab? PS, though you have mentioned that no imposed structure on D (covariance matrix of random effect), I am curious if it is possible to define the structure of D in this toolbox. Thanks for your kind attention and look forward to your reply soon. Regards, Ting

3 5

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