RF+ENH: in the case of a radius of zero only consider the node itself. This speeds things up tremendously in the case of just finding voxels intersecting with the surfaces
DOC: minor