RF: First pass at separating out the single voxel case from estimation.

Trying to break down the long function `compute_prf_estimate` and make it
easier to use (?).