自由手绘曲线弯折处角度计算不准确问题求助（附代码）

需求为计算自由手绘曲线（非直线/多段线，复杂度更高）各弯折位置的角度（单位：度），并区分角度大于45°与小于45°的位置。当前采用Ramer-Douglas-Peucker算法简化曲线，并用np.arccos计算连续线段间的角度，但得到的角度值不准确。

原代码片段：

....
largest_contour = max(filter_coontour, key=cv2.contourArea)

# Simplify the curve using the Ramer-Douglas-Peucker algorithm
epsilon = 0.001 * cv2.arcLength(largest_contour, True)
simplified_curve = cv2.approxPolyDP(largest_contour, epsilon, True)

# Get the coordinates of the points along the simplified curve
curve_points = simplified_curve.reshape(-1, 2)

for i in range(1, len(curve_points)):
     vector = curve_points[i] - curve_points[i - 1]
     if prev_vector is not None:
           dot = np.dot(vector, prev_vector)
           norm = np.linalg.norm(vector) * np.linalg.norm(prev_vector)       
           angle = np.arccos(dot / norm) * 180 / np.pi
      print("angle", angle)
....

问题根源 #

1. 角度逻辑偏差：当前计算的是相邻线段的最小夹角（0°-180°），但手绘曲线的弯折角度应是转向角（即曲线在该点的外角，更符合直观弯折感知），且代码未初始化prev_vector，会跳过第一个有效角度计算。
2. 曲线简化参数不合理：epsilon=0.001*弧长过小，会保留大量冗余点，导致计算出无意义的小角度；过大则会丢失真实弯折点。
3. 数值精度问题：dot/norm可能因浮点误差超出[-1,1]范围，导致np.arccos返回NaN，破坏计算准确性。

修正方案 #

1. 修复角度计算逻辑 #

使用np.arctan2计算向量方向角，通过方向角差得到转向角（考虑顺时针/逆时针转向），再转换为0-180°范围的弯折角度；同时初始化prev_vector，确保第一个有效线段对的角度被计算。

2. 优化曲线简化参数 #

调整epsilon至0.01-0.1*弧长区间，根据曲线复杂度灵活调整，平衡保留关键弯折点与去除冗余点的需求；可提前对轮廓做高斯平滑，减少噪声干扰。

3. 处理数值精度问题 #

对dot/norm做截断处理，限制在[-1,1]范围内，避免np.arccos报错。

修正后代码示例 #

import cv2
import numpy as np

# 修正原代码拼写错误：filter_coontour → filter_contour
largest_contour = max(filter_contour, key=cv2.contourArea)

# 可选：平滑轮廓减少噪声
# blurred_contour = cv2.GaussianBlur(largest_contour, (3,3), 0)

# 优化曲线简化参数
epsilon = 0.01 * cv2.arcLength(largest_contour, True)
simplified_curve = cv2.approxPolyDP(largest_contour, epsilon, True)
curve_points = simplified_curve.reshape(-1, 2)

angle_results = []
prev_vector = None

for i in range(1, len(curve_points)):
    curr_vector = curve_points[i] - curve_points[i-1]
    if prev_vector is not None:
        # 计算向量方向角
        prev_angle = np.arctan2(prev_vector[1], prev_vector[0])
        curr_angle = np.arctan2(curr_vector[1], curr_vector[0])
        # 计算转向角并转换为0-180°范围
        angle_diff = np.abs((curr_angle - prev_angle) * 180 / np.pi)
        bend_angle = angle_diff if angle_diff <= 180 else 360 - angle_diff
        
        # 区分角度类型
        angle_type = "大于45°" if bend_angle > 45 else "小于等于45°"
        angle_results.append({
            "point": curve_points[i],
            "bend_angle(°)": round(bend_angle, 2),
            "type": angle_type
        })
        print(f"点{curve_points[i]}处弯折角度: {round(bend_angle,2)}° ({angle_type})")
    prev_vector = curr_vector

# 批量输出结果
print("\n所有弯折角度汇总：")
for res in angle_results:
    print(res)

额外注意 #

• 若曲线为闭合轮廓，需额外处理最后一个点与第一个点的向量夹角；
• 可根据需求切换内角/外角计算，上述代码默认计算转向外角，更贴合手绘弯折的直观感受。

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