查找多边形内的最大矩形- Python

def calculate_area(rectangle):
    length = abs(rectangle[0][0] - rectangle[1][0])
    width = abs(rectangle[0][1] - rectangle[2][1])
    return length * width

def is_point_inside_polygon(point, polygon):
    x, y = point
    inside = False
    for i in range(len(polygon)):
        j = (i + 1) % len(polygon)
        if (polygon[i][1] > y) != (polygon[j][1] > y) and \
                x < (polygon[j][0] - polygon[i][0]) * (y - polygon[i][1]) / (polygon[j][1] - polygon[i][1]) + polygon[i][0]:
            inside = not inside
    return inside

def find_max_rectangle_in_polygon(polygon):
    max_area = 0
    max_rectangle = None
    for i in range(len(polygon)):
        for j in range(i + 1, len(polygon)):
            rectangle = [polygon[i], (polygon[i][0], polygon[j][1]), polygon[j], (polygon[j][0], polygon[i][1])]
            if is_point_inside_polygon(rectangle[1], polygon) and is_point_inside_polygon(rectangle[3], polygon):
                area = calculate_area(rectangle)
                if area > max_area:
                    max_area = area
                    max_rectangle = rectangle
    return max_rectangle

# 示例用法
polygon = [(0, 0), (0, 5), (5, 5), (5, 0)]
max_rectangle = find_max_rectangle_in_polygon(polygon)
print("最大矩形的四个角点坐标：", max_rectangle)
print("最大矩形的面积：", calculate_area(max_rectangle))