Hi Jean-Marc,
> I'm looking for a file where I could find t or z values of the totality of
> the pixels in a cluster?
> My aim is to obtain an average stat value in a cluster..
I wrote a little script once, using spm-functions, to do something like
this. It is rather crude but seems to work for me. It is attached.
> Thanks for yours answers.
Glad to be of service ;)
Best,
Marko
--
=====================================================================
Marko Wilke (Dr.med./M.D.)
[log in to unmask]
Universitäts-Kinderklinik University Children's Hospital
Abt. III (Neuropädiatrie) Dept. III (Pediatric neurology)
Hoppe-Seyler-Str. 1, D - 72076 Tübingen
Tel.: (+49) 07071 29-83416 Fax: (+49) 07071 29-5473
=====================================================================
% Little script to calculate the mean voxel intensity
% in an image after masking with another image (e.g.,
% a cluster defining my region of interest);
% implementation by Marko Wilke, USE AT YOUR OWN RISK!
% Get files (need to be in register!)
mask = spm_read_vols(spm_vol(spm_get(1,'*.img','Select Masking Image', pwd)));
contrast = spm_read_vols(spm_vol(spm_get(1,'*.img','Select Contrast Image to mask')));
% mask, clean masked data
masked_contrast = (mask>0).*contrast;
range = reshape(masked_contrast,1,[]);
range(find(range == 0)) = [];
range = sort(range);
% output (very crude)
disp([' Mean voxel value within the cluster is : ' num2str(mean(range))]);
disp([' number of voxels : ' num2str(size(range,2))]);
disp([' maximum value : ' num2str(range(size(range,2)))]);
disp([' minimum value : ' num2str(range(1))]);
% clean up, leave results
clear mask contrast masked_contrast
|