如何获得多边形质心与其外环之间的最大距离？

import math

# 多边形顶点坐标
polygon = [(0, 0), (0, 4), (4, 4), (4, 0)]

# 计算多边形质心
centroid_x = sum(x for x, y in polygon) / len(polygon)
centroid_y = sum(y for x, y in polygon) / len(polygon)
centroid = (centroid_x, centroid_y)

# 初始化最大距离为0
max_distance = 0

# 遍历多边形的边界
for i in range(len(polygon)):
    p1 = polygon[i]
    p2 = polygon[(i + 1) % len(polygon)]  # 下一个顶点

    # 计算当前边界上的点与质心之间的距离
    distance = math.sqrt((centroid_x - p1[0])**2 + (centroid_y - p1[1])**2)

    # 更新最大距离
    if distance > max_distance:
        max_distance = distance

# 输出最大距离
print("多边形质心与外环之间的最大距离为:", max_distance)