Freesurfer

freesurfer@nmr.mgh.harvard.edu

30132 discussions

nuisance factors
by Georgina Angelopoulou 30 Oct '19

30 Oct '19

External Email - Use Caution Dear FreeSurfer Developers, I am trying to design a matrix in order to run a group comparison analysis and I would like to verify if I am following the right steps. I have read a lot of things in the forum, but I got confused regarding how to treat nuisance factors. I have two groups (patients and controls) and I want to test if there is a significant difference in cortical thickness regressing out the effects of age and sex. My fsgd file is ControlMales ControlFemales PatientMales PatientFemales Again, based on what I have read so far, in order to use the age as a nuisance variable, I first tested the age slopes using the row values of age (meaning without centering them), is that correct? Since there is no significant difference in age slopes, that means that now I can run the analysis without including the age at all as a nuisance variable in the model? (what you usually refer as DOSS?) with the contrast matrix 0.5 0.5 -0.5 0.5 (meaning that I only control for the effect of sex)? Thank you in advance, Georgina

2 1

question about qdec 'Generate Stats Data Tables'
by 艺 30 Oct '19

30 Oct '19

External Email - Use Caution Hello, We had a problem when using QDEC for group analysis. After generating relevant files, I press 'Load Data Tabel' button to load my file 'qdec.tabel.dat'. But when I press the 'Generate Stats Data Tables' button, a popup window shows an error indicating that: --------------------------------------------------------------------------------------------------------------------- ERROR: In /usr/pubsw/packages/KWWidgets/CVS-vtk560/KWWidgets/vtKWTkUtilities.cxx,line230 vtKWQdecApp (0x27fbfce0): Script: vtKTemp2 GenetateStatsDataTables Returened Error on line 1: Uncaught exception: command failed: asegstats2table --common-segs --meas volume --tablefile /usr/local/freesurfer/subjects/qdec/stats_tables/aseg.volume.stats.dat --statsfile=aseg.stats --subjects PT1_01 PT1_03 PT1_04 PT1_06 PT1_11 PT1_12 PT1_14 PT1_27 PT1_32 PT1_33 PT1_34 PT1_47 PT1_55 PT1_05 PT1_08 PT1_09 PT1_16 PT1_18 PT1_19 PT1_20 PT1_21 PT1_23 PT1_24 PT1_26 PT1_38 PT1_45 PT1_48 PT1_52 PT1_75 PT1_81 PT1_87 Stack trace: Uncaught exception: command failed: asegstats2table --common-segs --meas volume --tablefile /usr/local/freesurfer/subjects/qdec/stats_tables/aseg.volume.stats.dat --statsfile=aseg.stats --subjects PT1_01 PT1_03 PT1_04 PT1_06 PT1_11 PT1_12 PT1_14 PT1_27 PT1_32 PT1_33 PT1_34 PT1_47 PT1_55 PT1_05 PT1_08 PT1_09 PT1_16 PT1_18 PT1_19 PT1_20 PT1_21 PT1_23 PT1_24 PT1_26 PT1_38 PT1_45 PT1_48 PT1_52 PT1_75 PT1_81 PT1_87 while executing "vtkTemp2 GenerateStatsDataTables" --------------------------------------------------------------------------------------------------------------------- Before that we were able to successfully ‘Generate Stats Data Tables’  (for our data but not all of the above data). We don't know if it is a server environment or something else. What‘s wrong? Please, could you advise me? Best regards,   Adam

2 1

Mapping coordinates from sphere file to ?h.white
by Oscar Esteban 30 Oct '19

30 Oct '19

External Email - Use Caution Dear experts, In the context of the BIDS transforms extension and the ".x5" format we chatted about in the workshop we made in March at MIT, we are addressing the following question. We have a sphere of fsLR coordinates in fsaverage space (e.g., https://github.com/Washington-University/HCPpipelines/blob/master/global/te… - but the question extends to the right hemisphere and also resolutions other than 59k). We want to expand this sphere to the fsaverage's white surface (lh.white, in this example). Since coordinates on the fsLR sphere don't match any of fsaverage{,3,4,5,6}'s coordinates, we wonder whether FreeSurfer stores somewhere/somehow the analytical transform of this mapping. In other words, how one would map arbitrary coordinates on the surface to a particular fsaverage's white (or midthickness, or pial, etc.) surface. Otherwise, I guess I would pick fsaverage's surfaces (i.e., 7th degree icosahedron), and interpolate the displacement vectors that start at each sample on the surface and end at the corresponding location of that vertex index of the destination surface (lh.white for the example above). Is there any tool that actually does the _inverse transform of the surface-to-sphere mapping_? Thanks very much. Regards, Oscar -- *___________________________* *Oscar Esteban, Ph.D.* Postdoctoral Fellow, Poldrack Lab Stanford University +1 (650) 733 33 82

2 1

recon-all -skullstrip issues
by Adam Rytina 30 Oct '19

30 Oct '19

External Email - Use Caution Hello all, I´m struggling with running "recon-all -skullstrip", in general with "recon-all". I added NIfTI series file to a new directory "QWE" (/usr/local/freesurfer/subjects/QWE) and running the skull stripping "mri_watershed /usr/local/freesurfer/subjects/QWE/ser.nii /usr/local/freesurfer/subjects/QWE/brainwithout.nii" without any proble. However, when trying to adjust watershed parameters the problem occurs. "recon-all -skullstrip -wsthresh 5 -clean-bm -subjid QWE". There are afterwards some errors occuring: "mv: cannot stat '/usr/local/freesurfer/subjects/QWE/mri/brainmask.mgz': No such file or directory", "mghRead(/usr/local/freesurfer/subjects/QWE/mri/nu.mgz, -1): could not open file mri_em_register: could not open input volume nu.mgz." I´m working with image series "ser.nii" and "brainwithout.nii", not with "brainmask.mgz" and "nu.mgz" in the tutorial. Please, where is the problem? Thanks a lot Adam [https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-anima…]<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> Bez virů. www.avast.com<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…>

2 3

mri_glmfit ill-conditioned matrix question
by Nillo, Ryan Michael R 30 Oct '19

30 Oct '19

Hello, I ran mri_glmfit with the following fsgd file: GroupDescriptorFile 1 Title CJD Class MM1_MV1 Class MV2_VV2 Variables frontal Input 7693 MV2_VV2 2.930667 Input 11985 MV2_VV2 2.410333 Input 13082 MV2_VV2 2.529333 Input 13425 MV2_VV2 2.593667 Input 13667 MM1_MV1 2.628667 Input 13900 MM1_MV1 2.675333 Input 14094 MM1_MV1 2.549333 Input 15566 MV2_VV2 2.666667 Input 15999 MM1_MV1 2.641 Input 16714 MV2_VV2 2.643 Input 16848 MV2_VV2 2.63 Input 16952 MV2_VV2 2.516 Input 7693a MV2_VV2 3.054 Input 11985a MV2_VV2 2.42 Input 13082a MV2_VV2 2.454 Input 13425a MV2_VV2 2.602 Input 13667a MM1_MV1 2.604333 Input 13900a MM1_MV1 2.601333 Input 14094a MM1_MV1 2.523333 Input 15566a MV2_VV2 2.713333 Input 15999a MM1_MV1 2.584333 Input 16714a MV2_VV2 2.673333 Input 16848a MV2_VV2 2.533333 Input 16952a MV2_VV2 2.486 My two contrast matrices are: one.mtx = [ 0 0 1 0 ] two.mtx = [ 0 0 0 1 ] When running the mri_glmfit command "mri_glmfit --glmdir new.frontal.lgi --fsgd new.frontal.fsgd.txt --fwhm 30 --surf fsaverage lh --y lh.pial_lgi.fwhm0.mgh --C one.mtx --C two.mtx”, I get the following output: INFO: gd2mtx_method is dods Saving design matrix to new.frontal.lgi/Xg.dat Computing normalized matrix Normalized matrix condition is 12820.7 Design matrix ------------------ 0.00000 1.00000 0.00000 2.93067; 0.00000 1.00000 0.00000 2.41033; 0.00000 1.00000 0.00000 2.52933; 0.00000 1.00000 0.00000 2.59367; 1.00000 0.00000 2.62867 0.00000; 1.00000 0.00000 2.67533 0.00000; 1.00000 0.00000 2.54933 0.00000; 0.00000 1.00000 0.00000 2.66667; 1.00000 0.00000 2.64100 0.00000; 0.00000 1.00000 0.00000 2.64300; 0.00000 1.00000 0.00000 2.63000; 0.00000 1.00000 0.00000 2.51600; 0.00000 1.00000 0.00000 3.05400; 0.00000 1.00000 0.00000 2.42000; 0.00000 1.00000 0.00000 2.45400; 0.00000 1.00000 0.00000 2.60200; 1.00000 0.00000 2.60433 0.00000; 1.00000 0.00000 2.60133 0.00000; 1.00000 0.00000 2.52333 0.00000; 0.00000 1.00000 0.00000 2.71333; 1.00000 0.00000 2.58433 0.00000; 0.00000 1.00000 0.00000 2.67333; 0.00000 1.00000 0.00000 2.53333; 0.00000 1.00000 0.00000 2.48600; -------------------------------- ERROR: matrix is ill-conditioned or badly scaled, condno = 12820.7 -------------------------------- Possible problem with experimental design: Check for duplicate entries and/or lack of range of continuous variables within a class. If you seek help with this problem, make sure to send: 1. Your command line: mri_glmfit --glmdir new.frontal.lgi --fsgd new.frontal.fsgd.txt --fwhm 30 --surf fsaverage lh --y lh.pial_lgi.fwhm0.mgh --C one.mtx --C two.mtx 2. The FSGD file (if using one) 3. And the design matrix above Attempting to diagnose further SumSq: Min=2.828427 (col 1), Max=10.485683 (col 4) ... could not determine the cause of the problem Can anyone give me some insight into what might be going wrong? Thanks so much, Ryan Michael Nillo Staff Research Associate I University of California San Francisco Department of Radiology and Biomedical Imaging External Email - Use Caution

2 5

[mri_cvs_register] Transform is zero
by Stefano Zappala 29 Oct '19

29 Oct '19

External Email - Use Caution Hi there, I am having problem reading the final_CVSmorph_toS131032017_Supine.m3z transform as output of the mri_cvs_register function. I used mri_convert final_CVSmorph_toS131032017_Supine.m3z final_CVSmorph_toS131032017_Supine.nii.gz to convert the transform file to nifti format but what I get it's just a map of zeros. I have checked any different source of error: * mri_convert: it works because if I convert /transforms/talairach.m3z to nifti the conversion is successful; * mri_cvs_register: it works because the two images are aligned very well and correspond perfectly; * final_CVSmorph_toS131032017_Supine: it is not empty because if I run mri_convert -at final_CVSmorph_toS131032017_Supine.m3z firstScan/mri/T1.mgz secondScan/mri/T1mine.mgz the result is morphed correctly. So the transform final_CVSmorph_toS131032017_Supine is actually not populated with zeros. So it could be a problem of type of data: the transform should be indeed populated with very small values, since the difference between the scans even if noticeable, it is very small. I don't know what to try anymore: I was wondering if you have any suggestion? In summary, the call to mri_cvs_register is successful giving me as output a transform in .m3z format which is not zero. When I try to convert it to nifti file, the transformation results full of zeros. 1) FreeSurfer version: freesurfer-Linux-centos6_x86_64-stable-pub-v6.0.0-2beb96c 2) Platform: Ubuntu 18.04.3 LTS 3) uname -a: Linux zappst-Q170M-C 5.0.0-29-generic #31~18.04.1-Ubuntu SMP Thu Sep 12 18:29:21 UTC 2019 x86_64 x86_64 x86_64 GNU/LinuxThu Aug 13 22:55:16 UTC 2015 x86_64 x86_64 x86_64 GNU/Linux 4) recon-all.log for only one of the two images registered: see attached 5) mri_cvs_register.log: see attached Many thanks, Stefano. Stefano Zappalá Stefano Zappalá PhD Student Myfyrwraig PhD School of Computer Science & Informatics Yr Ysgol Cyfrifiadureg a Gwybodeg Cardiff University Prifysgol Caerdydd 5 The Parade 5 The Parade Cardiff Caerdydd CF24 3AA CF24 3AA Tel: +44 (0)7<tel:+44%2029%202087%206070>943096049 Ffôn: +44 (0)79<tel:+44%2029%202087%206070>43096049

2 1

FW: Mapping coordinates from sphere file to ?h.white
by Glasser, Matthew 28 Oct '19

28 Oct '19

External Email - Use Caution From: Timothy Coalson <tsc5yc(a)mst.edu> Reply-To: Timothy Coalson <tsc5yc(a)mst.edu> Date: Monday, October 28, 2019 at 8:21 PM To: "Glasser, Matthew" <glasserm(a)wustl.edu> Subject: Re: FW: [Freesurfer] Mapping coordinates from sphere file to ?h.white That is not a sphere of fs_LR coordinates (in fact, "fs_LR coordinates" is not really a meaningful phrase). That sphere's per-vertex coordinates are specifically intended to spatially line up our fs_LR atlas data to fsaverage atlas data on the fsaverage atlas sphere, because this is the information needed to resample between the atlases. It has the same number of vertices as fs_LR surfaces because it needs to know the target location on the fsaverage atlas sphere for every fs_LR vertex. The inflation of an anatomical surface to a sphere is a tricky process, as the surface can easily fold back over itself (which would encode a non-invertible transformation). The path the vertices took to get to a sphere is neither recorded, nor particularly important to anyone not working on the implementation details of the inflation. The endpoint of the process, the spherical surface, is the only goal, the displacement vectors between them are not useful to surface operations (since the data is tied to the vertices, and all the vertices still exist in the sphere (just at new coordinates), it is trivial to use the same data "on" either the sphere or anatomical surfaces). These spherical surfaces are only used to do registration and resampling (which does not involve any kind of warp of anatomical coordinates), after which the data is used with anatomical surfaces (which always maintain the original shape and size even when resampled), because spatial processing should be done on brain-shaped geometry, not on spheres. It is unclear to me what you intend to use the surface you described for, and I don't know what the .x5 format is. The surface resampling weights are fairly easily derived from the correct pair of spheres, and resampling does not involve any kind of spatial displacement. Tim From: <freesurfer-bounces(a)nmr.mgh.harvard.edu<mailto:freesurfer-bounces@nmr.mgh.harvard.edu>> on behalf of Oscar Esteban <phd(a)oscaresteban.es<mailto:phd@oscaresteban.es>> Reply-To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>> Date: Monday, October 28, 2019 at 6:58 PM To: "freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>" <freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>> Subject: [Freesurfer] Mapping coordinates from sphere file to ?h.white External Email - Use Caution Dear experts, In the context of the BIDS transforms extension and the ".x5" format we chatted about in the workshop we made in March at MIT, we are addressing the following question. We have a sphere of fsLR coordinates in fsaverage space (e.g., https://github.com/Washington-University/HCPpipelines/blob/master/global/te… - but the question extends to the right hemisphere and also resolutions other than 59k). We want to expand this sphere to the fsaverage's white surface (lh.white, in this example). Since coordinates on the fsLR sphere don't match any of fsaverage{,3,4,5,6}'s coordinates, we wonder whether FreeSurfer stores somewhere/somehow the analytical transform of this mapping. In other words, how one would map arbitrary coordinates on the surface to a particular fsaverage's white (or midthickness, or pial, etc.) surface. Otherwise, I guess I would pick fsaverage's surfaces (i.e., 7th degree icosahedron), and interpolate the displacement vectors that start at each sample on the surface and end at the corresponding location of that vertex index of the destination surface (lh.white for the example above). Is there any tool that actually does the _inverse transform of the surface-to-sphere mapping_? Thanks very much. Regards, Oscar -- ___________________________ Oscar Esteban, Ph.D. Postdoctoral Fellow, Poldrack Lab Stanford University +1 (650) 733 33 82 ________________________________ The materials in this message are private and may contain Protected Healthcare Information or other information of a sensitive nature. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this email in error, please immediately notify the sender via telephone or return mail. ________________________________ The materials in this message are private and may contain Protected Healthcare Information or other information of a sensitive nature. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this email in error, please immediately notify the sender via telephone or return mail.

1 0

Re: [Freesurfer] username and password
by fsbuild 28 Oct '19

28 Oct '19

External Email - Use Caution OK, thanks for letting us know.   We will look into when/how that redirect started to happen.  Just to clarify though, you do not need to login to the Martinos center to download the freesurfer release binaries or source code, i.e., using that form was just to get on the freesurfer mail list. The Freesurfer 6.0  release binaries (for Linux and MacOS) can be downloaded from, https://surfer.nmr.mgh.harvard.edu/fswiki/DownloadAndInstall The latest development branch with the freesurfer source code can be obtained from GitHub, https://github.com/freesurfer/freesurfer - R. Rachel S October 2, 2019 at 15:36         External Email - Use Caution        Hi, I have the same problem. It directs me to the Martinos Centre Login, which never happened before. I've subscribed to the maillist. Best,Rachel _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer fsbuild October 1, 2019 at 16:08         External Email - Use Caution        Hello Takashi, Did you create a login to the list on this page, http://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer - and then receive the email notification, and then did you do the confirmation step? - R. _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

3 2

mri_cvs_register error while diffusion imaging processing
by Renew Andrade 28 Oct '19

28 Oct '19

External Email - Use Caution Dear FreeSurfer experts: I am trying to run Spatial Normalization but I cannot get the morph file because I get this error. I don’t know what to do with these libraries. Do I need to copy them from where they are in freesurfer to somewhere else in the Mac? sh-3.2# mri_cvs_register --mov IXI_25_Guys_T1 Using the CVS template as registration target... The outdir: /Volumes/FREESURFER/freesurfer_preprocessed/IXI_25_Guys_T1/cvs Log file is /Volumes/FREESURFER/freesurfer_preprocessed/IXI_25_Guys_T1/cvs/IXI_25_Guys_T1_to_cvs_avg35.mri_cvs_register.1910241722.log Thu Oct 24 17:22:55 CEST 2019 --mov IXI_25_Guys_T1 $Id: mri_cvs_register,v 1.48 2016/03/06 21:31:56 lzollei Exp $ Renes-Air.home Darwin Renes-Air.home 19.0.0 Darwin Kernel Version 19.0.0: Wed Sep 25 20:18:50 PDT 2019; root:xnu-6153.11.26~2/RELEASE_X86_64 x86_64 asegfname aseg Using m3z type morph files. dyld: lazy symbol binding failed: Symbol not found: ___emutls_get_address Referenced from: /Applications/freesurfer/bin/../lib/gcc/lib/libgomp.1.dylib Expected in: /usr/lib/libSystem.B.dylib dyld: Symbol not found: ___emutls_get_address Referenced from: /Applications/freesurfer/bin/../lib/gcc/lib/libgomp.1.dylib Expected in: /usr/lib/libSystem.B.dylib Abort Sincerely, Andrade.

2 1

recon-all at native resolution
by Alexopoulos, Dimitrios 28 Oct '19

28 Oct '19

External Email - Use Caution We have data acquired at 0.8mm isotropic and would like to maintain this resolution in the processing stream rather than FS resampling to 1mm. Is the command to do this 'recon-all -hires -all'? Are there any drawback to doing this? Jim ________________________________ The materials in this message are private and may contain Protected Healthcare Information or other information of a sensitive nature. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this email in error, please immediately notify the sender via telephone or return mail.

2 11

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