P=spm_select(Inf,'nifti','Select scans'); Nii=nifti(P); X=cat(4,Nii.dat); Then you can read out the value of voxel i, j, k by: x=squeeze(X(i,j,k,:)) Best regards, -John On 26 October 2011 13:39, Pavel Sokolov <[log in to unmask]> wrote: > Hello everybody! > I need help. > > I want to get intensity of BOLD-signal in preproccessed images with respect > to time ( I=f(t) ) for every voxel in the volume to analyze synchronously > activated voxels. > Tell please if any possibility to do this by spm or, if not, how to extract > the functions for every voxel to use its in other statistical sortware. > > Thanks a lot, > Pavel.