Yes, and it leads to a small difference...

> On Feb 5, 2014, at 7:26 PM, Eric89GXL <notifications@github.com> wrote:
> 
> In mne/label.py:
> 
> > +    elif label.subject is None:
> > +        pass
> > +    elif subject != label.subject:
> > +        err = ("The label specifies a different subject (%r) from the subject "
> > +               "parameter (%r)." % label.subject, subject)
> > +        raise ValueError(err)
> > +
> > +    # find the label coordinates on the spherical surface
> > +    surf_fname = '.'.join((label.hemi, 'sphere'))
> > +    surf_path = os.path.join(subjects_dir, subject, "surf", surf_fname)
> > +    surface = _read_surface_geom(surf_path, add_geom=True)
> > +    points = surface['rr'][label.vertices]
> > +
> > +    # find the axis vector
> > +    # center label coordinates
> > +    center = np.mean(points, axis=0)
> Basically I think all you need to do to increase numerical precision is add a conditional here, where for continuous you just do normal = center / norm(center), or (for fs) execute the next six lines. Right?
> 
> —
> Reply to this email directly or view it on GitHub.