Freesurfer

freesurfer@nmr.mgh.harvard.edu

17 participants
30099 discussions

Subfields on 1.5 Tesla
by Alan Francis 23 Apr '14

23 Apr '14

Hi Koen & FreeSurfers: I have received a manuscript to review for possible publication. The authors have used the subfields algorithm on 1.5T scans and obtained a parcellation with values. They have drawn some major conclusions on the basis of the findings. My understanding is that this method can only be done on 3T. Is the 1.5T results valid? Please advice. thanks, Alan Francis

1 0

Re: [Freesurfer] mri_convert Siemens .IMA error
by Douglas N Greve 22 Apr '14

22 Apr '14

Hi Paul, I was finally able to track down the bug. You can get a new version from here: ftp://surfer.nmr.mgh.harvard.edu/pub/dist/freesurfer/5.3.0-patch/mri_conver… copy it to $FREESURFER_HOME/bin/mri_convert doug On 12/08/2013 06:28 PM, prasser wrote: > Hi Doug, > > Sorry for the delay. Please find 1.IMA attached. > > Thanks, > Paul > > > ---- On Fri, 06 Dec 2013 13:15:29 -0800 *Douglas N Greve > <greve(a)nmr.mgh.harvard.edu>* wrote ---- > > > Can you upload a single dicom file to me? One should be enough > > On 12/06/2013 04:02 PM, Paul Rasser wrote: > >> > >> We are using Freesurfer5.3 on os x 10.8. > >> > > > > > > On 6 Dec 2013, at 2:19 pm, Douglas Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > > <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > > >> > >> This may be a known problem. What version of FS are you using and > >> what is your platform? > >> doug > >> > >> > >> > >> > >> On 12/5/13 3:26 PM, prasser wrote: > >>> Thanks, but the orientation shift appears as a skew. Also, we > don't > >>> have any fiducial markers for this group of scans so it would > be a > >>> little tricky to confirm the correct application of > >>> -out_orientation. I would rather let mri_convert handle this from > >>> the .IMAs. > >>> > >>> Is there any additional information I can provide from from > the .IMA > >>> header etc that could shed some light on why mri_convert is > behaving > >>> in this way? > >>> > >>> Thanks, > >>> Paul > >>> > >>> > >>> > >>> ---- On Wed, 04 Dec 2013 22:11:27 -0800 *Mark Plantz > >>> <markplantz2016(a)u.northwestern.edu > <mailto:markplantz2016@u.northwestern.edu>>* wrote ---- > >>> > >>> Not sure if this helps, but you can alter the output orientation > >>> using the flag "--out_orientation" followed by LIA (for > >>> example). Maybe the input file needs to be in a specific > >>> orientation before it can be converted properly. > >>> > >>> > >>> On Wed, Dec 4, 2013 at 9:56 PM, prasser <prasser(a)zoho.com > <mailto:prasser@zoho.com> > >>> <mailto:prasser@zoho.com>> wrote: > >>> > >>> Hi, > >>> > >>> I'm currently having problems converting .IMA files > >>> collected from a Siemens Skyra to .mgz format using > >>> mri_convert from Freesurfer5.3 on os x. > >>> > >>> Using: > >>> > >>> mri_convert -it siemens -ot mgz 1.IMA test.mgz > >>> > >>> I get the error message: > >>> > >>> does not contain a Siemens ASCII header has this file been > >>> anonymised? > >>> > >>> ERROR: cannot unpack mosiacs without ASCII header > >>> > >>> The .IMA files have been taken straight from the scanner > >>> console. > >>> > >>> If I try: > >>> > >>> mri_convert -it dicom -ot mgz 1.IMA test.mgz > >>> > >>> this creates the .mgz, however there is an orientation shift > >>> and a vertical stripe artefact appears in the sagittal view > >>> (please see attached) when compared to the original .IMA. > >>> > >>> Any suggestions on how to convert the Siemens .IMA files > >>> correctly? > >>> > >>> Thanks, > >>> > >>> Paul > >>> > >>> > >>> _______________________________________________ > >>> Freesurfer mailing list > >>> Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@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. > >>> > >>> > >>> _______________________________________________ > >>> Freesurfer mailing list > >>> Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@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. > >>> > >>> > >>> > >>> > >>> _______________________________________________ > >>> Freesurfer mailing list > >>> Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > >>> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >> > >> _______________________________________________ > >> Freesurfer mailing list > >> Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@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. > > -- > 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 > FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > Outgoing: > ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

2 1

Re: [Freesurfer] Extract vertex-wise cortical thickness, surface area and lgi values in ascii files?
by sabin khadka 22 Apr '14

22 Apr '14

Hi all, I ran mri_surf2surf and then mris_convert to get vertexwise thickness and surface area values in ascii files. I found that some vertices had negative values of surface area for few subjects. Should not surface area values be positive only? -Sabin On Thursday, April 17, 2014 10:30 AM, sabin khadka <sabinkhdk(a)yahoo.com> wrote: Ahh. Great. Thanks you very much for the information. I ran mri_surf2surf and then mris_convert to get thickness and surface area vertex wise values. I found that some vertices had negative values of surface area for few subjects. Should not surface area values be positive only? -Sabin On Thursday, April 17, 2014 10:10 AM, Bruce Fischl <fischl(a)nmr.mgh.harvard.edu> wrote: yes, you can use mri_surf2surf first to map them to fsaverage, then convert them to ascii On Thu, 17 Apr 2014, sabin khadka wrote: > Hi Bruce, > Thanks for the reply. I was wondering if I have to somehow resample or map > to fsaverage so as to have same # of vertices across all subjects. I tried > using mris_convert as you suggested above but I can see that different > subjects have different #of vertices. > > Thanks for help. > > -Sabin > On Thursday, April 17, 2014 9:29 AM, Bruce Fischl > <fischl(a)nmr.mgh.harvard.edu> wrote: > Hi Sabin > > mris_convert can convert binary "curvature" format files to ascii ones if > the output file has the extension asc. Something like: > > mris_convert -c lh.thickness lh.orig lh.thickness.asc > > this will give you a 5 column ascii file in the format > <vertex index> <x> <y> <z> <thickness> > > cheers > Bruce > > On Thu, 17 Apr 2014, > sabin khadka > wrote: > > > Hi all, > > I want to extract to extract thickness, area and lgi values vertex wise in > > ascii files for group of subjects. Could anyone suggest me on how to do > it. > > > > Thank you for your help. > > > > -Sabin > > > > > > > > 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. > > > > >

2 1

CUDA PC build
by Chiu, Bryan (PHTH) 22 Apr '14

22 Apr '14

Hi, Trying to build a CUDA enabled PC for the lab to do processing for a variety of neuroimaging programs including Freesurfer and FSL. Want to know input from others about the 'recommended' PC and what the 'standard' is based on the current PC prices and technology. Sample build is from NCIX Canada, if curious. SKU Description Price Amount 84500 ASUS Z87-A ATX LGA1150 Z87 DDR3 3PCI-E16 2PCI-E1 2PCI CrossFireX/SLI SATA3 USB3.0 HDMI Motherboard [Reg. $159.99] $154.99 $154.99 93608 ADATA Xpg V1.0 16GB DDR3-1600 CL9-9-9-24 2X8GB Dual Channel Memory Kit Black [Reg.$189.99] $149.99 $149.99 82953 ASUS DRW-24F1ST 24X SATA DVD Writer Black $24.98 $24.98 65701 Seagate Barracuda 1TB 7200RPM 64MB SATA 6Gbps 3.5in Internal Hard Drive - OEM [Reg. $73.44] $64.99 $64.99 76367 ASUS GeForce GTX 660 OC DirectCU 1020MHZ 2GB 6.0GHZ GDDR5 2xDVI HDMI DisplayPort PCI-E Video Card $229.99 $459.98 81327 Intel Core i5 4670K Unlocked Quad Core 3.4GHZ Processor LGA1150 Haswell 6MB Cache Retail [Reg. $269.98] $249.99 $249.99 38650 Corsair Professional HX850W 850W ATX 12V 70A 24PIN ATX Modular Power Supply Active PFC 140mm Fan $169.98 $169.98 67326 Antec Three Hundred Two Mid Tower Gaming Case 302 ATX 3X5.25 6X3.5INT No PS Front USB3.0 & Audio [Reg. $77.98] $69.99 $69.99 In short -- runs with 16 Gb RAM, 3.4 GHz Quad core, with 2 GTX 660s in SLI. Final quote comes out to $1506 CAD. This is my first time attempting to build a PC of this type, so I still have tons of questions regarding minor details. I think spec wise the computer is good, unless others have recommendations from previous experience. Most likely will be running a flavour of Linux on this rig. Any suggestions for what the 'standard' for this build would be? I'm thinking CentOS based on older webmail threads, but I would prefer Ubuntu if possible. I would also like it to be user friendly so the technologically challenged researchers can get the software up and running. Is there anything else I should think about before I order? I want to make sure I have all of my bases covered before ordering and getting this rig going. Regards, - Bryan Chiu

1 0

recon all problem
by Danielle Miller 22 Apr '14

22 Apr '14

Hi Freesurfer Experts, I seem to be having some difficulty with reconing a healthy young adult. I have tried running this on multiple different computers and multiple different versions (5.3 and 5.1) and all return the same weird recon (see attached image). I have used the same sequence and commands on another subject and the recon returned fine. This participant did not have excess motion and the orig.mgz file looks fine. Please let me know if there is any more information I could provide that would be of use. Commands used: 5.3: recon-all -all -FLAIR (FLAIR image) -subjid (SUBJID) -FLAIRpial -hippo-subfields recon-all -all -subjid (SUBJID) 5.1: recon-all -all -subjid (SUBJID) -hippo-subfields recon-all -all -subjid (SUBJID) Thanks, Danielle Miller -- Ph.D. Program in Behavioral Neuroscience Boston University School of Medicine L-815 72 E. Concord St Boston, MA 02118 VA Boston Healthcare System Jamaica Plain Memory Disorders Research Center 150 South Huntington Ave D11-103 Boston, MA 02130 OFFICE:(857) 364-2130

2 1

Re: [Freesurfer] viewing certain regions
by Douglas N Greve 22 Apr '14

22 Apr '14

HI Courtney, please post to the list. Thanks! Get the codes from $FREESURFER_HOME/FreeSurferColorLUT.txt. Also, you can use multiple numbers to match all the ROIs you want at once (if that is what you want). Eg, --match 2 26 17 doug On 04/22/2014 03:16 PM, Courtney Haswell wrote: > Thanks for the reply. I was able to make the hippocampus and amygdala masks easily, but was having some trouble for the other ones. > > Looking in the "aparc.annot.ctab" file, I see that caudalanteriorcingulate is labeled 2 and rostralanteriorcingulate is 26. So then is this the correct way to get the ROI extracted: > mri_binarize --i $SUBJDIR/aparc+aseg.mgz --o $ROIDIR/3_rostralACC.mgz --match 26 > mri_binarize --i $SUBJDIR/aparc+aseg.mgz --o $ROIDIR/3_caudalACC.mgz --match 2 > > > > I also did the caudal ACC (that is in aparc.annot) the same way but that ROI covers all of the white matter in the right hemisphere so I think I did something wrong. > mri_binarize --i $SUBJDIR/aparc+aseg.mgz --o $ROIDIR/3_caudalACC.mgz --match 2 > > mri_mask $SUBJDIR/aparc+aseg.mgz $ROIDIR/3_caudalACC.mgz $ROIDIR/acc_mask5.mgz > > > > Thank you, Courtney > > > > > Use mri_binarize with the --match option to create a mask of the > regions you want. Then use mri_mask applied to aparc+aseg to create a > new volume with only the regions you want > doug > On 04/22/2014 09:06 AM, Courtney Haswell wrote: > > Hi! > > > > I want to view only certain segmented regions on one of our subject's > > brains. It seems like it should be simple but I am struggling to find > > the steps. I want to display only hippocampus, amygdala, medial > > orbitofrontal, lateral orbitofrontal, sulcus/gyri anterior cingulate, > > subcallosal gyrus, and rACC and cACC. > > > > How do I extract these labels so that I can only show just these > > regions on the brain? I extracted stats for these areas, but want to > > visually show the locations. > > > > Thank you, > > Courtney -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

1 0

2 available spots in Boston FreeSurfer course
by Allison Stevens 22 Apr '14

22 Apr '14

We have had some last minute drop outs and there are now 2 spots available in the upcoming Boston FreeSurfer course (May 5-7th). If you would like to sign up for these spots, please register here: https://surfer.nmr.mgh.harvard.edu/cgi-bin/fsurfer/register.8263.cgi Note that credit card cannot be used to pay for the course. Also, this course will not have the diffusion analysis component (it will focus on fMRI analysis). Allison

2 3

viewing certain regions
by Courtney Haswell 22 Apr '14

22 Apr '14

Hi! I want to view only certain segmented regions on one of our subject's brains. It seems like it should be simple but I am struggling to find the steps. I want to display only hippocampus, amygdala, medial orbitofrontal, lateral orbitofrontal, sulcus/gyri anterior cingulate, subcallosal gyrus, and rACC and cACC. How do I extract these labels so that I can only show just these regions on the brain? I extracted stats for these areas, but want to visually show the locations. Thank you, Courtney

2 1

Re: [Freesurfer] RAS to MNI coordinate conversion function?
by Douglas N Greve 22 Apr '14

22 Apr '14

Hi John, please remember to post to the list. thanks! Yes, that should work. Use read read_surf.m to get the proper coordinates. doug On 04/20/2014 08:08 PM, John Pyles wrote: > Hi Doug, > > Thanks so much for the reply. > > I had looked through the wiki on coordinates but because I'm not a big > Freesurfer user was unsure which was the correct transformation. > > It seems like for what I need to do the example below is appropriate? > So in theory I could write a Matlab function using the calculation > below correct? > > Thanks again. > > Best, > John > > > ------------ > 2. I have an RAS point on the surface (tkrR tkrA tkrS) ("Vertex RAS" > from tksurfer) and want to compute the MNI305 RAS that corresponds to > this point: > > MNI305RAS = TalXFM*Norig*inv(Torig)*[tkrR tkrA tkrS 1]' > > TalXFM: subject/orig/transforms/talairach.xfm Norig: mri_info > --vox2ras orig.mgz Torig: mri_info --vox2ras-tkr orig.mgz > ------------ > > > ______________________________________ > > John Pyles, Ph.D. > Research Scientist > Center for the Neural Basis of Cognition > Department of Psychology > Carnegie Mellon University > email: jpyles(a)cmu.edu > phone: 206.552.0107 > ______________________________________ > > > On Fri, 18 Apr 2014, Douglas Greve wrote: > >> >> There is not a simple function, but you can look at our wiki page on >> coordinate systems. There are some examples there on >> tranforming from one coordinate to another >> doug >> >> >> On 4/18/14 2:40 PM, John Pyles wrote: >> Hi, >> >> I am not a heavy Freesurfer user and I haven't had good luck figuring >> this out through the FreeSurfer documentation and I >> thought an expert might have a fast answer. If anyone has an easy >> solution it would be very much appreciated. >> >> While I use FreeSurfer for automatic segmentation, I don't normally >> use if for analysis. However for an ECoG project I'm >> using parts of FreeSurfer to aid in electrode coregistration. >> >> Basically the issue is that I have electrode coordinates in what I >> believe are FreeSurfer RAS space and I need to convert >> them to MNI space. I have many of these electrode coordinates for >> each subject so I need a FreeSurfer or Matlab function to >> automate this process. Is there a simple way to do this? >> >> The *very* kludgy way I figured out how to do this "manually" is as >> follows (I'm sure this is going to seem very dumb to a >> FreeSurfer user, but I needed a quick solution that worked for a >> figure). >> 1) Enter the coordinate output from my electrode coregistration >> routines manually into the RAS field in Freeview. >> 2) Take the "brain volume" coordinates (all positive) in Freeview and >> enter them into the "Volume Coordinate" field in >> TkMedit. >> 3) Copy the TAL coordinates output in MkEdit when the brain volume >> coordinate is entered. >> 4) Convert the TAL coordinates to MNI in Matlab with tal2mni.m >> >> It seems like there should be a simple function to do this coordinate >> conversion using the transformation info from >> recon-all that I'm assuming both Freeview and TkMedit use, but I >> can't seem to find it? >> >> Any help would be appreciated. Thanks again! >> >> Best, >> John >> >> ______________________________________ >> >> John Pyles, Ph.D. >> Research Scientist >> Center for the Neural Basis of Cognition >> Department of Psychology >> Carnegie Mellon University >> ______________________________________ >> >> >> >> >> >> >> >> _______________________________________________ >> Freesurfer mailing list >> Freesurfer(a)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 Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

1 0

ROI analysis in Matlab
by Reza Rajimehr 22 Apr '14

22 Apr '14

Hi, For my dataset, it appears to be easier to do an ROI analysis in Matlab (instead of using funcroi commands). I just want to make sure the procedure described below is right: 1) I have a label/ROI on the right hemisphere surface of a subject (e.g. MT from the motion localizer). I convert it to a mask file (MT.mask.nii.gz). I load it to Matlab using MRIread. In an fMRI experiment, I have 3 conditions. For the right hemisphere vertices of the subject, I compute cespct.nii.gz for the contrasts of each condition vs. baseline (1_vs_0 , 2_vs_0 , 3_vs_0). Then I load the three cespct files to Matlab using MRIread. 3) For the vertices that are 1 in the MT.mask.nii.gz, I calculate the averaged cespct across those vertices, in each contrast. So I will have % signal change values for the three conditions. 4) I repeat this for all subjects and do statistics across subjects. Any comments? Reza

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