general purpose numerical utilities, eg for computational geometry, that are not available in numpy
inside = pnpoly(x, y, xyverts)
Return 1 if x,y is inside the polygon, 0 otherwise.
A point on the boundary may be treated as inside or outside. See pnpoly
mask = points_inside_poly(xypoints, xyverts)
Return a boolean ndarray, True for points inside the polygon.
A point on the boundary may be treated as inside or outside. See pnpoly