Freesurfer

freesurfer@nmr.mgh.harvard.edu

11 participants
30080 discussions

Re: ffx bug
by Doug Greve 13 May '08

13 May '08

There were a bunch of voxels with no significance because the ffx noise had essentially been multiplied by 100 (so 100 times more noise than there really was). make sense? xiaomin wrote: > Hi Doug, > > > > Thanks for fixing it. I will rerun the group analysis and let you > know if it works. However, it is still puzzling to me why the > significant histogram from the ffx without percent (as attached in my > last email) looks so weird. The maximum is too large. Also, the > shape of the significant ffx histogram should not be different for > percent vs. non-percent. > > > > Xiaomin > > > > From: Doug Greve [mailto:greve@nmr.mgh.harvard.edu] > Sent: Tuesday, May 13, 2008 5:13 PM > To: xiaomin > Cc: analysis-bugs(a)nmr.mgh.harvard.edu; 'Roger Tootell'; 'Daphne Holt'; > 'Brittany Cassidy'; Free Surfer > Subject: Re: ffx bug > > > > > Hi Xiaomin, I finally got around to looking at this, and this is a > bug. The "cespct" is not actually a percent, rather it is a fraction > of the mean functional. However, the "cesvarpct" is a percent, which > makes the group ffx analysis with percent fail. I have fixed this in > both dev and stable such that cespct is actually a percent, live on > next update. Let me know if this does not fix your problem. > > doug > > > > > xiaomin wrote: > >Hi, > > > >I am running the following command for group fixed factor and random factor > >analysis with and without percent parameter on. > > > >Here are the commands: > >1) for the no percent: > >isxconcat-sess -sf subjlist -df subjdir -analysis RotDepth -all-contrats > >-hemis -o . -subject fsaverage > >mri_glmfit --y rh.ces.nii --osgm --glmdir rh.ffx.osgm --nii --mask > >../rh.mash.nii --yffxvar rh.cesvar.nii --ffxdofdat ../ffxdof.dat --surf > >fsaverage rh > > > >2) for the percent: > >isxconcat-sess -sf subjlist -df subjdir -analysis RotDepth -all-contrats > >-percent -hemis -o . -subject fsaverage > >mri_glmfit --y rh.cespct.nii --osgm --glmdir pct.rh.ffx.osgm --nii --mask > >../rh.mash.nii --yffxvar rh.cesvarpct.nii --ffxdofdat ../ffxdof.dat --surf > >fsaverage rh > > > >However, the significance between percent and no percent is inconsistent. > >For testing if the data has problem, I run it with rfx design. The > >significance between percent and no percent is almost same. (see attached > >graph) Second, the significance with percent parameter for the fixed factor > >analysis is very small, which is abnormal. We need to use ffx significant > >map to define a ROI. > > > >The data is here: > > > >/autofs/space/yneuron_001/users/xiaomin/projects/FaceRotation/subjects/Rot_D > >epth. > > > > > >Xiaomin Yue > >Research Fellow > >NMR center, MGH > > > > > > > > > > ------------------------------------------------------------------------ > > > > > >-- > >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 > > > >In order to help us help you, please follow the steps in: > >surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > > > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 In order to help us help you, please follow the steps in: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting

1 0

Re: ffx bug
by Doug Greve 13 May '08

13 May '08

Hi Xiaomin, I finally got around to looking at this, and this is a bug. The "cespct" is not actually a percent, rather it is a fraction of the mean functional. However, the "cesvarpct" is a percent, which makes the group ffx analysis with percent fail. I have fixed this in both dev and stable such that cespct is actually a percent, live on next update. Let me know if this does not fix your problem. doug xiaomin wrote: >Hi, > >I am running the following command for group fixed factor and random factor >analysis with and without percent parameter on. > >Here are the commands: >1) for the no percent: >isxconcat-sess -sf subjlist -df subjdir -analysis RotDepth -all-contrats >-hemis -o . -subject fsaverage >mri_glmfit --y rh.ces.nii --osgm --glmdir rh.ffx.osgm --nii --mask >../rh.mash.nii --yffxvar rh.cesvar.nii --ffxdofdat ../ffxdof.dat --surf >fsaverage rh > >2) for the percent: >isxconcat-sess -sf subjlist -df subjdir -analysis RotDepth -all-contrats >-percent -hemis -o . -subject fsaverage >mri_glmfit --y rh.cespct.nii --osgm --glmdir pct.rh.ffx.osgm --nii --mask >../rh.mash.nii --yffxvar rh.cesvarpct.nii --ffxdofdat ../ffxdof.dat --surf >fsaverage rh > >However, the significance between percent and no percent is inconsistent. >For testing if the data has problem, I run it with rfx design. The >significance between percent and no percent is almost same. (see attached >graph) Second, the significance with percent parameter for the fixed factor >analysis is very small, which is abnormal. We need to use ffx significant >map to define a ROI. > >The data is here: > >/autofs/space/yneuron_001/users/xiaomin/projects/FaceRotation/subjects/Rot_D >epth. > > >Xiaomin Yue >Research Fellow >NMR center, MGH > > > > > > ------------------------------------------------------------------------ > -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 In order to help us help you, please follow the steps in: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting

1 0

imensions of an mgh file generated by
by Lars M. Rimol 13 May '08

13 May '08

Hi, I'm just wondering what the dimensions of an mgh file generated by mris_preproc with e.g. 300 subjects should be? Should it be 163842x1 or 163842x300 ? I would assume that it should contain a vector for each vertex, with the same length as number of subjects that were listed in the fsgd file used with mris_preproc. But when I run mris_preproc on 300 subjects I get a 163842x1 file (as read in Matlab). Thank you! -- yours, LMR

1 0

mgz to dicom or NiFTI
by Stoffers, D. 13 May '08

13 May '08

Dear all, Our imaging techs use a specific artifact correction which requires converting the DICOM images to mgz files. Does anybody know a good tool to convert the mgz files back to DICOM or, even better, to NiFTI? Any suggestions would be really appreciated! Cheers Diederick

3 2

Fwd: [Freesurfer] Re: mri_mask
by Daniel H Choi 12 May '08

12 May '08

Hi Bruce, The following is the output that my colleague was able to render after running the mri_mask again. Also attached are the brainmask.mgz and lh.stg.mask.mgz files too. Would you be able to point out what went wrong? Thanks. Dan [root@slr08 masks]# mri_info lh.stg.mask.mgz Volume information for lh.stg.mask.mgz type: MGH dimensions: 256 x 256 x 256 voxel sizes: 1.0000, 1.0000, 1.0000 type: INT (1) fov: 256.000 dof: 0 xstart: -128.0, xend: 128.0 ystart: -128.0, yend: 128.0 zstart: -128.0, zend: 128.0 TR: 2300.00 msec, TE: 3.36 msec, TI: 900.00 msec, flip angle: 9.00 degrees nframes: 1 ras xform present xform info: x_r = -1.0000, y_r = -0.0000, z_r = -0.0000, c_r = -0.4948 : x_a = -0.0000, y_a = -0.0000, z_a = 1.0000, c_a = 47.1750 : x_s = -0.0000, y_s = -1.0000, z_s = -0.0000, c_s = 19.0414 talairach xfm : /mnt/win/MRI_data/Trial_Runs/Freesurfer/E-08/mri/transforms/talairach.xfm Orientation : LIA Primary Slice Direction: coronal voxel to ras transform: -1.0000 -0.0000 -0.0000 127.5052 -0.0000 -0.0000 1.0000 -80.8250 -0.0000 -1.0000 -0.0000 147.0414 0.0000 0.0000 0.0000 1.0000 voxel-to-ras determinant -1 ras to voxel transform: -1.0000 -0.0000 0.0000 127.5052 0.0000 0.0000 -1.0000 147.0414 -0.0000 1.0000 0.0000 80.8250 0.0000 0.0000 0.0000 1.0000 Anil. On Sat, May 10, 2008 at 1:50 PM, Daniel H Choi <dhchoi(a)u.northwestern.edu> wrote: > Hey Anil, > > This is what Bruce sent me just in case you needed verbatim what he said. > Thanks. And i'll forward him your results as soon as I get them from you. > Laters > > Dan > > ---------- Forwarded message ---------- > From: Bruce Fischl <fischl(a)nmr.mgh.harvard.edu> > Date: Fri, May 9, 2008 at 8:21 PM > Subject: Re: [Freesurfer] Re: mri_mask > To: Daniel H Choi <dhchoi(a)u.northwestern.edu> > Cc: Krish Subramaniam <krish(a)nmr.mgh.harvard.edu> > > > Hi Daniel, > > I just tested this on some data and it seemed to work for me. Can you run > mri_info on the lh.STG.mgz volume and send us the results? And maybe you can > email me and Krish (ccd) the 2 volumes (brainmask.mgz and lh.STG.mgz) so we > can see if we can replicate it. > > cheers, > Bruce > > On Fri, 9 May 2008, Daniel H Choi wrote: > > Please note the command we used: >> >> mri_mask brainmask.mgz lh.stg.mask.mgz brainout.mgz >> >> and the result we got was: >> >> Writing masked volume to brainout.mgz...done. >> >> On Fri, May 9, 2008 at 2:39 PM, Daniel H Choi <dhchoi(a)u.northwestern.edu> >> wrote: >> >> Hi, >>> >>> I have tried applying a binary mask, specifically for STG, in conjunction >>> with the mri_mask command. We applied it to a brainmask volume, but we >>> did >>> not get the result we wanted. The output we recieved was simply the mask >>> that we used. We wish to produce an output that gives the STG volume >>> alone >>> from the brain, and not from the binary mask. Are we using the command >>> incorrectly? Thanks >>> >>> -- >>> Daniel Horim Choi >>> Northwestern University >>> Class of 2009 >>> dhchoi(a)u.northwestern.edu or >>> dchoi1024(a)gmail.com >>> >> >> >> >> >> >> > > > -- > Daniel Horim Choi > Northwestern University > Class of 2009 > dhchoi(a)u.northwestern.edu or > dchoi1024(a)gmail.com > -- Daniel Horim Choi Northwestern University Class of 2009 dhchoi(a)u.northwestern.edu or dchoi1024(a)gmail.com

4 4

mri_mask
by Daniel H Choi 09 May '08

09 May '08

Hi, I have tried applying a binary mask, specifically for STG, in conjunction with the mri_mask command. We applied it to a brainmask volume, but we did not get the result we wanted. The output we recieved was simply the mask that we used. We wish to produce an output that gives the STG volume alone from the brain, and not from the binary mask. Are we using the command incorrectly? Thanks -- Daniel Horim Choi Northwestern University Class of 2009 dhchoi(a)u.northwestern.edu or dchoi1024(a)gmail.com

1 1

Optseq2 and 2-way design?
by Elizabeth Reynolds Losin 09 May '08

09 May '08

Hello all, I am new to Optseq2 and I have a question about generating 2-way design orders using the program and about a question about temporal jitter. I am working on a design with 7 conditions. I have run the following design through optseq2 and gotten out reasonable looking orders. --ntp 160 --tr 3 --psdwin 0 21 --ev imi_touch 3 20 --ev imi_separate 3 20 --ev obs_touch 3 20 --ev obs_separate 3 20 --ev exe_touch 3 20 --ev exe_separate 3 20 --ev ident 3 20 - -o pilot_order_evrltd --nsearch 10000 --nkeep 3 1) The stimuli are hand signs presented for 3 seconds each and there are 20 repetitions for each of the conditions. The conditions that have seperate or touch in their labels involve presentation of the same set of hand signs in each condition; the different conditions just involve doing different things in response to seeing these hand signs, e.g. imitating or observing them, etc. For example, ev imi_touch, obs_touch and ev exe_touch involve seeing the same set of 20 signs and either observing, imitating or executing them. So, while I'm not interested in analyzing responses to the 20 different hands signs separately, I also don't want the same sign to appear in different conditions successively or even close in time, i.e. imitate and then observe the same sign close in time. I'm wondering if there is a way to tell optseq2 all of this information so that it will generate an order specifying not only my conditions but also randomize the specific repetition of the condition. Basically I want all of the considerations such as psdwin taken into account for condition order and only randomization taken into account for the order of trial repetitions. Somehow I need to let optseq know that the same repetition numbers in different conditions are linked. Is there any way to do this with Optseq2? If not, can you suggest some way I can achieve this, perhaps using optseq2 for condition order and then another program for trial order? Thanks! 2) The second question I have is regarding ISI jitter. I have heard that it is optimal in fast event-related designs that the stimulus duration NOT be a multiple of the TR, yet it seems that this is what optseq is looking for. It doesn't seem that optseq is generating jitter between successive stimuli but only between successive presentations of the same condition. Is this all that is necessary for deconvolution of the signal during post processing, or is a jittered ISI between successive stimuli also necessary (or desirable), and if so, how can this be achieved? One strategy I have heard of is moving the stimulus onset time around within a larger TR. For example shifting a 1.5 sec stimulus around within a 3 sec TR to create temporal jitter. Is this type of timing something optseq could generate? Thanks. Liz

1 0

unpacksdcmdir issue
by jake＠nmr.mgh.harvard.edu 09 May '08

09 May '08

Hello, I am getting an error saying that the path from the findsession command is not found. Here is the bugr information and following is the log from my terminal window; could anyone explain why the directory from findsession is missing? THANKS! FREESURFER_HOME: /usr/local/freesurfer/stable4 Build stamp: freesurfer-Linux-centos4-stable-v4.0.4-20080508 RedHat release: CentOS release 5 (Final) Kernel info: Linux 2.6.18-8.1.15.el5 i686 NMR Center info (/space/freesurfer exists): machine: striatum SUBJECTS_DIR: /space/amaebi/35/users/ablood/R01_Botox PWD: /space/amaebi/35/users/ablood/R01_Botox striatum:jake[81] cd /space/amaebi/35/users/ablood/R01_Botox striatum:jake[82] source /usr/local/freesurfer/nmr-std-env -------- freesurfer-Linux-centos4-stable-v4.0.4-20080508 -------- Setting up environment for FreeSurfer/FS-FAST (and FSL) FREESURFER_HOME /usr/local/freesurfer/stable4 FSFAST_HOME /usr/local/freesurfer/stable4/fsfast FSF_OUTPUT_FORMAT nii SUBJECTS_DIR /space/sake/3/users/inverse/subjects MNI_DIR /usr/local/freesurfer/stable4/mni FSL_DIR /usr/pubsw/packages/fsl/current [striatum:R01_Botox] (nmr-std-env) setenv SUBJECTS_DIR /space/amaebi/35/users/ablood/R01_Botox [striatum:R01_Botox] (nmr-std-env) findsession CD_RO1BTX_pat1_sess1 ======= SUBJECT: CD_RO1BTX_pat1_sess1 SUBJ ID: CD_RO1BTX_pat1_sess1 DATE : May 03, 2008 TIME : 19:27:12 XPRMNTR: PATH : /space/archive/160/siemens/TrioTim-35162-20080503-192630-556000 [striatum:R01_Botox] (nmr-std-env) unpacksdcmdir -src /space/archive/160/siemens/TrioTim-35162-20080503-192630-556000 -targ /space/amaebi/35/users/ablood/R01_Botox/CD_R01BTX_pat1_sess1 -scanonly CD_RO1BTX_pat1_sess1_info $Id: unpacksdcmdir,v 1.19.2.2 2008/03/11 19:56:38 nicks Exp $ /autofs/space/amaebi_035/users/ablood/R01_Botox -src /space/archive/160/siemens/TrioTim-35162-20080503-192630-556000 -targ /space/amaebi/35/users/ablood/R01_Botox/CD_R01BTX_pat1_sess1 -scanonly CD_RO1BTX_pat1_sess1_info Fri May 9 10:51:51 EDT 2008 mri_convert -all-info ProgramName: mri_convert ProgramArguments: -all-info ProgramVersion: $Name: stable4 $ TimeStamp: 2008/05/09-14:51:51-GMT BuildTimeStamp: May 8 2008 03:44:12 CVS: $Id: mri_convert.c,v 1.146.2.2 2008/03/02 02:53:20 nicks Exp $ User: jake Machine: striatum Platform: Linux PlatformVersion: 2.6.18-8.1.15.el5 CompilerName: GCC CompilerVersion: 30400 striatum Linux striatum 2.6.18-8.1.15.el5 #1 SMP Mon Oct 22 08:32:04 EDT 2007 i686 i686 i386 GNU/Linux ERROR: directory /space/archive/160/siemens/TrioTim-35162-20080503-192630-556000 does not exist [striatum:R01_Botox] (nmr-std-env) cd /space/archive/160/siemens/ /space/archive/160/siemens/: No such file or directory. [striatum:R01_Botox] (nmr-std-env)

2 1

Fixing Inaccuracies in White Matter Surfaces
by Pradeep Reddy Ramana 09 May '08

09 May '08

Hello Everybody, This is Pradeep Reddy, just started my PhD in Medical Image Analysis and hence a newbie to FreeSurfer. I've passed the auto-recon2 stage in the reconstruction procedure. Now I am trying to verify WM segmentation produced by the recon-all/auto-recon2, before I move on to further processing. The excercise given in the freesurfer tutorial ( http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/WhiteMatterEdits ) gives me one example on how to detect if there were a lesion and outlines how to fix it. But now when I observe the outputs I have, I am quite confused as to how to detect inaccuracies in the segmentation of white matter. To make it more specific, I show you a slice of wm.mgh I have ( wmConfusion1 attached). Please look at the small yellow circle on the right. If I just follow the exercise, this looks like a lesion, as the yellow line cuts into it. but when I move up to next few slices, this hole opens up ( close to the green line ) giving me indications that it may NOT be a hole ( keeping in mind the inherent 3D nature of the MR image ). Also I am not sure whether there is an error in the image ( wmConfusion5 ) attched is an error is WM segmentation. So my questions: 1. What are the things I have to look for, while trying to verify the WM segmentation ( after autorecon2 )? 2. Any thumb-rules laid out for this already? My googling was in fact in vain. 3. Any other tutorials/resources/books stating clearly what to verify, in different stages of reconstruction? I am sorry to write a long email. Please be kind to help me in this regard. Regards, Pradeep Reddy PhD Student, Simon Fraser Unviersity, Canada. -- There's more to life. Finding it is my job. -- raamana(a)gmail.com NOTICE my new contact Number in Canada: +1 604 720 2233

2 1

LGI of annotated cortex
by Martin Kavec 09 May '08

09 May '08

Hi, is it possible to obtain LGI values of cortices, similarly as curvature indices in ?h.aparc.stats? Could this value possibly be included in the ?h.aparc.stats files? Thanks in advance, Martin -- ********************************** Senior Clinical Research Associate MRI Unit of the Department of Radiology Erasme Hospital Lennik Street 808 B-1070 Brussels BELGIUM tel: +32-2-555-4325 fax: +32-2-555-3994 email: martin.kavec(a)gmail.com ********************************** ----------------------------------------------- Find a way, or make one!

4 10

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