曲线和面的交点

admin1

1. void findIntersectionWithSpline(AcDbSpline* pSpline, const AcGePlane& plane)
2. {
3.     if (!pSpline)
4.         return;
6.     // 获取Spline的几何表示
7.     AcGeSplineCurve3d splineGeom;
8.     pSpline->getSplineGeometry(splineGeom);
10.     // 交点数组
11.     AcGePoint3dArray intersectionPts;
13.     // 交点查找逻辑可能需要根据具体需求实现，这里提供一个简化的思路
14.     double paramStart = splineGeom.startParameter();
15.     double paramEnd = splineGeom.endParameter();
16.     double tol = 1e-6; // 容差值，根据需要调整
18.     // 数值方法求解交点，这里简化处理，实际情况可能需要更精细的迭代
19.     for(double t = paramStart; t <= paramEnd; t += (paramEnd - paramStart) / 100.0) // 简单分段处理
20.     {
21.         AcGePoint3d pt = splineGeom.evalPoint(t);
22.         if (plane.normal().isParallelTo(pt - plane.origin()) == Adesk::kFalse) // 判断点是否大致位于平面上
23.         {
24.             // 进一步检查点是否确实在平面上（考虑容差）
25.             if (plane.distanceTo(pt) <= tol)
26.             {
27.                 intersectionPts.append(pt);
28.             }
29.         }
30.     }
32.     // 处理交点数组intersectionPts
33.     // ...
34. }
38. static void findIntersectionByDichotomy(AcDbCurve* pCurve, const AcGePlane& plane, double tol = 1e-7)
39. {
40.         if (!pCurve)
41.                 return;
43.         // 获取Spline的几何表示
44.        
46.         double paramStart;
47.         pCurve->getStartParam(paramStart);
48.        
49.        
50.         double paramEnd ;
51.         pCurve->getEndParam(paramEnd);
53.         double midParam;
55.         while ((paramEnd - paramStart) > tol)
56.         {
57.                 midParam = (paramStart + paramEnd) / 2.0;
59.                 // 计算中点位置
60.                 AcGePoint3d midPoint;
61.                 pCurve->getPointAtParam(midParam,midPoint);
62.        
64.                 // 判断中点与平面的相对位置，调整搜索区间
65.                 double distToPlane = plane.distanceTo(midPoint);
66.                 if (distToPlane < tol)
67.                 {
68.                         // 点在平面上，记录交点并继续在两侧查找可能的其他交点
69.                         // 注意：这里直接记录了交点，实际应用中可能需要进一步逻辑处理连续交点的情况
70.                         // midPoint即为一交点
71.                         paramStart = midParam; // 继续在左侧搜索
72.                         paramEnd = midParam; // 也在右侧搜索
73.                 }
74.                 else {
75.                         // 使用向量点积判断方向，但需注意，这里逻辑需要明确化
77.                         AcGePoint3d staPoint,endPoint;
78.                         pCurve->getStartPoint(staPoint);
79.                         pCurve->getEndPoint(staPoint);
81.                         AcGeVector3d startVec = midPoint - staPoint;
82.                         AcGeVector3d endVec = midPoint - endPoint;
83.                         double dotStart = plane.normal().dotProduct(startVec);
84.                         double dotEnd = plane.normal().dotProduct(endVec);
86.                         // 如果一个点在平面上方，一个在下方，说明交点在区间内
87.                         if ((dotStart * dotEnd) < 0) {
88.                                 paramEnd = midParam; // 缩小区间在左侧
89.                         } else {
90.                                 paramStart = midParam; // 缩小区间在右侧
91.                         }
92.                 }
93.         }
94. }
102. static void findIntersectionWithSpline(AcDbCurve* pCurve, const AcGePlane& plane)
103. {
104.         if (!pCurve)
105.                 return;
107.         // 获取Spline的几何表示
109.         double paramStart;
110.         pCurve->getStartParam(paramStart);
113.         double paramEnd ;
114.         pCurve->getEndParam(paramEnd);
116.         // 交点数组
117.         AcGePoint3dArray intersectionPts;
119.         double tol = 1e-6; // 容差值，根据需要调整
121.         // 主循环，遍历初始参数区间
122.         while (paramStart <= paramEnd)
123.         {
124.                 // 简单分段处理，初步查找交点候选
125.                 for(double t = paramStart; t <= paramEnd; t += (paramEnd - paramStart) / 100.0)
126.                 {
127.                         AcGePoint3d pt ;
128.             pCurve->getPointAtParam(t,pt);
130.                         // 判断点是否接近平面上
131.                         if (std::abs(plane.distanceTo(pt)) <= tol)
132.                         {
133.                                 intersectionPts.append(pt);
134.                                 // 一旦发现交点，锁定当前区间两端点，进行更细粒度的搜索
135.                                 double newParamStart = max(paramStart, t - (paramEnd - paramStart) / 100.0);
136.                                 double newParamEnd =   min(paramEnd, t + (paramEnd - paramStart) / 100.0);
137.                                 // 递归细分并查找该区间内的交点
138.                                 findIntersectionSubdivided(pCurve, plane, newParamStart, newParamEnd, tol, intersectionPts);
139.                                 // 发现交点后，移动start参数避免重复查找
140.                                 paramStart = newParamEnd;
141.                                 // 退出循环，避免在当前粗略区间重复检查
142.                                 break;
143.                         }
144.                 }
145.                 // 缩小搜索区间，避免已查找到的交点区间重复检查
146.                 paramEnd = paramStart;
147.         }
149. }
153. // 辅助函数，用于递归细分区间查找交点
154. static void findIntersectionSubdivided(const AcDbCurve *splineGeom, const AcGePlane& plane, double startParam, double endParam, double tol, AcGePoint3dArray& intersectionPts)
155. {
156.         if (endParam - startParam <= tol)
157.                 return;
159.         for(double t = startParam; t <= endParam; t += (endParam - startParam) / 100.0)
160.         {
161.                  
163.                 AcGePoint3d pt ;
164.                 splineGeom->getPointAtParam(t,pt);
166.                 if (std::abs(plane.distanceTo(pt)) <= tol)
167.                 {
168.                         intersectionPts.append(pt);
169.                         // 继续细分找到的交点周围区域
170.                         findIntersectionSubdivided(splineGeom, plane, max(startParam, t - (endParam - startParam) / 100.0), min(endParam, t + (endParam - startParam) / 100.0), tol, intersectionPts);
171.                         break;
172.                 }
173.         }
174. }
184. static void findIntersectionByDichotomy(AcDbCurve* pCurve, const AcGePlane& plane, double tol = 1e-7)
185. {
186.         if (!pCurve)
187.                 return;
189.         // 获取Spline的几何表示
190.        
192.         double paramStart;
193.         pCurve->getStartParam(paramStart);
194.        
195.        
196.         double paramEnd ;
197.         pCurve->getEndParam(paramEnd);
199.         double midParam;
201.         while ((paramEnd - paramStart) > tol)
202.         {
203.                 midParam = (paramStart + paramEnd) / 2.0;
205.                 // 计算中点位置
206.                 AcGePoint3d midPoint;
207.                 pCurve->getPointAtParam(midParam,midPoint);
208.        
210.                 // 判断中点与平面的相对位置，调整搜索区间
211.                 double distToPlane = plane.distanceTo(midPoint);
212.                 if (distToPlane < tol)
213.                 {
214.                         // 点在平面上，记录交点并继续在两侧查找可能的其他交点
215.                         // 注意：这里直接记录了交点，实际应用中可能需要进一步逻辑处理连续交点的情况
216.                         // midPoint即为一交点
217.                         paramStart = midParam; // 继续在左侧搜索
218.                         paramEnd = midParam; // 也在右侧搜索
219.                 }
220.                 else {
221.                         // 使用向量点积判断方向，但需注意，这里逻辑需要明确化
223.                         AcGePoint3d staPoint,endPoint;
224.                         pCurve->getStartPoint(staPoint);
225.                         pCurve->getEndPoint(staPoint);
227.                         AcGeVector3d startVec = midPoint - staPoint;
228.                         AcGeVector3d endVec = midPoint - endPoint;
229.                         double dotStart = plane.normal().dotProduct(startVec);
230.                         double dotEnd = plane.normal().dotProduct(endVec);
232.                         // 如果一个点在平面上方，一个在下方，说明交点在区间内
233.                         if ((dotStart * dotEnd) < 0) {
234.                                 paramEnd = midParam; // 缩小区间在左侧
235.                         } else {
236.                                 paramStart = midParam; // 缩小区间在右侧
237.                         }
238.                 }
239.         }
240. }

复制代码