Pierre GM wrote: > Mathew, > Have you tried the solution that was suggested by Angus yesterday on the > numpy > mailing list ? > > http://projects.scipy.org/pipermail/scipy-user/2008-February/015418.html > >>>> import numpy as np >>>> import matplotlib.nxutils as nxutils >>>> polygon=np.array([(0,0),(0,1),(1,1),(1,0)]) >>>> points = np.array([(0.5,0.5),(0.4,1.5)]) >>>> nxutils.points_inside_poly(points, polygon) > array([1, 0], dtype=int32) > > Meaning that the first point is, the second is not.
It seems to me that the output should be boolean, so the line mask = (PyArrayObject *)PyArray_SimpleNew(1,dimensions,PyArray_INT); in nxutils.cpp would become mask = (PyArrayObject *)PyArray_SimpleNew(1,dimensions, NPY_BOOL); Can anyone think of anything this would break, or any disadvantages? Eric ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
