IntersectWith怎么排除多余的交点

你有种再说一遍

你有种再说一遍 发表于 2022-4-22 05:40
首先是内存有极限,这个极限来自于双精度浮点数,
而曲线求交点是经过多次牛顿迭代来逼近交点,当逼近的时 ...

那就不平移,概率出错,留点bug🤣

pmq

IntersectWith 投影平面 Plane 参数的值应该输入什么

你有种再说一遍

pmq 发表于 2024-9-20 09:26
IntersectWith 投影平面 Plane 参数的值应该输入什么

求交点函数的考虑了三维的.
投影的plan一般情况就是XOY平面的法向量vector.Z==(0,0,1)

wang2006zhi

2.     /// <summary>
   
3.     /// 获取两个曲线的交点集合（支持大坐标偏移优化）
   
4.     /// </summary>
   
5.     /// <param name="curve1">第一个曲线（线段或圆弧可延伸，多段线不可延伸）</param>
   
6.     /// <param name="curve2">第二个曲线（线段或圆弧可延伸，多段线不可延伸）</param>
   
7.     /// <param name="intersectOption">相交处理方式：
   
8.     /// <para>OnBothOperands - 两者都不延伸</para>
   
9.     /// <para>ExtendBoth - 两者都延伸</para>
   
10.     /// <para>ExtendThis - 仅延伸当前曲线(curve1)</para>
    
11.     /// <para>ExtendArgument - 仅延伸参数曲线(curve2)</para></param>
    
12.     /// <returns>输出的交点集合</returns>
    
13.     public static Point3dCollection GetInterPoint(
    
14.         this Curve curve1,
    
15.         Curve curve2,
    
16.         Intersect intersectOption = Intersect.OnBothOperands)
    
17.     {
    
18.         var intersectionPoints = new Point3dCollection();
    
19. 
    
20.         // 第一阶段：快速包围盒检查
    
21.         var box1 = curve1.GeometricExtents;
    
22.         var box2 = curve2.GeometricExtents;
    
23.         
    
24.         if (!box1.IsIntersect(box2))
    
25.             return intersectionPoints;
    
26. 
    
27.         var plan = new Plane(Point3d.Origin, Vector3d.ZAxis);
    
28.         // 第二阶段：大坐标偏移处理
    
29.         const double coordinateThreshold = 1e9; // 坐标偏移阈值
    
30. 
    
31.         // 计算合并包围盒中心
    
32.         var combinedBox = box1.CombinedBox(box2);
    
33. 
    
34.         var centroid = combinedBox.GetCenterPoint();
    
35. 
    
36.         // 判断是否需要坐标偏移
    
37.         bool needTransform =
    
38.             Math.Abs(centroid.X) > coordinateThreshold ||
    
39.             Math.Abs(centroid.Y) > coordinateThreshold;
    
40. 
    
41.         if (needTransform)
    
42.         {
    
43.             // 构建平移矩阵
    
44.             Vector3d offset = -centroid.GetAsVector();
    
45.             Matrix3d transform = Matrix3d.Displacement(offset);
    
46.             Matrix3d inverseTransform = Matrix3d.Displacement(offset.Negate());
    
47. 
    
48.             // 克隆并平移曲线
    
49.             using Curve transCurve1 = (Curve)curve1.Clone();
    
50.             using Curve transCurve2 = (Curve)curve2.Clone();
    
51.             transCurve1.TransformBy(transform);
    
52.             transCurve2.TransformBy(transform);
    
53. 
    
54.             // 在平移后坐标系求交
    
55.             var tempPoints = new Point3dCollection();
    
56.             transCurve1.IntersectWith(
    
57.                 transCurve2,
    
58.                 intersectOption,
    
59.                 plan,
    
60.                 tempPoints,
    
61.                 IntPtr.Zero,
    
62.                 IntPtr.Zero);
    
63. 
    
64.             // 坐标逆向转换
    
65.             foreach (Point3d pt in tempPoints)
    
66.             {
    
67.                 intersectionPoints.Add(pt.TransformBy(inverseTransform));
    
68.             }
    
69.         }
    
70.         else
    
71.         {
    
72.             // 直接求交
    
73.             curve1.IntersectWith(
    
74.                 curve2,
    
75.                 intersectOption,
    
76.                 plan,
    
77.                 intersectionPoints,
    
78.                 IntPtr.Zero,
    
79.                 IntPtr.Zero);
    
80.         }
    
81. 
    
82.         return intersectionPoints;
    
83.     }

    /// <summary><br />     /// 获取两个曲线的交点集合（支持大坐标偏移优化）<br />     /// </summary><br />     /// <param name="curve1">第一个曲线（线段或圆弧可延伸，多段线不可延伸）</param><br />     /// <param name="curve2">第二个曲线（线段或圆弧可延伸，多段线不可延伸）</param><br />     /// <param name="intersectOption">相交处理方式：<br />     /// <para>OnBothOperands - 两者都不延伸</para><br />     /// <para>ExtendBoth - 两者都延伸</para><br />     /// <para>ExtendThis - 仅延伸当前曲线(curve1)</para><br />     /// <para>ExtendArgument - 仅延伸参数曲线(curve2)</para></param><br />     /// <returns>输出的交点集合</returns><br />     public static Point3dCollection GetInterPoint(<br />         this Curve curve1,<br />         Curve curve2,<br />         Intersect intersectOption = Intersect.OnBothOperands)<br />     { <br />         var intersectionPoints = new Point3dCollection();<br /> <br />         // 第一阶段：快速包围盒检查<br />         var box1 = curve1.GeometricExtents;<br />         var box2 = curve2.GeometricExtents;<br />         <br />         if (!box1.IsIntersect(box2))<br />             return intersectionPoints;<br /> <br />         var plan = new Plane(Point3d.Origin, Vector3d.ZAxis);<br />         // 第二阶段：大坐标偏移处理<br />         const double coordinateThreshold = 1e9; // 坐标偏移阈值<br /> <br />         // 计算合并包围盒中心<br />         var combinedBox = box1.CombinedBox(box2);<br /> <br />         var centroid = combinedBox.GetCenterPoint();<br /> <br />         // 判断是否需要坐标偏移<br />         bool needTransform =<br />             Math.Abs(centroid.X) > coordinateThreshold ||<br />             Math.Abs(centroid.Y) > coordinateThreshold;<br /> <br />         if (needTransform)<br />         {<br />             // 构建平移矩阵<br />             Vector3d offset = -centroid.GetAsVector();<br />             Matrix3d transform = Matrix3d.Displacement(offset);<br />             Matrix3d inverseTransform = Matrix3d.Displacement(offset.Negate());<br /> <br />             // 克隆并平移曲线<br />             using Curve transCurve1 = (Curve)curve1.Clone();<br />             using Curve transCurve2 = (Curve)curve2.Clone();<br />             transCurve1.TransformBy(transform);<br />             transCurve2.TransformBy(transform);<br /> <br />             // 在平移后坐标系求交<br />             var tempPoints = new Point3dCollection();<br />             transCurve1.IntersectWith(<br />                 transCurve2,<br />                 intersectOption,<br />                 plan,<br />                 tempPoints,<br />                 IntPtr.Zero,<br />                 IntPtr.Zero);<br /> <br />             // 坐标逆向转换<br />             foreach (Point3d pt in tempPoints)<br />             {<br />                 intersectionPoints.Add(pt.TransformBy(inverseTransform));<br />             }<br />         }<br />         else<br />         {<br />             // 直接求交<br />             curve1.IntersectWith(<br />                 curve2,<br />                 intersectOption,<br />                 plan,<br />                 intersectionPoints,<br />                 IntPtr.Zero,<br />                 IntPtr.Zero);<br />         }<br /> <br />         return intersectionPoints;<br />     }