afni.nimh.nih.gov

Note that you can calculate whole brain correlation maps with 3dNetCorr; so, you would get a 503x503 matrix of those 'regions', and then you would get 503 whole brain correlation maps which is essentially your 503 x Nvoxels; and as Rick noted, Nvoxels might be really large.

--pt