SymPy多边形角度代码符号化报错处理及Lambda化方法

问题原因 #

当使用符号坐标创建SymPy多边形时，SymPy无法确定多边形的朝向（顺时针/逆时针），而Polygon.angles属性依赖朝向信息计算内角，因此抛出Can't determine orientation错误。

解决方案 #

绕过Polygon.angles属性，手动通过向量运算推导内角的符号表达式，同时处理多边形朝向的判断逻辑，再完成Lambda化。

代码实现 #

import sympy as sp

# 定义所有符号坐标变量
x0, y0, x1, y1, x2, y2, x3, y3 = sp.symbols('x0 y0 x1 y1 x2 y2 x3 y3')

# 定义多边形顶点
pt1 = sp.Point(x0, y0)
pt2 = sp.Point(x1, y1)
pt3 = sp.Point(x2, y2)

# 计算顶点pt2处的两个邻边向量（指向相邻顶点）
v_prev = pt1 - pt2
v_next = pt3 - pt2

# 计算向量点积与模长，用于推导夹角
dot_product = v_prev.dot(v_next)
mag_prev = v_prev.norm()
mag_next = v_next.norm()

# 计算向量夹角的余弦值与基础夹角
cos_theta = dot_product / (mag_prev * mag_next)
base_theta = sp.acos(cos_theta)

# 通过连续三点的叉积判断多边形朝向
cross_product = (pt2 - pt1).cross(pt3 - pt2)

# 根据朝向定义内角：逆时针取优角，顺时针取劣角，共线返回NaN
interior_angle = sp.Piecewise(
    (2 * sp.pi - base_theta, cross_product > 0),
    (base_theta, cross_product < 0),
    (sp.nan, cross_product == 0)
)

# 将符号表达式Lambda化为可计算数值的函数
angle_calculator = sp.lambdify(
    [x0, y0, x1, y1, x2, y2, x3, y3],
    interior_angle,
    modules='numpy'
)

# 测试原数值案例
print(angle_calculator(0, 0, 1, 1, 0, 1, 1, 0))  # 输出：5.49778714378214

说明 #

1. 向量运算：通过邻边向量的点积计算夹角余弦，再用反余弦得到基础夹角，不依赖SymPy的多边形朝向判断逻辑。
2. 朝向处理：利用连续三点的叉积符号区分多边形朝向，自动调整内角为优角或劣角，覆盖简单多边形与交叉多边形的情况。
3. Lambda化：将符号表达式转换为高效数值函数，保留原代码的参数结构，支持直接传入坐标数值计算。

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