Sounds good, though it's probably better to use example_func to register to the anatomical. If you've done registration inside FEAT you should find all the transforms you need in the .feat/reg directory. Cheers, Steve. On 17 May 2006, at 20:17, Martin M. Monti wrote: > Steve, > > my simple plan was the following: > > 1) Align filtered_func to the subject's anatomical via flirt > > 2) Register the mask (created in Standard space) to the individual > subject by using the inverse of highres2standard.mat , and use the > example_filtered_func as the reference. This way the mask is > registered to the functional (i.e. they both are registered to the > subject's anatomical) and the header of the mask is the same as the > functional. (at this point I'll use avwmeants to pull out > timecourses from that mask) > > [alternatively I could probably just apply to the 'standard' mask > the inverse of the example_func2standard.mat, but I have several > functionals and I'd like to concatenate them to then use once > avwmeants] > > Makes good sense? > > Cheers & thanks for the continuous support, > > Martin > > ----- Original Message ----- From: "Steve Smith" > <[log in to unmask]> > To: <[log in to unmask]> > Sent: Thursday, May 18, 2006 2:55 AM > Subject: Re: [FSL] filtered_func & Peri-Stimulus plot > > >> Hi, >> >> On 17 May 2006, at 18:13, Martin M. Monti wrote: >> >>> Mark, >>> >>> thank you, avwmeants works great! >>> Only one question: does the mask used have to be in the FOV of >>> the functional I'm pulling out the timecourses from? >>> >>> I.e. say I do: avwmeants -i filtered_func_data -o tseries.txt -m >>> Mask --showall >>> >>> Does 'Mask' have to have same FOV and voxel resolution as >>> filtered_func_data? >> >> Yes. >> >>> I presume I can use flirt on the mask with filtered_func_data as >>> a reference image to have them have these same header properties? >> >> Indeed, as long as you have an appropriate FLIRT transform file >> to tell it to apply to the mask. >> >> Cheers, Steve. >> >> >>> >>> Thank for the help, >>> >>> Martin >>> >>> ----- Original Message ----- From: "Mark Jenkinson" >>> <[log in to unmask]> >>> To: <[log in to unmask]> >>> Sent: Tuesday, May 16, 2006 6:36 AM >>> Subject: Re: [FSL] filtered_func & Peri-Stimulus plot >>> >>> >>>> Hi Martin, >>>> >>>> You can get individual time courses or averages using avwmeants >>>> (the output is in text form). This will let you process the >>>> values later >>>> in whatever way you want. >>>> >>>> You are correct about the filtered_func_data - it does not have >>>> registration applied to it (but is motion corrected). >>>> >>>> I'm not quite sure what is wrong with your Featquery, but your >>>> usage of avwroi is wrong, as you need to use voxel coordinates, >>>> not mm coordinates. All the voxel coordinates should be integers >>>> between 0 and N-1 (inclusive). Using -34 as a coordinate would be >>>> well outside the volume and hence not give you anything except >>>> zeros in the resulting output. >>>> >>>> All the best, >>>> Mark >>>> >>>> >>>> On 15 May 2006, at 19:42, Martin M. Monti wrote: >>>> >>>>> FSL masters, >>>>> >>>>> I have 2 ( somewhat unrelated questions) >>>>> >>>>> 1) the filtered_func_data: this is the data with all pre- >>>>> processing steps I asked for, EXCEPT for registration with the >>>>> High-res and MNI-Template? I'm trying a multi-voxel analysis >>>>> of my data and need to pull out time-courses for various >>>>> regions. Also, is there a way to pull out the time-courses >>>>> associated to EACH voxel within a mask? I know avwroi will >>>>> pull out one average from a given mask. Is there any pre- >>>>> compiled utility to pull out each and all the voxels from a >>>>> mask? >>>>> >>>>> 2) I tried featquery to get a peristimulus plot, but I get >>>>> more event plots than I actually have (i.e. 8 instead of 5!). >>>>> This is by selecting the filtered_func_data and clicking on >>>>> the PE/% signal change and peri-stimulus plot options (I use a >>>>> mask in MNI coordinated and use single-voxel coords). Indeed I >>>>> think I remember someone mentioning that FSL is left to >>>>> "guess" what an event is. Is there a way to instruct FSL as of >>>>> the timepoints that make an event? >>>>> Alternatively I was wondering whether I could just extract the >>>>> raw timecourse from the filtered_func_data using avwroi and >>>>> then manually create the plots (excel-style). Re this, one >>>>> thing I'm unclear about is, when the syntax requires size: is >>>>> that number of voxels or voxel size? I tried using 1 for all >>>>> dimensions (e.g. avwroi filtered_func_data TC_out -34 1 60 1 >>>>> 4 1) but i then couldn't open the image. By changing the 1s >>>>> to the voxel dimension I could successfully create a new >>>>> image. Is this correct? >>>>> >>>>> cheers >>>>> >>>>> Martin >>>>> >> >> >> --------------------------------------------------------------------- >> --- --- >> Stephen M. Smith, Professor of Biomedical Engineering >> Associate Director, Oxford University FMRIB Centre >> >> FMRIB, JR Hospital, Headington, Oxford OX3 9DU, UK >> +44 (0) 1865 222726 (fax 222717) >> [log in to unmask] http://www.fmrib.ox.ac.uk/~steve >> --------------------------------------------------------------------- >> --- --- ------------------------------------------------------------------------ --- Stephen M. Smith, Professor of Biomedical Engineering Associate Director, Oxford University FMRIB Centre FMRIB, JR Hospital, Headington, Oxford OX3 9DU, UK +44 (0) 1865 222726 (fax 222717) [log in to unmask] http://www.fmrib.ox.ac.uk/~steve ------------------------------------------------------------------------ ---