def get_points_inside_bbox(x, y, bbox):
x1, y1, x2, y2 = bbox.xmin, bbox.ymin, bbox.xmax, bbox.ymax<EOL>x_in = np.logical_and(x>x1, x<x2)<EOL>y_in = np.logical_and(y>y1, y<y2)<EOL>return np.asarray(np.nonzero(x_in & y_in)[<NUM_LIT:0>])<EOL>
Return the indices of points inside the given bbox.