Hi,
If you want to use Two-Thresholding (TT) - and I guess that is what you
described -, it is already implemented in SPM:
- publication:
http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6WNP-4X3PHYG-13F&_user=10&_rdoc=1&_fmt=&_orig=search&_sort=d&_docanchor=&view=c&_searchStrId=1084104426&_rerunOrigin=google&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=7008634d9068c6b3c3763d40ddb536f6
- download: SPM website/ Extensions: "TT Toolbox for SPM" (you can also
use it for SPM8)
Bests,
Auer, Tibor M.D. Ph.D.
Biomedizinische NMR Forschungs GmbH
am Max-Planck Institut für Biophysikalische Chemie
Am Fassberg 11
37077 Göttingen
Germany
Phone/Work: +49-(0)551-201-1725
Phone/Home: +49-(0)551-250-4576
Mobile: +49-(0)176-8012-7921
Mail: [log in to unmask]
> Ok, thank you very much! I think I will already succeed.
>
> On Sun, 8 Nov 2009 13:14:20 +0000, Reza Salimi <[log in to unmask]>
> wrote:
>
>>hi,
>>you can use *cluster* (to extract the clusters and their index numbers)
>>and fslmaths (to threshold the index volume and keep the one you want)
>>cheers
>>
>>
>>2009/11/8 Gabor Perlaki <[log in to unmask]>
>>
>>> This is exactly what I want. Do you know some command which can I use
>>> to
>>> keep count of the 2.3 clusters and split them into separate images?
>>>
>>> On Sat, 7 Nov 2009 08:37:22 +0000, Dave Flitney
>>> <[log in to unmask]>
>>> wrote:
>>>
>>> >Does something like this do what you want?
>>> >
>>> >First split your 2.3 clusters into separate images.
>>> >Multiply each 2.3 cluster image by the combined 3.1 clusters image and
>>> >discard 2.3 cluster image if result empty, i.e., no common voxels.
>>> >Re-combine all surviving 2.3 cluster images.
>>> >
>>> >Dave Flitney, IT Manager
>>> >University of Oxford, FMRIB Centre
>>> >JR Hospital, Oxford OX3 9DU
>>> >T: 01865 222713 F: 01865 222717
>>> >
>>> >On 7 Nov 2009, at 03:48, Gabor Perlaki <[log in to unmask]>
>>> >wrote:
>>> >
>>> >> I think the problem is not solvable with a simple multiplication of
>>> >> masks. I
>>> >> would like to keep the clusters from the thresholding 2.3 which
>>> >> contains at
>>> >> least one voxel which is active in thresholding 3.1 as well. But I
>>> >> would
>>> >> like only these clusters from thresholding 2.3. I don't need
>>> >> clusters which
>>> >> are active in thresholding 2.3, but don't contain any voxels which
>>> are
>>> >> active in thresholding 3.1. So I want to do a new map from clusters
>>> of
>>> >> thresholding 2.3, which have seed with z>3.1.
>>> >>
>>>
>>
>>
>>
>>--
>>G. Salimi-Khorshidi,
>>D.Phil. Student, Dept. of Clinical Neurology, University of Oxford.
>>[log in to unmask] http://www.fmrib.ox.ac.uk/~reza
>>FMRIB Centre, John Radcliffe Hospital,
>>Headington, Oxford, OX3 9DU
>>Tel: +44 (0) 1865 222466 Fax: +44 (0)1865 222717
>>
>
Tibor Auer, M.D. Ph.D.
Biomedizinische NMR Forschungs GmbH
am Max-Planck Institut für Biophysikalische Chemie
Am Fassberg 11
37077 Göttingen
Germany
Phone: +49-(0)551-201 1725
Mobile: +36-(06)20-986-9341
Mail: [log in to unmask]
|