Hi Cesar,
To be sure, do you want to correlate each voxel with its local neighborhood, radially? If so, that is what @radial_correlate does. It has fast, Gaussian methods, as well as a slow 3dLocalStat method. If there is an aspect where the correlation is restricted to good ones, this does not do it.
Does that do what you want?
- rick