Hi All,
I’ve added a function to complement nxutils.points_inside_poly (which takes a sequence of points but just one polygon) called point_inside_polys (which takes one point but a sequence of polygons). All of the polygons must have the same number of vertices though. I find it useful for finding out which poly conatins a point of interest to map back to the original array for pcolor, pcolormesh etc.
Patch is attached, comments welcome. Tested with 0.87.5 on OSX, w/numpy
Cheers,
Seth
nxutils.patch (4.04 KB)