Added ability to threshold by number of voxels in reduce - average_within_regions