def get_crop_coords(coords, bb):
    """
    Get all coordinates in coords which are within bounding box bb.

    Parameters
    ----------
    coords : numpy array of coordinates
    bb : bounding box tuple of ints structured as (x_min, x_max, y_min, y_max)

    Returns
    -------
    crop_coords : all coordinates in coords which are within bb
    """
    crop_coords = []
    for coord in coords:
        if (coord[0] >= bb[0]) and (coord[0] <= bb[1]):
            if (coord[1] >= bb[2]) and (coord[1] <= bb[3]):
                crop_coords.append(coord)
    return crop_coords