def intersection(bbox1, bbox2):
    """
    Return a bbox of the intersection of bbox1 and bbox2.

    """
    llx = max(bbox1[0], bbox2[0])
    lly = max(bbox1[1], bbox2[1])
    urx = min(bbox1[2], bbox2[2])
    ury = min(bbox1[3], bbox2[3])

    return llx, lly, urx, ury