Freesurfer

freesurfer@nmr.mgh.harvard.edu

4 participants
30091 discussions

LME and clusterwise correction
by Benjamín Garzón 12 Apr '13

12 Apr '13

Hi, I am doing a longitudinal analysis using the LME toolbox, and I would like to perform a clusterwise correction for multiple comparisons. Would it be correct to use mri_surfcluster on the resulting maps, or do you suggest any alternative otherwise ? I can only find functions implementing FDR correction in the toolbox. Thanks, Benjamin -- Benjamín Garzón, Ph.D. Karolinska Institutet Department of Neuroscience Retzius Väg 8 17177 Stockholm (SWEDEN)

2 1

bblabel values for DLPFC
by Anupa AV 12 Apr '13

12 Apr '13

Dear All, I 'd like to know the xmin, ymin and zmin for creating DLPFC label using bblabel string. I want to make a coronal cut at y=26 for seperating DLPFC from premotor cortex. I tried by giving a value of xmin=0, ymin=26 and zmin= 0. I'd like to know whether I'm on right track. Thanks in advance for your response.

2 2

Fwd: Sample size estimation from LME models
by nicolas guizard 11 Apr '13

11 Apr '13

Hi, I am new to Freesurfer and the Linear Mixed Effects (LME) Models packages. And I would like to compute sample size from volumetric measures for 2 groups. I followed the longitudinal example provided in the wiki ( http://surfer.nmr.mgh.harvard.edu/fswiki/LinearMixedEffectsModels#a.29Univa…) with the ADNI example. But I am not sure what to use for the "common random effect design matrix" Zi and "column" ZiCol with the lme_plannedSampleSize(Zi,ZiCol,Dhat,phisqhat,effsz,dr,pw,alpha,gr_pr) function: Here is what I do in matlab: > load('ADNI791_Hipp_and_Entorh.mat') > total_hipp_vol_stats = lme_fit_EM(X_Hipp,[1,2],Y(:,1)+Y(:,2),ni); > lme_plannedSampleSize(total_hipp_vol_stats.X,[1,2] ,total_hipp_vol_stats.Dhat,total_hipp_vol_stats.phisqhat,1,.2,.8,0.05) > ans = 1.0e+07 * 1.3584 0.0154 0.0154 0.0100 Is that right? Regards, Nicolas

2 1

"final" skull-stripped file?
by Marco Loggia, PhD 11 Apr '13

11 Apr '13

Dear all, what is the 'best' skull stripped file produced by freesurfer? From what I see, it would seem to be 'brainmask'. However this file -if I understand correctly- is produced in the very first stages of recon-all, and therefore does not take advantage of the detailed reconstruction of the brain 'borders' which is achieved in the subsequent stages of recon-all. For this reason, I have been masking my brainmask.mgz file with a binary file created using the aparc+aseg file (which labels only the brain), slightly smoothed to be conservative in the stripping. Or is there already a file skull-stripped this way? Thanks, Marco ______________________________ Marco L. Loggia, PhD Instructor in Radiology Harvard Medical School Massachusetts General Hospital 149 Thirteenth Street, Room 2301 Charlestown, MA 02129 Phone: (617) 643-7267 Fax: (617) 726-7422 marco(a)nmr.mgh.harvard.edu

2 1

Re: [Freesurfer] Transform SMOOTHWM to DTI B0
by Daniel Haehn 11 Apr '13

11 Apr '13

Thanks, Satra! This is helpful!! On Thu, Apr 11, 2013 at 12:32 PM, Satrajit Ghosh <satra(a)mit.edu> wrote: > hi all, > > here you go: > > https://gist.github.com/satra/1451812 > > this goes from trackvis streamlines to freesurfer surface space. you > would do the inverse of this transform to go the other way. > > cheers, > > satra > > On Thu, Apr 11, 2013 at 12:26 PM, Daniel Haehn < > daniel.haehn(a)childrens.harvard.edu> wrote: > >> Hi Satra et al, >> >> thanks a lot and it would be great if you could find some code! >> >> In the meantime, I tried the following based on information from Kiho >> Im: >> >> # we have to first remove transformation of 'vertex to RAS' embedded in >> smoothwm surface. >> mris_transform lh.smoothwm identity.xfm ./lh.smoothwm.clean >> >> # then use the flirt T1-to-b0 matrix converted from .mat to .xfm format >> to transform the surface >> mris_transform lh.smoothwm.clean T1-to-b0.xfm lh.b0.smoothwm >> >> And of course: it didn't work.. In the attached screenshot, the red >> surface is the original one which kind of fits the b0 image but not very >> well (it obviously fits the original T1 very well). The blue is the >> transformed surface (lh.b0.smoothwm) and is completely off. >> >> [image: Inline image 1] >> >> Do you know what I could be doing wrong? I also attached both >> tranforms. Maybe I just converted them wrong. >> >> Thanks!! >> Daniel >> >> >> On Thu, Apr 11, 2013 at 10:53 AM, Satrajit Ghosh <satra(a)mit.edu> wrote: >> >>> hi dan, >>> >>> if you use bbregister to match the DTI b0 to freesurfer space, then >>> you can transform the surface back (no guarantees on the tesselations >>> though i think). the only additional component here would be a tkr >>> transform that's based on the conformed space. i might even have some >>> python code lying around to do this. i'll take a look. >>> >>> cheers, >>> >>> satra >>> >>> On Thu, Apr 11, 2013 at 10:40 AM, Daniel Haehn < >>> daniel.haehn(a)childrens.harvard.edu> wrote: >>> >>>> Hi all, >>>> >>>> I want to transform the Freesurfer Smoothwm surfaces to DTI B0 space. >>>> >>>> I found >>>> http://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferTrackVisTransformsbut this describes transforming the B0 to T1 space and then match the >>>> surfaces. I don't want to change the DTI B0 image. >>>> >>>> What is the best way to do that? >>>> >>>> Thank you very much! >>>> Daniel >>>> >>>> -- >>>> Daniel Haehn >>>> FNNDSC / BCH >>>> +1.857.218.5140 >>>> >>>> _______________________________________________ >>>> 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. >>>> >>>> >>> >> >> >> -- >> Daniel Haehn >> FNNDSC / BCH >> +1.857.218.5140 >> > > -- Daniel Haehn FNNDSC / BCH +1.857.218.5140

1 0

DTI: problems with (high) FA values
by Rotem Saar 11 Apr '13

11 Apr '13

Hi Douge, In your last reply to my email (that contained questions regarding some really high FA values - see recent emails below), u asked me to write u the first value in register.dat.mincost - I'm copying the line I have: 0.280554 140.210312 156.795241 11.152485 Does this takes us one step further towards understanding the problem and its solution ? Thanks for your time, Rotem On 04/07/2013 06:48 AM, Rotem Saar wrote: > Hi Doug, > > Thanks for your answer. > I changed the last step as u suggested, but I think I need to clarify > my question: > Indeed values of the CC looks OK, but don't u think other values, like > the ventricles are too high ? *yes, the ventricles look too high* > > When I first run the script, I got all values and looked the the CC - > it looked fine. Then, I wanted to validate with some other regions, > just to show that I got the appropriate (low) FA values for regions I > don't expect to see high values in - and things looked odd (too high). > > I remember Bruce also writing me that the values seems a little high, > but we didn't further discuss. I read the following this link: > http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/MultiModalDtiIndividual > > "The CC has an average FA of about 0.75, gyral parcellations are about > 0.4, the left putamen is 0.27, and the ventricle is 0.2. This is as > expected because the CC is highly directional with no crossing fibers > so we would expect the CC to have the highest FA. The gyral white > matter is also directional but has fibers crossing in them, so one > expects the FA to be lower than CC. The gray matter (putamen) is still > lower. The ventricle has no fibers, so we expect it to have the lowest > FA". > > All my values are above 0.4.... > > Additionally, if in the last step I'm writing --seg > /usr/local/freesurfer/ subjects/Rotem_try/mri/wmparc.mgz > > why do I get 182 regions ? in my wmparc.stats I have only 70... The wmparc includes all the subcortical and cortical regions too, but only the WM regions are reported in the wmparc.stats file. *To only report the WM regions, add to mri_segstats --ctab $FREESURFER_HOME/**WMParcStatsLUT.txt* > > I'm probably doing something wrong, but I can't really point to the > problem, thus asking for your help. *I don't know either. If I had to guess, I'd say it is a registration problem. What is the first value in register.dat.mincost?* doug > > Attached is the table again, > > Thanks > Rotem > > > The values look about right in the table. Your pipeline looks ok, > thought the last step uses fa_FOLDER-NAME.mask.nii instead of the > output > of mri_vol2vol (fa_FOLDER-NAME.nii). > doug > > > > On 03/17/2013 02:16 AM, Rotem Saar wrote: > > > > Hi all, > > I run into somthing that seems odd to me and wanted to consult - > > I run the following script for getting the FA values from my DTI > > scans: > > > > 1) dt_recon --i > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/I00001.dcm --s > > FOLDER-NAME --o > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI --b > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/bvals.bval > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/bvecs.bvec > > 2) tkregister2 --s fsaverage --surf white --reg > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa-tal.nii.reg > > --mov /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa-tal.nii > > 3) tkmedit FOLDER-NAME orig.mgz -aux brain.mgz -seg wmparc.mgz > > -reg /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/register.dat > > -overlay /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa.nii > > -fthresh 0.2 -fmax 1 > > 4) mri_vol2vol --mov > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa.nii --targ > > /usr/local/freesurfer/subjects/FOLDER-NAME/mri/orig.mgz --s > > FOLDER-NAME --interp nearest --o > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa_FOLDER-NAME.nii > > --reg > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/register.dat > > 5) tkregister2 --mov > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa_FOLDER-NAME.nii > > --targ /usr/local/freesurfer/subjects/FOLDER-NAME/mri/orig.mgz > > --reg > > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa_FOLDER-NAME.nii.reg > > 6) mri_segstats --seg > > /usr/local/freesurfer/subjects/FOLDER-NAME/mri/wmparc.mgz --ctab > > $FREESURFER_HOME/FreeSurferColorLUT.txt --i > > > /usr/local/freesurfer/subjects/FOLDER-NAME/DTI/fa_FOLDER-NAME.mask.nii > > --sum > > > /usr/local/freesurfer/subjects/FOLDER-NAME/stats/all_stats_fa_FOLDER-NAME > > > > > > I got a table with all the FA values (see attached), for each > > segment, but I suspect a problem: I think that the values > are too > > high (I set the threshold to 0.2-1), am I right ? > > Can u guide me regarding what I can do to solve the problem ? in > > addition I attached a figure of the corpus-callosum, in > which I'm > > interested - i't seems that the green area in not big enough, is > > this fixable ? > > > > Thanks ! > > > > Rotem > > > > > _______________________________________________ > 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: www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

2 1

mri_aparc2aseg issue
by Paul Beach 11 Apr '13

11 Apr '13

Hi Freesurfers, I'm currently trying to take manually made ROIs (label files) that have been incorporated into new annotation files (one for each hemisphere - rh.aug_aparc.annot & lh.aug_aparc.annot) and turn them into a new (augmented) annotation into a aseg file (an aug_aparc+aseg.mgz file). To do so I've been using the mri_aparc2aseg command and, at face value, it works well. However, some of my ROIs (specifically my insular ROIs) are being either shifted or lost in the transition from annotation to segmentation. Note that each hemisphere's new annotation files look perfect on tksurfer (which is where they were originally drawn). It's only when I try and make them fill the full cortical ribbon with the mri_aparc2aseg command that I get this issue. I piloted everything out on averaged subject and it all worked perfectly. However, when I attempted things on individual subject data I received the problem described above. I've attached a picture illustrating the discrepancy between the an averaged subject and a single subject. I have also attached my script code. The overall goal is to create a gmroi_volume which will allow me to determine an aux number of each manual ROI so as to export each one to AFNI space. I'd *really* appreciate help on how to fix this. Thanks, Paul

2 2

Re: [Freesurfer] Transform SMOOTHWM to DTI B0
by Daniel Haehn 11 Apr '13

11 Apr '13

Hi Satra et al, thanks a lot and it would be great if you could find some code! In the meantime, I tried the following based on information from Kiho Im: # we have to first remove transformation of 'vertex to RAS' embedded in smoothwm surface. mris_transform lh.smoothwm identity.xfm ./lh.smoothwm.clean # then use the flirt T1-to-b0 matrix converted from .mat to .xfm format to transform the surface mris_transform lh.smoothwm.clean T1-to-b0.xfm lh.b0.smoothwm And of course: it didn't work.. In the attached screenshot, the red surface is the original one which kind of fits the b0 image but not very well (it obviously fits the original T1 very well). The blue is the transformed surface (lh.b0.smoothwm) and is completely off. [image: Inline image 1] Do you know what I could be doing wrong? I also attached both tranforms. Maybe I just converted them wrong. Thanks!! Daniel On Thu, Apr 11, 2013 at 10:53 AM, Satrajit Ghosh <satra(a)mit.edu> wrote: > hi dan, > > if you use bbregister to match the DTI b0 to freesurfer space, then you > can transform the surface back (no guarantees on the tesselations though i > think). the only additional component here would be a tkr transform that's > based on the conformed space. i might even have some python code lying > around to do this. i'll take a look. > > cheers, > > satra > > On Thu, Apr 11, 2013 at 10:40 AM, Daniel Haehn < > daniel.haehn(a)childrens.harvard.edu> wrote: > >> Hi all, >> >> I want to transform the Freesurfer Smoothwm surfaces to DTI B0 space. >> >> I found >> http://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferTrackVisTransformsbut this describes transforming the B0 to T1 space and then match the >> surfaces. I don't want to change the DTI B0 image. >> >> What is the best way to do that? >> >> Thank you very much! >> Daniel >> >> -- >> Daniel Haehn >> FNNDSC / BCH >> +1.857.218.5140 >> >> _______________________________________________ >> 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. >> >> > -- Daniel Haehn FNNDSC / BCH +1.857.218.5140

2 1

Transform SMOOTHWM to DTI B0
by Daniel Haehn 11 Apr '13

11 Apr '13

Hi all, I want to transform the Freesurfer Smoothwm surfaces to DTI B0 space. I found http://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferTrackVisTransforms but this describes transforming the B0 to T1 space and then match the surfaces. I don't want to change the DTI B0 image. What is the best way to do that? Thank you very much! Daniel -- Daniel Haehn FNNDSC / BCH +1.857.218.5140

3 2

mgz to dicom
by Varghese Chikku 11 Apr '13

11 Apr '13

Hi All, Is there any tool,which can convert mgz files to Dicom? In Thanks chikku

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