Freesurfer

freesurfer@nmr.mgh.harvard.edu

19 participants
30088 discussions

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

view labels in native space
by Nasiriavanaki, Zahra 28 Oct '19

28 Oct '19

Hi Freesurfer experts I have a few labels that correspond to activation clusters in fsaverage space for an individual subject. I used mri_label2label to transfer them to native space of the subject. export SUBJECTS_DIR=/autofs/space/oprah_001/users/zn025/looming_7T/SUBJECTS_DIR mri_label2label --srclabel ./first_levels/$subj/label/aparc/pos/lh-0001.label --trglabel lh-new.label --srcsubject fsaverage --trgsubject ${subj} --regmethod surface --hemi lh When I overlay the native label on top of the fsaverage significance map (tksurfer fsaverage lh inflated -ov sig.nii.gz), it looks like small dots all over the map. Is this because I'm viewing the native label in fsaverage space? Is there a way to view these labels in this subject's native space instead of fsaverage, something like tksurfer native(instead of fsaverage) lh inflated -ov sig.nii.gz? Thank you very much in advance Mona Zahra (Mona) Nasiriavanaki Postdoctoral Research Fellow Martinos Center for Biomedical Imaging Massachusetts General Hospital 149 13th Street, 149-2615 Charlestown, MA, USA, 02129

3 4

Re: [Freesurfer] Error in glmfit-sim cache
by Renew Andrade 28 Oct '19

28 Oct '19

External Email - Use Caution I don’t know what do you mean by not knowing what it means. I said that I have tried all the values the error pinpoints to. That is I have tried mri_glmfit-sim --glmdir lh.IXI.glmdir --cache 4 abs —cap 0.05 --2spaces with "1.3", "2.0", "3.0", “4.0”. And it doesn’t work. I know it doesn’t seem to be a software issue but I cannot run this functionality of the software only for the moment everything else works fine. I have two laptops one Mac and one BMAX I have done recon-all in one and I have done glmfit in the MAC. But this I hope is only temporary because makes no sense that a software works partially. Sincerely, Andrade. ERROR: > thresh = 4,0 must be 1.3, 2.0, 3.0, 4.0 I don't know what this refers to On 10/18/19 4:45 PM, Renew Andrade wrote: > External Email - Use Caution > > I have tried all the values it pinpoints to but none seem to work. > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer <https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer> _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer <https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer>

2 1

Vertex and voxel-wise linear regression in Freesurfer
by Larry Lai 28 Oct '19

28 Oct '19

External Email - Use Caution Hello FreeSurfer Users I'm attempting to explore the correlations between the surface-based and voxel-based values and adjusted for nuisance covariate in FreeSurfer. *Correlational study:* A. Cortical thickness and VBM results B. VBM_1 and VBM_2 results *Platform:* 1) freesurfer-Linux-centos6_x86_64-stable-pub-v6.0.0-2beb96c 2) Ubuntu 16.0.4 Thanks for considering my request. Best regards -- *Ying-Liang (**Larry)** Lai, Ph.D. Student* Research Assistant *Ph.D. Program for Neural Regenerative Medicine* College of Medical Science and Technology Taipei Medical University and National Health Research Institutes

2 1

b0 to T1’s co-registration without b0 inhommogeneities
by Daniel Callow 28 Oct '19

28 Oct '19

External Email - Use Caution Hello, I am trying to do a within-subject co-registration between b0 volumes and T1's for my subjects. This is an old dataset that does not have anything for correcting b0 inhomogeneities. The registration using --dti and normal parameters for bbregister led to "okay" within-subject registration. I am wondering if modifying the --dof (like 12) option for bbregister or applying some type of linear/nonlinear registration might improve the registration? Best, -- *Daniel Callow* *PhD Student, Neuroscience and Cognitive Science* Exercise for Brain Health Lab University of Maryland, College Park *ddcc2442(a)gmail.com <ddcc2442(a)gmail.com>* 443-254-6298

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