Hi Koushik,
if you look into the archive you will find several code .. I use this
one (already posted but I can remember who)
P=spm_get(1,'*.img','Select statistic image');
V=spm_vol(P);
[x,y,z] = ndgrid(1:V.dim(1),1:V.dim(2),0);
for i=1:V.dim(3),
z = z + i;
tmp = spm_sample_vol(V,x,y,z,0);
msk = find(tmp~=0 & finite(tmp));
if ~isempty(msk),
tmp = tmp(msk);
xyz1=[x(msk)'; y(msk)'; z(msk)'; ones(1,length(msk))];
xyzt=V.mat(1:3,:)*xyz1;
for j=1:length(tmp),
fprintf('%.4g %.4g %.4g\t%g\n',xyzt(1,j),xyzt(2,j),xyzt(3,j),tmp(j));
end;
end;
end;
you have all the voxel in the MNI and the corresponding value
hope it helps
best
- cyril
>hi SPMers,
> Is there a way I can list the activated voxels and the
>activations in a t-map without a design matrix? Is there any
>extension or code that I can use to get a list of the
>activated voxels.
>
>Thank You
>Regards
>Koushik Athreya Govindarajan
>
>
>***************************************
>
>Graduate Research Assistant
>WSU MR Research Facility
>Harper University Hospital
>3990 John R
>Detroit, MI - 48201.
>phone(off): 313-993-8821
> (home): 313-832-7149
> (cell): 313-2472336
>Email: [log in to unmask]
>
>***************************************
>
>
>
|