area = (b[:, 2] - b[:, 0] + 1) * (b[:, 3] - b[:, 1] + 1)