I know, but scipy norm is missing an axis argument...


> On Feb 5, 2014, at 7:26 PM, Eric89GXL <notifications@github.com> wrote:
> 
> In mne/label.py:
> 
> > +    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)
> > +    centered_points = points - center
> > +    distance = np.linalg.norm(centered_points, axis=1)
> BTW @agramfort is against using np.linalg.
> 
> —
> Reply to this email directly or view it on GitHub.