def _find_box_cartesian(latlon, boxes):
    """Return the box which encloses the coordinates."""
    from shapely.geometry import Point
    p = Point(latlon[::-1])
    for box in boxes:
        poly = box['poly']
        if p.within(poly):
            return box
