计算外部点对简单多边形的视角角度及中心指向

核心思路 #

直接取多边形顶点的最小/最大极角会因角度跨0点（如从3π/2到π/2）导致区间判断错误，必须结合多边形的顺/逆时针顶点顺序来确定正确的视角区间。

具体步骤 #

1. 计算顶点相对点A的极角 #

对每个多边形顶点pᵢ，计算向量vec = (pᵢ.x - A.x, pᵢ.y - A.y)，用atan2(vec.y, vec.x)得到极角θᵢ（建议统一转换为[0, 2π]范围，避免负角度的跨区间判断混乱）。

2. 判断多边形的顺逆时针方向 #

通过计算多边形的有向面积确定顺序：

def get_polygon_winding_order(points):
    n = len(points)
    area = 0.0
    for i in range(n):
        x1, y1 = points[i]
        x2, y2 = points[(i+1)%n]
        area += (x1 * y2) - (x2 * y1)
    return "counterclockwise" if area > 0 else "clockwise"

• 面积为正：逆时针（CCW）
• 面积为负：顺时针（CW）

3. 确定正确的视角区间 #

按多边形顶点顺序遍历相邻顶点对(pᵢ, pᵢ₊₁)（pₙ₊₁ = p₁），对应极角(θᵢ, θᵢ₊₁)：

• 若多边形为逆时针且点A在外部：
  • 当θᵢ₊₁ > θᵢ：覆盖区间为[θᵢ, θᵢ₊₁]
  • 当θᵢ₊₁ < θᵢ：覆盖区间为[θᵢ, 2π] ∪ [0, θᵢ₊₁]
• 若多边形为顺时针且点A在外部：
  • 当θᵢ₊₁ < θᵢ：覆盖区间为[θᵢ₊₁, θᵢ]
  • 当θᵢ₊₁ > θᵢ：覆盖区间为[θᵢ₊₁, 2π] ∪ [0, θᵢ]
    （注：非自交多边形外点的覆盖视角必然是单一连续区间，无需处理多区间情况）

4. 计算总视角角度Θ #

• 区间不跨0：Θ = θ_end - θ_start
• 区间跨0：Θ = (2π - θ_start) + θ_end

5. 计算视角中心方向 #

由于角度是环形量，直接取区间端点算术平均会出错（如区间3π/2 ~ π/2的算术平均为π，但实际中心是0），正确方法：

• 不跨0：中心角度为(θ_start + θ_end) / 2
• 跨0：中心角度为(θ_start + θ_end + 2π) / 2 % 2π（转换为[0,2π]范围）
  也可通过将区间端点转为单位向量求和后再取极角，结果一致。

关键注意事项 #

• 极角计算必须统一范围，避免跨区间判断混乱
• 务必结合多边形顺逆时针顺序，否则会得到完全相反的视角区间
• 点A在多边形外的前提保证了覆盖视角是单一连续区间，无需考虑多段分割

内容的提问来源于stack exchange，提问作者Earlh21