用Bulge保持多段线圆弧连续性

在CAD多段线分割操作中，利用Bulge属性保持圆弧几何连续性的核心在于：精确计算分割点对应的凸度值，并正确更新分割后两条新多段线的顶点序列和凸度数组。这确保了分割点前后的圆弧段在几何上仍然是原圆弧的一部分，从而维持了曲率的连续性。

一、 Bulge属性的几何意义与连续性原理

多段线的Bulge（凸度）定义了两个相邻顶点之间圆弧段的几何信息。对于顶点i到i+1的段：

Bulge = tan(θ/4)，其中θ是该圆弧段对应的包含角（圆心角）。
Bulge的符号表示圆弧方向：正值为逆时针，负值为顺时针。
Bulge的绝对值决定了圆弧的弯曲程度：0表示直线，绝对值越大，圆弧越弯曲。

几何连续性要求 ：当在一个圆弧段上的点P处分割时，原圆弧段A->B被拆分为两个子圆弧段A->P和P->B。这两个子圆弧必须与原圆弧同心、同半径 ，只是包含角发生了变化。因此，需要为A->P和P->B分别计算新的、正确的Bulge值。

二、分割算法步骤与关键计算

假设我们要在由顶点V_i和V_{i+1}定义的圆弧段上一点P_split（该点位于圆弧上）进行分割。以下是保证连续性的具体步骤：

定位与参数计算：
- 首先需确认P_split位于多段线的第i个段（V_i到V_{i+1}）上，并计算出该点在圆弧上的参数t （ 0 <= t <= 1），表示从起点V_i到P_split的弧长占总弧长的比例。这通常通过求解点与圆弧的最短距离或进行弧长参数化得到。
- 获取该段的凸度值bulge_original。
计算原圆弧几何参数：
- 根据V_i、V_{i+1}和bulge_original，计算出圆弧的圆心C 、半径R 、起始角α_start和包含角θ_total 。包含角θ_total可由凸度直接求得：θ_total = 4 * atan(|bulge_original|)，方向由bulge_original的符号决定。
计算分割点参数与新凸度：
- 计算分割点P_split对应的圆心角θ_split = θ_total * t。
- 第一段（V_i -> P_split）的凸度 ：bulge1 = tan(θ_split / 4)。其符号与bulge_original相同。
- 第二段（P_split -> V_{i+1}）的凸度 ：bulge2 = tan((θ_total - θ_split) / 4)。其符号与bulge_original相同。
- 特殊情况 ：如果分割点恰好非常接近端点（t ≈ 0 或 t ≈ 1），则其中一个凸度会接近0（直线），另一个接近原凸度。需要处理数值精度问题。
重构两条新多段线：
- 第一条多段线 ：顶点序列为 [V_0, V_1, ..., V_i, P_split]，凸度序列为原bulge[0]到bulge[i-1]，加上新计算的bulge1。
- 第二条多段线 ：顶点序列为 [P_split, V_{i+1}, ..., V_n]，凸度序列为新计算的bulge2，后接原bulge[i+1]到bulge[n-2]。
- 需注意闭合多段线(Closed)属性的处理：如果原多段线是闭合的，分割后可能需要调整或创建两条独立的多段线。

三、核心代码实现

以下C#（基于.NET及类似AutoCAD API环境）伪代码演示了关键计算步骤：

csharp 复制代码

/// <summary>
/// 计算在给定圆弧段上分割点处的新凸度值，以保持几何连续性。
/// </summary>
/// <param name="startPoint">圆弧段起点 (V_i)</param>
/// <param name="endPoint">圆弧段终点 (V_{i+1})</param>
/// <param name="bulgeOriginal">原凸度值</param>
/// <param name="splitParam">分割点参数t (0到1之间，表示从起点到分割点的弧长比例)</param>
/// <param name="bulge1">输出：第一段（起点->分割点）的新凸度</param>
/// <param name="bulge2">输出：第二段（分割点->终点）的新凸度</param>
public static void CalculateBulgesAtSplit(
    Point2d startPoint, Point2d endPoint, double bulgeOriginal,
    double splitParam, out double bulge1, out double bulge2)
{
    // 1. 计算原圆弧的总包含角（弧度）
    double totalIncludedAngle = 4.0 * Math.Atan(Math.Abs(bulgeOriginal));
    // 保留符号以确定方向
    double sign = Math.Sign(bulgeOriginal);

    // 2. 计算分割点前后的包含角
    double angleToSplit = totalIncludedAngle * splitParam;
    double angleFromSplit = totalIncludedAngle - angleToSplit;

    // 3. 计算新凸度值
    bulge1 = sign * Math.Tan(angleToSplit / 4.0);
    bulge2 = sign * Math.Tan(angleFromSplit / 4.0);

    // 4. 处理数值精度边界情况：当角度极小时，凸度接近0（直线）
    const double tolerance = 1e-10;
    if (Math.Abs(angleToSplit) < tolerance) bulge1 = 0.0;
    if (Math.Abs(angleFromSplit) < tolerance) bulge2 = 0.0;
}

/// <summary>
/// 分割多段线，并保持圆弧连续性。
/// </summary>
/// <param name="pline">待分割的多段线</param>
/// <param name="splitPoint">位于多段线上的分割点</param>
/// <returns>包含两条新多段线的元组</returns>
public static (Polyline, Polyline) SplitPolylineWithBulgeContinuity(Polyline pline, Point3d splitPoint)
{
    // 步骤1: 找到分割点所在的段索引和参数t
    int segmentIndex;
    double paramOnSegment; // 该段上的参数t
    FindSegmentAndParameter(pline, splitPoint, out segmentIndex, out paramOnSegment);

    // 获取该段的起点、终点和原凸度
    Point2d start2d = pline.GetPoint2dAt(segmentIndex);
    Point2d end2d = pline.GetPoint2dAt(segmentIndex + 1);
    double originalBulge = pline.GetBulgeAt(segmentIndex);

    // 步骤2: 计算分割点处的新凸度
    double bulge1, bulge2;
    CalculateBulgesAtSplit(start2d, end2d, originalBulge, paramOnSegment, out bulge1, out bulge2);

    // 步骤3: 创建两条新的多段线
    Polyline pline1 = new Polyline();
    Polyline pline2 = new Polyline();

    // 第一条线：从顶点0到分割点
    for (int i = 0; i <= segmentIndex; i++)
    {
        pline1.AddVertexAt(pline1.NumberOfVertices, pline.GetPoint2dAt(i), pline.GetBulgeAt(i), 0, 0);
    }
    // 在第一条线的末尾添加分割点，并设置凸度bulge1
    pline1.AddVertexAt(pline1.NumberOfVertices, new Point2d(splitPoint.X, splitPoint.Y), bulge1, 0, 0);

    // 第二条线：从分割点到最后一个顶点
    pline2.AddVertexAt(0, new Point2d(splitPoint.X, splitPoint.Y), bulge2, 0, 0);
    for (int i = segmentIndex + 1; i < pline.NumberOfVertices; i++)
    {
        // 注意凸度索引的偏移：原多段线中，顶点i对应的凸度索引是i（对于最后一点，凸度无意义或为0）
        // 当i == segmentIndex + 1时，应该使用原凸度数组中的下一个凸度（如果有）
        double bulgeForPline2 = (i == segmentIndex + 1) ? pline.GetBulgeAt(i) : pline.GetBulgeAt(i - 1);
        pline2.AddVertexAt(pline2.NumberOfVertices, pline.GetPoint2dAt(i), bulgeForPline2, 0, 0);
    }

    // 处理闭合属性：如果原多段线闭合，分割后通常变为两条开放多段线
    // 但根据需求，可能需要特殊处理最后一段与第一段的连接
    // pline1.Closed = false;
    // pline2.Closed = false;

    return (pline1, pline2);
}

// 辅助函数：查找点所在的段及参数（需要实现点在线上的精确投影或参数计算）
private static void FindSegmentAndParameter(Polyline pline, Point3d point, out int segmentIndex, out double param)
{
    segmentIndex = -1;
    param = 0.0;
    // 实现逻辑：遍历多段线的每个段，计算点到该段（直线或圆弧）的最近点及参数t。
    // 可使用API如`GetClosestPointTo`或自行几何计算。
    // 这是一个简化示例，实际实现需处理几何计算和容差。
    double minDist = double.MaxValue;
    for (int i = 0; i < pline.NumberOfVertices; i++)
    {
        Curve seg = pline.GetSegmentAt(i);
        if (seg != null)
        {
            Point3d closest = seg.GetClosestPointTo(point, false);
            double dist = closest.DistanceTo(point);
            if (dist < minDist)
            {
                minDist = dist;
                segmentIndex = i;
                // 获取参数化值（假设曲线有GetParameterAtPoint方法或类似功能）
                // param = seg.GetParameterAtPoint(closest); // 伪代码，需根据实际API调整
            }
        }
    }
}

四、关键注意事项与验证

精度处理：
- 在计算凸度时，特别是当splitParam接近0或1时，tan(θ/4)的计算可能产生数值误差。必须设置合理的容差（如 1e-10），将极小的凸度视为0（直线）。
- 分割点P_split必须精确位于圆弧上。通常通过GetClosestPointTo方法获取线上最近点作为分割点，但需确保该点确实在线段参数范围内。
验证连续性：
- 分割后，应验证新生成的两段圆弧是否与原始圆弧同心同半径。可以通过检查P_split点是否同时满足两个新圆弧的方程，或计算分割点处两个新圆弧的切线方向是否一致来进行验证。
处理多段线属性：
- 分割操作不仅影响顶点和凸度，还可能影响图层、颜色、线型、线宽等属性。新生成的多段线应继承原多段线的这些属性。
- 对于闭合多段线，分割逻辑更复杂，因为最后一个顶点与第一个顶点相连。通常，分割闭合多段线会将其"打开"，生成两条开放的多段线，需要妥善处理最后一段与第一段之间的隐含连接。
性能考虑：
- 在批量处理时，应避免在循环中频繁创建和销毁临时几何对象。可以复用计算出的圆心、半径等中间结果。

总结：通过严格遵循Bulge的几何定义，并基于分割点参数t按比例分配原圆弧的包含角θ，可以准确计算出分割后两段子圆弧的凸度值。这是确保CAD多段线分割后圆弧几何连续性最根本和精确的方法。所有操作必须建立在精确的几何计算和稳健的数值处理基础上。

用Bulge保持多段线圆弧连续性

一、 Bulge属性的几何意义与连续性原理

二、 分割算法步骤与关键计算

三、 核心代码实现

四、 关键注意事项与验证

参考来源

二、分割算法步骤与关键计算

三、核心代码实现

四、关键注意事项与验证