Freesurfer May 2013

freesurfer@nmr.mgh.harvard.edu

174 participants
299 discussions

nuintensitycor and skullstrip problems
by LAOUCHEDI MAKHLOUF 28 May '13

28 May '13

hi i send you the rawavg.mgz of the subject with the massive artefact. i inform you that i proceeded differently by splitting the recon-all -all into three steps (as in freesurfer wiki) : (1) recon-all -autorecon1, (2) recon-all -autorecon2 (3) recon-all -autorecon3. i just finished with the first step which still didn't work for some subjects so i deleted then in my SUBJECTS_DIR then i run "recon-all -i" to regenerate them. after that i re-run "recon-all -autorecon1" and it worked even for the subject with the artefacted rawavg.mgz i am sending you. i noticed also that some extra cerebral structures like dura is not fully removed after the skullstripping step, and i reviewed some of my previous datasets i processed with freesurfer and i noticed that the cortical ribbon (containing the cortical parcellation) fits well the gray matter. so is a wise idea to not bother a lot about these "some" extra cerebral structures left and just go forward with the rest of the recon-all pipeline (but of course i'll take a look a posteriori). Thanks

2 1

minor bug report
by Joseph Dien 28 May '13

28 May '13

when I run: isxconcat-sess -sf sessidlist.dat -analysis CPA.sm05.lh -all-contrasts -o secondLevel I get the error messages: secondLevel/CPA.sm05.mni305/log/isxconcat.log stat: illegal option -- c usage: stat [-FlLnqrsx] [-f format] [-t timefmt] [file ...] stat: illegal option -- c usage: stat [-FlLnqrsx] [-f format] [-t timefmt] [file ...] stat: illegal option -- c usage: stat [-FlLnqrsx] [-f format] [-t timefmt] [file …] ... I'm assuming this isn't a serious problem as it just seems to affect what goes into the log file but thought you might want to know. This is on OS X 10.8.2 with Freesurfer 5.2.0. Cheers! Joe -------------------------------------------------------------------------------- Joseph Dien, Senior Research Scientist University of Maryland E-mail: jdien07(a)mac.com Phone: 301-226-8848 Fax: 301-226-8811 http://joedien.com//

3 3

fix talairach transform
by LAOUCHEDI MAKHLOUF 28 May '13

28 May '13

Hi i have a problem when running "recon-all -all", so i decided to split the process into three staps : "recon-all -autorecon1", ""recon-all -autorecon2" and "recon-all -autorecon3". i have just finished with the first step (recon-all -autorecon1) and looked for misalignment in talairach transform for all subjects using "tkregister2 --mgz --s ctr1 --fstal" (i can't use "--surf orig " because i haven't the surfaces yet). so i used the small red cross in the coronal, sagittal and axial views to align as much as possible. my goal is to extract the cortical labels (cortical parcellation). is it necessary to have a precise talairach transform just for this purpose ? Thank you

2 1

nuintensitycor and skullstrip problems
by LAOUCHEDI MAKHLOUF 28 May '13

28 May '13

hi i send you further information of the -nuintensitycor and skullstrip problem,about three control subjects. after running "recon_all -i", i decided to process my data step by step, i run "recon_all -autorecon1" for all of them and it ended with errors. i picked up three cases of the recon-all.log file : subject 1: #@# Nu Intensity Correction Mon May 27 18:08:32 CEST 2013 mri_nu_correct.mni --i orig.mgz --o nu.mgz --uchar transforms/talairach.xfm --n 2 /home_cluster/mlaouche/freesurf_makh/ctr3/mri /home_cluster/amesse/software/freesurfer5.1/bin/mri_nu_correct.mni --i orig.mgz --o nu.mgz --uchar transforms/talairach.xfm --n 2 nIters 2 $Id: mri_nu_correct.mni,v 1.18 2011/03/02 20:16:39 nicks Exp $ Linux cn09 2.6.32-358.2.1.el6.x86_64 #1 SMP Wed Mar 13 00:26:49 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux Mon May 27 18:08:36 CEST 2013 Program nu_correct, built from: Package MNI N3, version 1.10, compiled by nicks@minerva (x86_64-unknown-linux-gnu) on 2010-02-20 at 17:32:37 tmpdir is ./tmp.mri_nu_correct.mni.12456 /home_cluster/mlaouche/freesurf_makh/ctr3/mri mri_convert orig.mgz ./tmp.mri_nu_correct.mni.12456/nu0.mnc -odt float mri_convert orig.mgz ./tmp.mri_nu_correct.mni.12456/nu0.mnc -odt float $Id: mri_convert.c,v 1.179.2.2 2011/05/16 20:53:47 greve Exp $ reading from orig.mgz... TR=7.05, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) changing data type from uchar to float (noscale = 0)... writing to ./tmp.mri_nu_correct.mni.12456/nu0.mnc... -------------------------------------------------------- Iteration 1 Mon May 27 18:23:46 CEST 2013 nu_correct -clobber ./tmp.mri_nu_correct.mni.12456/nu0.mnc ./tmp.mri_nu_correct.mni.12456/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.12456/0/ ./tmp.mri_nu_correct.mni.12456/ is a directory but not writeable [mlaouche@cn09:/home_cluster/mlaouche/freesurf_makh/ctr3/mri/] [2013-05-27 18:23:56] running: /home_cluster/amesse/software/freesurfer5.1/mni/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir ./tmp.mri_nu_correct.mni.12456/0/ ./tmp.mri_nu_correct.mni.12456/nu0.mnc ./tmp.mri_nu_correct.mni.12456/nu1.imp Assertion failed at line 827 in file templates/CachedArray.cc nu_estimate_np_and_em: crashed while running volume_stats (termination status=256) nu_correct: crashed while running nu_estimate_np_and_em (termination status=256) ERROR: nu_correct Linux cn09 2.6.32-358.2.1.el6.x86_64 #1 SMP Wed Mar 13 00:26:49 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux recon-all -s ctr3 exited with ERRORS at Mon May 27 18:47:13 CEST 2013 To report a problem, see http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting subject 2: #@# Nu Intensity Correction Mon May 27 18:08:21 CEST 2013 mri_nu_correct.mni --i orig.mgz --o nu.mgz --uchar transforms/talairach.xfm --n 2 /home_cluster/mlaouche/freesurf_makh/ctr5/mri /home_cluster/amesse/software/freesurfer5.1/bin/mri_nu_correct.mni --i orig.mgz --o nu.mgz --uchar transforms/talairach.xfm --n 2 nIters 2 $Id: mri_nu_correct.mni,v 1.18 2011/03/02 20:16:39 nicks Exp $ Linux cn09 2.6.32-358.2.1.el6.x86_64 #1 SMP Wed Mar 13 00:26:49 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux Mon May 27 18:08:31 CEST 2013 Program nu_correct, built from: Package MNI N3, version 1.10, compiled by nicks@minerva (x86_64-unknown-linux-gnu) on 2010-02-20 at 17:32:37 tmpdir is ./tmp.mri_nu_correct.mni.12252 /home_cluster/mlaouche/freesurf_makh/ctr5/mri mri_convert orig.mgz ./tmp.mri_nu_correct.mni.12252/nu0.mnc -odt float mri_convert orig.mgz ./tmp.mri_nu_correct.mni.12252/nu0.mnc -odt float $Id: mri_convert.c,v 1.179.2.2 2011/05/16 20:53:47 greve Exp $ reading from orig.mgz... TR=7.07, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) changing data type from uchar to float (noscale = 0)... writing to ./tmp.mri_nu_correct.mni.12252/nu0.mnc... tee: ./mri_nu_correct.mni.log: Input/output error -------------------------------------------------------- Iteration 1 Mon May 27 18:23:44 CEST 2013 nu_correct -clobber ./tmp.mri_nu_correct.mni.12252/nu0.mnc ./tmp.mri_nu_correct.mni.12252/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.12252/0/ ./tmp.mri_nu_correct.mni.12252/ is a directory but not writeable [mlaouche@cn09:/home_cluster/mlaouche/freesurf_makh/ctr5/mri/] [2013-05-27 18:23:56] running: /home_cluster/amesse/software/freesurfer5.1/mni/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir ./tmp.mri_nu_correct.mni.12252/0/ ./tmp.mri_nu_correct.mni.12252/nu0.mnc ./tmp.mri_nu_correct.mni.12252/nu1.imp Computing ranges for each class.: ............................................. Processing:.................................................................Done Fitting splines: .............................................................. Computing ranges for each class.: ............................................. Processing:.................................................................Done Computing ranges for each class.: ............................................. Processing:.................................................................Done Processing:.................................................................Done Processing:.................................................................Done Processing:.................................................................Done Processing:.................................................................Done Linux cn09 2.6.32-358.2.1.el6.x86_64 #1 SMP Wed Mar 13 00:26:49 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux recon-all -s ctr5 exited with ERRORS at Mon May 27 19:13:37 CEST 2013 To report a problem, see http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting subject 3: #@# Nu Intensity Correction Mon May 27 18:08:34 CEST 2013 mri_nu_correct.mni --i orig.mgz --o nu.mgz --uchar transforms/talairach.xfm --n 2 /home_cluster/mlaouche/freesurf_makh/ctr6/mri /home_cluster/amesse/software/freesurfer5.1/bin/mri_nu_correct.mni --i orig.mgz --o nu.mgz --uchar transforms/talairach.xfm --n 2 nIters 2 $Id: mri_nu_correct.mni,v 1.18 2011/03/02 20:16:39 nicks Exp $ Linux cn09 2.6.32-358.2.1.el6.x86_64 #1 SMP Wed Mar 13 00:26:49 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux Mon May 27 18:08:38 CEST 2013 Program nu_correct, built from: Package MNI N3, version 1.10, compiled by nicks@minerva (x86_64-unknown-linux-gnu) on 2010-02-20 at 17:32:37 mkdir: cannot create directory `./tmp.mri_nu_correct.mni.12552': Input/output error tee: ./mri_nu_correct.mni.log: Input/output error tmpdir is ./tmp.mri_nu_correct.mni.12552 Linux cn09 2.6.32-358.2.1.el6.x86_64 #1 SMP Wed Mar 13 00:26:49 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux

2 1

Re: [Freesurfer] rerun interrupted processes
by Bruce Fischl 28 May '13

28 May '13

Hi Valeria can you address this type of question to the list so that others can answer and learn from the answers? Yes, you can run recon-all -make all -s <subject id> and it should pick up from where it left off Bruce On Tue, 28 May 2013, Valeria Contarino wrote: > Dear Bruce, > I have a new question for you. > I was running 30 subjects when processes were interrupted because of no more space left on device. > Do I have the chance for starting again the processes from where they stopped? > > Thank you. > > Regards, > Valeria > >

1 0

recon-all with FLAIR error no FSL installed
by Popescu, Veronica 28 May '13

28 May '13

Dear all, I am running FreeSurfer 5.3 with FLAIR refinement. I do not have FSL installed. I do get results but there's an error message in the log file (attached) about the FSL installation missing. How bad is it? Do I have to rerun everything with FSL installed? And if so, is there any version that is preferable? Thank you very much in advance! Best wishes, Veronica

3 6

Re: [Freesurfer] possible to use freesurfer segmentations to get mean values for anatomic structures from other sequences
by Salil Soman 27 May '13

27 May '13

Dear Doug, Thank you for your email. Following up on our previous discussions: 1) *bbregister fail *- this was happening because of error from fslregister. The variable fsvol was translating to "ERROR". In fslregister th value for fsvol is set to brainmask, along with numerous other parameters being initialized. By adding the line "set fsvol = brain mask" (without quotes) at the beginning of the section "#Convert reference to analyze" the CentOS 6 version I am running on ubuntu 12.04 (with SGE), now works. 2) when you suggested I run bbregister on the lowb.nii file, was it because I had suggested I wanted to register to a 2D ADC map? If I want to instead use a 3D volume (e.g. a 3D ASL), would I use the actual volume I want to get value from instead of lowb.nii? I ended up running 3) for the mri_vol2vol step - I had to change targ to fstarg, and I left out the interp flag so it would do its default trilinear: mri_vol2vol --mov Vol_906.nii --fstarg T1.mgz --o Vol_906_Reg_Bold.nii --no-save-reg --reg register_bold.dat --inv 4) I would like to get the segmentations from aparc+aseg.mgz as the areas where I get the values summed. So I entered mri_segstats --seg mri/aparc+aseg.mgz --ctab $FREESURFER_HOME/FreeSurferColorLUT.txt --i Vol_906_Reg_Bold.nii --sum Vol_aparc_aseg.stats but then I received this error: ERROR: dimension mismatch between input volume and seg input 224 224 170 seg 256 256 256 * * *So, in summary, while I've fixed bbregister to work on my system, I am still having trouble normalizing my 3D non T1 volume data into the FS space to be able to use the labels from aparc+aseg.mgz to get statistics for each anatomic area.* Any ideas on how to proceed? Thanks, Sal On Sat, May 25, 2013 at 11:51 AM, Douglas Greve <greve(a)nmr.mgh.harvard.edu>wrote: > > On 5/25/13 1:39 PM, Salil Soman wrote: > > Thank you Doug. > > Just so I am clear - is the anatomical the nifti T1.mgz or is it a > different file. From your email I gather should do the following: > > *1) Use bbregister to register the lowb image to the anatomical. This > creates a registration matrix.* > > bbregister -s SUBJECTNAME --mov lowb.nii --initfsl --reg register.dat > > That is right, but add --t2 (since it is t2 weighted). > > > *2) Use mri_vol2vol and the registration to map the ADC map to the > anatomical* > > mri_vol2vol --mov lowb.nii --targ $vol --inv --interp nearest --o $vol2diff --reg $outdir/register.dat --no-save-reg > > Use the ADC as the moveable (not lowb, but use the lowb for > bbregister). The target volume should be the T1.mgz (or any conformed > volume). The output will be the adc in the anatomical space (not sure why > you call it vol2diff). > > *3) use mri_segstats to extract the values* > > *mri_segstats* --seg $SUBJECTS_DIR/SUBJECTNAME/mri/wmparc.mgz --ctab $FREESURFER_HOME/FreeSurferColorLUT.txt --i lowb.nii --sum fa.stats > > The input would be the adc volume mapped the anatomical space. I would > use WMParcStatsLUT.txt or else you'll get a billion different areas not > represented in wmparc. You may also want to add "--seg-erode 1" to erode > the segmentations by a voxel. This helps to prevent partial voluming. > > doug > > > > Where lowb.nii is the other MRI tissue sequence I am analyzing (e.g. > ADC), $vol2diff is the output image of the registration, and fa.stats will > be my output stats table? > > I suspect there is part of the syntax for these tools I do not > understand. Also, what input would $vol be? > > Best wishes, > > Sal > > > On Sat, May 25, 2013 at 10:06 AM, Douglas Greve <greve(a)nmr.mgh.harvard.edu > > wrote: > >> >> Hi Sal, yes you can. Use bbregister to register the lowb image to the >> anatomical. This creates a registration matrix. Use mri_vol2vol and the >> registration to map the ADC map to the anatomical, then use mri_segstats to >> extract the values >> doug >> >> ps. Please post questions to the FS list instead of us directly so that >> others can benefit and it gets archived. thanks! >> >> >> >> On 5/25/13 1:03 PM, Salil Soman wrote: >> >>> Dear Dr. Greve, >>> >>> Is it possible to register a nonstructural MR sequences (e.g. an ADC >>> map) with the aseg.mgz file (or aparc+aseg.mgz file, etc) and then, using >>> the segmentation from the aseg (or aparc+aseg) file to output a mean ADC >>> value for each anatomic area segmented? >>> >>> Thank you for your time and consideration. >>> >>> Best wishes, >>> >>> Sal >>> >>> Salil Soman, MD, MS >>> >> >> >> >> 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 6

white matter volume per lobe...
by Gonzalo Rojas Costa 27 May '13

27 May '13

Hi: I used the following commands to get the cortical volume per each lobe: mri_annotation2label --subject OAS1_0001_MR1 --hemi rh --lobesStrict lobulos --seg lobulos --stat lobulos mri_annotation2label --subject OAS1_0001_MR1 --hemi lh --lobesStrict lobulos --seg lobulos --stat lobulos mris_anatomical_stats -a lobulos -f OAS1_0001_MR1_lobulos_rh -b OAS1_0001_MR1 rh mris_anatomical_stats -a lobulos -f OAS1_0001_MR1_lobulos_lh -b OAS1_0001_MR1 lh Which commands (sintaxis) must I use to compute the white matter volume per each lobe ?... I tested the mri_aparc2aseg and mri_segstats commans with different options, but I couldn't get the volume that I need... Sincerely, -- Gonzalo Rojas Costa Laboratory for Advanced Medical Image Processing Department of Radiology Clínica las Condes Lo Fontecilla 441, Las Condes, Santiago, Chile. Tel: 56-2-2105170 Cel: 56-9-97771785 www.clc.cl

2 1

mri_glmfit files
by Joshua deSouza 27 May '13

27 May '13

Hello all, I would just like to know what is the significance of the files that are produced after a glmfit is done. I was trying to find a write up online but I was not able to. They include : - cnr.mgh - F.mgh - gamma.mgh - gammavar.mgh - sig.mgh Thank you

2 1

Re: [Freesurfer] FMRI surface transformation for vertex-wise connectivity analysis
by Bruce Fischl 27 May '13

27 May '13

you would change fsaverage5 to fsaverage4 (fewer vertices) or fsaverage6 (more vertices) On Thu, 16 May 2013, Ludovico Minati wrote: > Hi Thomas > > many thanks. So, in the example you have sent me, how should I go > about if I wished to vary the granularity? Are there different meshes > to reshape or should I generate them somehow? > > Ludovico > > > On 16 May 2013 14:47, Thomas Yeo <ythomas(a)csail.mit.edu> wrote: >> Hi Ludovico, >> >> The mri_vol2surf example I sent you directly resample to fsaverage5, >> which is around 4mm spacing, resulting in 18715 vertices for both >> hemispheres combined after excluding the medial wall in the surfaces. >> >> I am also including the freesurfer mailing list in the email as >> suggested by Bruce. >> >> --Thomas >> >> On Thu, May 16, 2013 at 8:25 PM, Bruce Fischl >> <fischl(a)nmr.mgh.harvard.edu> wrote: >>> Hi Ludovico >>> >>> I wouldn't downsample before calling mri_surf2surf. We have tools for >>> uniformly parcellating the cortex (mris_make_face_parcellation) that can be >>> used for this purpose after sampling the timecourses to the cortex. >>> >>> >>> cheers >>> Bruce >>> On Thu, 16 May 2013, Ludovico Minati wrote: >>> >>>> Hi Bruce, >>>> >>>> many thanks. Yes sorry I shall use the list from now on. But let me >>>> just ask you a clarification as we are already on these questions >>>> >>>> For our connectivity analyses, at present we parcellated in 3D volume >>>> space the cortex into a pre-set number of roughly equal volume >>>> parcels, ranging in number between around 500 to 10,000. I would need >>>> to translate this approach to the surface space in Freesurfer, i.e. be >>>> able to smooth on the cortical surface and undersample it. What is the >>>> best approach that can be used? >>>> >>>> Along this like: you wrote around 150k vertices/hemisphere. For our >>>> purposes that's an enormous number, much larger than the number of EPI >>>> voxels. How can we mesh downsample before calling mri_vol2surf? >>>> >>>> Again thanks for your kind help and patience >>>> >>>> Ludovico >>>> >>>> >>>> On 16 May 2013 14:09, Bruce Fischl <fischl(a)nmr.mgh.harvard.edu> wrote: >>>>> >>>>> Hi Ludovico >>>>> >>>>> can we move this discussion to the FreeSurfer list in the future so >>>>> others >>>>> can answer? >>>>> >>>>> 1. Hard to say as sometimes one works better than the other, you'll just >>>>> have to try. The transform should be rigid (6DOF) since the scalings are >>>>> read from the header, why would you want affine? I think Doug has worked >>>>> on >>>>> this and you can specify it in the bbregister command line, but I've >>>>> never >>>>> used it. >>>>> >>>>> 2. About 150K vertices/hemisphere. They are distributed initially at 1mm >>>>> spacing (the surface of the voxels), and cannot be tuned. There are some >>>>> tools around for mesh downsampling if you want, and this is even easier >>>>> if >>>>> you go to fsaverage space (with mri_surf2surf) where we have >>>>> representations >>>>> based on recursive icosahedral supersamplings. >>>>> >>>>> 3. This is due to a bug in the initial nifti format that didn't allow >>>>> any >>>>> one dimension greater than 64K. If you output in .mgz it should be >>>>> nvertices >>>>> x 1 x 1 x ntimepoints. >>>>> >>>>> 4. Search for "coordinate systems" on our wiki. This is documented >>>>> extensively. >>>>> >>>>> there are also matlab i/o functions for reading/writing all our formats >>>>> (e.g. MRIread.m) >>>>> >>>>> cheers >>>>> Bruce >>>>> >>>>> >>>>> On Thu, 16 May 2013, Ludovico Minati wrote: >>>>> >>>>>> Hi Thomas, >>>>>> >>>>>> thanks. A few quick questions, which hopefully will take little of >>>>>> your time but save me a lot of hassle: >>>>>> >>>>>> 1) bbregister. Do you advise using FLIRT or SPM for initialization? Is >>>>>> this a rigid-body transformation, or an affine one? Any way I can >>>>>> control the DOF? If I use SPM shall I feed the mean realigned volume >>>>>> here, as I expect? >>>>>> >>>>>> 2) recon-all. What is the typical number of vertices for the cortex? >>>>>> Any way I can control it, i.e. sample more or less finely? This would >>>>>> be particularly relevant for our task. >>>>>> >>>>>> 3) mri_vol2surf. Perhaps naively I was expecting an output file simply >>>>>> consisting of a list of values, i.e. one per vertex (assuming a single >>>>>> input volume that is). But the documentation says: "The output will be >>>>>> a data set with Nv/R colums, 1 row, R slices, and Nf frames..." What >>>>>> is the reshaping factor R for? How do I deal with this, i.e. how do I >>>>>> simply get a value per vertex? Or am I missing something here? Is >>>>>> mri_vol2surf averaging over neighbouring vertices? >>>>>> >>>>>> 4) Where do I read the Euclinean co-ordinates of the vertices from? >>>>>> >>>>>> Cheers >>>>>> >>>>>> Ludovico >>>>>> >>>>>> On 16 May 2013 04:30, Thomas Yeo <ythomas(a)csail.mit.edu> wrote: >>>>>>> >>>>>>> >>>>>>> Hi Ludovico, >>>>>>> >>>>>>> Just remember that you need to run "recon-all" on your T1 anatomical >>>>>>> first. >>>>>>> >>>>>>> For bbregister, the command goes something like this: bbregister >>>>>>> --bold --s $anat_s --init-fsl --mov $input --reg $reg. I find the >>>>>>> bbregister cost function value a good indication of whether the >>>>>>> registration has good well. Below 0.5 is good. Below 0.6 is probably >>>>>>> acceptable. >>>>>>> >>>>>>> For mri_vol2surf, the command goes something like this: mri_vol2surf >>>>>>> --mov $f --reg $reg --hemi $hemi --projfrac 0.5 --trgsubject >>>>>>> fsaverage5 --o $output --reshape --interp trilinear >>>>>>> >>>>>>> Cheers, >>>>>>> Thomas >>>>>>> >>>>>>> On Wed, May 15, 2013 at 8:57 PM, Bruce Fischl >>>>>>> <fischl(a)nmr.mgh.harvard.edu> wrote: >>>>>>>> >>>>>>>> >>>>>>>> Hi Ludovico >>>>>>>> >>>>>>>> you can get vertex-level timecourse by using bbregister to register >>>>>>>> the >>>>>>>> EPI >>>>>>>> data to the surface, then mri_vol2surf to sample onto the surface. >>>>>>>> >>>>>>>> cheers >>>>>>>> Bruce >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> On Wed, 15 May 2013, Ludovico Minati wrote: >>>>>>>> >>>>>>>>> Hello Thomas, >>>>>>>>> >>>>>>>>> yes sure. Our area of interest is early dementia diagnosis, and we >>>>>>>>> are >>>>>>>>> exploring the predictive value of graph topological indices. As you >>>>>>>>> know >>>>>>>>> there are major volume changes in dementia, and from the literature I >>>>>>>>> have >>>>>>>>> convinced myself we may be better off resolving this confound through >>>>>>>>> working on the surface, rather than playing around with ad-hoc >>>>>>>>> atlases >>>>>>>>> for >>>>>>>>> volumetric normalization. >>>>>>>>> >>>>>>>>> I have got plenty of signal processing code to perform the analyses >>>>>>>>> from >>>>>>>>> BOLD timeseries; what I am looking for now is a way to apply >>>>>>>>> Freesurfer >>>>>>>>> normalization to obtain vertex-level time-courses. >>>>>>>>> >>>>>>>>> Hence, I am so far only interested in solving the following problem >>>>>>>>> with >>>>>>>>> Freesurfer: we have as input an individual T1 scan and a series of >>>>>>>>> EPI >>>>>>>>> volumes, we want to preprocess/coregister/normalize them with >>>>>>>>> Freesurfer and >>>>>>>>> extract the BOLD timecourses for each vertex. No fMRI analysis in >>>>>>>>> Freesurfer, no group analysis. >>>>>>>>> >>>>>>>>> Hopefully it is not too complex? I am dwarfed by the range of >>>>>>>>> commands >>>>>>>>> available! >>>>>>>>> >>>>>>>>> Thank you very much in advance >>>>>>>>> >>>>>>>>> Ludovico >>>>>>>>> >>>>>>>>> -----Original Message----- >>>>>>>>> From: Thomas Yeo <ythomas(a)csail.mit.edu> >>>>>>>>> Sender: yeoyeo02(a)gmail.com >>>>>>>>> Date: Wed, 15 May 2013 10:07:00 >>>>>>>>> To: Bruce Fischl<fischl(a)nmr.mgh.harvard.edu> >>>>>>>>> Cc: lminati(a)ieee.org<lminati(a)ieee.org>; Dennis >>>>>>>>> Chan<D.Chan(a)bsms.ac.uk> >>>>>>>>> Subject: Re: FMRI surface transformation for vertex-wise connectivity >>>>>>>>> analysis >>>>>>>>> >>>>>>>>> Hi Ludovico, >>>>>>>>> >>>>>>>>> Can you give me the context of what you need? >>>>>>>>> >>>>>>>>> Cheers, >>>>>>>>> Thomas >>>>>>>>> >>>>>>>>> On Wed, May 15, 2013 at 9:52 AM, Bruce Fischl >>>>>>>>> <fischl(a)nmr.mgh.harvard.edu> wrote: >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> Hi Ludovico >>>>>>>>>> Thomas Yeo (ccd) has a set of scripts for this purpose >>>>>>>>>> Cheers >>>>>>>>>> Bruce >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> On May 14, 2013, at 8:50 PM, "Ludovico Minati" <lminati(a)ieee.org> >>>>>>>>>> wrote: >>>>>>>>>> >>>>>>>>>>> Ludovico >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> 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. >>>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>> >>>>>> >>>>>> >>>>>> >>>>> >>>> >>>> >>>> >>> > > >

4 3

← Newer
1
...
6
7
8
9
10
11
12
...
30
Older →

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 May 2013