waveforms = fast_vol2mat(MRIread('waveform.nii.gz'));

will give me a v x time matrix of values from my functional values.

In the case of fsaverage, I believe I can expect somewhere around 160K rows, 1 per vertex, for each hemisphere.

I was wondering how I might be able to determine the neighbors (i.e., vertices connected by an edge) for any vertex v?

Thanks,

Chris